A Study on Automated Lineament Extraction with Respect to Spatial Resolution of Digital Elevation Model

수치표고모형 공간해상도에 따른 선구조 자동 추출 연구

Abstract

The lineament is a linear or curved terrain element to discriminate adjacent geological structures in each other. It has been widely used for analysis of geology, mineral exploration, natural disasters, and earthquake, etc. In the past, the lineament has been extracted using cartographic map or field survey. However, it is possible to extract more efficiently the lineament for a very wide area thanks to development of remote sensing technique. Remotely sensed observation by aircraft, satellite, or digital elevation model (DEM) has been used for visual recognition for manual lineament extraction. Automatic approaches using computer science have been proposed to extract lineament more objectively. In this study, we evaluate the characteristics of lineament which is automatically extracted with respect to difference of spatial resolution of DEM. We utilized two types of DEM: one is Shuttle Radar Topography Mission (SRTM) with spatial resolution of about 90 m (3 arc sec), and the other is the latest world DEM of TerraSAR-X add-on for Global DEM with 12 m spatial resolution. In addition, a global DEM was resampled to produce a DEM with a spatial resolution of 30 m (1 arc sec). The shaded relief map was constructed considering various sun elevation and solar azimuth angle. In order to extract lineament automatically, we used the LINE module in PCI Geomatica software. We found that predominant direction of the extracted lineament is about $N15-25^{\circ}E$ (NNE), regardless of spatial resolution of DEM. However, more fine and detailed lineament were extracted using higher spatial resolution of DEM. The result shows that the lineament density is proportional to the spatial resolution of DEM. Thus, the DEM with appropriate spatial resolution should be selected according to the purpose of the study.

선구조는 인접한 지형구조 사이의 구분이 명확한 직선 또는 만곡의 지형요소로서 일반지질, 광물탐사, 자연재해, 지구조 분석 등에 널리 활용된다. 과거에는 현장조사 혹은 지도를 이용하여 선구조를 추출하였으나, 현재는 원격탐사 기술의 발달로 인하여 넓은 지역에 대한 선구조를 효율적으로 추출할 수 있게 되었다. 선구조 추출을 위해서 항공기 혹은 인공위성 원격탐사 영상 혹은 지형표고모형의 육안판독 방법 이외에 보다 객관적인 결과 도출을 위한 자동화 방법이 개발되었다. 본 연구에서는 지형표고모형의 공간해상도 차이에 따라 자동 추출된 선구조의 특성을 분석, 평가하고자 한다. 연구 지역은 대한민국 경상도 언양, 모량 주변 지역으로 약 90 m (3초) 공간해상도를 갖는 Shuttle Radar Topography Mission(SRTM) Digital Elevation Model(DEM)과 12 m 공간해상도를 갖는 TerraSAR-X add-on for Global Digital Elevation Measurement(Global DEM) 자료를 사용하였다. 또한 Global DEM을 재배열(resampling)하여 30 m (1초) 간격의 공간해상도를 갖는 지형고도모형을 제작하였다. 다양한 각도의 태양고도와 태양방위각을 고려한 음영기복도를 제작하였으며, 선구조 자동 추출을 위해 PCI Geomatica 소프트웨어의 LINE 모듈을 이용하였다. 수치표고모형 공간해상도에 상관없이 선구조의 최빈값은 $N15-25^{\circ}E$로 북북동(NNE)의 방향성을 보였다. 그러나 공간해상도가 좋을수록 보다 많고 세밀한 선구조가 추출되었다. 본 연구 결과를 통해 선구조 밀도는 수치표고모형의 공간해상도에 비례함을 알 수 있었으며 연구 목적에 따라 적절한 공간해상도를 갖는 수치표고모형이 선택되어야 함을 알 수 있었다.

 

1. 서론

