KR101973726B1 - 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치 - Google Patents

모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치 Download PDF

Info

Publication number
KR101973726B1
KR101973726B1 KR1020180170868A KR20180170868A KR101973726B1 KR 101973726 B1 KR101973726 B1 KR 101973726B1 KR 1020180170868 A KR1020180170868 A KR 1020180170868A KR 20180170868 A KR20180170868 A KR 20180170868A KR 101973726 B1 KR101973726 B1 KR 101973726B1
Authority
KR
South Korea
Prior art keywords
numerical
model
elevation model
minimum
morphological filtering
Prior art date
Application number
KR1020180170868A
Other languages
English (en)
Inventor
윤준희
김영석
김태훈
Original Assignee
한국건설기술연구원
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 한국건설기술연구원 filed Critical 한국건설기술연구원
Priority to KR1020180170868A priority Critical patent/KR101973726B1/ko
Application granted granted Critical
Publication of KR101973726B1 publication Critical patent/KR101973726B1/ko

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
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
    • G06T5/006—
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/20—Image enhancement or restoration using local operators
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/80—Geometric correction

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Remote Sensing (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Image Processing (AREA)

Abstract

모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치가 제공된다. 전자장치의 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법은, 전자장치가, 대상지역이 경사 지역을 포함하는 경우, 격자 기반의 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하는 단계와, 생성된 임시 수치표고모델에 최대 모폴로지컬 필터링을 적용하여 경사 지역에 해당하는 수치표고모델을 생성하는 단계를 포함한다.

Description

모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치{Apparatus and method for generating a DEM from a DAM of a target area using morphlogical filter}
본 발명은 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치에 관한 것으로서, 보다 상세하게는, 수치표면모델에 적어도 2회의 모폴로지컬 필터링을 연속 적용하여 왜곡을 최소화할 수 있는 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치에 관한 것이다.
수치표고모델(DEM: Digital Elevation Model)은 수치고도모델이라고도 불린다. 일반적으로 수치표고모델은 지표면위의 모든 인공구조물(빌딩, 탑, 전선 등)과 나무, 식생들이 제거된 bare-ground의 높이를 격자 형태로 표현한 것을 의미한다.
유사한 의미로 수치표면모델(DSM: Digital Surface Model)이 있으며, 수치표면모델은 지표면에 인공구조물과 나무, 식생들이 포함된 높이를 격자 형태로 표현한 것을 의미한다.
수치표면모델은 삼차원 도시 모형 구현, 통신, 도시계획 등으로 많이 활용되며, 수치표고모델은 건설현장 업무량 산정을 위한 토공량 계산, 수리학적 모델링, 지형 안정성 분석, 토양 매핑, 지형의 고도 차로 인해 발생하는 영상자료의 기하학적 왜곡 등에 많이 활용된다.
한편, 최근 항공라이다 측량, 드론 측량의 발전으로 인하여 수치표면모델의 제작이 쉬워지고 있다. 항공라이다의 경우, 레이저 스캐너를 통해 지상의 3차원 포인트 클라우드를 취득하고 이를 격자형태로 만들어 수치표면모델로 사용한다. 드론측량의 경우, 항공사진측량 기법에 의해 연속 촬영된 복수의 영상을 이용하여 지상의 3차원 포인트 클라우드를 만들고 이를 격자 형태로 만들어 수치표면모델로 활용한다.
지상의 포인트클라우드 데이터는 빌딩, 전선, 탑, 나무, 식생의 높이들을 모두 감안한 값이 나오기 때문에 바로 수치표면모델을 만들 수 있으며, 이 데이터로부터 Non-ground 데이터를 제외하면 수치표고모델이 된다.
지상 포인트클라우드 데이터에서 non-ground 데이터를 제외하는 방법 중 일반적인 방법은 모폴로지컬 필터를 활용하는 방식이며 많은 연구 그룹들에 의해서 시도되어 왔다(Killian et al., 1996, Hug, 1997, Morgan and Tempfli, 2000, Morgan and Habib, 2002). 모폴로지컬 필터를 활용하는 방식은 수치표고모델의 격자에 써치 윈도우를 만들고 그 써치 윈도우 내에서 가장 낮은 높이 값을 검사하는 것으로서, 전체 지역에 적용하면 수치표고모델이 만들어 진다.
그러나, 기존의 방식으로 수치표고모델을 생성하는 경우, 모폴로지컬 필터는 경사가 있는 지역에 적용하면, 실제 지면의 고도보다 낮은 고도가 설정되는 심각한 왜곡이 발생한다. 또한, 도심지역에서는 지하철 환풍구 밑의 지역에 의해 반사된 라이다 데이터로 인해 실제 지면 높이를 왜곡하는 단점이 있다.
국내 등록특허 제10-1008394호(2011.01.07. 등록)
전술한 문제점을 해결하기 위하여 본 발명이 이루고자 하는 기술적 과제는, 수치표면모델에 적어도 2회의 모폴로지컬 필터링을 적용하여 수치표고모델 생성시 발생하는 왜곡을 최소화할 수 있는 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치를 제시하는 데 있다.
본 발명의 해결과제는 이상에서 언급된 것들에 한정되지 않으며, 언급되지 아니한 다른 해결과제들은 아래의 기재로부터 당업자에게 명확하게 이해될 수 있을 것이다.
전술한 기술적 과제를 해결하기 위한 수단으로서, 본 발명의 실시 예에 따르면, 전자장치의 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법은, (A) 상기 전자장치가, 상기 대상지역이 경사 지역을 포함하는 경우, 격자 기반의 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하는 단계; 및 (B) 상기 전자장치가, 상기 (A) 단계에서 생성된 임시 수치표고모델에 최대 모폴로지컬 필터링을 적용하여 상기 경사 지역을 포함하는 대상지역에 해당하는 수치표고모델을 생성하는 단계;를 포함한다.
상기 (A) 단계는, 상기 전자장치가, 상기 격자 기반의 수치표면모델에 사전에 정해진 크기의 써치 윈도우를 적용하고, 상기 써치 윈도우 내에서 높이가 가장 낮은 격자의 높이를 상기 써치 윈도우의 중심점의 최소 지면 높이로 정한 후, 상기 써치 윈도우의 중심점과 상기 정해진 최소 지면 높이를 이용하여 상기 임시 수치표고모델을 생성한다.
상기 (B) 단계는, 상기 전자장치가, 상기 임시 수치표고모델에 써치 윈도우를 적용하고, 상기 써치 윈도우 내에서 높이가 가장 높은 격자의 높이를 상기 써치 윈도우의 중심점의 최종 지면 높이로 정한 후, 상기 써치 윈도우의 중심점과 상기 정해진 최종 지면 높이를 이용하여 상기 경사 지역을 포함하는 대상지역에 해당하는 수치표고모델을 생성한다.
상기 (A) 단계 이전 및 상기 (B) 단계 이후 중 어느 하나에, (C) 상기 전자장치가, 상기 대상지역이 평지인 경우, 상기 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 상기 평지에 해당하는 수치표고모델을 생성하는 단계;를 더 포함한다.
상기 (A) 단계 이전 및 상기 (B) 단계 이후 중 어느 하나에, (D) 상기 전자장치가, 상기 대상지역이 환풍구를 포함하는 경우, 상기 수치표면모델에 최대 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하는 단계; 및 (E) 상기 전자장치가, 상기 (D) 단계에서 생성된 임시 수치표고모델에 최소 모폴로지컬 필터링을 적용하여 상기 환풍구를 포함하는 대상지역에 해당하는 수치표고모델을 생성하는 단계;를 더 포함한다.
한편, 본 발명의 다른 실시 예에 따르면, 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 전자장치는, 상기 대상지역이 경사 지역을 포함하는 경우, 격자 기반의 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하는 최소 필터링부; 및 상기 최소 필터링부에서 생성된 임시 수치표고모델에 최대 모폴로지컬 필터링을 적용하여 상기 경사 지역을 포함하는 대상지역에 해당하는 수치표고모델을 생성하는 최대 필터링부;를 포함한다.
상기 최소 필터링부는, 상기 격자 기반의 수치표면모델에 사전에 정해진 크기의 써치 윈도우를 적용하고, 상기 써치 윈도우 내에서 높이가 가장 낮은 격자의 높이를 상기 써치 윈도우의 중심점의 최소 지면 높이로 정한 후, 상기 써치 윈도우의 중심점과 상기 정해진 최소 지면 높이를 이용하여 상기 임시 수치표고모델을 생성한다.
상기 최대 필터링부는, 상기 임시 수치표고모델에 써치 윈도우를 적용하고, 상기 써치 윈도우 내에서 높이가 가장 높은 격자의 높이를 상기 써치 윈도우의 중심점의 최종 지면 높이로 정한 후, 상기 써치 윈도우의 중심점과 상기 정해진 최종 지면 높이를 이용하여 상기 경사 지역을 포함하는 대상지역에 해당하는 수치표고모델을 생성한다.
상기 대상지역이 평지인 경우, 상기 최소 필터링부는, 상기 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 상기 평지에 해당하는 수치표고모델을 생성한다.
상기 대상지역이 환풍구를 포함하는 경우, 상기 최대 필터링부는, 상기 수치표면모델에 최대 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하고, 상기 최소 필터링부는, 상기 최대 필터링부에서 생성된 임시 수치표고모델에 최소 모폴로지컬 필터링을 적용하여 상기 환풍구를 포함하는 대상지역에 해당하는 수치표고모델을 생성한다.
본 발명에 따르면, 수치표면모델의 특성에 따라 적어도 2회의 최소 모폴로지컬 필터링과 최대 모폴로지컬 필터링을 적용함으로써 수치표고모델 생성 시 발생하는 왜곡을 최소화할 수 있다.
또한, 본 발명에 따르면, 대상지역이 경사 지역 또는 환풍구 지역 또는 평지 지역인지에 따라 적응적으로 필터링을 적용함으로써 보다 고품질의 수치표고모델을 생성할 수 있다.
본 발명의 효과는 이상에서 언급된 것들에 한정되지 않으며, 언급되지 아니한 다른 효과들은 아래의 기재로부터 당업자에게 명확하게 이해될 수 있을 것이다.
도 1은 본 발명의 실시 예에 따른 전자장치(1000)의 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법을 도시한 흐름도,
도 2의 (a)는 실제 표면의 프로파일을 도시한 도면,
도 2의 (b)는 실제 표면에 대해 생성되는 수치표고모델의 이상적인 프로파일을 도시한 도면
도 3의 (a)는 실제 표면에 대해 생성된 수치표면모델의 프로파일을 도시한 도면,
도 3의 (b)는 S120단계의 최소 모폴로지컬 필터 적용 결과를 도시한 도면,
도 4는 도 3의 (b)와 같은 임시 수치표고모델에 최대 모폴로지컬 필터 적용 결과를 도시한 도면,
도 5는 S140단계에 의해 생성된 수치표고모형과 도 2에 도시된 원지형을 비교한 도면,
도 6은 수치표면모델과 모폴로지컬 필터링에 사용할 써치 윈도우의 일 예를 도시한 도면,
도 7은 수치표면모델을 격자 형태로 표시하고, 써치 윈도우를 적용하여 모폴로지컬 필터링을 수행하는 동작을 설명하기 위한 도면,
도 8은 도 1의 S120단계 또는 S170단계 또는 S195단계를 보다 자세히 설명하기 위한 흐름도,
도 9는 도 1의 S130단계 또는 S160단계를 보다 자세히 설명하기 위한 흐름도,
도 10은 본 발명의 일 실시 예에 따른 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 전자장치를 도시한 블록도,
도 11은 본 발명의 다른 실시 예에 따른 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 전자장치를 도시한 블록도,
도 12는 본 발명의 실시 예에 따라 최소 모폴로지컬 필터링 및 최대 모폴로지컬 필터링을 적용한 결과를 보여주는 그래프, 그리고,
도 13은 본 발명의 실시 예에 따라 환풍구 지역에 최대 모폴로지컬 필터링 및 최소 모폴로지컬 필터링을 적용한 결과를 설명하기 위한 그래프이다.
본 발명의 실시를 위한 구체적인 내용을 설명하기에 앞서, 본 명세서 및 청구범위에 사용된 용어나 단어는 발명자가 그 자신의 발명을 가장 최선의 방법으로 설명하기 위해 용어의 개념을 적절하게 정의할 수 있다는 원칙에 입각하여 본 발명의 기술적 사항에 부합하는 의미와 개념으로 해석되어야 할 것이다.
또한, 본 발명에 관련된 공지 기능 및 그 구성에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우에는, 그 구체적인 설명을 생략하였음에 유의해야 할 것이다.
이하, 본 발명에서 실시하고자 하는 구체적인 기술내용에 대해 첨부도면을 참조하여 상세하게 설명하기로 한다.
도 10에 도시된 전자장치(1000)의 각각의 구성은 기능 및 논리적으로 분리될 수 있음을 나타내는 것이며, 반드시 각각의 구성이 별도의 물리적 장치로 구분되거나 별도의 코드로 생성됨을 의미하는 것은 아님을 본 발명의 기술분야의 평균적 전문가는 용이하게 추론할 수 있을 것이다.
상기 전자장치(1000)는 소정의 데이터 프로세싱 장치에 설치되어 본 발명의 기술적 사상을 구현할 수 있다.
또한, 본 발명의 실시 예에 따른 상기 전자장치(1100)는 마이크로 프로세서, 메모리, FPGA(Field Programmable Gate Array), ASIC(Application Specific Integrated Circuit) 등 프로그램의 설치 및 실행이 가능한 전자 장치를 사용하여 구현될 수 있다.
도 1은 본 발명의 실시 예에 따른 전자장치(1000)의 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법을 도시한 흐름도이다.
도 1에 도시된 본 발명의 실시 예에 따른 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하는 동작을 수행하는 전자장치(1000)는 도 10 또는 도 11을 참조하여 설명할 전자장치(1000, 1100) 중 하나일 수 있다.
도 1을 참조하면, 전자장치(1000)는 수치표고모델(DEM)로 생성하기 위한 대상지역의 격자 기반의 수치표면모델(DSM)을 입력받는다(S100).
입력된 수치표면모델에 해당하는 대상지역이 경사 지역을 포함하면(S110), 전자장치(1000)는 입력된 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성할 수 있다(S120).
S120단계에서, 전자장치(1000)는 수치표면모델에 사전에 정해진 크기의 써치 윈도우를 적용하고, 써치 윈도우 내에서 높이가 가장 낮은 격자의 높이를 써치 윈도우의 중심점의 최소 지면 높이로 정한 후, 써치 윈도우의 중심점과 정해진 최소 지면 높이를 이용하여 임시 수치표고모델을 생성할 수 있다.
도 2의 (a)는 실제 표면의 프로파일을 도시한 도면이고, 도 2의 (b)는 실제 표면에 대해 생성되는 수치표고모델의 이상적인 프로파일을 도시한 도면이다. 도 2를 참조하면, 이상적인 수치표고모델의 프로파일(200b)은 실제 표면(200a)에서 나무나 건물 등의 객체만 제거되고 표면의 윤곽은 동일함을 알 수 있다.
도 3의 (a)는 실제 표면에 대해 생성된 수치표면모델의 프로파일을 도시한 도면이고, 도 3의 (b)는 S120단계의 최소 모폴로지컬 필터 적용 결과를 도시한 도면이다.
도 3에 도시된 것처럼 수치표면모델(300a)에 최소 모폴로지컬 필터에 의한 필터링을 적용하면, 경사 지역에는 고도 왜곡이 심하게 발생한다. 즉, 도 3의 (b)에 도시된 임시 수치표고모델(300b)처럼, 검사점(A)에서의 이상적인 수치표고모델의 값은 a 고도가 되어야 하지만, 써치 윈도우를 사용함으로써 써치 윈도우 내 최소값을 갖는 B 위치에서의 값(b 고도)이 검사점 A의 수치표고모델 결과값으로 되어, 왜곡이 발생하였다.
이러한 왜곡을 보상하기 위하여, 전자장치(1000)는 S120단계에서 생성된 임시 수치표고모델(300b)에 최대 모폴로지컬 필터링을 적용하여 도 4와 같은 수치표고모델(300c)을 생성할 수 있다(S130, S140).
도 4는 도 3의 (b)와 같은 임시 수치표고모델(300b)에 최대 모폴로지컬 필터 적용 결과를 도시한 도면이다.
도 4를 참조하면, S130단계에서, 전자장치(1000)는 임시 수치표고모델(300b)에 써치 윈도우를 적용하고, 써치 윈도우 내에서 높이가 가장 높은 격자의 높이를 써치 윈도우의 중심점의 최종 지면 높이로 정한 후, S140단계에서, 써치 윈도우의 중심점과 정해진 최종 지면 높이를 취합하여 경사 지역을 포함하는 대상지역에 해당하는 수치표고모델(300c)을 생성할 수 있다.
즉, 도 4에 도시된 것처럼, 임시 수치표고모델(300b)에 최대 모폴로지컬 필터링을 적용함으로써, 최소 모폴로지컬 필터링에 의해 왜곡된 고도가 보상(화살표 방향으로 이동)되었음을 알 수 있다.
도 5는 S140단계에 의해 생성된 수치표고모형과 도 2에 도시된 원지형을 비교한 도면이다.
도 5를 참조하면, 최소 모폴로지컬 필터와 최대 모폴로지컬 필터를 적용하여 생성된 수치표고모형(300c)과 원지형(200b) 사이의 왜곡은 최소화되었음을 알 수 있다.
한편, S100단계에서 입력받은 수치표면모델에 해당하는 대상지역이 지하로 형성된 환풍구를 포함하면(S150), 전자장치(1000)는 입력된 수치표면모델에 최대 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성할 수 있다(S160). 지하철 부근과 같이 땅 속에 형성되는 환풍구는 환풍구의 바닥에 의해 반사된 라이다 데이터로 인해 실제 지면 높이가 왜곡되어 표시되므로, 전자장치(1000)는 S120단계와는 반대로 최대 모폴로지컬 필터링을 먼저 적용한다.
그리고, 전자장치(1000)는 S160단계에서 생성된 임시 수치표고모델에 최소 모폴로지컬 필터링을 적용하여 환풍구를 포함하는 대상지역에 해당하는 수치표고모델을 생성할 수 있다(S170, S180).
또한, S100단계에서 입력받은 수치표면모델에 해당하는 대상지역이 평지인 경우(S190), 전자장치(1000)는 입력된 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 평지에 해당하는 수치표고모델을 생성할 수 있다(S195, S200).
이하에서는 도 6 내지 도 9를 참조하여 본 발명의 실시 예에 따라 최소 모폴로지컬 필터링을 적용하는 방식 및 최대 모폴로지컬 필터링을 적용하는 방식에 대해 자세히 설명한다.
도 6은 수치표면모델과 모폴로지컬 필터링에 사용할 써치 윈도우의 일 예를 도시한 도면이다.
써치 윈도우는 도 6에 도시된 화살표 방향으로 이동하면서 모폴로지컬 필터링에 사용되고, 우측까지 이동하면 다시 한 행 아래로 이동하여 좌측부터 화살표 방향으로 이동한다.
도 7은 수치표면모델을 격자 형태로 표시하고, 써치 윈도우를 적용하여 모폴로지컬 필터링을 수행하는 동작을 설명하기 위한 도면이다.
도 7을 참조하면, 써치 윈도우의 크기는 대상지역에 위치하는 평지를 기준으로 볼록한 형태의 객체(예를 들어, 건물, 탑, 전선, 나무, 식생 등)들의 최대 평면 넓이보다 크도록 설정된다. 이는 S120단계, S130단계, S160단계, S170단계 및 S195단계에서, 최소 지면 높이 또는 최대 지면 높이로서 객체의 높이가 정해지는 것을 방지하고 지면의 높이가 정해지도록 하기 위함이다.
써치 윈도우가 1~p의 격자로 이루어지며, (l, s)를 검사점이라고 하면, (k×k) 크기의 정사각형 써치 윈도우 내 격자 높이의 집합은 다음 [수학식 1]과 같이 표현된다.
Figure 112018131132824-pat00001
S120단계에서 최소 모폴로지컬 필터링을 수행한 결과값(Clsmin)은 [수학식 2]와 같이 표현될 수 있다.
Figure 112018131132824-pat00002
[수학식 2]에서, Clsmin은 검사점 (l, s)에서 검사점 (l, s)와 써치 윈도우의 중심점을 일치시켰을 때, 써치 윈도우 내에서 가장 낮은 격자의 높이값(고도)을 의미한다.
또한, S130단계에서 최대 모폴로지컬 필터링을 수행한 결과값(Clsmax)은 [수학식 3]과 같이 표현될 수 있다.
Figure 112018131132824-pat00003
[수학식 3]에서, Clsmax는 검사점 (l, s)에서 검사점 (l, s)와 써치 윈도우의 중심점을 일치시켰을 때, 써치 윈도우 내에서 가장 높은 격자의 높이값(고도)을 의미한다.
이러한 최소 모폴로지컬 필터링과 최대 모폴로지컬 필터링은 i(또는 l)는 1부터 m까지, j(또는 s)는 1부터 n까지 수행된다.
도 8은 도 1의 S120단계 또는 S170단계 또는 S195단계를 보다 자세히 설명하기 위한 흐름도이다.
도 1 내지 도 8을 참조하면, 전자장치(1000)는 수치표고모델을 생성하기 위한 대상지역이 결정되면, 대상지역의 수치표면모델(300a)에 적용할 모폴로지컬 필터의 써치 윈도우의 크기를 결정한다(S800). S800단계에서, 전자장치(1000)는 대상지역에 위치하는 가장 큰 객체의 최대 평면 넓이보다 크도록 설정할 수 있다.
전자장치(1000)는 초기 검사점의 위치로서 l=1, s=1을 정한다(S810). S810단계에서 정해지는 검사점의 위치는 예를 들어 도 7에 도시된 DSM의 최상단 좌측에 위치한 격자에 해당한다.
전자장치(1000)는 S810단계에서 정해진 검사점 (l, s)에 써치 윈도우의 중심점이 위치하도록 써치 윈도우의 중심점을 정렬한다(S820). 검사점 (l, s)=(1, 1)이고, 써치 윈도우의 중심점은 도 7의 경우 (4, 4)이므로, 전자장치(1000)는 (4,4)의 위치가 검사점 (1, 1)에 위치하도록 써치 윈도우를 정렬한다.
전자장치(1000)는 써치 윈도우의 중심점이 정렬되면, (l, s) 위치에서 최소 모폴로지컬 필터링을 적용하여 결과값을 산출하고, 산출된 결과값을 (l, s)=(1, 1)에 매핑저장한다(S830). 즉, S830단계에서, 전자장치(1000)는 써치 윈도우의 중심점이 (1, 1)인 경우, 최소 모폴로지컬 필터를 이용하여 써치 윈도우에 속하는 격자들의 높이를 수치표면모델로부터 확인하고, 확인된 격자들의 높이 중 가장 작은 높이를 결과값으로서 출력할 수 있다.
그리고, 전자장치(1000)는 s를 1 증가하여 검사점의 위치를 열 방향으로 1칸 이동하고(S840), s의 위치가 n에 도달하였는지 확인한다(S850).
s=n이 아니면(S850-No), 전자장치(1000)는 S820단계로 진입하여, 검사점의 위치를 (l, s)=(1, 2)로 이동시키고, 써치 윈도우의 중심점이 검사점의 위치에 오도록 정렬하여 S830단계를 수행한다.
이러한 동작의 반복에 의해 s=n에 도달하면(S850-Yes), 즉, l행에 있는 모든 격자들(모든 검사점들)에 대해 최소 모폴로지컬 필터링이 수행되면, 전자장치(1000)는 l을 1 증가하여 검사점의 위치를 행 방향으로 1줄 이동((l, s)=(2, 1))하고(S860), l의 위치가 m에 도달하였는지 확인한다(S870).
l=m이 아니면(S870-No), 전자장치(1000)는 S820단계로 진입하여 써치 윈도우의 중심점이 검사점 (l, s)=(2, 1)의 위치에 오도록 정렬하여 S830단계를 수행한다.
이러한 동작의 반복에 의해 l=m에 도달하고, s=n에 도달하면(S870-Yes), 즉, 수치표면모델에 있는 모든 격자들에 대해 최소 모폴로지컬 필터링이 수행되면, 전자장치(1000)는 최소 모폴로지컬 필터링의 결과값과 해당하는 중심점의 위치를 매핑하고, 고도를 이용하여 임시 수치표고모델(300b)을 생성한다(S880). S880단계에서 생성되는 임시 수치표고모델은 도
도 9는 도 1의 S130단계 또는 S160단계를 보다 자세히 설명하기 위한 흐름도이다.
도 1 내지 도 9를 참조하면, 전자장치(1000)는 임시 수치표고모델(300b)이 생성되면(S130), 초기 검사점의 위치로서 l=1, s=1을 정한다(S900).
전자장치(1000)는 S900단계에서 정해진 검사점 (l, s)에 써치 윈도우의 중심점이 위치하도록 써치 윈도우의 중심점을 정렬한다(S910).
전자장치(1000)는 써치 윈도우의 중심점이 정렬되면, (l, s) 위치에서 최대 모폴로지컬 필터링을 적용하여 결과값을 산출하고, 산출된 결과값을 (l, s)=(1, 1)에 매핑저장한다(S920). 즉, S920단계에서, 전자장치(1000)는 써치 윈도우의 중심점이 (1, 1)인 경우, 최소 모폴로지컬 필터를 이용하여 써치 윈도우에 속하는 격자들의 높이를 수치표면모델로부터 확인하고, 확인된 격자들의 높이 중 가장 큰 높이를 결과값으로서 출력할 수 있다.
그리고, 전자장치(1000)는 s를 1 증가하여 검사점의 위치를 열 방향으로 1칸 이동하고(S930), s의 위치가 n에 도달하였는지 확인한다(S940).
s=n이 아니면(S940-No), 전자장치(1000)는 S910단계로 진입하여, 검사점의 위치를 (l, s)=(1, 2)로 이동시키고, 써치 윈도우의 중심점이 검사점의 위치에 오도록 정렬하여 S920단계를 수행한다.
이러한 동작의 반복에 의해 s=n에 도달하면(S940-Yes), 즉, l행에 있는 모든 격자들(모든 검사점들)에 대해 최대 모폴로지컬 필터링이 수행되면, 전자장치(1000)는 l을 1 증가하여 검사점의 위치를 행 방향으로 1줄 이동((l, s)=(2, 1))하고(S950, l의 위치가 m에 도달하였는지 확인한다(S960).
l=m이 아니면(S960-No), 전자장치(1000)는 S910단계로 진입하여 써치 윈도우의 중심점이 검사점 (l, s)=(2, 1)의 위치에 오도록 정렬하여 S920단계를 수행한다.
이러한 동작의 반복에 의해 l=m에 도달하고, s=n에 도달하면(S960-Yes), 즉, 임시 수치표고모델에 있는 모든 격자들에 대해 최대 모폴로지컬 필터링이 수행되면, 전자장치(1000)는 최대 모폴로지컬 필터링의 결과값과 해당하는 중심점의 위치를 매핑하고, 고도를 이용하여 최종 수치표고모델(300c)을 생성한다(S970).
도 10은 본 발명의 일 실시 예에 따른 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 전자장치(1000)를 도시한 블록도이다.
도 10을 참조하면, 본 발명의 실시 예에 따른 전자장치(1000)는 입력부(1010), 최소 모폴로지컬 필터링부(1020), 최대 모폴로지컬 필터링부(1030) 및 저장부(1040)를 포함한다.
입력부(1010)는 수치표고모델로 생성하기 위한 대상지역의 수치표면모델을 입력받는다. 입력부(1010)는 대상지역이 언덕이나 협곡과 같은 경사 지역을 포함하면, 입력된 대상지역의 수치표면모델을 최소 모폴로지컬 필터링부(1020)로 전달한다.
대상지역이 평지인지 또는 경사 지역을 포함하는지 또는 환풍구를 포함하는지의 여부는 입력부(1010)가 프로세서를 통해 직접 판단하거나 별도의 프로세서에 의해 판단될 수 있다.
최소 모폴로지컬 필터링부(1020)는 입력부(1010)로부터 입력된 수치표면모델에 모든 검사점마다 최소 모폴로지컬 필터링을 적용하여 임시 수치표고모델(300b)을 생성하고, 생성된 임시 수치표고모델을 최대 모폴로지컬 필터링부(1030)로 전달할 수 있다.
최소 모폴로지컬 필터링부(1020)는 최소 모폴로지컬 필터를 이용하여, 수치표면모델에 써치 윈도우를 적용하고, 써치 윈도우 내에서 높이가 가장 낮은 격자의 높이를 써치 윈도우의 중심점의 최소 지면 높이로 정한 후, 써치 윈도우의 중심점과 정해진 최소 지면 높이를 이용하여 임시 수치표고모델을 생성할 수 있다. 최소 모폴로지컬 필터링부(1020)의 동작은 도 8을 참조하여 상세히 설명하고 있다.
최대 모폴로지컬 필터링부(1030)는 임시 수치표고모델(300b)에서 발생하는 왜곡을 보상하기 위하여, 임시 수치표고모델(300b)에 최대 모폴로지컬 필터링을 적용하여 수치표고모델(300c)을 생성할 수 있다.
최대 모폴로지컬 필터링부(1030)는 최대 모폴로지컬 필터를 이용하여, 임시 수치표고모델(300b)에 써치 윈도우를 적용하고, 써치 윈도우 내에서 높이가 가장 높은 격자의 높이를 써치 윈도우의 중심점의 최종 지면 높이로 정한 후, 써치 윈도우의 중심점과 정해진 최종 지면 높이를 취합하여 수치표고모델(300c)을 생성할 수 있다.
이로써, 임시 수치표고모델(300b)에서 발생한 왜곡이 보상된, 경사 지역을 포함하는 대상지역에 해당하는 수치표고모델(300c)이 생성되어 저장부(1040)에 저장된다. 최대 모폴로지컬 필터링부(1030)의 동작은 도 9를 참조하여 상세히 설명하고 있다.
또한, 입력부(1010)는 입력된 수치표면모델에 해당하는 대상지역이 지하에 형성된 환풍구를 포함하면, 수치표면모델을 최대 모폴로지컬 필터링부(1030)로 전달한다.
최대 모폴로지컬 필터링부(1030)는 입력부(1010)로부터 입력된 환풍구를 포함하는 수치표면모델에 모든 검사점마다 최대 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하고, 생성된 임시 수치표고모델을 최소 모폴로지컬 필터링부(1020)로 전달할 수 있다.
최소 모폴로지컬 필터링부(1020)는 최대 모폴로지컬 필터링부(1030)로부터 입력되는 임시 수치표고모델에 최소 모폴로지컬 필터링을 적용하여 환풍구를 포함하는 대상지역에 해당하는 수치표고모델을 생성할 수 있다.
또한, 입력부(1010)는 입력받은 수치표면모델에 해당하는 대상지역이 평지인 경우, 입력된 대상지역의 수치표면모델을 최소 모폴로지컬 필터링부(1020)로 전달한다.
최소 모폴로지컬 필터링부(1020)는 입력부(1010)로부터 입력된 평지에 해당하는 수치표면모델에 모든 검사점마다 최소 모폴로지컬 필터링을 적용하여 수치표고모델을 생성할 수 있다.
도 11은 본 발명의 다른 실시 예에 따른 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 전자장치(1100)를 도시한 블록도이다.
도 11을 참조하면, 본 발명의 실시 예에 따른 전자장치(1100)는 버스(1110), 사용자 인터페이스부(User Interface, UI부)(1120), 표시부(1130), DSM 저장부(1140), DEM 저장부(1150), 제어부(1160), 최소 모폴로지컬 필터링부(1170) 및 최대 모폴로지컬 필터링부(1180)를 포함한다.
버스(1110)는 각 구성요소들(1120~1180)을 서로 연결하고, 제어 메시지 및/또는 데이터와 같은 각종 신호를 전달하는 회로를 포함할 수 있다.
UI부(1120)는 사용자로부터 입력된 명령 또는 데이터를 전자 장치(1100)의 다른 구성요소(들)에 전달할 수 있는 인터페이스의 역할을 할 수 있다. 예를 들어, UI부(1120)는 사용자로부터 대상지역의 수치표면모델로부터 수치표고모델을 생성하는 필터링 모드를 선택하기 위한 모드 선택 명령을 수신하고, 대상지역 지정 명령을 수신할 수 있다. 필터링 모드는, 자동 모드, 경사 모드, 환풍구 모드 및 평지 모드를 포함한다.
자동 모드는 제어부(1160)가 대상지역의 고도를 기반으로 경사 지역인지, 환풍구 지역인지 또는 평지인지 자동으로 판별하여, 최소 모폴로지컬 필터링과 최대 모폴로지컬 필터링을 적응적으로 적용하는 모드이다.
경사 모드는 최소 모폴로지컬 필터링과 최대 모폴로지컬 필터링을 순차적으로 적용하는 모드이고, 환풍구 모드는 최대 모폴로지컬 필터링과 최소 모폴로지컬 필터링을 순차적으로 적용하는 모드이고, 평지 모드는 최소 모폴로지컬 필터링을 적용하는 모드이다.
경사 모드, 환풍구 모드 및 평지 모드는 수동 모드로서, 사용자가 직접 대상지역에 적용할 모드, 즉, 필터링 순서를 선택하는 모드이다. 이 때, 사용자는 모드와 모드를 적용할 대상지역을 함께 지정할 수도 있고, 대상지역이 두 개 이상의 지역(예를 들어, 경사 지역+평지, 경사 지역+환풍구, 평지+환풍구, 경사 지역+환풍구+평지 등)을 포함하는 경우, 모드만 지정할 수도 있다.
표시부(1130)는 필터링할 대상지역, 대상지역의 수치표면모델, 임시 수치표고모델, 최종으로 생성되는 수치표고모델을 표시한다.
DSM 저장부(1140)는 수치표고모델을 생성하기 위한 대상지역의 수치표면모델을 저장한다.
DEM 저장부(1150)는 최소 모폴로지컬 필터링부(1170) 또는 최대 모폴로지컬 필터링부(1180)에 의해 생성되는 수치표고모델을 저장한다.
제어부(1160)는 전자장치(1100)의 전반적인 동작을 제어하며, 프로세서로 구현될 수 있다.
자동 모드가 선택되면, 제어부(1160)는 DSM 저장부(1140)에 저장된 대상지역의 수치표면모델을 읽어와 이웃하는 격자들 간의 고도차를 계산하고, 계산된 고도차와 임계값을 비교하여 대상지역에 적용할 모드를 판단할 수 있다.
예를 들어, 계산된 고도차가 제1임계값보다 크며, 제1임계값보다 큰 고도차가 대상지역에서 N개 이상 연속 또는 이격 산출되면, 제어부(1160)는 대상지역에 경사 모드를 적용할 수 있다. 또한, 계산된 고도차가 제2임계값보다 크며, 제2임계값보다 큰 고도차가 대상지역에서 M개 이상 산출되면, 제어부(1160)는 대상지역에 환풍구 모드를 적용할 수 있다. 또한, 계산된 고도차가 제3임계값보다 작으며, 제3임계값보다 작은 고도차가 대상지역에서 L개 이상 산출되면, 제어부(1160)는 대상지역에 평지 모드를 적용할 수 있다.
여기서, 제1임계값은 제2임계값보다 작으며, 경사 지역은 완만한 경사도 포함하므로, 제1임계값은 예를 들어, 10도보다 작게 설정될 수 있다. 또한, 환풍구는 대부분 평지 기준으로 직각에 가깝게 형성되므로 80도보다 크게 설정될 수 있다. 또한, 평지는 고도차가 거의 없으므로 제3임계값은 1도로 설정될 수 있다. 1도와 10도와 80도는 발명의 이해를 위해 기재한 값으로서, 제1 내지 제3임계값은 이에 한정되지 않으며, 사용자가 대상지역의 특성을 고려하여 변경할 수 있다.
수동 모드로서 경사 모드가 선택되면, 제어부(1160)는 DSM 저장부(1140)에 저장된 수치표면모델이 최소 모폴로지컬 필터링부(1170)로 입력되도록 하고, 최소 모폴로지컬 필터링과 최대 모폴로지컬 필터링이 순차적으로 수행되도록 처리한다.
환풍구 모드가 선택되면, 제어부(1160)는 저장된 수치표면모델이 최대 모폴로지컬 필터링부(1180)로 입력되도록 하고, 최대 모폴로지컬 필터링과 최소 모폴로지컬 필터링이 순차적으로 수행되도록 처리한다.
평지 모드가 선택되면, 제어부(1160)는 저장된 수치표면모델이 최소 모폴로지컬 필터링부(1170)로 입력되도록 하여, 최소 모폴로지컬 필터링이 수행되도록 처리한다.
경사 모드가 선택되거나, 자동 모드에 따라 제어부(1160)가 판단한 결과에 의해 DSM 저장부(1140)에 저장된 대상지역의 수치표면모델이 입력되면, 최소 모폴로지컬 필터링부(1170)는 입력된 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하고, 생성된 임시 수치표고모델을 최대 모폴로지컬 필터링부(1180)로 전달할 수 있다.
최대 모폴로지컬 필터링부(1180)는 임시 수치표고모델에 최대 모폴로지컬 필터링을 적용하여 경사 지역을 포함하는 대상지역에 대한 수치표고모델을 생성할 수 있다.
또한, 환풍구 모드가 선택되거나, 자동 모드에 따라 제어부(1160)가 판단한 결과에 의해 DSM 저장부(1140)에 저장된 대상지역의 수치표면모델이 입력되면, 최대 모폴로지컬 필터링부(1180)는 입력된 수치표면모델에 최대 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하고, 생성된 임시 수치표고모델을 최소 모폴로지컬 필터링부(1170)로 전달할 수 있다.
최소 모폴로지컬 필터링부(1170)는 임시 수치표고모델에 최소 모폴로지컬 필터링을 적용하여 환풍구 지역을 포함하는 대상지역에 대한 수치표고모델을 생성할 수 있다.
또한, 평지 모드가 선택되거나, 자동 모드에 따라 제어부(1160)가 판단한 결과에 의해 DSM 저장부(1140)에 저장된 대상지역의 수치표면모델이 입력되면, 최소 모폴로지컬 필터링부(1170)는 입력된 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 평지에 대한 수치표고모델을 생성할 수 있다.
도 12는 본 발명의 실시 예에 따라 최소 모폴로지컬 필터링 및 최대 모폴로지컬 필터링을 적용한 결과를 보여주는 그래프이다.
도 12를 참조하면, 샌프란시스코 수치표면모델에 본 발명을 적용한 결과, 최소 모폴로지컬 필터링과 최대 모폴로지컬 필터링을 순차적으로 적용하여 생성한 수치표고모델은 이상적인 수치표고모델에 비해 거의 왜곡되지 않은 결과를 제공한다.
도 13은 본 발명의 실시 예에 따라 환풍구 지역에 최대 모폴로지컬 필터링 및 최소 모폴로지컬 필터링을 적용한 결과를 설명하기 위한 그래프이다.
도 13의 (a)는 대상지역의 수치표면모델, 즉, 원본 DSM으로서, 빨강색 원은 환풍구가 위치하는 지역을 의미한다. 도 13의 (b)는 (a)의 수치표면모델을 3D로 보여준 것이고, (c)는 (b)에 도시된 수치표면모델을 아래에서 본 것이고, (d)는 (a)의 수치표면모델을 최소 모폴로지컬 필터링한 결과로서 환풍구가 남아 있음을 알 수 있다. 도 13의 (e)는 (a)의 수치표면모델에 최대 모폴로지컬 필터링한 결과이고, (f)는 (e)의 결과를 3D로 보여준 것으로서, 환풍구가 제거되었음을 확인할 수 있다.
한편, 상술한 본 발명의 다른 실시 예에서, 제어부(1160)는 대상지역이 경사 지역, 환풍구 지역 및 평지 지역 중 적어도 두 개를 포함하면, 대상지역을 해당하는 지역으로 구획하고, 구획된 대상지역의 수치표면모델이 적응적으로 최소 모폴로지컬 필터링부(1170) 또는 최대 모폴로지컬 필터링부(1180)로 먼저 입력되도록 처리할 수도 있다.
즉, 제어부(1160)는, 대상지역 중 경사 지역에 해당하는 부분은 최소 모폴로지컬 필터링부(1170)에서 먼저 필터링된 후 최대 모폴로지컬 필터링부(1180)에서 필터링되도록 처리한다. 또한, 제어부(1160)는 대상지역 중 환풍구 지역에 해당하는 부분은 최대 모폴로지컬 필터링부(1180)에서 먼저 필터링된 후 최소 모폴로지컬 필터링부(1170)에서 필터링되도록 처리한다.
한편, 본 발명에 따른 전자장치의 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 방법은 이를 구현하기 위한 명령어들의 프로그램이 유형적으로 구현됨으로써, 컴퓨터를 통해 판독될 수 있는 기록매체에 포함되어 제공될 수도 있음은 통상의 기술자가 쉽게 이해할 수 있다.
따라서, 본 발명은 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 방법을 구현하기 위하여 상기 전자장치를 제어하는 컴퓨터 상에서 수행되는 컴퓨터 판독 가능한 기록매체에 저장된 프로그램을 함께 제공한다.
한편, 이상으로 본 발명의 기술적 사상을 예시하기 위한 바람직한 실시 예와 관련하여 설명하고 도시하였지만, 본 발명은 이와 같이 도시되고 설명된 그대로의 구성 및 작용에만 국한되는 것이 아니며, 기술적 사상의 범주를 일탈함이 없이 본 발명에 대해 다수의 변경 및 수정 가능함을 당업자들은 잘 이해할 수 있을 것이다. 따라서, 그러한 모든 적절한 변경 및 수정과 균등물들도 본 발명의 범위에 속하는 것으로 간주하여야 할 것이다. 따라서, 본 발명의 진정한 기술적 보호 범위는 첨부된 등록청구범위의 기술적 사상에 의해 정해져야 할 것이다.
1000, 1100: 전자장치
1010: 입력부
1020, 1170: 최소 모폴로지컬 필터링부
1030, 1180: 최대 모폴로지컬 필터링부
1040: 저장부

Claims (10)

  1. 전자장치의 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법에 있어서,
    (A) 상기 전자장치가, 상기 대상지역이 경사 지역을 포함하는 경우, 격자 기반의 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하는 단계;
    (B) 상기 전자장치가, 상기 (A) 단계에서 생성된 임시 수치표고모델 중 경사 지역에 발생하는 고도 왜곡이 보상되도록 상기 (A) 단계에서 생성된 임시 수치표고모델에 최대 모폴로지컬 필터링을 적용하여 상기 경사 지역을 포함하는 대상지역에 해당하는 수치표고모델을 생성하는 단계;
    (C) 상기 전자장치가, 상기 대상지역이 환풍구를 포함하는 경우, 상기 수치표면모델에 최대 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하는 단계; 및
    (D) 상기 전자장치가, 상기 (C) 단계에서 생성된 임시 수치표고모델에 최소 모폴로지컬 필터링을 적용하여 상기 환풍구를 포함하는 대상지역에 해당하는 수치표고모델을 생성하는 단계;를 포함하고,
    상기 (A) 단계는,
    상기 전자장치가, 상기 격자 기반의 수치표면모델에 사전에 정해진 크기의 써치 윈도우를 적용하고, 상기 써치 윈도우 내에서 높이가 가장 낮은 격자의 높이를 상기 써치 윈도우의 중심점의 최소 지면 높이로 정한 후, 상기 써치 윈도우의 중심점과 상기 정해진 최소 지면 높이를 이용하여 상기 임시 수치표고모델을 생성하며,
    상기 최소 지면 높이로서 상기 대상지역에 위치하는 볼록한 형태의 객체들 중 가장 큰 객체의 높이가 정해지는 것을 방지하기 위하여, 상기 써치 윈도우의 크기는 상기 가장 큰 객체의 최대 평면 넓이보다 크도록 정해지는 것을 특징으로 하는 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 방법.
  2. 삭제
  3. 제1항에 있어서,
    상기 (B) 단계는,
    상기 전자장치가, 상기 임시 수치표고모델에 써치 윈도우를 적용하고, 상기 써치 윈도우 내에서 높이가 가장 높은 격자의 높이를 상기 써치 윈도우의 중심점의 최종 지면 높이로 정한 후, 상기 써치 윈도우의 중심점과 상기 정해진 최종 지면 높이를 이용하여 상기 경사 지역을 포함하는 대상지역에 해당하는 수치표고모델을 생성하는 것을 특징으로 하는 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 방법.
  4. 제1항에 있어서,
    상기 (A) 단계 이전 및 상기 (B) 단계 이후 중 어느 하나에,
    (E) 상기 전자장치가, 상기 대상지역이 평지인 경우, 상기 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 상기 평지에 해당하는 수치표고모델을 생성하는 단계;를 더 포함하는 것을 특징으로 하는 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 방법.
  5. 삭제
  6. 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 전자장치에 있어서,
    상기 대상지역이 경사 지역을 포함하는 경우, 격자 기반의 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하는 최소 필터링부; 및
    상기 최소 필터링부에서 생성된 임시 수치표고모델 중 경사 지역에 발생하는 고도 왜곡이 보상되도록 상기 임시 수치표고모델에 최대 모폴로지컬 필터링을 적용하여 상기 경사 지역을 포함하는 대상지역에 해당하는 수치표고모델을 생성하는 최대 필터링부;를 포함하고,
    상기 대상지역이 환풍구를 포함하는 경우, 상기 최대 필터링부는, 상기 수치표면모델에 최대 모폴로지컬 필터링을 적용하여 임시 수치표고모델을 생성하고, 상기 최소 필터링부는, 상기 최대 필터링부에서 생성된 임시 수치표고모델에 최소 모폴로지컬 필터링을 적용하여 상기 환풍구를 포함하는 대상지역에 해당하는 수치표고모델을 생성하며,
    상기 최소 필터링부는,
    상기 격자 기반의 수치표면모델에 사전에 정해진 크기의 써치 윈도우를 적용하고, 상기 써치 윈도우 내에서 높이가 가장 낮은 격자의 높이를 상기 써치 윈도우의 중심점의 최소 지면 높이로 정한 후, 상기 써치 윈도우의 중심점과 상기 정해진 최소 지면 높이를 이용하여 임시 수치표고모델을 생성하며,
    여기서, 상기 최소 지면 높이로서 상기 대상지역에 위치하는 볼록한 형태의 객체들 중 가장 큰 객체의 높이가 정해지는 것을 방지하기 위하여, 상기 써치 윈도우의 크기는 상기 가장 큰 객체의 최대 평면 넓이보다 크도록 정해지는 것을 특징으로 하는 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 전자장치.
  7. 삭제
  8. 제6항에 있어서,
    상기 최대 필터링부는,
    상기 임시 수치표고모델에 써치 윈도우를 적용하고, 상기 써치 윈도우 내에서 높이가 가장 높은 격자의 높이를 상기 써치 윈도우의 중심점의 최종 지면 높이로 정한 후, 상기 써치 윈도우의 중심점과 상기 정해진 최종 지면 높이를 이용하여 상기 경사 지역을 포함하는 대상지역에 해당하는 수치표고모델을 생성하는 것을 특징으로 하는 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 전자장치.
  9. 제6항에 있어서,
    상기 대상지역이 평지인 경우,
    상기 최소 필터링부는, 상기 수치표면모델에 최소 모폴로지컬 필터링을 적용하여 상기 평지에 해당하는 수치표고모델을 생성하는 것을 특징으로 하는 모폴로지컬 필터링을 이용하여 수치표면모델로부터 수치표고모델을 생성하기 위한 전자장치.
  10. 삭제
KR1020180170868A 2018-12-27 2018-12-27 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치 KR101973726B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020180170868A KR101973726B1 (ko) 2018-12-27 2018-12-27 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020180170868A KR101973726B1 (ko) 2018-12-27 2018-12-27 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치

Publications (1)

Publication Number Publication Date
KR101973726B1 true KR101973726B1 (ko) 2019-04-30

Family

ID=66285655

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020180170868A KR101973726B1 (ko) 2018-12-27 2018-12-27 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치

Country Status (1)

Country Link
KR (1) KR101973726B1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112184900A (zh) * 2019-07-04 2021-01-05 北京四维图新科技股份有限公司 高程数据的确定方法、装置及存储介质

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20070026814A (ko) * 2004-07-15 2007-03-08 해리스 코포레이션 지형적 지점으로부터 3-차원 등록에 대한 맨 지구 디지털고도 모델 추출
KR100927642B1 (ko) * 2009-05-14 2009-11-20 한국 한의학 연구원 영상 분할 히스토그램 분석을 이용한 얼굴 윤곽선 검출 방법
KR101006729B1 (ko) * 2010-07-23 2011-01-10 (주)동광지엔티 수치 표고모델 제작 방법 및 시스템
KR101008394B1 (ko) 2010-10-12 2011-01-14 (주)동광지엔티 수치표면모형의 경계선을 이용한 수치표고모형의 건물영역 평탄화 장치
KR101193115B1 (ko) * 2011-10-07 2012-10-19 한국항공우주산업 주식회사 3ｄ 전자 지도 시스템
KR20140049100A (ko) * 2012-10-11 2014-04-25 (주)아이비리더스 항공정보 현시 시스템 및 그 방법
KR101389126B1 (ko) * 2013-08-07 2014-04-25 (주) 지오씨엔아이 라이다 자료와 정사영상을 이용한 제방법선 매핑 방법
KR20140142381A (ko) * 2013-05-28 2014-12-12 삼성테크윈 주식회사 단일영상 내의 안개 제거 방법 및 장치
KR101781009B1 (ko) * 2016-08-31 2017-10-23 노틸러스효성 주식회사 낡은 지폐 식별 방법
KR101795547B1 (ko) * 2011-04-01 2017-11-09 에스케이 주식회사 디지털 고도모델 생성방법 및 장치
KR20180131932A (ko) * 2017-06-01 2018-12-11 충남대학교산학협력단 드론과 공간정보를 이용한 하천지형정보 생성 방법

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20070026814A (ko) * 2004-07-15 2007-03-08 해리스 코포레이션 지형적 지점으로부터 3-차원 등록에 대한 맨 지구 디지털고도 모델 추출
KR100927642B1 (ko) * 2009-05-14 2009-11-20 한국 한의학 연구원 영상 분할 히스토그램 분석을 이용한 얼굴 윤곽선 검출 방법
KR101006729B1 (ko) * 2010-07-23 2011-01-10 (주)동광지엔티 수치 표고모델 제작 방법 및 시스템
KR101008394B1 (ko) 2010-10-12 2011-01-14 (주)동광지엔티 수치표면모형의 경계선을 이용한 수치표고모형의 건물영역 평탄화 장치
KR101795547B1 (ko) * 2011-04-01 2017-11-09 에스케이 주식회사 디지털 고도모델 생성방법 및 장치
KR101193115B1 (ko) * 2011-10-07 2012-10-19 한국항공우주산업 주식회사 3ｄ 전자 지도 시스템
KR20140049100A (ko) * 2012-10-11 2014-04-25 (주)아이비리더스 항공정보 현시 시스템 및 그 방법
KR20140142381A (ko) * 2013-05-28 2014-12-12 삼성테크윈 주식회사 단일영상 내의 안개 제거 방법 및 장치
KR101389126B1 (ko) * 2013-08-07 2014-04-25 (주) 지오씨엔아이 라이다 자료와 정사영상을 이용한 제방법선 매핑 방법
KR101781009B1 (ko) * 2016-08-31 2017-10-23 노틸러스효성 주식회사 낡은 지폐 식별 방법
KR20180131932A (ko) * 2017-06-01 2018-12-11 충남대학교산학협력단 드론과 공간정보를 이용한 하천지형정보 생성 방법

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Junhee Youn, et al., "AUTOMATIC EXTRACTING DEM FROM DSM WITH CONSECUTIVE MORPHOLOGICAL FILTERING", GLOBAL JOURNAL OF ENGINEERING SCIENCE AND RESEARCHES(2018.11.) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112184900A (zh) * 2019-07-04 2021-01-05 北京四维图新科技股份有限公司 高程数据的确定方法、装置及存储介质
CN112184900B (zh) * 2019-07-04 2024-03-19 北京四维图新科技股份有限公司 高程数据的确定方法、装置及存储介质

Similar Documents

Publication Publication Date Title
US11144758B2 (en) System and method for object detection and classification in aerial imagery
Sithole et al. Filtering of laser altimetry data using a slope adaptive filter
Elmqvist Ground surface estimation from airborne laser scanner data using active shape models
KR102415767B1 (ko) 지물 높이별 컬러 화상 생성 장치 및 지물 높이별 컬러 화상 생성 프로그램
KR102415768B1 (ko) 지물 지반 높이별 컬러 화상 생성 장치 및 지물 높이별 컬러 화상 생성 프로그램
KR100967838B1 (ko) 항공 라이다 정보와 디지털 항공 사진 정보를 이용하여 3차원 지리 정보를 생성하는 방법 및 그 시스템
US10650583B2 (en) Image processing device, image processing method, and image processing program
JP6534296B2 (ja) ３次元モデル生成装置、３次元モデル生成方法、及びプログラム
JP6007312B2 (ja) 画像処理装置、画像処理方法及び画像処理プログラム
Masaharu et al. A filtering method of airborne laser scanner data for complex terrain
KR101973726B1 (ko) 모폴로지컬 필터링을 이용하여 대상지역의 수치표면모델로부터 수치표고모델을 생성하기 위한 방법 및 전자장치
JP7170575B2 (ja) 地すべり面推定装置、及び地すべり面推定方法
KR102567800B1 (ko) 드론을 이용한 3차원 지도 제작 방법
CN110399807B (zh) 检测地面障碍物的方法、装置、可读存储介质及电子设备
KR102491674B1 (ko) 가상시점 영상을 생성하는 방법 및 장치
CN109255809A (zh) 一种光场图像深度估计方法及装置
KR101021013B1 (ko) 경계강화필터링과 표고정보를 이용하여 3차원 지리정보를 생성하는 시스템
KR102014897B1 (ko) 항공측량을 위한 3차원 경로 비행 방법 및 장치
CN116438581A (zh) 三维点群高密度化装置、三维点群高密度化方法及程序
WO2018159079A1 (ja) 地形情報処理装置、地形情報処理方法、およびプログラム
JP2020021368A (ja) 画像解析システム、画像解析方法及び画像解析プログラム
KR102564218B1 (ko) 드론 영상과 인공지능을 이용한 침수 범람 예측 시스템
Ismail et al. DEM derived from photogrammetric generated DSM using morphological filter
KR102550233B1 (ko) 수치건물지면모델 생성 방법 및 장치, 컴퓨터 판독 가능한 기록 매체 및 컴퓨터 프로그램
KR101114904B1 (ko) 도화원도와 항공 레이저 측량 데이터를 이용한 도시공간정보 구축 시스템 및 그 방법

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant