KR101006977B1 - 수치지도 제작 시 수치데이터를 보상하는 방법 - Google Patents

Abstract

본 발명에서 수치지도의 정확성을 높이기 위한 수치데이터를 보상하는 방법을 개시한다.
본 발명에 따른 방법은, 디지털 항공 수치지도를 보정하기 위한 방법에 있어서, a) 좌표 정보를 갖는 디지털 항공 사진을 소정 거리로 구획하는 단계; b) 상기 디지털 항공 사진의 구획 위치에 대응하는 좌표를 추출하는 단계; c) 해당 좌표 및 구획거리에 대응하는 지구체의 원주 거리를 자오선 및 횡 곡률 반경(Rm, Rn)에 따라 산출하는 단계; d) 상기 디지털 항공 사진에서의 구획 거리와, 상기 자오선 및 횡 곡률 반경(Rm, Rn)에 따라 산출된 거리로부터 각각의 차이 값(D, D')를 산출하는 단계; e) 상기 차이 값(D, D')에 근거하여 현재의 구획 좌표를 중심으로 상기 디지털 항공 사진을 보정 하는 단계; 및 f) 상기 구획 위치에 대응하는 좌표별로 상기 c) 단계 내지 e) 단계를 반복하여 상기 디지털 항공 사진을 보정하여 수치지도를 완성하는 단계;로 이루어진다.
따라서, 본 발명은 지구체의 곡률에 대응하는 수치지도 보정을 통해, 실질적인 지상 실측 데이터에 근접할 수 있어 데이터의 신뢰성을 높이는 효과가 있다.

Description

수치지도 제작 시 수치데이터를 보상하는 방법{METHOD FOR SUPPLEMENTING MAP DATA DURING A PRODUCTIVE DIGITAL MAP}

본 발명은 실거리 수치지도 제작에 관한 것으로, 더욱 상세하게는 GPS 및 INS를 이용한 수치지도 제작 시 실질적인 거리 산출이 가능하도록 평면 지형 데이터를 곡면 지형 데이터로 변환할 수 있는 수치지도 제작 시 수치데이터를 보상하는 방법에 관한 것이다.

일반적으로, 인공위성 자동 위치측정 시스템(GPS : Global Positioning System)이란 범지구적 위치 결정 시스템이라 할 수 있다. 인공위성 자동 위치측정 시스템을 이용하면 지상, 해상, 공중 등 지구상의 어느 곳에서나 시간 제약 없이 인공위성에서 발신하는 정보를 수신하여 정지 또는 이동하는 물체의 절대적 위치를 측정할 수 있다. GPS를 이용한 측량은 시공간의 제약 없이 높은 정확도의 위치정보 취득이 가능하므로 수치지도 제작에 있어 필수적인 요소라 할 수 있다.

한편, 수치지도는 지형 요소에 관한 정보들을 계층화된 레이어별로 구분하여 데이터베이스화하고, 전산화를 통해 이용 목적에 따라 지도를 자유로이 작성할 수 있으므로 지형공간정보 시스템뿐만 아니라 여러 분야에 다양한 응용이 가능하다. 또한, 수치지도는 입출력, 편집, 수정된 수치지도 자료를 그대로 이용할 수 있고, 변환과 검색이 가능하며, 건물, 도로 하천, 지형 등의 구조화된 데이터베이스 구축이 가능하기 때문에 다른 종료의 자료와 조합시켜 각종 형식의 지도를 만들 수 있다. 그러므로 구조화된 수치지도는 지리정보 시스템의 가장 기본적인 자료이기 때문에 매우 중요하다.

항공사진을 이용한 수치지도 제작은 일반적으로 항공 삼각 측량 방법을 이용하여 제작된다. 항공 삼각 측량은 촬영노선을 따라 60%의 중복도를 갖도록 사진을 촬영하고, 사진에서는 실제 지상을 측량한 점의 위치를 사진에서 찾아내어 사진에 실제 좌표를 등록할 수 있다. 또한 등록된 지상좌표를 이용하여, 수학적으로 사진의 중심위치와 회전각을 나타내는 외부표정 요소를 구할 수 있는데, 이러한 일련의 과정을 항공 삼각 측량이라 한다.