선구조는 인접한 표면 및 천부의 지형요소와 지형구조 사이의 구분이 명확한 직선 또는 만곡의 지형요소로서 지질, 탐사, 자연재해, 도로건설 등 다양한 분야에 널리 활용되고 있다. 특히 선구조는 지구조 분석, 단층면 해석 등 지질학 연구 및 지형 활용 분야에 있어 유용한 도구로 인식되어 왔다(Hobbs, 1904; 1912). 위성영상을 통해 추출 가능한 선구조는 크게 단층, 습곡 또는 파쇄대에서 주로 나타나는 지질학적 선구조와 기복, 색조의 대비로 생성되는 지형학적 선구조가 있다(Sabins, 1996). 일반적으로 선구조는 급격한 지형학적 변화, 식생의 종류와 높이의 변화, 식생의 정렬, 토양의 색상변화, 지표면의 침강 또는 융기로 인한 계곡 등의 구조에서 확인가능하다. 이들은 대부분 단층 또는 절리 등 지구조 운동으로 인해 생성된 지질학적 구조로 인해 만들어진 결과물이다(O’ Leary et al., 1976). 이러한 선구조의 분석을 위해 과거에는 현장조사 혹은 지형도에 기반한 방법이 이용되었으나, 현재는 항공기와 인공위성 등의 발달로 인하여 넓은 지역에 대한 원격탐사 자료를 쉽게 취득할 수 있어 보다 효과적으로 넓은 범위의 선구조를 추출할 수 있다.

원격탐사 자료에서 나타나는 선구조는 주로 영상에 나타나는 화소의 음영 차이를 이용하여 확인 가능하며, 육안 분석을 통한 전통적인 방법과 이웃한 화소간 음영 차이 분석 또는 방향성 추출 알고리즘을 이용한 선구조 추출 방법이 사용되어 왔다. 육안 분석을 통한 선구조 분석은 주관적인 해석이 포함될 수 있으며 넓은 지역에 대한 육안 분석을 수행하는데 한계점이 나타나, 자동화된 선구조 추출 연구가 주목을 받게 되었으며 다양한 알고리즘의 한계점 도출과 보완에 대한 다양한 연구가 수행되고 있다. 국내의 경우, Kim et al.(1996; 1997)는 경상분지 의성지역, 옥천대 남서경계부 등에서 원격탐사자료와 수치표고모형(DEM)을 이용하여 선구조의 방향성을 분석하였다. 해당 연구에서는 선구조 추출시 한계점으로 지적되었던 관측 방향에 따른 선택적 증감 효과를 보정하기 위해 음영기복도를 사용하였다. 또한 Hough 변환 등 기존의 선구조 추출 알고리즘이 국내에 적합하지 않다는 한계가 지적되어 Won et al.(1998)는 Dynamic Segment Tracing Algorithm(DSTA), Automatic Lineament Extraction by Hierarchical Hough Transform(ALEHHT), Automatic Lineament Extraction by Generalized Hough Transform(ALEGHT) 알고리즘을 사용하여 국내에서의 선구조 추출 기술을 향상시켰다. Hong et al.(2000)는 영상레이더에서 나타나는 지질학적 선구조의 특성과 방향편차에 대한 분석 연구가 수행된 바 있다. Lee and Moon(2002)는 Landsat TM, JERS-1 Synthetic Aperture Radr(SR)와 지형고도모형 자료를 바탕으로 Hough 변환과 다중 해상도 방법을 사용하여 다양한 규모에서 선구조를 추출하였고, 다중 해상도 방법이 국내에서 보다 효과적임을 보고하였다.

해외에서는 보다 넓은 지역에 대한 알고리즘 및 선구조 추출 프로그램에 기반한 선구조 추출 연구가 수행되어 왔다. Hung et al.(2005)은 ASTER와 Landsat ETM 자료를 이용하여 선구조를 추출하였으며 PCI Geomatica의 LINE 모듈과 ArcGIS^TM 등 다양한 소프트웨어가 활용되었다. 또한, 방향성과 선구조 밀도 분석 등 통계분석을 통해 결과를 점검하는 일련의 과정을 통해 자동적이고 디지털화된 선구조 분석법을 개발하였다.

