• 한국환경농학회지
• /
• 제27권2호
• /
• pp.120-132
• /
• 2008

직접 조사가 힘든 비접근 지역인 북한 전역을 대상으로 DEM을 이용하여 표고 및 경사별 분포 현황을 분석하였고, Landsat TM 위성영상을 이용하여 논, 밭, 산림, 나지, 초지, 물, 간척지, 염전, 건물. 주거지, 기타10개의 분류 항목에 대한 토지피복도를 작성하였다. DEM을 이용한 지형분석 결과 개마고원이 위치한 량강도를 중심으로 동쪽 지역의 표고가 1,000 m 이상으로 높게 나타났고, 평안남도와 황해남도 지역이 낮게 나타났다. 산악지로 구분되는 심한 경사인 E 등급이 전체 면적 대비 38.2%로 가장 넓게 분포하는 것으로 나타났다. 편평한 A 경사는 주로 북한 서해안 지역에 넓게 분포하고 E 경사는 동북부 산악지형에서 높은 비율을 보이고 있어 북한의 전형적인 동고서저의 지형특성을 잘 반영하였다. 위성영상을 이용하여 분류한 북한의 토지피복 항목을 살펴보면 전체 면적 중 산림이 69.6%로 가장 넓게 분포하고 있고, 밭이 15.7%, 나지가 6.6%, 논이 4.2%, 하천과 저수지 등을 포함한 물이 1.6%, 초지가 1.1%, 도시와 주거지가 0.9%인 것으로 나타났다. 행정구역별 지표면 피복을 살펴보면 황해남도와 평안남도 등 서쪽에 위치한 해안가 저위평탄지에 주로 논이 넓게 분포하는 것으로 나타났다. 밭의 분포는 논과 같이 서쪽 지역에 많이 분포하는 경향이었으나 북동쪽에 위치한 함경도와 자강도 및 량강도에도 비교적 고르게 분포하는 것으로 나타났다. 경사등급별로 농경지의 분포를 살펴보면, $0{\sim}2%$인 A 경사에 약 80% 이상 논이 분포하고 있고, 반면 밭은 A, B, C, D, E 등급에 비교적 고르게 분포하고 있는 것으로 나타났다. 농업기후지대별 토지피복 현황을 살펴보면, 논과 밭은 북부 평야지대와 북부 서해안지대에 전체의 약 79%와 45%가 분포하였고 산림은 비교적 모든 농업기후지대별로 고르게 분포하였다. 위성영상을 이용한 원격탐사 기술은 접근이 힘든 지역에 대한 농업기반 및 농경지 정보를 주기적으로 파악할 수 있고, 넓은 지역에 대한 정보 수집이 가능한 장점이 있어, 3년$\sim$5년 주기로 영상분류를 통한 토지피복도를 작성하여 토지이용 및 분류에 대한 시간적 공간적인 변화를 분석한다면 농경지와 산림에 대한 이용 현황 자료를 제공할 수 있고 앞으로의 이용계획 수립에 효율적으로 사용될 수 있을 것으로 생각된다.

• 대한원격탐사학회지
• /
• 제33권5_2호
• /
• pp.699-708
• /
• 2017