이와 같이 항공 삼각 측량을 이용한 수치지도 제작 과정은, 먼저 항공기에 장착된 카메라를 이용하여 일정 지역의 항공사진을 촬영한 후, 지상 측량을 통해 지상 기준점의 좌표를 구하고, 인접 사진에서 동일점을 추출하는 공액점 추출과정을 통해 여러 장의 사진을 연결하고, 지상기준점의 위치를 항공사진에서 찾아내어 실제 측량 좌표를 입력함으로써 촬영된 모든 항공사진이 좌표를 가질 수 있다. 이러한 과정을 통해 얻어진 외부 표정을 요소를 이용하여 항공사진의 수치지도를 제작한다.

전술된 방법에 기인하는 수치지도 제작 방법을 이하의 선행기술문헌을 통해 살펴 보면 다음과 같다.

대한민국 등록특허 10-0473950호 "수치사진측량 기술을 이용한 항공사진 촬영면적이 기록된 수치지도의 제작방법"에서 수치사진측량기술을 통해 얻어지는 정사영상 및 표정요소들을 이용하여 항공사진 촬영위치를 점사상(point feature)이 아닌 면 사상(polygon feature)으로 표시함으로써 명확하게 항공사진 촬영면적을 제공하기 위한 수치지도의 제작방법을 개시하고 있다.

이를 첨부된 예시도면에 의거 상세히 설명하면 다음과 같다.

먼저, 도 1은 종래 수치사진측량시스템과 항공사진 촬영면적이 기록된 수치지도 제작과의 관계도이다. 도시된 바와 같이, 수치지도의 제작에는 수치사진측량 기술을 통해 얻어지는 정사영상 및 표정요소 데이터가 이용된다. 이를 상세히 설명하면, 항공사진을 촬영(S211)하고 촬영된 항공사진 필름은 필름현상기를 통해 현상하며 이를 자동독취기를 통해 항공사진영상으로 변환(S212)한 후, 수치사진측량 시스템에 입력된다. 입력된 항공사진영상은 수치사진측량 시스템 상에서 내부표정(S221)하는 과정 및 GPS(global positioning system) 측량(S25a) 결과를 지상기준점 좌표로 이용하여 외부표정(S222)하는 과정과 정사편위 수정 과정을 거쳐 항공사진의 정사영상(S231)으로 제작된다.

상기와 같이 제작된 정사영상과 수치사진측량 시스템을 이용하여 항공사진 촬영면적이 기록된 수치지도가 제작된다. 상기 단계에서 얻어진 정사영상으로부터 촬영면적을 추출한다. 촬영면적의 추출은 (i) 대상파일(정사영상)을 선택하고 영상수치값을 검색(S311~S313)하는 단계와; (ii) 검색영역을 그룹화하고 폴리곤으로 변환(S321~S324)하는 단계와; (iii) 촬영면적 폴리곤을 단순화(S331~S332)한 다음 저장(S341)하는 단계로 이루어진다. 구체적으로, 촬영영역을 추출한 정사영상 대상파일을 선택한 후, 영상의 화소값 검색을 통해 항공사진의 촬영영역을 검색한다.

이때, 촬영영역 바깥쪽 화소들의 수치값은 일반적으로 '0' 또는 '0'에 근사한 양(+)의 값으로 화상강조가 이루어진 영상으로부터 화소의 수치값을 이용해 촬영면적을 효과적으로 검색할 수 있다. 검색결과를 확인한 후, 촬영영역이 충분히 검색되지 않았다면 수치를 +1씩 증가시켜 다시 검색하여 촬영면적이 최대한 명확히 추출되도록 한다. 검색이 완료된 촬영영역 화소들을 그룹화하여 촬영영역과 촬영영역 이외의 지역으로 구분하고 레스터(raster) 형태의 자료를 벡터(vector) 형태의 자료로 변환한다. 이때, 검색과정에서 발생한 소수의 화소 그룹들이 작은 폴리곤들로 변환된다.

실질적인 촬영영역은 일반적으로 변환된 폴리곤들 중 가장 큰 폴리곤을 가리키므로 최대면적 폴리곤을 제외한 기타 폴리곤들을 제거하고 나면 촬영면적을 가리키는 단 하나의 폴리곤만 남겨지게 된다. 결과를 확인하여 이상이 없으면 폴리곤 변환 단계를 완료한다. 추출된 폴리곤은 대개 작은 티클 화소들로 인해 요철이 있는 복잡한 형태로 나타나 있어 불필요한 버틱스(vertix)를 많이 포함하고 있으므로, 폴리곤을 구성하는 개개의 선분들 중 특정 길이(화소 길이)보다 작은 선분들을 제거하여 형태가 일반화된 촬영면적을 구한다. 상기와 같은 단순화는 정사영상 개개의 형태에 따라 적용에 차이를 줄 수 있으므로, 단순화의 결과를 확인하여 이상이 있다면 제거할 길이를 조정하여 다시 단순화를 시행하고, 이상이 없다면 해당 정사영상의 촬영면적을 폴리곤으로 저장한다.