원격탐사 자료에 기반한 초기 선구조 추출 연구에서는 큰 선구조의 분포나 방향성 검증이 주요 관심사였으며 주로 중/저 공간해상도의 영상자료를 이용하였다. 그러나 이와 같은 낮은 공간해상도의 원격탐사 영상으로부터 추출된 선구조는 정확성과 정밀성을 요구하는 지질공학이나 광상 조사에는 이용하기 부적합하다는 한계가 지적되었다(Kocal et al., 2007). 따라서 1 m의 고해상 공간분해능을 가지는 다중 스펙트럼 IKONOSPAN-sharpened 이미지를 바탕으로 한 선구조 추출연구가 제안되었다. Abdullah et al.(2010)은 지형고도모형 자료에서 선구조를 추출할 때 보다 명확한 결과 추출을 위해 다양한 방향의 태양방위각을 고려한 음영기복도를 사용하는 방법을 제시되었다. Kusák and Krbcová(2017)은 SRTM 자료를 바탕으로 ArcGIS^TM와 PCI Geomatica의 LINE 모듈을 이용한 넓은 지역에서의 선형 추출 방법을 제시하였으며 특히, RADI, GTHR, LTHR, FTHR, ATHR, DTHR 등 6가지 변수를 조정하여 이에 따른 선구조의 경향성을 파악하고자 하였다.

본 연구에서는 소프트웨어에 기반한 자동화된 선구조 추출 연구의 국내 지역에 대한 적용 가능성을 평가하는데 그 목적이 있다. 또한 다양한 공간해상도를 갖는 지형고도모형의 자료 특성에 따라 나타나는 선구조 추출 결과를 비교하고자 한다. 나아가 본 연구를 통해 원격탐사 자료 및 소프트웨어에 기반한 자동화된 선구조 추출 자료가 야외지질조사 전 기초 참고자료로 활용 가능한지에 대해 고찰해보고자 한다.

2. 연구 지역

본 연구에서는 선구조 추출 분석을 위해 대한민국 동남부 경상권에 위치한 언양, 모량을 포함한 지역을 선정하였다(Fig. 1). 한국지질자원연구원에서 발행한 1:50000 지질 도폭상 모량, 언양, 자인, 유천도폭에 해당하며, 경산, 경주, 밀양, 언양 지역을 부분적으로 포함한다. 해당 면적은 약 40 km × 40 km에 달하며, 지리좌표계상 약 위도 35°30′28″N ~ 35°50′28″N, 경도 128deg;52&prime2″E ~ 129°12′2″E 에 위치한다. 연구지역에서 서쪽에서부터 자인, 밀양, 모량, 양산단층이 순서대로 분포하고 있으며(Korea Institute of Geoscience And Mineral Resources, 1998), 대부분 산악지대로 이루어져 있으며, 단층은 주요 방향은 북북동으로 나타나고 있다(Choi, 2003; Han et al., 2009).

Fig. 1. Landsat-8 optic color composite image of the study area (courtesy of USGS). The red inlet represent the area which is used for automatic lineament extraction.

Fig. 2. Main fault lines overlaid with shaded relief map of Global DEM in the study area (courtesy of DLR;Choi, 2003; Han et al., 2009, Lee and Moon, 2003).

3. 연구 방법

본 선구조 추출 연구는 수치표고모형을 이용한 음영기복도 제작, 선구조 추출, 선구조 수정, 통계 분석 순으로 진행하였다. 전체적인 자료 처리 순서도는 Fig. 3과 같다. 입력 자료는 두 가지 형태의 지형고도모형을 수집하였으며, 약 90 m(3초) 공간해상도를 갖는 Shuttle Radar Topography Mission(SRTM) Digital Elevation Model(DEM)과 12 m 공간해상도를 갖는 TerraSAR-X add-on for Global Digital Elevation Measurement(Global DEM) 자료를 사용하였다. 또한 중간 공간해상도를 갖는 약 1초 간격의 SRTM DEM을 수집하였으나 본 연구 지역에 해당하는 자료가 존재하지 않아 30 m(1초) 간격의 공간해상도를 갖는 지형고도모형을 제작하기 위하여 12 m 공간해상도를 갖는 Global DEM을 재배열(resampling)하였다.

Fig. 3. Flow chart used in this study. The software which used in this study is presented at each processing step.

1) 음영기복도 제작