본 연구에서는 재배 방법, 토양 특성 등의 정보를 상세하게 수집하기 어려운 지역단위의 콩 작황을 작물생육 모델을 이용하여 예측하는 방법을 개발하고자 하였다. 작물 생육 모델은 DSSAT에 포함된 CROPGRO-Soybean 모델을 이용하였고, 미국의 주요 콩 생산지역인 Illinois주를 연구 사례지역으로 선택하였다. CROPGRO-Soybean 모델을 이용하여 Illinois주의 콩 수량을 예측하기 위한 첫 단계로 다양한 성숙군에 속하는 국내외 품종들을 수집하여 서울대학교농장($37.27^{\circ}N$, $126.99^{\circ}E$)에서 2년동안 파종기 실험을 하여 성숙군(maturity group) I~VI까지의 성숙군별 대표 품종모수(genetic coefficients)를 추정하였다. 대표 품종모수는 각 성숙군 내에 포함되어 있는 품종들의 발육을 매우 정확하게 추정하였다. $10km{\times}10km$ 격자 단위의 기상자료를 바탕으로 성숙군(3), 파종시기(3), 관개여부(2) 등을 조합하여 18가지 조건으로 2000년에서 2011년까지 수량을 각각 모의 하였다. 성숙군과 파종시기는 Illinois주를 위도에 따라 3등분하여 각각 다르게 설정하였다. 관개 및 무관개 조건으로 구분하여 격자 별 모의결과로부터 Illinois주 전체 평균 모의수량을 구하여 연도 별 통계 수량과 비교한 결과 두 경우 모두 실제 수량과 큰 차이를 보일 뿐만 아니라 연차에 따른 수량 변동과 증가 경향을 반영하지 못하였다. 이러한 한계를 극복하고자 처리 별 격자 별로 모의된 수량을 수량을 18개 모의 조건 별로 평균하여 구한 9개 농업지구의 연도별 수량을 독립 변수, 농업지구의 연도별 통계수량을 종속 변수로 하는 중회귀 모델을 구축하였다. 18개 모의 조건 별 수량 외에 품종 개량, 재배 기술 발전 등에 따른 수량의 연차적 변화경향을 반영하기 위하여 연도를 독립변수로 추가하였으며, 중회귀모델은 농업지구와 연도별 수량 변이를 비교적 잘 예측($R^2=0.61$, n=108)하였다. 중회귀 모델로 추정한 9개 농업지구의 연도별 수량을 농업지구별 재배 면적으로 가중 평균한 Illinois의 연도별 추정수량은 통계수량에 매우 근사하였다($R^2=0.80$). 뿐만 아니라 모델 구축 대상연도가 아니고 가뭄으로 수량이 크게 감소한 2012년의 예측 수량은 $3006kg\;ha^{-1}$로 통계수량 $2890kg\;ha^{-1}$과 $116kg\;ha^{-1}$의 근사한 차이를 보였다.

• 대한원격탐사학회지
• /
• 제33권5_1호
• /
• pp.571-585
• /
• 2017

적색경계밴드(red-edge band)가 식물의 생물리적 특성과 밀접한 관계를 가지고 있다고 알려지고 있으며, 이에 따라 최근 적색경계밴드를 포함한 위성영상센서가 증가하고 있다. 본 연구는 향후 농림업중형위성에 적색경계밴드 탑재를 계획하고 있는 점을 감안하여, 적색경계밴드와 관련된 연구 현황과 활용 가치를 분석하고자 한다. 수관울폐도가 높은 우리나라 산림의 엽면적지수(Leaf Area Index, LAI) 추정에 있어서 적색경계밴드의 효과를 분석하였고, 더 나아가 LAI 추정을 위한 최적의 파장폭과 파장영역을 도출하고자 하였다. LAI가 5 이상인 갈참나무와 리기다소나무를 대상으로 4월부터 10월까지 시계열 분광반사 측정자료를 이용하여 LAI와의 상관관계를 분석하였다. 분광반사측정자료에서 5개의 파장폭(10 nm, 20 nm, 30 nm, 40 nm, 50 nm)과 71개의 중심파장(680 nm부터 750 nm까지 1 nm 간격)을 달리하여 모두 355개의 적색경계밴드를 모의 생성했다. 적색경계밴드를 기반으로 하는 두 개의 분광지수 NDRE(normalized difference red-edge index)와 CIRE(chlorophyll index red-edge)를 산출하여 LAI와 상관관계를 분석하였다. 적색경계밴드 기반의 분광지수인 NDRE 및 CIRE는 수관울폐도가 높은 갈참나무와 리기다소나무의 LAI와 높은 상관관계를 얻을 수 있었다. 이는 수관울폐도가 높은 국내 산림에서 일반적으로 사용되는 NDVI가 LAI와의 상관관계가 낮게 나타났던 한계를 해결할 수 있는 가능성을 보여주었다. 10 nm부터 50 nm까지 적색경계밴드의 파장폭 효과는 산림의 LAI와 관계에서 큰 차이를 보이지 않았다. LAI와 최대 상관관계를 보이는 적색경계밴드의 중심파장은 갈참나무에서는 720 nm 부근, 그리고 리기다소나무에서는 710 nm 주변으로 나타났다. 우리나라 농작물 및 산림의 식생정보 획득과 모니터링을 위한 최적의 적색경계밴드의 파장폭과 파장영역을 결정하기 위해서는 다른 생물리적인자(엽록소, 질소, 수분함량, 생체량 등)와의 관계도 충분히 고려하여야 한다.

• 대한원격탐사학회지
• /
• 제36권2_2호
• /
• pp.351-363
• /
• 2020