수치지도를 제작하기 위하여, 상기 단계에서 수집된 폴리곤 자료를 수치지도로 저장하여 항공사진 촬영면적이 기록된 수치지도를 제작한다. 사업수행 지역을 포함하는 촬영은 일반적으로 종방향으로 60%, 횡방향으로 30%의 중복영역을 가진 다수의 항공사진들로 구성된다. 따라서, 항공사진 촬영면적이 기록된 수치지도는 상기의 과정을 통해 수집된 촬영면적 폴리곤들의 집합으로 구성되게 된다. 첫 번째로 추출한 항공사진 촬영면적 폴리곤에서부터 n번째로 추출한 항공사진 촬영면적 폴리곤까지를 통합하여 항공사진 촬영면적이 기록된 수치지도를 제작한다. 이때, 해당 폴리곤들이 모두 포함된 것을 검사하여 이상이 없다면 작업을 종료한다.

이와 같이, 종래 수치지도는 수치사진측량기술을 통해 얻어지는 정사영상 및 표정요소들을 이용하여 항공사진 촬영위치를 면 사상(polygon feature)으로 표시함으로써, 촬영면적, 인접 촬영사진과의 종중복도 및 횡중복도 면적을 시각적으로 명확하게 나타낼 수 있도록 하고 있다. 그러나, 항공사진의 촬영위치를 면사상으로 표시할 경우에도, 지구체의 곡률로 인해 사진의 왜곡이 발생될 수밖에 없으며, 곡률에 대한 영상 보정이 수행되지 못하고 있다. 즉, 수치지도를 위한 항공사진은 국소적인 부분에 한하여 이루어지기 때문에, 지구체의 곡률에 대한 오차를 무시할 수 있어 수치지도와 지면상에서의 실질적인 이동 거리에 오차가 발생한다.

본 발명은 이와 같은 문제점을 해결하기 위해 창출된 것으로, 본 발명의 목적은 항공사진을 기반으로 제작된 수치지도를 좌표에 대응하는 지구체의 곡률정보를 가미시켜 실측 거리와 유사한 수치지도를 생성할 수 있는 수치지도 제작 시 수치데이터를 보상하는 방법을 제공함에 있다.

본 발명의 다른 목적은, 항공사진 촬영 지점에 대한 GPS 좌표 정보와, 이에 대응하는 지구체의 법선 각도 정보를 토대로 항공사진을 보정함으로써, 실질적인 실측 거리에 근접한 수치지도를 생성할 수 있는 수치지도 제작 시 수치데이터를 보상하는 방법을 제공함에 있다.

상기 목적을 달성하기 위한 본 발명의 관점에 따른 수치지도 제작 시 수치데이터를 보상하는 방법은, 디지털 항공 수치지도를 보정하기 위한 방법에 있어서, a) 좌표 정보를 갖는 디지털 항공 사진을 소정 거리로 구획하는 단계; b) 상기 디지털 항공 사진의 구획 위치에 대응하는 좌표를 추출하는 단계; c) 해당 좌표 및 구획거리에 대응하는 지구체의 원주 거리를 자오선 및 횡 곡률 반경(Rm, Rn)에 따라 산출하는 단계; d) 상기 디지털 항공 사진에서의 구획 거리와, 상기 자오선 및 횡 곡률 반경(Rm, Rn)에 따라 산출된 거리로부터 각각의 차이 값(D, D')를 산출하는 단계; e) 상기 차이 값(D, D')에 근거하여 현재의 구획 좌표를 중심으로 상기 디지털 항공 사진을 보정 하는 단계; 및 f) 상기 구획 위치에 대응하는 좌표별로 상기 c) 단계 내지 e) 단계를 반복하여 상기 디지털 항공 사진을 보정하여 수치지도를 완성하는 단계;로 이루어진 것을 특징으로 한다.

본 발명의 바람직한 실시 예에 따른 e) 단계는, e-1) 상기 구획거리를 N개로 분할하는 단계; e-2) 상기 오차(D)를 상기 N개로 나눈 값을 해당 구획점으로부터 1/N씩 분할된 구획 거리별로 누적시켜 오차를 보정하는 단계;로 이루어진 것을 특징으로 한다.