음영기복도는 특정 태양방위각(sun azimuth angle) 및 고도(sun elevation angle)에 대한 지형고도모형 자료의 음영 등을 사용하여 3차원 지형을 2차원의 평면 위에 도시한 지도이다. 직접 손으로 만드는 기법, 프로그램을 사용하는 기법 등 여러 방법이 있으나(Chang et al., 2003), 본 연구에서는 ENVI^® 프로그램의지형학적 모델링(topography modeling)모듈을 이용하여 수집된 지형고도모형 자료에 대한 음영기복도를 제작하였다. 지형고도모형은 식생, 빌딩 등의 인공 구조물을 제외한 노출 지형(bare earth)을 표현하는 것으로(Yoon et al., 2008), 지형에 대한 고도값을 표현하고 있으며, 일정한 격자를 가진 grid 형태 혹은 TIN(Triangulated Irregular Network) 등 여러 형태가 있다. 음영기복도 제작시 태양고도와 태양방위각에 따라 동일 지역의 음영기복을 다르게 표현할 수 있다. 태양고도는 지평선에서 태양까지의 수직 각도이며, 태양방위각은 정북으로부터 시계방향으로 태양까지의 수평 각도를 의미한다. Kim et al.(1997)는 태양의 고도 및 방위각을 각각 45°, 135°NE로, Abdullah et al.(2010)는 태양의 고도를 30°, 방위각을 0°, 45°, 90°, 135°를 사용하였으며, Kusák and Krbcová(2017)는 태양의 고도를 30°, 방위각을 0°, 45°, 90°, 135°로 지정하여 선구조 추출을 실시하였다. 본 연구에서는 음영기복도 제작과 관련한 선행 연구를 참고하여 태양고도는 45°, 태양방위각은 0°, 45°, 90°, 135°로 지정하여 세 종류의 지형고도모형의 공간 해상도에 따라 총 12개의 음영기복도를 제작하였다.

2) 선구조 추출

위성 및 항공기에서 얻어진 원격탐사 자료에서의 선구조는 주로 색조의 차이에 의해 구별된다. 선구조를 추출하는 방법은 크게 사람에 의한 수동 방법과 컴퓨터 알고리즘에 기반한 자동 방법이 있다. 수동 방법은 소벨(Sobel), 라플라시안(Laplacian)과 같은 이미지 경계 탐지(Image Edge Detection)를 사용하여 경계 향상(edge enhancement)을 시킨 후 직접 선구조를 찾아내는 방법이고(Hung et al., 2005), 자동 방법은 선구조 추출 프로그램을 제작 혹은 활용하여 선구조를 자동으로 추출하는 것이다. 자동화된 선구조 추출은 여러 프로그램으로 제작 수행되어 왔으나 본 연구에서는 PCI Geomatica에서 제공하는 LINE 모듈을 이용하여 선구조 추출을 수행하였다. LINE 모듈에는 선구조 추출을 위한 6가지 변수가 존재하며 연구지역에 적합한 각 변수의 조정을 통해 선구조를 효과적으로 추출할 수 있다.

선행 연구(Kusák and Krbcová, 2017)에서 제시된 6가지 변수를 연구 지역에 그대로 적용해보았을 때 선구조가 제대로 추출되지 않는 문제점이 발생하여 연구 지역에 맞는 변수를 새로 조정할 필요가 있었다. 또한, 6가지 변수를 조정할 때 보다 객관적인 기준을 찾기 위해 한국지질자원연구원의 1:50,000 지질도 상에 기재된 지질학적 선구조 중 연구 지역에 가 뚜렷하게 나타나는 밀양단층과 유사하게 선구조를 추출하는 기으로 찾기 위한 변수를 조정하였다. Table 1은 PCI Geomatica LINE 모듈에서 입력값으로 사용되는 6개의 변수들에 대한 설명이며, 연구 지역에 적합한 각 공간해상도 별 지형고도모형에 따른 6가지 변숫값을 Table 2에 제시하였다.

Table 1. Definition of six parameters in PCI Geomatica LINE Module

Table 2. Six parameters of PCI Geomatica LINE module with respect to spatial resolution of digital elevation model in this study

3) 선구조 수정

동일 지역에서 태양방위각에 따라 선구조 추출 양상이 크게 다르게 나타나는 부분이 있었으며, 지형고도모형 자료의 경계 부분에서는 경계조건에서의 화소값의 손실로 인하여 선구조로 잘못 추출되는 경우도 있었다. 따라서 선구조 수정은 자동화 기법에 의해 얻어진 선구조에 대해 분석자의 주관성이 개입될 수 있으나, 프로그램의 한계를 보완하여 더욱 정확한 선구조 추출하기 위해 필수적이다.