연안 충적층 지역은 토양침식, 지반침하 등 지질 재해와 태풍, 홍수 해양 재해에 노출되어 인명 또는 재산 피해의 위험에 취약하다. 지반 침하는 지표 물질이 지하로 이동하면서 지반의 점진적인 또는 급격한 침강이 발생하는 현상으로 지속적인 감시가 요구된다. 영상레이더 위상간섭기법 (Radar Interferometry)은 마이크로파 영역에서 관측된 위상정보를 이용하여 지형 변위를 정밀하게 관측할 수 있는 기술이다. Small BAseline Subset(SBAS) 기법은 최소 20장 이상의 영상레이더 자료를 사용하여 대상 지역에 대한 시계열 지표 변위를 정밀하게 분석할 수 있어 지반침하 감시에 유용하다. X- 또는 C-밴드에 비해 장파장인 L-밴드 영상레이더 자료는 긴밀도 유지에 보다 유리하여 지구과학용으로 적합하다. 하지만 L-밴드 ALOS-2 PALSAR-2은 전지구 관측 프로그램을 운영하고 있어 시계열 분석을 하기에 영상 획득이 충분하지 않을 수 있다. 관심 지역인 부산의 경우, Stripmap 영상은 11장으로서 SBAS 시계열 분석 기법을 적용하기에는 부족한 촬영 수이다. 일반적으로 동일한 모드의 영상 간 위상간섭기법의 적용이 가능하지만, 비슷한 관측 기하에서 관측한 Stripmap과 ScanSAR 영상을 이용한 위상간섭기법의 적용이 성공적으로 수행된 바 있다. 해당 지역의 ScanSAR 영상은 18장으로 Stripmap-ScanSAR 위상간섭기법을 적용하면 SBAS 기법을 이용한 향상된 시계열 분석이 가능하게 된다. 본 연구에서는 Gamma 소프트웨어를 이용하여 ALOS-2 PALSAR-2로부터 획득된 L-band 영상 자료에 대한 위상간섭기법 적용 가능여부를 평가하였다. Stripmap-ScanSAR 위상간섭기법 적용을 위해 이종 모드 관측 영상 사이의 chirp bandwidth와 pulse repetition frequency (PRF)의 차이를 고려한 영상레이더 자료 처리를 수행하였으며, radar carrier frequency의 차이 보정과 common band filtering 적용 여부에 따라 발생하는 위상간섭도의 품질을 분석하였다. Radar carrier frequency 차이 보정에 따른 위상간섭도의 긴밀도 변화는 크게 나타나지 않았으나, common band filtering으로 인해, 위상간섭도에서 azimuth 방향으로 주기적인 잡음과 전체적인 긴밀도 저하가 발생하였다. 따라서 Stripmap-ScanSAR 위상간섭기법을 적용하는 경우 두 관측 모드의 range와 azimuth 방향의 밴드 폭의 특성에 따라 주의 깊은 자료 처리가 요구된다.

• 대한원격탐사학회지
• /
• 제31권5호
• /
• pp.385-394
• /
• 2015

소나무재선충병은 우리나라와 중국을 포함한 동아시아 지역뿐만 아니라 유럽 지역의 소나무림에도 막대한 피해를 주고 있다. 소나무재선충에 의한 피해는 임분 내에서 다발적으로 발현되고 급진적으로 진행되고 있으나 기존의 현장조사 방법은 광범위한 지역에 대한 피해목 탐지에 한계가 있다. 본 연구에서는 시계열 초분광 항공영상을 이용하여 소나무재선충병 피해목을 추출하고, 추출된 자료를 이용하여 소나무재선충병 확산 특성을 분석하고자 하였다. 6월, 9월, 10월에 1 m 공간해상도의 초분광 항공사진을 취득하였다. 9월 영상의 Normalized Difference Vegetation Index(NDVI)와 Vegetation Index green(VIgreen)을 이용하여 소나무재선충 피해목을 추출하였다. 추출된 피해목을 잎이 있는 고사목과 잎이 없는 고사목으로 구분하였으며, 6월, 9월, 10월 영상의 분광반사값을 비교하여 9월과 10월에 새로 발생한 피해목을 추출하였다. 추출된 피해목의 시계열 분포를 토대로 공간통계분석 기법을 활용하여 초분광 항공사진 촬영지역 내 소나무재선충병 피해목의 확산 특성을 분석하였다. 그 결과, 대상지 내 소나무재선충병 피해목은 총 2,262본이 추출되었으며, 피해시기에 따라 잎 있는 고사목(작년 피해목)은 604본, 9월과 10월에 새로 발생한 피해목은 각각 300본, 101본으로 분류되었다. 구축된 자료를 이용한 공간 분포형 분석 결과, 작년 피해목과 당년 피해목 모두 집중분포 형태로 나타났으며, 최근거리 분석 결과, 당해년도 고사목의 약 80%는 전년도 고사목 주변 60 m 이내에서 발생하는 것으로 나타났다.

