KR100473959B1 - Method for Amending Digital Map Using Rational Polynomial Coefficient Data of High Accurate - Google Patents
Method for Amending Digital Map Using Rational Polynomial Coefficient Data of High Accurate Download PDFInfo
- Publication number
- KR100473959B1 KR100473959B1 KR1020040103101A KR20040103101A KR100473959B1 KR 100473959 B1 KR100473959 B1 KR 100473959B1 KR 1020040103101 A KR1020040103101 A KR 1020040103101A KR 20040103101 A KR20040103101 A KR 20040103101A KR 100473959 B1 KR100473959 B1 KR 100473959B1
- Authority
- KR
- South Korea
- Prior art keywords
- data
- image
- digital
- rpc
- reference point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Classifications
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
- G09B29/00—Maps; Plans; Charts; Diagrams, e.g. route diagram
- G09B29/003—Maps
- G09B29/004—Map manufacture or repair; Tear or ink or water resistant maps; Long-life maps
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Mathematical Physics (AREA)
- Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Educational Technology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Processing Or Creating Images (AREA)
- Image Processing (AREA)
Abstract
본 발명은 고정확도의 RPC(rational polynomial coefficient) 데이터를 이용하여 수치지도를 수정제작하는 방법에 관한 것으로, 더욱 상세하게는, 고해상도 위성영상 및 단일 기준점을 이용하여 지형정보 추출시 필요한 기준점의 수를 대폭 줄이고 RPC 데이터의 정확도를 향상시켜, 수치지도의 수정제작시 소요되는 경비와 시간을 최소화할 수 있는 수치지도의 수정제작 방법에 관한 것이다. The present invention relates to a method of modifying a digital map using high accuracy rational polynomial coefficient (RPC) data. More specifically, the number of reference points required to extract topographic information using a high resolution satellite image and a single reference point is described. The present invention relates to a method of modifying a digital map that can greatly reduce the cost and time required to modify the digital map by greatly reducing and improving the accuracy of the RPC data.
본 발명에 따르면, 고해상도 위성영상과 단일 기준점만을 이용하여 기 구축되어진 수치지도의 수정제작이 가능하도록 함으로써, 수치지도 수정제작 작업 공정을 단순화하고, 작업시간을 단축시켜, 수치지도 수정제작의 효율성 및 채산성(採算性)을 향상시킬 수 있다.According to the present invention, by using a high-resolution satellite image and a single reference point to be able to modify the pre-built digital map, it is possible to simplify the process of creating a digital map correction, shorten the work time, efficiency and efficiency of digital map correction Profitability can be improved.
Description
본 발명은 고정확도의 RPC(rational polynomial coefficient) 데이터를 이용하여 수치지도를 수정제작하는 방법에 관한 것으로, 더욱 상세하게는, 고해상도 위성영상 및 단일 기준점을 이용하여 지형정보 추출시 필요한 기준점의 수를 대폭 줄이고 RPC 데이터의 정확도를 향상시켜, 수치지도의 수정제작시 소요되는 경비와 시간을 최소화할 수 있는 수치지도의 수정제작 방법에 관한 것이다. The present invention relates to a method of modifying a digital map using high accuracy rational polynomial coefficient (RPC) data. More specifically, the number of reference points required to extract topographic information using a high resolution satellite image and a single reference point is described. The present invention relates to a method of modifying a digital map that can greatly reduce the cost and time required to modify the digital map by greatly reducing and improving the accuracy of the RPC data.
수치지도를 제작하는 방법으로는 지상에서 측량기기를 이용하여 제작하는 방법, 항공 또는 보다 높은 위치에서 비행기 또는 위성을 이용하여 촬영한 사진 또는 영상과 지상측량 자료를 결합하여 제작하는 방법, 또는 기본 지형도를 수정하여 제작하는 방법이 있다. 이중 항공사진 또는 위성영상을 이용하여 제작하는 방법은 지상에서 측량장비를 이용하여 제작하는 방법에 비하여 경제적이므로 광대한 지역 지도 제작에 일반적으로 사용된다. 항공사진을 이용한 수치지도 제작은, 항공사진 촬영, 사진제작, 평면기준점측량, 표석매설측량을 통해 촬영된 항공사진과, 항공사진에서 판독이 분명하고 항공사진과 지상의 물체를 연결할 지상기준점들을 취득한 후, 사진기준점측량, 해석도화, 도화원도 출력을 통해 수치도화자료 제작 및 출력을 한다. 수치도화자료의 정확도를 검사하기 위해 출력된 도화원도를 가지고 현지조사를 실시하여 도로 및 시설들의 폭과 위치를 파악하고 사무실로 돌아와 정위치 편집 단계를 거쳐 도화과정에서 발생한 오류를 수정하고 도면제작편집 과정을 통해 지도도식에 맞추어 수정편집을 수행하여 수치지도 제작을 완료하게 된다. Digital maps can be created using survey instruments on the ground, combined with ground survey data and photographs or images taken using airplanes or satellites from aerial or higher locations, or basic topographic maps. There is a way to modify the production. The production of aerial photographs or satellite imagery is generally economical compared to the production of surveying equipment on the ground, and is therefore commonly used for making large area maps. Numerical map production using aerial photographs acquires aerial photographs taken through aerial photographing, photographing, plane reference point surveying, and surface laying surveying, and ground reference points that can be clearly read from aerial photographs and connect aerial photographs with ground objects. After that, digital drawing data is produced and output through photographic reference surveying, analysis drawing, and drawing drawing. In order to check the accuracy of the digital drawing data, the field research is conducted with the printed drawing diagram to identify the width and location of the roads and facilities, and return to the office to correct the errors during the drawing process. Through the process, correction and editing are performed according to the map diagram to complete the digital map production.
수치지도의 수정제작은 상기와 같은 과정을 통해 제작되어진 수치지도가 일정 기간이 지나 토지이용의 변화로 지형 및 지물이 변화하였을 때 변화된 지역을 파악 및 수정제작하여 지도의 최신성을 유지하는 작업을 말한다. 수치지도 수정제작 작업이 수행되는 부분은 수치지도 전체를 볼때 극히 일부이나, 일부의 수정을 위해서는 상기의 수치지도 제작과 동일한 단계를 수행해야 하므로 이로 인해 사업 목적에 비해 작업량의 부담이 많은 단점이 있다.Modification of digital map is to maintain the latest map by identifying and revising the changed area when the digital map produced through the above process changes terrain and features due to land use change after a certain period of time. Say. Part of the digital map correction production work is very partial when viewing the whole digital map, but it is necessary to perform the same steps as the above-mentioned digital map production for the partial correction, which is a burden of the work load compared to the business purpose. .
한편, 지구관측 위성은 크게, AVHRR이나 Landsat TM과 같이 과학적인 목적으로 발사된 위성과, SPOT, IRS, IKONOS, QUICKBIRD 등과 같이 상업적인 목적을 위한 인공위성으로 구분할 수 있다. 또한, 인공위성이 지구를 관측하는데 사용하는 파장대에 따라서 Landsat TM이나 SPOT, IRS, IKONOS와 같은 광학 영상과, RADARSAT이나 ERS-1, Seasat과 같은 레이다(Radar) 영상으로 구분할 수 있다.Earth observation satellites can be broadly classified into satellites launched for scientific purposes such as AVHRR and Landsat ™, and satellites for commercial purposes such as SPOT, IRS, IKONOS, and QUICKBIRD. In addition, according to the wavelength band that the satellite uses to observe the earth, it can be classified into optical images such as LandsatTM, SPOT, IRS, IKONOS, and radar images such as RADARSAT, ERS-1, and Seasat.
이러한 인공위성들은 일반적으로, 인공위성을 발사하는 목적에 따라서 그리고, 인공위성의 주된 관측의 대상이 무엇인지에 따라서 서로 다른 공간해상력, 관측주기, 분광해상력 및 파장대를 가지게 되는데, 이러한 이유로 각각의 인공위성에서 얻어지는 영상들은 영상마다 독특한 특징을 갖는다. These satellites generally have different spatial resolutions, observation periods, spectral resolutions, and wavelength bands, depending on the purpose of launching the satellite and what the primary object of the satellite is being observed. They have unique features for each image.
최근, 고해상도 위성영상 기술의 발달은 1m급 해상도의 영상정보를 상업적 목적으로 제공할 수 있도록 하여, 고해상도 위성영상의 항공사진 측량 및 수치지도 수정제작 분야에서의 활용방안에 대한 관심이 매우 높아지고 있다. 고해상도 위성영상은 항공사진 촬영, 사진제작 비용과 비교하여 가격이 경제적이고, 단영상의 촬영면적이 1/5,000 축척의 사진(약 1㎢) 또는 1/20,000 축척의 사진(약 17 ㎢)에 비해 100㎢ 이상의 넓은 면적을 포함한다.Recently, the development of high-resolution satellite image technology has been able to provide image information of 1m resolution for commercial purposes, the interest in the application of high-resolution satellite imagery in the field of aerial surveying and digital map correction production has been very high. High resolution satellite images are more cost-effective than aerial photographing and photo production costs, and the area of the single image is 1 / 5,000 scale (about 1㎢) or 1 / 20,000 scale (about 17㎢). It covers a large area of over 100 square meters.
일반적으로, 대표적인 상업적 고해상도 위성영상인 'IKONOS' 또는 'QUICKBIRD' 영상의 경우, SPOT, KOMPSAT와 같은 위성영상에서 제공하는 위성의 항법자료 또는 센서정보를 보안을 이유로 공개하지 않는 대신, RPC(rational polynomial coefficients)라는 매개변수를 제공하고 있다. 영상의 최종 사용자는 RPC를 이용하여 영상에 나타난 대상의 위치정보를 획득할 수 있으며, 일련의 사진측량학적 처리를 할 수 있다. 그러나, 제공되는 RPC의 정확도가 고해상도 위성영상의 해상도에 비하여 현저히 낮은 정확도로 제공되고 있기 때문에 항공사진측량 공정 및 수치지도 제작공정에서 정확도가 확보되지 않아 실질적으로 사용하기 어려운 단점이 있다. In general, 'IKONOS' or 'QUICKBIRD' images, which are representative commercial high resolution satellite images, do not disclose satellite navigation data or sensor information provided by satellite images such as SPOT and KOMPSAT for security reasons. provide a parameter called coefficients. The end user of the image can acquire the location information of the object shown in the image using the RPC, and can perform a series of photogrammetric processing. However, since the accuracy of the provided RPC is provided at a significantly lower accuracy than the resolution of the high resolution satellite image, it is difficult to use practically because the accuracy is not secured in the aerial survey process and the digital map manufacturing process.
현재 판매되고 있는 고해상도 위성영상인 IKONOS 위성영상의 등급별 정확도를 보면 Geo Product는 15m, Refernece Product는 25m, Pro Product는 10m, Precision Product는 4m, PrecisionPlus Product는 2m의 정확도를 제공하지만, 상 이상의 경우 안보적인 이유와 관련되어 국가에 따라 구매가 제한되며 Precision 이상은 미국을 제외하고는 구매가 불가능하다. 따라서 Geo와 Reference를 기준으로 본다면 IKONOS 본래의 영상해상력인 1m에 비해 현저히 낮은 정확도로 RPC 데이터가 제공되고 있다. 낮은 정확도의 RPC 데이터는 수치지도 수정제작시 경비와 시간을 증가시켜 작업능률을 저하시키게 된다. 따라서, 수치지도 수정제작에 사용 가능하도록 RPC 데이터의 품질을 갱신해야할 필요성이 있다. The accuracy of IKONOS satellite imagery, which is currently on sale, is 15m in Geo Product, 25m in Refernece Product, 10m in Pro Product, 4m in Precision Product, and 2m in PrecisionPlus Product. Purchasing is restricted by country for certain reasons and Precision is not available for purchase except in the United States. Therefore, based on Geo and Reference, RPC data is provided with a significantly lower accuracy than 1m, which is the original image resolution of IKONOS. Low-accuracy RPC data can increase the cost and time of digital map revisions, reducing work efficiency. Therefore, there is a need to update the quality of RPC data so that it can be used for digital map correction production.
본 발명은 이와 같은 문제점을 해결하기 위해서 창출된 것으로, 본 발명의 목적은 수치지도를 수정제작하는데 있어 경비와 시간을 최소화하고 효율성을 높일 수 있는 고정확도의 RPC 데이터를 이용한 수치지도의 수정제작 방법을 제공하는데 있다. The present invention was created to solve the above problems, and an object of the present invention is to modify a digital map using high accuracy RPC data to minimize the cost and time and improve the efficiency of digital map production. To provide.
상기 목적을 달성하기 위하여, 본 발명은 RPC(rational polynomial coefficient) 데이터를 이용하여 수치지도를 수정제작하는 방법에 있어서, (a) 기 구축된 자료 또는 위성영상 신규촬영을 통해 고해상도 위성영상 및 RPC 데이터를 취득하는 단계; (b) 상기 고해상도 위성영상으로부터 영상을 대표하는 단일 기준점을 선점한 후 고정밀 3차원 위치 결정이 가능한 GPS를 이용하여 단일 기준점을 측량하고 상대측위 방식에 의해 단일 기준점 좌표를 취득하는 단계; (c) 상기 취득된 단일 기준점 좌표을 이용하여 정확도가 향상된 RPC 데이터를 갱신하는 단계; (d) 수치표고모형(DEM)을 입력한 후, 입력된 수치표고모형과 상기 (c) 단계의 갱신된 RPC 데이터를 이용하여 고해상도 위성영상 지형을 보정한 수치영상지도를 제작하는 단계; (e) 상기 제작된 수치영상지도와 기존의 수치지도를 중첩한 후, 비교하여 변화된 지역을 추출하고 묘사하는 단계; (f) 변화된 지역의 묘사 자료를 출력하여 현지조사를 실시하고, 변화된 지역을 측량하는 단계; 및 (g) 상기 현지조사 결과와 측량자료를 이용하여 기존 수치지도의 정위치를 편집하고 도면편집하는 단계를 포함하는 것을 특징으로 하고,In order to achieve the above object, the present invention provides a method of modifying the digital map using the rational polynomial coefficient (RPC) data, (a) high resolution satellite image and RPC data through the newly constructed data or satellite image Obtaining a; (b) preempting a single reference point representative of the image from the high-resolution satellite image, surveying a single reference point using GPS capable of high-precision three-dimensional positioning, and obtaining a single reference point coordinate by a relative positioning method; (c) updating the RPC data with improved accuracy using the acquired single reference point coordinates; (d) inputting a digital elevation model (DEM) to produce a digital image map correcting the high resolution satellite image topography using the input digital elevation model and the updated RPC data of step (c); (e) overlapping the produced digital image map with an existing digital map, and comparing and extracting and describing the changed region; (f) output field description data of the changed area to conduct a field survey and survey the changed area; And (g) editing and drawing-editing the exact position of the existing numerical map using the field survey results and the survey data.
상기 (c) 단계는 (i) 단일 기준점 좌표로부터 RPC 데이터 갱신모형을 유도하는 단계; (ii) 추가된 기준점 및 적용영상에 따라 편위수정모형을 선택하고, 선택된 편위수정모형과 상기 (i) 단계의 RPC 데이터 갱신모형을 이용하여 최소제곱조정을 통해 편위수정 매개변수를 계산하는 단계; (iii) 상기 편위수정 매개변수를 정규좌표로 변환하는 단계; 및 (iv) 상기 변환된 정규좌표를 이용하여 기존의 낮은 정확도를 가지고 있던 RPC 데이터를 고정확도의 RPC 데이터로 갱신하는 단계로 구성되는 것을 특징으로 하는 고정확도의 RPC 데이터를 이용한 수치지도의 수정제작 방법을 제공한다.Step (c) includes (i) deriving an RPC data update model from a single reference point coordinate; (ii) selecting the deviation correction model according to the added reference point and the applied image, and calculating the deviation correction parameter through the least square adjustment using the selected deviation correction model and the RPC data update model of step (i); (iii) converting the deviation correction parameter to normal coordinates; And (iv) updating the existing low accuracy RPC data with the high accuracy RPC data using the converted normal coordinates. Provide a method.
본 발명에 있어서, 상기 (d) 단계의 수치표고모형은 갱신된 RPC 데이터에 의해 수치표고모형의 각 그리드(grid) 점의 3차원 좌표에 해당하는 영상좌표를 결정하고, 상기 영상좌표에 해당하는 영상내의 수치값(digital value)을 보강하여 취득되는 것을 특징으로 할 수 있다.In the present invention, the numerical elevation model of step (d) determines the image coordinates corresponding to the three-dimensional coordinates of each grid point of the numerical elevation model by the updated RPC data, and corresponds to the image coordinates. The digital value may be obtained by reinforcing a digital value in the image.
이하 첨부된 도면을 참조하여 본 발명을 보다 상세히 설명한다. Hereinafter, the present invention will be described in more detail with reference to the accompanying drawings.
도 2는 본 발명에 따른 고정확도의 RPC 데이터를 이용한 수치지도의 수정제작 흐름도이다. 도 2를 참조하면, 먼저, 기 구축된 자료 또는 위성영상 신규촬영을 통해 원하는 고해상도 위성영상 및 RPC 데이터를 취득한다. 수치지도 갱신을 요구하는 지역의 고품질의 영상이 매건 준비되어 있는 것이 아니므로, 촬영지역의 위치 및 조건을 요청해 기 구축된 자료에서 구하거나 위성영상 신규촬영을 통해 원하는 고해상도 위성영상을 획득한다. 이때 RPC 데이터의 취득도 함께 이루어진다.2 is a flowchart of a modified map of a digital map using high accuracy RPC data according to the present invention. Referring to FIG. 2, first, desired high resolution satellite image and RPC data are acquired through pre-built data or new satellite image recording. Since high quality images of regions requiring digital map updates are not prepared, the locations and conditions of the photographing regions are requested, and the desired high resolution satellite images can be obtained from existing data or through new satellite images. At this time, the RPC data is also acquired.
상기 획득된 고해상도 위성영상으로부터 영상을 대표하는 단일 기준점을 선점한 후 고정밀 3차원 위치 결정이 가능한 GPS를 이용하여 단일 기준점을 측량하여 단일 기준점 좌표를 취득한다. 단일 기준점 좌표의 취득은 국토지리정보원에 의해 높은 정확도의 좌표가 관리되고 수평 및 표고 기준점들을 기지의 점으로 하여, 미지의 점인 위성영상에서 사용할 단일 기준점을 결정하는 상대측위(differential GPS) 방식을 사용하여 이루어진다. 상기와 같은 과정을 통하여 단일 기준점 좌표로 영상기준점 좌표와 지상기준점 좌표를 수득하게 된다. After preempting a single reference point representing an image from the obtained high-resolution satellite image, a single reference point coordinate is acquired by surveying a single reference point using GPS capable of high-precision three-dimensional positioning. Acquisition of single reference point coordinates uses a differential GPS method that determines the single reference point to use in satellite imagery, which is managed by the National Geographic Information Institute, with high accuracy coordinates and the horizontal and elevation reference points as known points. It is done by Through the above process, the image reference point coordinates and the ground reference point coordinates are obtained as a single reference point coordinate.
상기 단일 기준점을 이용하여 수득된 고해상도 영상기준점 좌표와 지상기준점 좌표를 이용하여 정확도가 향상된 RPC 데이터를 갱신한다. The RPC data having improved accuracy is updated by using the high resolution image reference point coordinates and the ground reference point coordinates obtained using the single reference point.
정확도가 향상된 RPC 데이터의 갱신은 Updating RPC data with improved accuracy
(i) 단일 기준점 좌표로부터 RPC 데이터 갱신모형을 유도하는 단계; (i) deriving an RPC data update model from a single reference point coordinate;
(ii) 추가된 기준점 및 적용영상에 따라 편위수정모형을 선택하고, 선택된 편위수정모형과 상기 (i) 단계의 RPC 데이터 갱신모형을 이용하여 최소제곱조정을 통해 편위수정 매개변수를 계산하는 단계; (ii) selecting the deviation correction model according to the added reference point and the applied image, and calculating the deviation correction parameter through the least square adjustment using the selected deviation correction model and the RPC data update model of step (i);
(iii) 상기 편위수정 매개변수를 정규좌표로 변환하는 단계; 및 (iii) converting the deviation correction parameter to normal coordinates; And
(iv) 상기 변환된 정규좌표를 이용하여 기존의 낮은 정확도를 가지고 있던 RPC 데이터를 고정확도의 RPC 데이터로 갱신하는 단계로 이루어진다.(iv) updating the existing low accuracy RPC data with the high accuracy RPC data using the converted normal coordinates.
도 4는 상기의 정확도가 향상된 RPC 데이터의 갱신과정을 나타낸 흐름도이다. 먼저, RPC 데이터 갱신모형(S421)은 상기 단일 기준점을 이용하여 수득된 영상기준점 좌표(r, c)(S411)와 지상기준점 좌표(ψ, λ, h)(S412)를 이용하여 하기 수학식 2를 통해 유도된다. RPC 모형은 하기 수학식 1에서 보는 바와 같이 지상좌표와 영상좌표 간의 변환관계를 나타내며, 현재 사용되고 있는 RPC 모형은 전방변환방식과 후방변환방식으로 대별되나 일반적으로 활용되는 형태가 전방변환방식이다 (수학식 1). 전방변환방식의 정확도가 후방변환방식의 정확도보다 우수하므로, 본 발명에서는 전방변환방식을 기본모형으로 채택한다.4 is a flowchart illustrating a process of updating RPC data with improved accuracy. First, the RPC data update model S421 uses the image reference point coordinates ( r , c ) (S411) and the ground reference point coordinates (ψ, λ, h ) (S412) obtained using the single reference point. Guided through. As shown in Equation 1 below, the RPC model represents a transformation relationship between the ground coordinates and the image coordinates, and the currently used RPC model is roughly classified into a forward transform method and a backward transform method, but a commonly used form is a forward transform method. Equation 1). Since the accuracy of the forward conversion method is superior to that of the backward conversion method, the front conversion method is adopted as the basic model in the present invention.
여기서, r은 영상의 행 좌표이며, c는 영상의 열 좌표이고, ψ는 WGS-84 타원체 기준의 위도, λ는 WGS-84 타원체 기준의 경도, h는 고도값을 의미한다.Where r is the row coordinate of the image, c is the column coordinate of the image, ψ is the latitude of the WGS-84 ellipsoid, λ is the longitude of the WGS-84 ellipsoid, and h is the altitude value.
상기 수학식 1을 통해 지상의 좌표에 해당하는 영상 내의 위치를 결정하게 되며, 반대로 입체영상에서의 공액점을 관측하게 되면 지상의 3차원 좌표를 최소제곱조정을 통해 결정할 수 있다. 본 발명에서 낮은 정확도로 제공되는 RPC 데이터를 갱신하기 위한 갱신모형은 하기 수학식 2를 통해 유도된다.Through Equation 1, the position in the image corresponding to the coordinates of the ground is determined. On the contrary, when observing a conjugate point in the stereoscopic image, the three-dimensional coordinates of the ground may be determined through the least squares adjustment. An update model for updating RPC data provided with low accuracy in the present invention is derived through Equation 2 below.
여기서, R i 와 C i 는 사전에 관측된 지상기준점의 행과 열좌표이며, △R i 와 △C i 는 편위수정모형이고, r i (ψ, λ, h)와 c i (ψ, λ, h)는 지상기준점의 지상좌표를 기준으로 제공되는 상기 RPC 모형에 의해 계산된 영상좌표이며, row s 와 column s 는 제공된 영상의 스케일 벡터이고, row o 와 column o 는 제공된 영상의 오프셋 벡터이며, 와 는 편위수정모형에 의해 발생 가능한 랜덤오차를 의미한다.Where R i and C i are the row and column coordinates of the ground reference point previously observed, △ R i and Δ C i are the deviation-corrected models, and r i (ψ, λ, h ) and c i (ψ, λ , h ) is the image coordinate calculated by the RPC model provided based on the ground coordinates of the ground reference point, row s and column s are the scale vectors of the provided image, and row o and column o are the offset vectors of the provided image. , Wow Denotes a random error that can be generated by the deviation correction model.
상기 편위수정모형은 다음의 4가지 경우와 같다.The deviation correction model is the same as the following four cases.
여기서, α 1 , α 2 , α 3 , β 1 , β 2 , β 3 는 편위수정을 위한 매개변수를 의미한다. 상기 4가지 경우의 편위수정모형은 각각 추가된 기준점 및 적용영상의 특성에 따라 다르게 사용할 수 있다.Here, α 1 , α 2 , α 3 , β 1 , β 2 , β 3 mean parameters for deviation correction. The four cases of deviation correction model can be used differently according to the characteristics of the added reference point and the applied image, respectively.
상기와 같은 RPC 데이터 갱신모형 및 선택된 편위수정모형을 통해 추가된 지상기준점으로부터 최종적으로 편위수정모형의 매개변수를 계산하게 되며, 계산과정은 최소제곱조정이론에 따른다. Finally, the parameters of the deviation correction model are calculated from the ground reference point added through the RPC data update model and the selected deviation correction model, and the calculation process follows the least squares adjustment theory.
최소제곱 조정을 통해 편위수정모형의 매개변수가 계산된 이후, 편위수정모형의 매개변수를 하기 수학식 3을 이용하여 정규좌표로 변환하며, 최종적으로 하기 수학식 4를 통해 RPC 데이터를 갱신한다. 모든 RPC 모형의 분자항은 보정된 편위수정모형에 의해 갱신이 이루어진다.After the parameters of the deviation correction model are calculated through the least squares adjustment, the parameters of the deviation correction model are converted into regular coordinates using Equation 3 below, and finally, the RPC data is updated through Equation 4 below. The molecular terms of all RPC models are updated by the corrected deviation correction model.
상기와 같이 정확도가 향상된 RPC 데이터를 갱신한 다음, 수치표고모형과 갱신된 RPC 데이터를 이용하여 수치영상지도를 제작한다. 일반적으로 수치영상지도는 해당영상에 대응하는 수치표고모형을 이용하여 제작되며, 수치표고모형은 수치사진측량 S/W와 기존 수치지도의 표고관련 자료 및 외부 자료를 이용하여 제작된다. 본 발명에서 수치표고모형은 수치표고모형의 각 그리드(grid) 점의 3차원 좌표에 해당하는 영상좌표를 상기 갱신된 RPC 데이터에 의해 결정하고, 이에 해당하는 영상내의 수치값(digital value)을 보강하여 취득된다. 상기의 방식으로 수치표고모형의 모든 그리드 점에 해당하는 수치값을 계산하고, 이를 취합하여 고해상도 위성영상에서 존재하는 편위를 수정한 수치영상지도를 제작한다.After updating the RPC data with improved accuracy as described above, a digital image map is produced using the digital elevation model and the updated RPC data. In general, a digital image map is produced by using a digital elevation model corresponding to the corresponding image, and a digital elevation model is produced by using digital photogrammetry software and elevation data and external data of existing digital maps. In the present invention, the numerical elevation model determines the image coordinates corresponding to the three-dimensional coordinates of each grid point of the numerical elevation model by the updated RPC data, and reinforces the digital values in the corresponding images. Is obtained. In this way, the numerical values corresponding to all grid points of the digital elevation model are calculated, and the numerical values are combined to produce a digital image map that corrects the deviations in the high resolution satellite image.
상기와 같이 제작된 수치영상지도를 이용하여 수치지도의 수정제작 작업을 수행한다. 상기의 수치영상지도와 기존의 수치지도를 중첩한 다음 비교하여 변화된 지형지물을 파악하고 파악된 지형지물에 대한 묘사를 행한다. 지형지물에 대한 묘사가 끝나면, 묘사된 자료를 출력하여 직접 현지조사 및 측량을 통해 묘사의 정확성을 재고하고, 위성영상 촬영시기 이후 변화된 자료를 수집하여 자료의 최신성을 유지하게 된다. 상기 현지조사 및 측량성과는 실내에서의 작업을 통해 변화된 지물 갱신 및 정위치 편집과정과 최종적으로 도면편집하는 과정을 거치게 되며, 이로써 본 발명에 따른 수치지도 수정제작 작업이 완료하게 된다.Modification and production of the digital map is performed using the digital image map produced as described above. The digital image map and the existing digital map are superimposed and then compared to identify the changed feature and a description of the identified feature is performed. When the description of the feature is finished, the depicted data is output and the accuracy of the description is reconsidered through direct field surveys and surveys, and the data that have changed since the time of satellite imaging was collected to maintain the latest data. The field survey and surveying results undergo a process of updating a feature and editing an in-situ location and finally editing a drawing through indoor work, thereby completing the digital map correction manufacturing work according to the present invention.
이상으로 본 발명 내용의 특정한 부분을 상세히 기술하였는 바, 당업계의 통상의 지식을 가진 자에게 있어서, 이러한 구체적 기술은 단지 바람직한 실시양태일 뿐이며, 이에 의해 본 발명의 범위가 제한되는 것은 아니다. 본 발명의 실질적인 범위는 첨부된 청구항들과 그것들의 등가물에 의하여 정의된다고 할 것이다.As described above in detail specific parts of the present invention, those skilled in the art, these specific descriptions are only preferred embodiments, and the scope of the present invention is not limited thereby. It is intended that the scope of the invention be defined by the claims appended hereto and their equivalents.
이상에서 상세히 설명한 바와 같이, 본 발명에 따르면 수치지도를 수정제작하는데 있어, 고해상도 위성영상 및 단일 기준점만을 이용하여 기 구축되어진 수치지도의 수정제작이 가능하도록 함으로써, 수치지도의 수정제작시 소요되는 경비와 시간을 최소화하고 효율성을 높일 수 있으며, 또한, 수치지도 수정제작 작업 공정을 단순화하여 수치지도 수정제작의 채산성(採算性)을 향상시킬 수 있다.As described in detail above, according to the present invention, in the manufacture of digital maps, by using the high-resolution satellite image and a single reference point, it is possible to modify the digital map that is already built, the cost required for the production of digital maps In addition, it is possible to minimize the time, increase the efficiency, and to improve the profitability of digital map correction by simplifying the process of digital map correction.
도 1은 기존의 항공사진 측량기술을 이용한 수치지도 제작 및 수정제작의 흐름도를 나타낸 것이다.Figure 1 shows a flow chart of the digital map production and correction using the conventional aerial surveying technology.
도 2는 본 발명의 일실시예에 따른 수치지도의 수정제작 흐름도를 나타낸 것이다.2 is a flowchart of a modified production of a digital map according to an embodiment of the present invention.
도 3은 본 발명의 일실시예에 따른 단일 기준점을 이용한 RPC 데이터의 정확도 향상 및 수치지도 수정제작을 위한 시스템 구성도를 나타낸 것이다.3 is a diagram illustrating a system configuration for improving accuracy of a RPC data using a single reference point and manufacturing a digital map correction according to an embodiment of the present invention.
도 4는 본 발명의 일실시예에 따른 단일 기준점을 이용한 RPC 데이터의 정확도 향상을 위한 처리 흐름도를 나타낸 것이다.4 is a flowchart illustrating a process for improving accuracy of RPC data using a single reference point according to an embodiment of the present invention.
Claims (2)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020040103101A KR100473959B1 (en) | 2004-12-08 | 2004-12-08 | Method for Amending Digital Map Using Rational Polynomial Coefficient Data of High Accurate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020040103101A KR100473959B1 (en) | 2004-12-08 | 2004-12-08 | Method for Amending Digital Map Using Rational Polynomial Coefficient Data of High Accurate |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| KR100473959B1 true KR100473959B1 (en) | 2005-03-14 |
Family
ID=37228385
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| KR1020040103101A Expired - Lifetime KR100473959B1 (en) | 2004-12-08 | 2004-12-08 | Method for Amending Digital Map Using Rational Polynomial Coefficient Data of High Accurate |
Country Status (1)
| Country | Link |
|---|---|
| KR (1) | KR100473959B1 (en) |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100572119B1 (en) | 2005-10-26 | 2006-04-18 | 한진정보통신(주) | Method of detecting feature changes using digital map and aerial laser survey data |
| KR100797518B1 (en) | 2007-06-13 | 2008-02-01 | 대한민국 | Forest image reading and numerical clinical drawing system and its method |
| KR100850043B1 (en) | 2008-01-18 | 2008-08-05 | 대한민국 | Numerical clinical drawing production management system |
| KR101102415B1 (en) * | 2011-05-04 | 2012-01-05 | 주식회사 우리강산시스템 | Editing system and editing method for security area in aviation / satellite image |
| KR101166505B1 (en) | 2012-01-26 | 2012-07-20 | 주식회사 우리강산시스템 | Digital map data upgrade system |
| KR101166508B1 (en) | 2012-01-26 | 2012-07-20 | 주식회사 우리강산시스템 | Digital map update system |
| KR101166510B1 (en) | 2012-01-26 | 2012-07-20 | 주식회사 우리강산시스템 | Image drawing management system |
| KR101200523B1 (en) | 2012-05-10 | 2012-11-13 | 주식회사 첨단공간정보 | Digital Map Update System |
| KR101213608B1 (en) | 2012-05-31 | 2012-12-18 | 주식회사 첨단공간정보 | Digital map accuracy improvement system |
| KR101223126B1 (en) | 2012-10-10 | 2013-01-17 | 주식회사 고원항공정보 | Precision digital map making system |
| KR101224046B1 (en) * | 2010-08-13 | 2013-01-24 | 대한민국 | System for reading forest image and fabricating detailed digital forest map using data back-up type |
| KR101295747B1 (en) * | 2012-04-25 | 2013-08-20 | 서울시립대학교 산학협력단 | System of automatic geometric correction using rational polynomial cofficient and method thereof |
| KR101313871B1 (en) | 2013-07-19 | 2013-10-01 | (주)올포랜드 | Surveying system for seawater road by sound exploration |
| KR101954708B1 (en) * | 2018-10-12 | 2019-03-06 | 새한항업(주) | System editting spatial imagery by revising draw images of image per altitude |
| 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 |
| KR20190053404A (en) * | 2017-11-10 | 2019-05-20 | 한국항공우주연구원 | Method and apparatus for constructing high precision 3D GIS data using lattice type GCPs |
| KR102058626B1 (en) | 2017-11-29 | 2019-12-23 | (주)엔젤스윙 | Method and Device for visualizing the change in topography |
| KR102072229B1 (en) * | 2019-10-15 | 2020-01-31 | 주식회사 동신지티아이 | System For Synthesizing Drawing Images |
| US11657437B2 (en) | 2017-11-29 | 2023-05-23 | Angelswing Inc | Method and apparatus for providing drone data by matching user with provider |
| CN116625354A (en) * | 2023-07-21 | 2023-08-22 | 山东省国土测绘院 | High-precision topographic map generation method and system based on multi-source mapping data |
| WO2024066586A1 (en) * | 2022-09-27 | 2024-04-04 | 杭州海康威视系统技术有限公司 | Terrain data updating method and apparatus, device, and storage medium |
-
2004
- 2004-12-08 KR KR1020040103101A patent/KR100473959B1/en not_active Expired - Lifetime
Cited By (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100572119B1 (en) | 2005-10-26 | 2006-04-18 | 한진정보통신(주) | Method of detecting feature changes using digital map and aerial laser survey data |
| KR100797518B1 (en) | 2007-06-13 | 2008-02-01 | 대한민국 | Forest image reading and numerical clinical drawing system and its method |
| KR100850043B1 (en) | 2008-01-18 | 2008-08-05 | 대한민국 | Numerical clinical drawing production management system |
| KR101224046B1 (en) * | 2010-08-13 | 2013-01-24 | 대한민국 | System for reading forest image and fabricating detailed digital forest map using data back-up type |
| KR101102415B1 (en) * | 2011-05-04 | 2012-01-05 | 주식회사 우리강산시스템 | Editing system and editing method for security area in aviation / satellite image |
| KR101166510B1 (en) | 2012-01-26 | 2012-07-20 | 주식회사 우리강산시스템 | Image drawing management system |
| KR101166508B1 (en) | 2012-01-26 | 2012-07-20 | 주식회사 우리강산시스템 | Digital map update system |
| KR101166505B1 (en) | 2012-01-26 | 2012-07-20 | 주식회사 우리강산시스템 | Digital map data upgrade system |
| US9378585B2 (en) | 2012-04-25 | 2016-06-28 | University Of Seoul Industry Cooperation Foundation | System and method for automatic geometric correction using RPC |
| KR101295747B1 (en) * | 2012-04-25 | 2013-08-20 | 서울시립대학교 산학협력단 | System of automatic geometric correction using rational polynomial cofficient and method thereof |
| WO2013162207A1 (en) * | 2012-04-25 | 2013-10-31 | 서울시립대학교 산학협력단 | System and method for automatic geometric correction using rpc |
| KR101200523B1 (en) | 2012-05-10 | 2012-11-13 | 주식회사 첨단공간정보 | Digital Map Update System |
| KR101213608B1 (en) | 2012-05-31 | 2012-12-18 | 주식회사 첨단공간정보 | Digital map accuracy improvement system |
| KR101223126B1 (en) | 2012-10-10 | 2013-01-17 | 주식회사 고원항공정보 | Precision digital map making system |
| KR101313871B1 (en) | 2013-07-19 | 2013-10-01 | (주)올포랜드 | Surveying system for seawater road by sound exploration |
| 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 |
| KR101965965B1 (en) | 2017-09-05 | 2019-08-13 | 순천대학교 산학협력단 | A method of automatic geometric correction of digital elevation model made from satellite images and provided rpc |
| KR20190053404A (en) * | 2017-11-10 | 2019-05-20 | 한국항공우주연구원 | Method and apparatus for constructing high precision 3D GIS data using lattice type GCPs |
| KR102013625B1 (en) | 2017-11-10 | 2019-10-21 | 한국항공우주연구원 | Method and apparatus for constructing high precision 3D GIS data using lattice type GCPs |
| KR102058626B1 (en) | 2017-11-29 | 2019-12-23 | (주)엔젤스윙 | Method and Device for visualizing the change in topography |
| US11657437B2 (en) | 2017-11-29 | 2023-05-23 | Angelswing Inc | Method and apparatus for providing drone data by matching user with provider |
| KR101954708B1 (en) * | 2018-10-12 | 2019-03-06 | 새한항업(주) | System editting spatial imagery by revising draw images of image per altitude |
| KR102072229B1 (en) * | 2019-10-15 | 2020-01-31 | 주식회사 동신지티아이 | System For Synthesizing Drawing Images |
| WO2024066586A1 (en) * | 2022-09-27 | 2024-04-04 | 杭州海康威视系统技术有限公司 | Terrain data updating method and apparatus, device, and storage medium |
| CN116625354A (en) * | 2023-07-21 | 2023-08-22 | 山东省国土测绘院 | High-precision topographic map generation method and system based on multi-source mapping data |
| CN116625354B (en) * | 2023-07-21 | 2023-09-29 | 山东省国土测绘院 | High-precision topographic map generation method and system based on multi-source mapping data |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR100473959B1 (en) | Method for Amending Digital Map Using Rational Polynomial Coefficient Data of High Accurate | |
| Hu et al. | Understanding the rational function model: methods and applications | |
| KR101965965B1 (en) | A method of automatic geometric correction of digital elevation model made from satellite images and provided rpc | |
| Granshaw | Photogrammetric terminology | |
| Tao et al. | Photogrammetric exploitation of IKONOS imagery for mapping applications | |
| CN102735216B (en) | CCD stereoscopic camera three-line imagery data adjustment processing method | |
| Teo et al. | DEM-aided block adjustment for satellite images with weak convergence geometry | |
| Toutin et al. | 3D models for high resolution images: examples with QuickBird, IKONOS and EROS | |
| JP2008524684A (en) | How to process images using automatic georeferencing of images obtained from pairs of images acquired in the same focal plane | |
| KR20040055510A (en) | Ikonos imagery rpc data update method using additional gcp | |
| CN110986888A (en) | Aerial photography integrated method | |
| Zhang et al. | Block adjustment for satellite imagery based on the strip constraint | |
| Gonçalves | Automatic orientation and mosaicking of archived aerial photography using structure from motion | |
| WO2024207677A1 (en) | Method and system for rapid updating of topographic map of mining area | |
| Heisig et al. | Re-engineering the past: countrywide geo-referencing of archival aerial imagery | |
| Jain | How photogrammetric software works: a perspective based on UAV’s exterior orientation parameters | |
| Cogliati et al. | Extraction of dems and orthoimages from archive aerial imagery to support project planning in civil engineering | |
| Lee et al. | Orthorectification of KOMPSAT optical images using various ground reference data and accuracy assessment | |
| KR100473958B1 (en) | Method for Updating Rational Polynomial Coefficient Data of High Accurate Using GPS Single Control Point | |
| Kim et al. | Estimation and improvement in the geolocation accuracy of rational polynomial coefficients with minimum GCPs using KOMPSAT-3A | |
| Gonçalves et al. | Use of the SRTM DEM as a geo-referencing tool by elevation matching | |
| Fraser et al. | Precise georefrencing of long strips of ALOS imagery | |
| Dowman | Encoding and validating data from maps and images | |
| Sepehry Amin et al. | True orthophoto generation using google earth imagery and comparison to UAV orthophoto | |
| Smith | Topographic mapping |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A201 | Request for examination | ||
| A302 | Request for accelerated examination | ||
| PA0109 | Patent application |
Patent event code: PA01091R01D Comment text: Patent Application Patent event date: 20041208 |
|
| PA0201 | Request for examination | ||
| PA0302 | Request for accelerated examination |
Patent event date: 20041208 Patent event code: PA03022R01D Comment text: Request for Accelerated Examination |
|
| E701 | Decision to grant or registration of patent right | ||
| PE0701 | Decision of registration |
Patent event code: PE07011S01D Comment text: Decision to Grant Registration Patent event date: 20050219 |
|
| GRNT | Written decision to grant | ||
| PR0701 | Registration of establishment |
Comment text: Registration of Establishment Patent event date: 20050221 Patent event code: PR07011E01D |
|
| PR1002 | Payment of registration fee |
Payment date: 20050222 End annual number: 3 Start annual number: 1 |
|
| PG1601 | Publication of registration | ||
| PR1001 | Payment of annual fee |
Payment date: 20080219 Start annual number: 4 End annual number: 4 |
|
| PR1001 | Payment of annual fee |
Payment date: 20090121 Start annual number: 5 End annual number: 5 |
|
| PR1001 | Payment of annual fee |
Payment date: 20100218 Start annual number: 6 End annual number: 6 |
|
| PR1001 | Payment of annual fee |
Payment date: 20110218 Start annual number: 7 End annual number: 7 |
|
| PR1001 | Payment of annual fee |
Payment date: 20120221 Start annual number: 8 End annual number: 8 |
|
| FPAY | Annual fee payment |
Payment date: 20130221 Year of fee payment: 9 |
|
| PR1001 | Payment of annual fee |
Payment date: 20130221 Start annual number: 9 End annual number: 9 |
|
| FPAY | Annual fee payment |
Payment date: 20140221 Year of fee payment: 10 |
|
| PR1001 | Payment of annual fee |
Payment date: 20140221 Start annual number: 10 End annual number: 10 |
|
| FPAY | Annual fee payment |
Payment date: 20150223 Year of fee payment: 11 |
|
| PR1001 | Payment of annual fee |
Payment date: 20150223 Start annual number: 11 End annual number: 11 |
|
| FPAY | Annual fee payment |
Payment date: 20160222 Year of fee payment: 12 |
|
| PR1001 | Payment of annual fee |
Payment date: 20160222 Start annual number: 12 End annual number: 12 |
|
| FPAY | Annual fee payment |
Payment date: 20170221 Year of fee payment: 13 |
|
| PR1001 | Payment of annual fee |
Payment date: 20170221 Start annual number: 13 End annual number: 13 |
|
| FPAY | Annual fee payment |
Payment date: 20180409 Year of fee payment: 14 |
|
| PR1001 | Payment of annual fee |
Payment date: 20180409 Start annual number: 14 End annual number: 14 |
|
| FPAY | Annual fee payment |
Payment date: 20190220 Year of fee payment: 15 |
|
| PR1001 | Payment of annual fee |
Payment date: 20190220 Start annual number: 15 End annual number: 15 |
|
| FPAY | Annual fee payment |
Payment date: 20200205 Year of fee payment: 16 |
|
| PR1001 | Payment of annual fee |
Payment date: 20200205 Start annual number: 16 End annual number: 16 |
|
| PR1001 | Payment of annual fee |
Payment date: 20210209 Start annual number: 17 End annual number: 17 |
|
| PR1001 | Payment of annual fee |
Payment date: 20220221 Start annual number: 18 End annual number: 18 |
|
| PR1001 | Payment of annual fee |
Payment date: 20230221 Start annual number: 19 End annual number: 19 |
|
| PR1001 | Payment of annual fee |
Payment date: 20240221 Start annual number: 20 End annual number: 20 |
|
| PC1801 | Expiration of term |
Termination date: 20250608 Termination category: Expiration of duration |