선구조 수정 작업은 ArcGIS^TM 소프트웨어를 이용하여 수행하였다. 태양방위각에 따라 선구조 추출 양상에 차이가 있었기 때문에 모든 태양방위각을 고려하기 위하여 네 개의 선구조 지도를 중첩하여 도시하였다. 네 개의 선구조를 중첩한 결과, 동일하게 겹쳐 있거나, 너무 가까이 붙어있거나, 교차 되어있는 것 등의 중복 선구조가 발견되었다. 선구조 수정 시 최대한 객관성을 보장하기 위하여 가장 높은 공간 해상도의 음영기복도를 배경으로 하여 중복되는 선구조 중 가장 적합한 하나의 선구조를 제외하고 나머지 선구조를 수작업으로 제거하였다. 제거 후 프로그램의 data management tools의 병합(merge) 기능을 통해 하나의 통합된 선구조 지도를 제작하였다.

4) 통계 분석

선구조 수정 후, 선구조 지도에 대해 선구조의 밀도와 방향성 검증 등의 통계 분석을 하였다. 우선, 소프트웨어에 기반한 통계 처리를 위해 벡터화(vectorization)를 실시하여 각 선구조의 첫 점과 끝 점을(x1, y1),(x2, y2) 좌표로 표시하였다. 우선 공간해상도에 따른 선구조 추출 분포를 확인하기 위해 선구조 밀도 분석을 하였다. 밀도 분석은 ArcGIS^TM 프로그램의 spatial analyst tools에 있는 line density 모듈을 사용하여 선구조 밀도 지도를 제작하는 방식으로 진행하였다. 선구조 밀도의 단위는 [선구조의 개수/km²]로 선구조가 적게 분포한 곳은 연두색, 많이 분포한 곳은 붉은색으로 나타내었다. 또한, 공간해상도에 따른 밀도 분포를 비교하기 위해 선구조가 가장 많 추출된 Global DEM에 밀도 범위의 최댓값을 맞추었다. 또한 주요 선구조의 향성을 검증하기 위해 RockWorks^TM 프로그램을 사용하여 지형고도모형과 태양방위각에 따른 rose diagram을 제작하였다.

4. 연구 결과

1) 음영기복도

Fig. 4, 5, 6에서 나타나듯이 태양방위각에 따라 0°에서는 180°방향, 45°에서는 225°방향, 90°에서는 270°방향, 135°에서는 315°방향 등 태양 빛의 반대방향으로 음영이 생기는 것을 확인할 수 있었으며, 지정한 태양방위각에 따라 음영기복도에서 각각 다른 색조 경계가 나타났다. 또한, 지형고도모형의 공간해상도가 증가할수록 보다 자세한 지형 구조에 대한 음영이 표현되었다.

Fig. 4. Shaded relief of SRTM at sun azimuth angle of (a) 0°, (b) 45°, (c) 90°, (d) 135°.

Fig. 5. Shaded relief of 1-arc Global DEM at sun azimuth angle of (a) 0°, (b) 45°, (c) 90°, (d) 135°.

Fig. 6. Shaded relief of Global DEM at sun azimuth angle of (a) 0°, (b) 45°, (c) 90°, (d) 135°.

2) 선구조 추출

PCI Geomatica의 LINE 모듈에서 RADI, GTHR, LTHR, FTHR, ATHR, DTHR의 6가지 변수를 조정하여 선구조 추출을 수행하였다. 6가지 변수 모두 선구조 추출에 영향을 끼쳤으나, 그 중에서도 RADI, GTHR, LTHR이 선구조 추출에 가장 큰 영향을 끼쳤다. 4개의 태양방위각에 따라 선구조의 추출 경향이 다르게 나타났으며, 공간해상도가 증가할수록 보다 많은 선구조가 추출되는 경향이 있음을 알 수 있었다(Fig. 7, 8, 9). 따라서 DEM을 이용한 선구조 추출 시 다양한 태양 방위각을 고려하는 것이 필수적임을 알 수 있다.