• 대한원격탐사학회지
• /
• 제34권6_2호
• /
• pp.1165-1178
• /
• 2018

남극 빙붕의 붕괴 및 흐름속도의 변화는 빙상에 대한 지지력을 약화시킬 수 있어 해수면 상승에 잠재적인 원인이 될 수 있다. 이 연구에서는 2016년 4월 대규모 붕괴가 발생한 동남극 난센 빙붕에 대해 Landsat-7 Enhanced Thematic Mapper Plus(ETM+) 및 Landsat-8 Operational Land Imager(OLI) 영상을 이용하여 2000년부터 2017년까지의 연간 흐름속도 변화를 분석하였다. 흐름속도 산출을 위해 Landsat의 청색, 녹색, 적색, 근적외선, 전정색 및 첫 번째 주성분 영상 등 총 6개 영상에 orientation correlation 기법을 적용하고, 각각의 변위 산출 결과를 융합하는 다중분광 영상정합 기법을 사용하였다. Landsat 다중분광 영상정합은 난센 빙붕에서 전정색 단일 밴드 영상정합을 사용하는 경우보다 최소 14% 더 넓은 영역에 대해 신뢰할 수 있는 흐름속도를 산출하였고, Global Positioning System(GPS)로 관측된 흐름속도와 비교한 결과 ${\pm}2.1m\;a^{-1}$의 매우 작은 오차를 가지는 것으로 분석되었다. 난센 빙붕에서 2000-2017년 사이에 가장 급격한 흐름속도 증가를 나타낸 곳은 Drygalski 빙하설과 인접한 영역이었으며, 빙붕의 중앙 유선을 따라 측정된 흐름속도는 빙붕 전면(ice front)에 rift가 발달하기 전인 2010년까지 거의 변화가 없었다(${\sim}228m\;a^{-1}$). Rift가 발달하기 시작한 2011-2012년에 rift 상류에서 흐름속도의 가속화가 관측되었으나(${\sim}255m\;a^{-1}$), 이는 2010년에 비해 약 11% 빨라진 것에 불과하였다. 난센 빙붕의 rift가 완전히 발달한 2014년부터 rift 상류의 흐름속도는 다소 감소한 상태(${\sim}225m\;a^{-1}$)로 안정화 되었다. 이는 rift의 발달 및 빙붕 전면의 붕괴가 난센 빙붕의 흐름속도에 거의 영향을 주지 않았음을 의미한다.

• 대한원격탐사학회지
• /
• 제39권3호
• /
• pp.283-295
• /
• 2023

센서 및 위성 기술의 발전에 따라 전세계적으로 다양한 고해상도 다중대역 위성영상이 활용 가능해지고 있다. 다중대역 센서가 가지는 파장에 기인한 고유한 반사, 투과, 산란 특성에 따라 다중대역 위성영상은 지구 관측에 대한 다양한 상호보완적 지표정보를 제공한다. 특히, short-wave infrared (SWIR) 대역은 긴 파장으로 인해 가시광 대역에 비해 Rayleigh 산란에 적게 영향을 받으며, 이로 인해 특정 대기입자를 투과할 수 있다는 특징을 지닌다. 산불, 폭발 등에 의해 발생된 짙은 연기는 가시광 대역의 영상의 가시성을 저하시키고 일부 지역에 대한 지표를 차폐시키는데, SWIR 대역은 이러한 연기에 의해 가려진 지역에 대한 지표정보를 추가로 제공해주기도 한다. 본 연구에서는 이러한 SWIR 대역과 가시광 대역의 영상 정보를 융합하는 다중해상도 변환 기반의 영상 융합 기법을 제안하였다. 제안된 융합 기법의 목적은 상호보완적 관계에 있는 가시광 대역에서의 고해상도 세부적 배경정보와 SWIR 대역에서의 연기 지역에 대한 지표정보를 모두 내포하고 있는 단일 영상을 생성하는 것이다. 이를 위해 본 연구에서는 라플라시안(Laplacian) 피라미드 기반의 다중해상도 변환 기법을 가시광-SWIR 영상 융합에 적용하였다. 다중해상도 변환 기법은 영상 융합에 널리 활용되는 대표적인 영상분해 기반의 방법론으로, 각각의 원 영상을 다양한 스케일로 분해하여 융합하는 기법이다. 또한, 본 연구는 다중해상도 변환 기법에 haze-guided weight map을 융합한 방법론을 제안하였다. Haze-guided weight map은 SWIR 대역이 연기와 같은 특정 대기입자를 투과하여 지표에 대한 정보를 제공해줄 수 있다는 사전지식에 기반하여 제안된 알고리즘으로 다중해상도로 분해된 두 영상을 융합하는 기준이 되는 가중치 지도로써 활용되었다. 제안된 방법론은 가시광 및 SWIR 대역을 포함하고 있는 고해상도 다중대역 위성영상인 Worldview-3 위성영상을 활용하여 검증되었다. 실험 데이터는 주변 산불로 인해 연기가 발생하여 제한된 가시성을 지닌 연기 지역을 포함하고 있으며, 제안된 방법론의 투과 특성을 검증하기 위해 선정되었다. 제안된 기법에 대한 실험결과는 영상 품질 평가 지표를 활용한 정량평가 및 시각평가를 통해 분석되었으며, 결과분석을 통해 연기 지역에 대한 지표정보를 내포하는 SWIR 대역의 밝은 특징값과 가시광 대역 내의 고해상도 정보가 손실없이 최종 융합 영상에 내포됨을 확인할 수 있었다.