본 발명에서 제시되는 수치지도 제작 시 수치데이터를 보상하는 방법은, 항공사진을 기반으로 작성된 수치지도를 지상에서의 관측 데이터에 근거한 보정 절차를 줄일 수 있어 수치지도 제작에 따른 비용을 절감할 수 있는 효과가 있다. 또한 본 발명은 지구체의 곡률에 대응하는 수치지도 보정을 통해, 실질적인 지상 실측 데이터에 근접할 수 있어 데이터의 신뢰성을 높이는 효과가 있다.

도 1은 종래 수치지도 제작방법을 나타낸 도면이다.
도 2는 본 발명에 따른 수치지도 보정을 위한 원리를 설명하는 도면이다.
도 3은 본 발명에 따른 수치지도 보정 방법을 설명하기 위한 플로우챠트이다.
도 4 및 도 5는 본 발명의 보정원리를 부연설명하기 위한 도면이다.

이하, 본 발명의 바람직한 실시 예를 첨부된 예시도면에 의거 상세히 설명하면 다음과 같다.

도 2a, 2b는 본 발명에 따른 원리를 설명하기 위한 도면이다. 먼저, 도시된 바와 같이, 항공사진의 촬영위치에서 임의의 면적을 촬영할 경우, 도 2b와 같이 지구체 곡률에 따른 곡면은 평면화된 사진정보로 변환된다. 이 과정에서 사진정보에서 도시되는 임의의 두 지점 간 거리와, 지구체에서의 동일한 두 지점 간 거리는 차이가 발생한다. 따라서, 다수의 항공사진을 중첩 및 연결하여 특정 지역에 대한 항공사진을 제작할 경우, 그 오차는 더욱 커지게 된다. 본 발명에서는 이러한 오차를 줄이기 위해 도 2c와 같이 이미 촬영된 영상을 촬영 좌표에 따라 지구체와 동일한 곡률을 갖고 항공사진을 보정 하는 것이다.

여기서, 곡률이라 함은 곡면이나 곡선상 각 점의 굽은 정도를 표시하는 것으로, 본 발명의 실시 예에서는 자오선 곡률반경(Radius of Curvature of the Meridian : Rm)과 횡 곡률반경(Radius of Curvature of the Prime Vertical : Rn)을 기반으로 이루어진다.

상기 자오선 곡률반경(Rm)은 아래의 수학식 1;

이고,

횡 곡률반경(Rn)은 아래의 수학식 2;

으로 표현된다. 여기서, a : 장반경, e : 이심율(離心率 : eccentricity), φ : 임의 위도 값을 나타낸다.

이하, 본 발명에서는 상기 자오선 곡률반경(Rm)을 Y축으로 인지하고, 상기 횡 곡률반경(Rn)은 X축으로 인지할 것이며, 항공사진의 X축 방향을 자오선 곡률반경에 따라 연장하고, 항공사진의 Y축 방향을 횡 곡률반경에 따라 연장하여 항공사진을 보정할 것이다. 그리고, 수치지도의 보정은 본 발명의 실시 예에 따라 1km를 기준으로 보정할 것이나, 오차를 줄이기 위한 보정 단위는 그 이하가 될 수 있음은 당연할 것이다.

도 3은 본 발명에 따른 수치지도의 보정 방법을 설명하기 위한 플로우챠트이고, 더불어 첨부되는 도면은 본 발명의 원리를 구체적으로 설명하기 위한 도면이다.

먼저, 도시된 바와 같이 S301 단계에서 디지털 항공 수치사진을 데이터베이스화하고, 수치사진을 소정 단위의 정방형으로 구획한다. 구획단위는 관리자의 의도에 따라 달리 될 수 있으나, 본 발명의 실시 예에서는 가로 세로 1km 단위로 설정한다. 구획단위는 수치사진을 보정하기 위한 보정의 기준점으로, 1km 단위로 좌표를 설정하고, 각 좌표를 중심으로 보정한다.