Fig. 7. Lineament extraction pattern of SRTM by shadow relief at sun azimuth angle of (a) 0°, (b) 45°, (c) 90°, (d) 135°.

Fig. 8. Lineament extraction pattern of 1-arc Global DEM by shadow relief at sun azimuth angle of (a) 0°, (b) 45°, (c) 90°, (d) 135°.

Fig. 9. Lineament extraction pattern of Global DEM by shadow relief at sun azimuth angle of (a) 0°, (b) 45°, (c) 90°, (d) 135°.

3) 선구조 수정

선구조 추출 후 ArcGIS^TM 프로그램을 사용하여 0°, 45°, 90°, 135°의 모든 태양방위각의 선구조 지도를 병합한 후에 중복되거나 오류가 있는 선구조를 제거하는 작업을 수행하여 최종적인 선구조 지도를 제작하였다. 선구조 수정 작업을 통해 Table 3과 같이 선구조의 개수가 감소됨을 확인할 수 있었다. 최종 선구조의 수는 높은 공간해상도를 가질수록 많이 추출됨을 알 수 있었으며, 12 m 해상도를 갖는 Global DEM의 경우 90 m 공간해상도를 가지는 SRTM DEM에 비해 약 4배 가량 많은 선구조가 추출됨을 알 수 있었다. 흥미로운 점은 30 m의 공간해상도를 갖는 재배열된 Global DEM의 경우 SRTMDEM에 비해 다소 선구조의 개수가 증가하였으나 증가 비율은 Global DEM에 비해 매우 낮음을 알 수 있었다.

Table 3. Number of lineaments according to spatial resolution of DEM

Fig. 10. Final lineament map by DEM:(a) SRTM,(b) 1-arc Global DEM,(c) Global DEM.

4) 통계 분석

(1) 선구조 밀도 분석

선구조 밀도 분석은 ArcGIS^TM 프로그램을 통해 수행하였으며, 격자의 크기와 선구조 밀도의 범주는 공간해상도차이에 따른 동일한 선구조 밀도의 비교를 위하여 제일 낮은 공간해상도와 선구조 밀도를 보이는 SRTM DEM 을 기준에 맞춰 도시하였다. 선구조 밀도의 범위는 0∼5.589 개/km²이며, 선구조가 가장 많이 분포하는 곳은 붉은색으로 표시되며 km² 당 약 6개가 추출되었고, 적은 곳은 초록색으로 km² 당 0개로 나타난다. Fig. 11~12에서와 같이 공간 해상도가 증가할수록 전체적인 선구조 밀도가 증가하는 경향을 보이고 있음을 알 수 있다.

Fig. 11. Lineament density map by DEM with lineament: (a) SRTM, (b) 1-arc Global DEM, (c) Global DEM.

Fig. 12. Lineament density map by DEM without lineament: (a) SRTM, (b) 1-rc Global DEM, (c) Global DEM.

(2) 선구조 방향성 검증

지형표고모형의 공간해상도 차에 따라 태양방위각 각 범위 내의 선구조 최대 빈도가 다르게 나타났기 때문에, 비교 검증을 위해 각 태양방위각에 대해 SRTM 은 90, 1-arc Global DEM은 100, Global DEM은 350을 최대 빈도로 설정하여 도시하였다.

선구조 수정 전에는 추출된 선구조의 전반적인 방향성과 최빈값의 각도 범위가 태양방위각에 따라 다르게 나타나는 것을 확인할 수 있었다. 선구조 방향성을 추출은 5도 단위의 범위를 가지고 추출하였다(Fig. 13~16). 그러나 0°, 45°, 90°, 135°의 모든 태양방위각을 고려한 수정 후 선구조 지도에서는 최빈값의 각도 범위가 SRTM은 N15-20°E, 1-arc Global DEM은 N20-25°E, Global DEM은 N15-20°E로 매우 유사하게 나타났다. 또한 모든 지형고도모형 자료의 선구조를 종합한 결과 연구지역 선구조의 주 방향성은 N15-25°E로 나타났다. 본 연구지역과 정확히 동일한 지역에 대해 선구조 추출 연구가 수행된 사례가 없어 직접 비교는 어려웠지만 양산단층, 밀양단층과 주변 단층의 주 방향인 NNE와 같은 방향성을 보이고 있었다(Choi, 2003; Han et al., 2009; Lee et al., 2003).