• 대한원격탐사학회지
• /
• 제39권2호
• /
• pp.169-181
• /
• 2023

최근 지구온난화의 주원인인 이산화탄소(carbon dioxide, CO₂)의 배출량을 줄이기 위하여 한국은 탄소 배출량 감축목표와 탄소 중립을 선언하였으며, 이에 따른 지역별 배출량과 대기 중 CO₂ 농도의 정확한 평가가 중요해지고 있다. 본 연구에서는 Orbiting Carbon Observatory-2 위성자료와 CO₂ 배출량 자료를 활용하여 위성기반 대기 중 CO₂ 농도와 배출량의 시공간적 차이를 확인하고, 이러한 차이를 식생 성장에 따른 광합성 반응지수인 태양유도 엽록소 형광(solar-induced fluorescence, SIF)을 이용하여 설명하고자 하였다. 2014년부터 2018년까지 한국 지역에서 환경부 온실가스종합정보센터(Greenhouse Gas Inventory and Research Center, GIR) 및 Emissions Database for Global Atmospheric Research (EDGAR) 배출량은 지속적으로 증가하였지만, 위성에서 관측된 CO₂ 농도는 2018년에 전년 대비 감소하는 것으로 나타났다. 지역적으로 살펴보면 경기도, 충청북도는 2018년에 GIR, EDGAR 배출량이 증가하였지만 CO₂ 농도는 감소하였다. 또한, 배출량과 위성관측 CO₂ 농도의 상관성분석에서 서울과 강원도 지역에서 각각 0.22 (GIR), 0.16 (EDGAR)으로 낮은 상관성을 보였다. 대기 중 CO₂ 농도는 SIF와 지역별로 상이한 상관관계를 보였는데, 5~9월의 CO₂-SIF 상관성분석에서 서울과 경기지역은 -0.26의 음의 상관계수를, 충청북도와 강원도는 0.46의 양의 상관계수를 보이며 CO₂ 흡수와 대기 중 농도의 관계성이 지역별로 차이가 있음을 밝혔다. 따라서 대기 중 CO₂ 농도와 배출량 사이의 관계성을 분석함에 있어 CO₂ 흡수 과정에 대한 고려가 필요하다는 것을 시사한다.

• 대한원격탐사학회지
• /
• 제37권5_1호
• /
• pp.1029-1042
• /
• 2021