즉, 도 4에 도시된 바와 같이 항공사진(P)은 지구체(E)의 굴곡면을 평면화한 것으로, 임의의 위치 또는 최초 기준점(A)은 지구체에서 a점에 대응하고, 상기 A점으로부터 구획된 위치(A')는 지구체에서 a'점에 대응한다. 물론, AA' 거리와 aa' 거리는 상이할 것이다. 한편, 도 4에서 도시된 항공사진(P)는 X축과 Y축으로 그 방향이 존재하는데, 이는 자오선을 따라 정의되는 위도가 Y축이고, X축은 경도를 나타낸다. 도면에서 도시된 방향은 설명의 편의상 X축과 Y축을 통념에서 벗어나 설명하고 있다.

디지털 항공 수치사진을 기반으로 소정 단위별로 구획되면, S303 단계에서 디지털 항공 수치사진의 구획 위치에 대응하는 좌표를 추출한다. 상기 좌표는 항공사진 촬영시 INS, GPS 장비에 의해 정확한 좌표 정보가 저장되는 것으로, 항공기 촬영 당시 항공기의 운행 경로, 항공기의 운행 각도, 위치 등을 기반으로 생성된 좌표이다. 수치지도에서의 좌표 정보는 다수의 장비를 이용함으로써 그 정확도가 매우 높으며, 이에 대한 구체적인 설명은 생략한다.

이후, S305 단계로 진입하여 해당 좌표 및 기 설정된 구획 거리에 대응하는 지구체 원주 거리 즉, 실질적인 거리를 산출한다. 이는 항공사진에서의 거리와 타원 구조를 갖는 지구체에서의 실질적인 거리가 다르기 때문에, 항공 사진에서의 1km 거리가 실질적으로 1km를 넘게 된다. 더욱이, 1km는 매우 작은 거리이나 작은 오차가 중첩될 경우 많은 절차를 거쳐 보정해야 한다. 본 발명은 이러한 문제를 해결하기 위한 것으로, 오차의 범위를 최소화하기 위해 1km의 구획 거리를 기준으로 보정하는 것이다.

기본적으로 지구체의 원주를 인지하기 위해서는 도 4에 도시된 aa'에 해당하는 원주는 진원에 의한 원주로 정의한다. 그리고, a'a'' 또한 그 원주는 진원에 의해 정의된다. 그러나, aa''에 해당하는 원주는 진원이 아닌 타원으로 간주하는 것으로, 이는 a, a', a'' 각 지점에 대한 좌표를 중심으로 각 위도 또는 경도에서의 새로운 곡률을 적용하기 때문이다.

항공사진(P)에서 A'A'는 계산의 편의상 1km라 할 때, AA'는 0.5km이고, 수학식 1에서 전술된 자오선 곡률에 따른 지구체의 반경을 Rm이라 할 때, 자오선 좌표 a 지점에서의 자오선 반경은 Rma라 하면, 이때 aa'의 실질적인 거리(T)는 수학식 3;

로 정의된다. 여기서 α는 항공사진(P)에서 임의의 위치(A)를 중심으로 1km에 해당하는 지구체의 중심각도이다. 이를 근거로, 수학식 4와 같이;

해당 구획거리에 대응하는 지구체의 중심각도이다.

따라서, 상기 원주(T)는 수학식 5;

이 된다. 결국, 항공사진(P)에서의 1km 거리가 실질적인 거리로서 원주(T)에 해당하는 것으로, 거리 오차(D)는 (T-1)이 된다. 그러나, 이러한 오차(D)는 중심점으로부터 멀어져 가는 위치에서 오차의 적용이 용이한 반면, 중심점에 근접한 위치에서는 오차(D)를 적용하기에 또 다른 오차를 발생시키게 된다. 즉, a' 지점에서 오차(D)를 적용하는 것은 바람직하나 a 지점에서는 오차(D)의 의미가 없게 된다.

S307 단계에서는 이러한 오차를 적용하기 위한 과정으로, 기 설정된 구획거리를 N개로 분할한 후, 상기 오차(D)를 1/N씩 누적시켜 보정한다. 즉, 도 5에 도시된 도면과 같이 구획거리를 10개로 분할할 경우, 항공사진(P)를 양방향(Y방향, -Y방향)으로 각각 5개씩 N1, N2, ...N5로 분할한 후, N1에서 상기 오차(D)의 10%만을 적용하고, N2에서 상기 오차(D)의 20%를 적용하며, N5에서는 상기 오차(D)를 100% 적용하는 것이다. 여기서, N의 개 수 즉, 분할 개 수가 많을수록 오차의 폭은 줄어들 것이며, 본 발명에서는 계산의 편의상 10개로 분할하고 있으나, 관리자의 의도에 따라 정밀한 수치지도를 보정 하고자 할 경우에는 그 분할 개 수를 충분히 늘릴 수 있음은 당연할 것이다.

전술된 수치지도의 보정은 한 방향 즉, 자오선 방향으로 보정 하였으나 X축 방향으로서 횡 방향 보정이 필요할 것이다. S309 단계는 Y축 이외에 X축을 보정하는 것으로, 수학식 2에서 서술된 횡 곡률반경(Rn)을 근거로 상기 과정을 반복한다.

예컨대, a 지점에서의 횡 곡률반경(Rn)에 따른 지구체의 곡면 거리(T')는 수학식 6;

과 같이 나타내며, 오차(D')는 T'- 1로 나타낸다. 물론, 이는 구획거리를 1km를 기준으로 이루어지며, 기준거리(k)에 따라 오차(D')는 T'-k로 표현될 것이다. 이와 같이 횡 방향 곡률에 따른 오차(D')에 대한 실거리와의 차이 값이 산출되면, 이를 구획 거리 내의 분할 개 수에 따라 1/N' 등분하여 전술한 바와 같이, 오차 값을 누적시켜 횡축 방향의 곡률 보정을 수행한다. 이와 같이, X축 방향 및 Y축 방향에 대한 보정으로서, 항공사진을 자오선 곡률 및 횡 곡률에 따라 보정한 후, S311 단계와 같이 보정된 디지털 항공 수치 사진을 데이터베이스화한다. 따라서, 본 발명에서는 항공사진에 의한 수치 지도를 정확하게 보정하여 지도의 신뢰성을 높이게 되는 것이다.

전술된 바와 같이, 본 발명에 따른 수치지도 제작 시 수치데이터를 보상하는 방법은, 항공사진을 기반으로 제작된 수치지도를 좌표에 대응하는 지구체의 곡률정보를 가미시켜 실측 거리와 유사한 수치지도를 생성하기 때문에, 지면에서의 두 지점 간 거리와, 수치지도상에서의 동일한 두 지점 간 거리가 동일하게 산출되어 정확한 거리 측정이 가능하게 된다. 따라서 본 발명에 따른 수치지도는 그 신뢰성을 높이고, 활용도를 증대시킬 수 있기 때문에 산업적 이용가치가 충분히 높을 것으로 판단된다.

P : 항공사진 E : 지구체
Rm : 자오선 곡률 반경 Rn : 횡 곡률 반경
D, D' : 오차 N : 구획 내 분할 개 수
T : 지구체 상에서의 실질적인 거리

Claims (4)

1. 수치지도 제작 시 수치데이터를 보상하는 방법에 있어서,
   a) 좌표 정보를 갖는 디지털 항공 사진을 소정 거리로 구획하는 단계;
   b) 상기 디지털 항공 사진의 구획 위치에 대응하는 좌표를 추출하는 단계;
   c) 해당 좌표 및 구획거리에 대응하는 지구체의 원주 거리를 자오선 및 횡 곡률 반경(Rm, Rn)에 따라 산출하는 단계;
   d) 상기 디지털 항공 사진에서의 구획 거리와, 상기 자오선 및 횡 곡률 반경(Rm, Rn)에 따라 산출된 거리로부터 각각의 차이 값(D, D')를 산출하는 단계;
   e) 상기 차이 값(D, D')에 근거하여 현재의 구획 좌표를 중심으로 상기 디지털 항공 사진을 보정 하는 단계; 및
   f) 상기 구획 위치에 대응하는 좌표별로 상기 c) 단계 내지 e) 단계를 반복하여 상기 디지털 항공 사진을 보정하여 수치지도를 완성하는 단계;로 이루어진 것을 특징으로 하는 수치지도 제작 시 수치데이터를 보상하는 방법.
2. 제 1 항에 있어서,
   상기 자오선 곡률반경(Rm) 및 횡 곡률반경(Rn)은 아래의 수학식 7
   [수학식 7]
   

   

   
   <a : 장반경, e : 이심율(離心率 : eccentricity), φ : 임의 위도 값>
   에 근거한 것을 특징으로 하는 수치지도 제작 시 수치데이터를 보상하는 방법.
3. 제 1 항에 있어서 상기 e) 단계는,
   e-1) 상기 구획거리를 N개로 분할하는 단계;
   e-2) 상기 오차(D)를 상기 N개로 나눈 값을 해당 구획점으로부터 1/N씩 분할된 구획 거리별로 누적시켜 오차를 보정하는 단계;로 이루어진 것을 특징으로 하는 수치지도 제작 시 수치데이터를 보상하는 방법.
4. 제 1 항에 있어서,
   상기 구획 거리는 자오선 방향으로 1km, 횡 방향으로 1km인 것을 특징으로 하는 수치지도 제작 시 수치데이터를 보상하는 방법.
   

Images