화산 플리니안 및 서브플리니안 분출 표면 질감 분석 방법

초록

본 논문은 Landsat 5 TM·Landsat 7 ETM+ 위성 영상을 이용해 화산 주변 지표면 질감 변화를 정량적으로 분석하는 새로운 방법론을 제시한다. 공분산 행렬과 웨이브릿 변환을 결합한 알고리즘을 적용하고, 회색조 이미지 유사도 지표로 결과를 검증한다. 라스카르와 차이텐 화산을 사례로 기후·식생 차이에 따른 적용 가능성을 확인하였다.

상세 분석

이 연구는 화산 활동에 따른 표면 변형을 조기에 탐지하기 위한 이미지 기반 텍스처 분석 프레임워크를 구축한다는 점에서 의미가 크다. 먼저, Landsat 5 TM과 Landsat 7 ETM+에서 획득한 30 m 해상도 다중 스펙트럼 영상을 전처리하여 대기 보정·지형 보정을 수행하고, 동일 지역에 대한 시계열 데이터를 정렬한다. 텍스처 특징 추출 단계에서는 회색공동행렬(Gray Level Co‑occurrence Matrix, GLCM)을 이용해 에너지, 대조도, 상관관계, 엔트로피 등 4가지 2차 통계량을 계산한다. GLCM은 방향(0°, 45°, 90°, 135°)과 거리(1~3 픽셀) 파라미터를 다변화하여 다중 스케일 텍스처 정보를 확보한다.

다음으로, 웨이브릿 변환은 2‑차원 이산 웨이브릿(DWT)을 적용해 이미지의 저주파(approximation)와 고주파(detail) 계수를 분리한다. 여기서는 Haar와 Daubechies‑4 두 종류의 웨이브릿을 실험적으로 비교했으며, 특히 고주파 계수의 에너지 분포가 토양·식생 변화에 민감하게 반응함을 확인하였다. 웨이브릿 계수와 GLCM 통계량을 결합한 다중 피처 벡터는 차원 축소를 위해 주성분 분석(PCA)을 적용한 뒤, 유사도 평가지표인 구조적 유사도 지수(SSIM)와 평균 제곱 오차(MSE)로 정량화한다.

라스카르(23 S)와 차이텐(42 S) 화산은 각각 건조·반건조 기후와 습윤·고식생 기후를 대표한다. 두 지역 모두 플리니안·서브플리니안 유형의 폭발을 겪었으나, 침식·식생 복구 속도가 크게 달라 텍스처 변동 패턴이 다르게 나타났다. 라스카르에서는 화산재가 빠르게 노출되어 GLCM 대조도와 웨이브릿 고주파 에너지가 급격히 상승했으며, 이후 토양 재형성 단계에서 점진적인 감소 추세를 보였다. 반면 차이텐에서는 풍부한 식생이 초기 텍스처 변화를 완충했으며, 텍스처 지표의 변화가 더 완만하고 장기적인 추세를 나타냈다. 이러한 차이는 제안된 알고리즘이 기후·식생 조건에 관계없이 일관된 변화를 감지할 수 있음을 시사한다.

마지막으로, 파라미터 최적화 과정에서는 GLCM 거리 파라미터를 2픽셀, 방향을 평균화, 웨이브릿 레벨을 2단계로 설정했을 때 SSIM 점수가 가장 높았다. 또한, 알고리즘의 실행 시간은 평균 3.2 초(한 이미지당)로, 실시간 모니터링 시스템에 적용 가능함을 입증하였다. 전체적으로 본 연구는 위성 영상 기반 텍스처 분석을 통해 화산 주변 지표면 변화를 정량화하고, 위험도 평가 및 조기 경보 시스템에 활용할 수 있는 실용적인 방법론을 제공한다.