Fig. 13. Rose diagram of lineament at sun azimuth angle of (a) 0° (N=221), (b) 45° (N=228), (c) 90° (N=220), (d) 135° (N=211) using SRTM.

Fig. 14. Rose diagram of lineament at sun azimuth angle of (a) 0° (N=221), (b) 45° (N=228), (c) 90° (N=220), (d) 135° (N=211) using 1-arc Global DEM.

Fig. 15. Rose diagram of lineament at sun azimuth angle of (a) 0° (N=1,043), (b) 45° (N=1,033), (c) 90° (N=959), (d) 135° (N=1,013) using Global DEM.

Fig. 16. Rose diagram of lineament at all of sun azimuth angle of (a) SRTM (N=486), (b) 1-arc Global DEM (N=569), (c) Global DEM (N=2,007).

5. 결론 및 토의

본 논문에서는 지형고도모형 자료의 공간해상도 차이에 따른 자동화된 선구조 추출 경향을 분석하였다. 선구조 추출 결과 태양방위각의 변화에 따라 음영기복도의 형태와 선구조의 추출 양상 모두 다르게 추출됨을 확인할 수 있었다. 최종적으로 모든 태양방위각을 고려하여 선구조를 병합한 경우, 공간해상도가 증가함에 따라 선구조 추출 빈도가 증가함을 확인하였으며, 동시에 선구조 밀도 또한 증가하는 양상을 보였다. 고해상도 지형고도모형을 사용할수록 많은 수의 선구조가 추출되는 것을 확인하였다. 그러나 30 m의 공간해상도를 갖는 재배열된 Global DEM의 경우 90 m의 SRTM DEM에 비해 다소 선구조의 개수가 증가하였으나 증가 비율은 Global DEM에 비해 매우 낮았다.

본 논문에서는 소프트웨어를 이용한 선구조 추출이 성공적으로 추출됨을 확인하였다. 최종 선구조 추출을 위해서 일정 부분 수동 보정작업이 필요하였으나 거의 모든 부분에서 자동화가 가능함을 확인하였다. PCI Geomatica의 LINE 모듈을 사용한 선구조 추출 시 해당 연구 지역 지형 구조에 알맞은 6가지 변수를 찾는 작업이 필요함을 확인하였으며 연구 지역에 맞는 변수를 조정하여 선구조 추출 작업을 수행하였다. 6가지 변수가 모두 선구조 추출에 영향이 있음을 확인하였으며, 본 연구의 경우 RADI, GTHR, LTHR의 값이 선구조 추출에 가장 큰 영향을 미침을 알 수 있었다.

SRTM DEM 자료에서 선구조를 추출할 때, 태양방위각이 0°에서 135°로 변화함에 따라 선구조 최빈값의 범위는 N05-10°E에서 N25-30°E로 변화하는 양상을 보였다. 1-arc Global DEM과 Global DEM 또한 태양방위각에 따라 최빈값의 범위가 북에서 북동 방향으로 점진적으로 달라지는 것을 확인하였다. 그러나 선구조 수정 후 모든 태양방위각을 고려한 지형고도모형 자료의 방향성 검증시, 선구조 최빈값의 범위는 N15-25°E로 북북동(NNE)의 방향성을 보였다.

소프트웨어를 이용한 자동화된 선구조 추출은 선구조에 대한 야외조사 전 보조자료로서 활용 가능성이 매우 크다. 그러나 연구 목적 및 지역의 크기에 따라 적절한 공간해상도의 지형고도모형 자료의 선택이 필요함을 알 수 있었으며, 연구 지역 특성을 고려하고, 선구조 추출에 필요한 알맞은 변숫값을 찾는 작업이 필요하며, 태양방위각에 따른 방향 편차도 함께 고려해야 함을 알 수 있었다. 또한 소프트웨어를 이용한 선구조 추출에는 비지질학적 선구조를 포함하고 있을 가능성이 있으므로, 보다 정밀한 수정 작업이 요구될 것으로 판단된다.