KR20150065302A - Method deciding 3-dimensional position of landsat imagery by Image Matching - Google Patents
Method deciding 3-dimensional position of landsat imagery by Image Matching Download PDFInfo
- Publication number
- KR20150065302A KR20150065302A KR1020130150460A KR20130150460A KR20150065302A KR 20150065302 A KR20150065302 A KR 20150065302A KR 1020130150460 A KR1020130150460 A KR 1020130150460A KR 20130150460 A KR20130150460 A KR 20130150460A KR 20150065302 A KR20150065302 A KR 20150065302A
- Authority
- KR
- South Korea
- Prior art keywords
- image
- matching
- satellite
- modeling
- point
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/32—Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
Description
The present invention processes a satellite image through a GCP (Ground Control Point) and a DEM (Digital Elevation Model) using a reference point and a digital map for three-dimensional modeling of a satellite image, . However, in the present invention, a satellite image is processed using an RPC (Rational Polynomial Coefficient) coefficient including a satellite attitude and position information according to the characteristics and types of the satellite image, and in order to produce it as a stereoscopic image, And then performs image matting to perform 3D modeling of the satellite image using the DEM obtained through the stereo image processing, thereby producing an orthoimage image.
Recently, high-resolution satellite imagery has emerged and the rapid development of GIS technology has enabled various types of terrain information to be acquired and managed systematically. High-resolution optical satellites such as IKONOS, QuickBird, SPOT and KOMSAT The use of RADARSAT, which is a geographical information based data, is used for city management, various construction management, resource management, and environmental management. It has become necessary to construct information system for various purposes by linking information data with GIS, and integration of various information obtained by using low-resolution resolution satellite image and high resolution image data and integration of multiple sensor data can be utilized Various researches and inventions are in progress to make it very diverse. .
Due to the development of such high-resolution observing satellites, it has rapidly grown to the level of technology capable of replacing aerial photographing as well as local surveying. Digital orthoimage images produced using high-resolution satellites are the spatial information that can represent the real world most accurately and are used to produce various map images and user-centered thematic maps. Even if you have not been able to see the situation in your area of interest.
In the case of satellite image processing technology, development of tools and module technology for image processing centering on developed countries have been developed in a large part, and processing technology has been progressing with the development of payload and sensor technology. The most important aspect of processing technology for utilizing future satellite images is the management and utilization of vast amount of construction information, the automation of processing, and the development and acquisition of accuracy verification technology.
When the technology that can improve the utilization of satellite images is developed as described above, not only the development of technology such as ortho image production but also the improvement of the utilization of the national security dimension by improving the image acquisition and processing technology for the non-access areas Will be.
Conventional image modeling techniques are general modeling and physical modeling. In general modeling, DLT (Direct Linear Transformation) does not need an initial value to perform modeling, and calculation process is simple using modified form of collinear conditional expression, but it is less accurate than physical modeling. Modeling using the other polynomial, like the DLT, does not need an initial value to calculate the model. Depending on the situation, several types such as 1, 2, and so on are used and the calculation process is simple. However, . We used a method of modeling the external facial expressions for each line to be scanned using a rigorous adjustment method, which is a modeling technique using satellite orbital elements, with a polynomial over time, and determining a large number of external facial expressions with a small number of ground reference points. In this case, the first to third order formulas are used for the external facial expression element based on the collinear condition expression. However, if the information provided in the auxiliary data is accurate, it is possible to obtain very accurate results. However, There is a drawback in that it takes a lot of time and cost for image processing.
Recently, the processing of satellite image satellite image is using RFM image processing technique using RPC data. This method compensates the disadvantages of conventional general modeling method and physical modeling method. There is no need, and depending on the situation, various forms such as 1, 2, etc. are used. The calculation process is simple and complicated according to the order, and it is more accurate than the general modeling, but it has a disadvantage of saving the time and cost for the image processing although the accuracy is somewhat lower than the strict method. can do.
Previously, for the 3D modeling of the satellite images, a lot of cost and reference data were required to process the satellite image through the GCP and DEM extracted using the reference point and the digital map. However, in the present invention, the image is processed by the RFM method using the RPC coefficient including the satellite attitude and position information according to the characteristics and types of the satellite image, and in order to produce it as a stereoscopic image, Dimensional modeling of the satellite image using the DEM obtained by extracting the same points among the stereo images using the extracted points.
In order to solve this problem, the present invention processes a satellite image through an RFM method using RPC data, extracts the same points from the processed left and right stereoscopic images, and performs three-dimensional modeling of the satellite image by using image matching I want to.
In order to achieve the above object, the present invention adopts an RFM modeling method using improved RPC data rather than the existing modeling technique, extracts DEMs of overlapping regions through image matching using the same point extraction technique of two images A method of performing three-dimensional modeling, the method comprising: (a) acquiring a stereo satellite image and RPC data through satellite imaging; (b) Image modeling by RFM method using RPC data for image modeling; (c) extracting an identical point between stereoscopic images; (d) performing image matching of two images based on the same point extraction result; (e) extracting a DEM (Digital Elevation Model) from the Stereo image subjected to image matching; (f) A method for three-dimensional positioning through orthoimage production using the modeled satellite image and DEM data.
According to the present invention, it is possible to further increase the usability of the satellite image, which is being utilized in the military and civilian fields, because the cost and time for orthoimage production due to the reduction of the processing time and cost of the satellite image are shortened .
According to the present invention, it is possible to perform three-dimensional modeling without reference data constructed in advance. In particular, when the image matching technique proposed in the present invention is used, it is possible to improve the accuracy of the DEM, And the productivity is higher than that of the conventional method, thereby realizing high quality and low cost.
FIG. 1 shows an image processing process that requires GCP and DEM as a step of performing 3D modeling of a satellite image.
2 is a flowchart illustrating an overall image processing process proposed in the present invention.
FIG. 3 is a flow chart showing the sequence of image matching proposed in the present invention.
4 is a conceptual diagram of a patch and a search area in correlation matching.
5 is a conceptual diagram showing the position of the maximum correlation coefficient in the calculation of the peripheral pixel accuracy.
6 shows a hierarchical structure of an image.
FIG. 7 is a conceptual diagram of a stereoscopic image photographing for assisting an understanding of a stereoscopic image.
Hereinafter, the present invention will be described in more detail with reference to the drawings.
While the present invention has been described in connection with certain embodiments, it is obvious that the invention is not limited to the disclosed embodiments, but, on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the invention. It is to be understood, however, that the invention is not intended to be limited to the particular forms disclosed, but on the contrary, is intended to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention.
In the drawings, the same reference numerals are used for the same reference numerals, and in particular, the numerals of the tens and the digits of the digits, the digits of the tens, the digits of the digits and the alphabets are the same, Members referred to by reference numerals can be identified as members corresponding to these standards.
In the drawings, the components are expressed by exaggeratingly larger (or thicker) or smaller (or thinner) in size or thickness in consideration of the convenience of understanding, etc. However, It should not be.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The singular expressions include plural expressions unless the context clearly dictates otherwise. In the present application, the term " comprising " or " consisting of ", or the like, refers to the presence of a feature, a number, a step, an operation, an element, a component, But do not preclude the presence or addition of one or more other features, integers, steps, operations, components, parts, or combinations thereof.
Unless defined otherwise, all terms used herein, including technical or scientific terms, have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Terms such as those defined in commonly used dictionaries are to be interpreted as having a meaning consistent with the contextual meaning of the related art and are to be interpreted as either ideal or overly formal in the sense of the present application Do not.
Image Matching in (d) is a step of finding a matching point in a stereoscopic image (i) applying an image-based correlation matching technique; (Ii) applying a matching point approximation equation of surrounding pixel accuracy to improve the precision of the generated matching point; (Iii) applying a hierarchical structure in order to solve a problem that occurs when a correlation matching technique of (i) is applied to a region lacking a characteristic change, a linear boundary region, or a matching pattern error occurrence region; (Iv) correcting the matching error using the parallax information each time the matching of each layer is completed, in order to prevent a matching error of the upper layer from propagating to the lower layer. The accuracy of the satellite image can be improved by extracting the DEM data with improved accuracy by performing the 3D positioning using the extracted DEM data.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS The present invention will be described in detail with reference to the accompanying drawings.
2 is a flowchart of three-dimensional modeling using an image matching technique according to the present invention. Referring to FIG. 2, first, a satellite image and RPC data are acquired. In this case, when acquiring a satellite image, a stereoscopic image must be acquired so that stereoscopic vision can be performed. Stereo images are satellite images of the same area or object taken at different positions, and the height value can be calculated by using the satellite images.
Using the acquired satellite image and RPC data, geometric modeling of the image is performed using the RFM technique. RFM is a mathematical sensor modeling defined as one of four real-time geometry models with universal real-time image geometry model as an image conversion standard in OGC. RFM is a polynomial proportional model in which the denominator and the numerator are polynomials of higher order functions. The left and right satellite images acquired using this model are processed.
By using the satellite image and RPC data acquired from the above, the left and right satellite images obtained by geometric modeling by RFM processing method are obtained.
The image matching of the modeled satellite image is performed as shown in FIG.
(I) The correlation matching method is a step of applying the real-time correlation matching method. The correlation matching method calculates the correlation coefficient with the reference point of the left image at all positions included in the search area calculated as the reference point, To find the position that it has. As shown in FIG. 4, a patch f of pxp size is set around a pixel of the reference image (left image). For all the points in the search region, a patch g of the same size as the reference image patch is set to calculate the correlation coefficient r, and the point having the greatest correlation coefficient value is selected as the matching point. The correlation coefficient r is calculated by the following equation.
(Ii) a step of applying a matching point approximation formula of the neighboring pixel accuracy to improve the accuracy of the generated matching point, and as described in the previous step, by selecting a point having the largest correlation coefficient value in the search area, The matching point can be calculated. Here, by fitting the correlation coefficient to the surrounding pixels of the matching point having the maximum correlation coefficient as a quadratic function, it is possible to calculate the matching point of the peripheral pixel unit more precisely.
(Iii) a step of applying a hierarchical structure to solve a matching error in a region with insufficient characteristic change, a linear boundary region, and a repeated pattern occurrence region, which occurs when the correlation matching technique of (i) is applied, For example, in the case of an image structure having four layers as shown in FIG. 6, the lowest image is a raw image, and the width and width of the image are reduced to 1/2 every time one layer is ascended. In stereo matching of hierarchical structure, matching is performed from the upper layer and the result is applied to the lower layer. The position of the matching point in the lower layer is predicted using the upper layer matching point. When the image reduction ratio between the layers is 1: 2 as shown in FIG. 6, when the matching point (x_m, y_m) of the right image with respect to the pixel (x, y) (2p_x, 2p_y) for the reference image points (2x, 2y), (2x + 1,2y), (2x, 2y + 1) and do. The equation is expressed as follows.
(Iv) correcting the matching error using the parallax information each time the matching of each layer is completed, in order to prevent a matching error of the upper layer from propagating to the lower layer. Way
As described above, by using the RFM method using the RPC data, which is easy to perform the geometric modeling and has a certain accuracy, it is possible to perform geometric correction of the satellite image while reducing time and cost for modeling, We use a step - by - step image matching technique to improve the accuracy of DEM data acquisition for modeling.
Although the techniques that have been proposed in the basic image processing proposed in the present invention have been proposed through many studies, the present invention proposes additional techniques that can effectively use the techniques and improve the accuracy in the image matching step, And suggests ways to improve accuracy.
The background of the present invention can be said to be the presentation of an optimal solution for each step for three-dimensional positioning.
While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it is to be understood that the invention is not limited to the disclosed exemplary embodiments. Should be interpreted as belonging to the scope.
Claims (2)
(b) Image modeling by RFM method using RPC data for image modeling
(c) extracting the same point between stereoscopic images
(d) Image Matching of two images based on the same point extraction result
(e) Extracting a DEM (Digital Elevation Model) from the Stereo image subjected to Image Matching
(f) positioning the 3D image using three-dimensional image processing using the modeled satellite image and DEM data.
Image Matching in (d) is a step of finding a matching point in a stereoscopic image,
(I) Applying the real-time correlation-matching technique
(Ii) applying a matching point approximation equation of surrounding pixel accuracy to improve the precision of the generated matching point;
(Iii) Applying a hierarchical structure to solve the problem of matching error in (a), (b), (c), (c)
(Iv) correcting the matching error using the parallax information each time the matching of each layer is completed, in order to prevent a matching error of the upper layer from propagating to the lower layer. Way
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020130150460A KR20150065302A (en) | 2013-12-05 | 2013-12-05 | Method deciding 3-dimensional position of landsat imagery by Image Matching |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020130150460A KR20150065302A (en) | 2013-12-05 | 2013-12-05 | Method deciding 3-dimensional position of landsat imagery by Image Matching |
Publications (1)
Publication Number | Publication Date |
---|---|
KR20150065302A true KR20150065302A (en) | 2015-06-15 |
Family
ID=53504244
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020130150460A KR20150065302A (en) | 2013-12-05 | 2013-12-05 | Method deciding 3-dimensional position of landsat imagery by Image Matching |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR20150065302A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101668006B1 (en) * | 2015-09-08 | 2016-10-20 | 한국항공우주연구원 | Satellite Based Method and System for Constructing 3D GIS Data |
KR101692278B1 (en) | 2016-01-18 | 2017-01-03 | 최쌍임 | Jig for manufacturing glass bottles |
CN108846436A (en) * | 2018-06-13 | 2018-11-20 | 武汉朗视软件有限公司 | A kind of more view stereoscopic matching process |
KR20190026452A (en) * | 2017-09-05 | 2019-03-13 | 순천대학교 산학협력단 | A method of automatic geometric correction of digital elevation model made from satellite images and provided rpc |
CN112184546A (en) * | 2020-06-10 | 2021-01-05 | 中国人民解放军32023部队 | Satellite remote sensing image data processing method |
CN117889831A (en) * | 2024-03-13 | 2024-04-16 | 成都本原星通科技有限公司 | Terminal positioning method based on low-orbit satellite image matching |
-
2013
- 2013-12-05 KR KR1020130150460A patent/KR20150065302A/en not_active Application Discontinuation
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101668006B1 (en) * | 2015-09-08 | 2016-10-20 | 한국항공우주연구원 | Satellite Based Method and System for Constructing 3D GIS Data |
KR101692278B1 (en) | 2016-01-18 | 2017-01-03 | 최쌍임 | Jig for manufacturing glass bottles |
KR20190026452A (en) * | 2017-09-05 | 2019-03-13 | 순천대학교 산학협력단 | A method of automatic geometric correction of digital elevation model made from satellite images and provided rpc |
CN108846436A (en) * | 2018-06-13 | 2018-11-20 | 武汉朗视软件有限公司 | A kind of more view stereoscopic matching process |
CN112184546A (en) * | 2020-06-10 | 2021-01-05 | 中国人民解放军32023部队 | Satellite remote sensing image data processing method |
CN112184546B (en) * | 2020-06-10 | 2024-03-15 | 中国人民解放军32023部队 | Satellite remote sensing image data processing method |
CN117889831A (en) * | 2024-03-13 | 2024-04-16 | 成都本原星通科技有限公司 | Terminal positioning method based on low-orbit satellite image matching |
CN117889831B (en) * | 2024-03-13 | 2024-06-04 | 成都本原星通科技有限公司 | Terminal positioning method based on low-orbit satellite image matching |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111062905B (en) | Infrared and visible light fusion method based on saliency map enhancement | |
CN111080724B (en) | Fusion method of infrared light and visible light | |
Wu et al. | Integration of aerial oblique imagery and terrestrial imagery for optimized 3D modeling in urban areas | |
US11830222B2 (en) | Bi-level optimization-based infrared and visible light fusion method | |
CN101604018B (en) | Method and system for processing high-definition remote sensing image data | |
KR20150065302A (en) | Method deciding 3-dimensional position of landsat imagery by Image Matching | |
CN104935909B (en) | Multi-image super-resolution method based on depth information | |
CN105389774A (en) | Method and device for aligning images | |
CN110211169B (en) | Reconstruction method of narrow baseline parallax based on multi-scale super-pixel and phase correlation | |
Peña-Villasenín et al. | 3-D modeling of historic façades using SFM photogrammetry metric documentation of different building types of a historic center | |
CN110462685B (en) | Three-dimensional model reconstruction method and system | |
CN106023147B (en) | The method and device of DSM in a kind of rapidly extracting linear array remote sensing image based on GPU | |
CN105005964A (en) | Video sequence image based method for rapidly generating panorama of geographic scene | |
CN110738731A (en) | 3D reconstruction method and system for binocular vision | |
CN111582022A (en) | Fusion method and system of mobile video and geographic scene and electronic equipment | |
Oh | Novel approach to epipolar resampling of HRSI and satellite stereo imagery-based georeferencing of aerial images | |
CN103679680A (en) | Stereo matching method and system | |
Tatar et al. | Stereo rectification of pushbroom satellite images by robustly estimating the fundamental matrix | |
Agisoft | Metashape python reference | |
Parmehr et al. | Automatic parameter selection for intensity-based registration of imagery to LiDAR data | |
CN117291808B (en) | Light field image super-resolution processing method based on stream prior and polar bias compensation | |
CN108615221B (en) | Light field angle super-resolution method and device based on shearing two-dimensional polar line plan | |
CN112002007B (en) | Model acquisition method and device based on air-ground image, equipment and storage medium | |
CN103914807A (en) | Non-locality image super-resolution method and system for zoom scale compensation | |
Zhang et al. | A CNN-based subpixel level DSM generation approach via single image super-resolution |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WITN | Withdrawal due to no request for examination |