차세대 중형위성 사업의 일환으로 농지 및 산림에서의 원격 탐사를 위하여 농림위성 (차세대 중형위성 4호)이 발사 예정에 있다. 위성 영상에서 식생의 정량적인 정보를 얻기 위해서는 대기보정을 통한 지표 반사도 취득이 선행되어야 하므로 농림위성을 위한 대기보정 기술 개발은 불가피할 것으로 생각된다. 특히 대기에서의 흡수와 산란 특성은 파장에 따라 다르게 나타나므로 농림위성 파장 영역을 고려한 대기보정 파라미터 민감도 분석이 필요하다. 또한, 농림위성은 5개 채널(Blue, Green, Red, Red edge, Near-infrared)을 보유하고있어 대기보정 주요 파라미터인 AOD (Aerosol optical depth)와 WV (Water vapor)를 직접 산출하기 어려우므로 이를 외부에서 제공할 수 있는 방안을 마련할 필요가 있다. 따라서, 본 연구에서는 농림위성과 유사한 사양을 가진 Sentinel-2 위성 영상을 이용하여 주요 파라미터인 AOD, WV, O₃ 민감도 분석을 수행하고, 파라미터 제공을 위해 천리안 2A (GK2A; GEO-KOMPSAT-2A) 정지궤도 복합위성의 산출물을 이용하여 대기보정 파라미터로서의 활용 가능성을 살펴보았다. 민감도 분석 결과는 AOD가 가장 중요한 파라미터임을 보여주었으며, 근적외선 채널보다는 가시광 채널에서 더 큰 민감도를 가지는 것으로 나타났다. 특히 Blue 채널에서 AOD의 20%의 변화는 지표 반사도에서 약 100%의 오차율을 야기하므로 정확한 지표 반사도 취득을 위해서는 높은 신뢰성을 가진 AOD가 필요할 것으로 생각된다. GK2A AOD 산출물을 이용한 대기보정 결과는 토지피복별 분류 가능성을 이용하여 Sentienl-2 L2A 자료와 비교한 결과, 두 모델별 분류 가능성은 유사하였으나, 파장대가 짧은 영역일수록 GK2A AOD 산출물을 적용한 대기보정 결과가 Sentinel-2 L2A보다 높게 나타났다. 이를 통해 GK2A에서 제공되는 산출물이 향후 농림위성 대기보정 파라미터로서 충분히 활용될 수 있을 것으로 판단된다. 본 연구의 결과는 추후 농림위성 발사 후 대기보정에 참고 자료로서 활용될 수 있을 것으로 기대된다.

• 대한원격탐사학회지
• /
• 제37권3호
• /
• pp.531-542
• /
• 2021

노천 채광을 수행하는 광산은 지표 변화와 환경 교란을 발생시킬 수 있기 때문에 지속적인 모니터링이 필요하다. 노천 광산은 채광 작업장에 식생이 거의 분포하지 않아 InSAR 긴밀도 영상을 이용한 모니터링이 가능하다. 본 연구는 최근 개발된 InSAR 긴밀도 영상 기반의 Normalized Difference Activity Index(NDAI)를 적용하여 광산에서 발생하는 활동을 분석하였다. 3월5일청년광산은 2008년 이후 본격적으로 개발이 확장된 북한의 광산이다. 3월5일청년광산을 촬영한 12일 간격의 Sentinel-1 SAR 영상을 이용하여 획득된 InSAR 긴밀도 영상으로 NDAI 분석을 진행하였다. 우선 2000년부터 약 14년간 발생한 75.24 m의 고도 하강 지역과 약 9.85 m의 고도 상승 지역을 채광 작업장 및 광미 적치장으로 정의하였다. 이후 NDAI 영상을 이용하여 기간별 활동 분석을 진행하기 위해 전체 기간의 평균 영상, 1년 단위의 평균 영상, 및 4개월 단위의 평균 영상을 제작하였다. 2017년부터 2019년까지 광산 활동은 평균적으로 채광 작업장의 중심에서 비교적 활발하였다. 보다 자세한 광산의 활동 변화를 확인하기 위해 시간 간격을 좁혀 1년간의 활동을 알아보고자 하였다. 2017년은 지진파 자료의 정보와 NDAI 영상을 이용하여 인공 지진의 발생 시점과 그 전후에 대하여 RGB 합성 영상을 제작하고 채광 작업장의 활동 변화를 분석하였다. 2017년 4월 30일 발생한 대규모 발파 이후 채광 작업장의 서쪽에서 활발한 활동이 감지되었다. 9월 30일의 두 차례의 발파 이후에는 채광 작업장의 크기가 확장된 것으로 추정된다. 2018년 및 2019년의 활동 변화는 4개월 단위의 시간 평균 영상을 RGB 영상으로 합성하여 분석하였다. 연도별 활동을 분석한 결과, 2018년은 채광 작업장의 북동쪽에서 활발하게 활동하는 영역을 찾을 수 있었으며, 2019년은 광미 적치장에서 확장에 따른 특징적인 활동이 확인되었다. NDAI를 이용한 시계열 분석으로 광학 영상으로는 확인하기 어려운 노천 광산의 무작위적인 지표 변화를 탐지할 수 있었다. 특히 현장 자료를 획득할 수 없는 지역의 광산 활동을 원격 탐사를 이용하여 효과적으로 수행할 수 있었다.