본 연구의 목적은 (1) MCD12Q1 자료를 이용하여 시계열 마스크 맵 작성하여 (2) 시계열 마스크 맵의 연차별 논벼 재배면적 변화를 추출하고, (3) MYD13Q1 NDVI와 단수와의 상관계수를 비교하여 (4) 적용성을 검토하는 것이다. 이를 위해 2002년부터 2019년까지의 연도별 MCD12Q1 PFT 자료를 수집하여 시계열 마스크 맵을 작성하고 고정형 마스크 맵과 비교하였다. 그 결과, 시계열 마스크 맵에 의한 논벼 재배면적이 실제 논벼 재배면적의 변동 특성을 잘 반영하고 있는 것으로 나타났으며, 고정형 마스크 맵과 비교하여 NDVI와 논벼 단수와의 상관계수도 높게 나타나 시계열 마스크 맵이 논벼 고유의 식생지수 추출 및 작황 모니터링, 단수 추정의 정확도를 높이는 방법임을 확인하였다.
본 연구의 목적은 (1) MCD12Q1 자료를 이용하여 시계열 마스크 맵 작성하여 (2) 시계열 마스크 맵의 연차별 논벼 재배면적 변화를 추출하고, (3) MYD13Q1 NDVI와 단수와의 상관계수를 비교하여 (4) 적용성을 검토하는 것이다. 이를 위해 2002년부터 2019년까지의 연도별 MCD12Q1 PFT 자료를 수집하여 시계열 마스크 맵을 작성하고 고정형 마스크 맵과 비교하였다. 그 결과, 시계열 마스크 맵에 의한 논벼 재배면적이 실제 논벼 재배면적의 변동 특성을 잘 반영하고 있는 것으로 나타났으며, 고정형 마스크 맵과 비교하여 NDVI와 논벼 단수와의 상관계수도 높게 나타나 시계열 마스크 맵이 논벼 고유의 식생지수 추출 및 작황 모니터링, 단수 추정의 정확도를 높이는 방법임을 확인하였다.
The objectives of this study were (1) to develop a temporal mask map using MCD12Q1 data, and (2) to extract the annual variations in paddy, (3) to investigate the correlation analysis between MYD13Q1 NDVI and rice yield, and (4) to review its applicability. For these purposes, the temporal mask map ...
The objectives of this study were (1) to develop a temporal mask map using MCD12Q1 data, and (2) to extract the annual variations in paddy, (3) to investigate the correlation analysis between MYD13Q1 NDVI and rice yield, and (4) to review its applicability. For these purposes, the temporal mask map was created using annual MCD12Q1 PFT data from 2002 to 2019, and compared with the fixed mask map. As a result, it found that the temporal mask map well reflected the variations of the paddy area. In addition, the correlation coefficient between NDVI and rice yield was also high significant as compared to the fixed mask map. Therefore, the temporal mask map will be useful for NDVI extraction, crop monitoring, and estimation of rice yield.
The objectives of this study were (1) to develop a temporal mask map using MCD12Q1 data, and (2) to extract the annual variations in paddy, (3) to investigate the correlation analysis between MYD13Q1 NDVI and rice yield, and (4) to review its applicability. For these purposes, the temporal mask map was created using annual MCD12Q1 PFT data from 2002 to 2019, and compared with the fixed mask map. As a result, it found that the temporal mask map well reflected the variations of the paddy area. In addition, the correlation coefficient between NDVI and rice yield was also high significant as compared to the fixed mask map. Therefore, the temporal mask map will be useful for NDVI extraction, crop monitoring, and estimation of rice yield.
* AI 자동 식별 결과로 적합하지 않은 문장이 있을 수 있으니, 이용에 유의하시기 바랍니다.
문제 정의
그러나 해당 자료들의 갱신주기가 3~10년으로 시계열 자료 수집이 불가능하고 위성영상 기반의 토지피복도 제작에 많은 노력이 필요한 관계로 대부분 최근 자료 기반의 고정형 마스크 맵(Fixed Mask map; FM)을 사용하고 있으며, 이는 매년 변화하는 농경지의 현황을 반영하지 못하고 있다는 점에서 정밀한 작황 분석 및 단수 추정에 문제점으로 지적되어 왔다. 따라서 본 연구에서는 MCD12Q1 기반의 1년 단위 시계열 마스크 맵(Temporal Mask map; TM)을 작성하고 이를 이용하여 추출한 정규화 식생지수(Normalized Difference Vegetation Index; NDVI)와 단수와의 상관분석을 통해 시계열 마스크 맵이 논벼 식생지수와 단수와의 관계에 미치는 영향을 검토하였다.
(2012)은 국내 벼의 경우 8월 21일의 Aqua MODIS NDVI와 단수와의 상관계수가 가장 높다고 보고한 바 있다. 따라서 본 연구에서는 선행연구 결과와 우리나라에서 가장 많이 재배되고 있는 중만 생종 벼의 대표 출수기가 8월 20일이라는 점을 고려하여 8월 21일(DOY 233)의 NDVI와 논벼 단수를 대상으로 마스크 맵 변화에 따른 상관계수 변화를 살펴보았다.
본 연구에서는 MCD12Q1 기반의 연차별 시계열 마스크 맵을 작성하고 이를 이용하여 추출한 식생지수와 단수와의 상관분석을 통해 시계열 마스크 맵이 논벼 식생지수와 단수와의 관계에 미치는 영향을 검토하였다. 그 결과, 시계열 마스크 맵에 의한 논벼 재배면적이 실제 논벼 재배면적의 변동 특성을 잘 반영하고 있는 것으로 나타났으며, 마스크 맵에 따른 논벼 NDVI의 연차별 변화도 NDVITM이 NDVIFM 보다 적게 나타나 매년 논벼 재배현황의 변화를 잘 반영하고 있는 것으로 판단되었다.
제안 방법
고정형 마스크 맵과의 비교를 위하여 2018년 기준의 스마트 팜 맵을 대상으로 논에 해당하는 영역만을 추출하여 벡터(Vector) 자료인 상태로 접합(Merge)한 후 90 m 격자의 래스터(Raster) 영상으로 변환하여 전국을 대상으로 스마트 팜 맵 기반의 고정형 마스크 맵을 작성하였다. 시계열 마스크 맵은 MCD12Q1 PFT와 스마트팜 맵을 중첩한 후, PFT의 분류항목 중 ‘Cereal croplands’와 스마트 팜 맵의 논이 공통으로 해당하는 영역만을 추출하여 연도별로 작성하였다.
, 2012). 그 중 MCD12Q1 product는 MODIS Terra 및 Aqua 반사율을 감독분류하여 International Geosphere Biosphere Programme (IGBP), University of MarylanD (UMD), Leaf Area Index/fraction of Photosynthetically Active Radiation (LAI/fPAR), Net Primary Production (NPP), Plant Functional Type (PFT) 5가지 유형의 토지피복도 형태로 제공하고 있으며, 본 연구에서는 식물의 기능에 따라 12개의 항목으로 분류된 PFT 유형을 이용하였다(Table 1).
또한, 자료의 중첩을 위하여 Sinusoidal 투영법을 USGS의 MRT (MODIS Reprojection Tool)를 이용하여 공간해상도 500 m의 UTM WGS-84 좌표로 재투영 (Reproject) 하였다(Fig. 1(a)).
시계열 마스크 맵은 MCD12Q1 PFT와 스마트팜 맵을 중첩한 후, PFT의 분류항목 중 ‘Cereal croplands’와 스마트 팜 맵의 논이 공통으로 해당하는 영역만을 추출하여 연도별로 작성하였다.
대상 데이터
2002년부터 2019년까지 논벼 생육기간인 5월부터 10월까지의 NDVI를 대상으로 시계열 마스크 맵(TM)과 고정형 마스크 맵(FM)을 MYD13Q1 NDVI 자료에 적용하여 마스크 맵에 따른 논벼 시계열 NDVI 변화를 나타내면 Fig. 6과 같다. 그림을 살펴보면 마스크 맵에 따른 NDVI의 연차별 증감 추이는 유사하게 나타났으나 연도별 최고치는 시계열 마스크 맵(TM)에 의한 NDVI (이하 NDVITM)가 고정형 마스크 맵(FM)에 의한 NDVI (이하 NDVIFM)와 비교하여 상대적으로 낮은 것으로 나타난 반면에 연도별 최저치는 NDVIFM가 낮게 나타났다.
논벼단수통계자료는통계청의국가통계포털(http://kosis.kr/)에서 배포하는 2002년∼2019년의 논벼 생산량 자료 중 미곡생산량의 정곡과 백미(92.9%)를 대상으로 단위면적(10 a)당 수량을 수집하였다.
MYD13Q1 product는 구름의 영향을 제거하기 위하여 매일 취득한 영상을 16일씩 중첩하여 조건별로 최대인 값을 계산하는 MVC (Maximum value composite) 방법으로 산정한 16일 단위의 NDVI 및 EVI (Enhanced Vegetation Index; EVI) 자료를 250 m, 1 km, 25 km 세 종류의 공간해상도로 제공하고 있다. 본 연구에서는 250 m 공간해상도의 NDVI 자료를 대상으로 우리나라를 포함하는 h27v04, h27v05, h28v05 세 개의 granule 자료를 2002년부터 2019년까지 수집하여 영상접합(Mosaic)하였다. 접합된 영상은 공간해상도 250 m의 UTM WGS84 좌표로 재투영(Reproject) 하였으며, 영상자료가 데이터 값으로 표현되도록 보정계수를 적용하여 분석에 사용하였다(Fig.
본 연구에서는 마스크 맵에 따른 논벼 NDVI와 단수와의 상관관계를 비교하기 위하여 2002년∼2019년의 18년 동안의 MODIS Product (National Aeronautics and SpaceAdministration;NASA)와농작물생산량통계(KOrea Statistical Information Service; KOSIS) 자료를 수집하였으며, 우리나라의 논벼를 대상으로 산출한 식생지수의 대표성 및 균질성 확보를 위하여 제주도 및 도서 지역을 제외한 내륙 지역만을 대상지역으로 선정하였다.
본 연구에서는 250 m 공간해상도의 NDVI 자료를 대상으로 우리나라를 포함하는 h27v04, h27v05, h28v05 세 개의 granule 자료를 2002년부터 2019년까지 수집하여 영상접합(Mosaic)하였다. 접합된 영상은 공간해상도 250 m의 UTM WGS84 좌표로 재투영(Reproject) 하였으며, 영상자료가 데이터 값으로 표현되도록 보정계수를 적용하여 분석에 사용하였다(Fig. 1(b)).
데이터처리
9%)를 대상으로 단위면적(10 a)당 수량을 수집하였다. 또한, 논벼 NDVI와 단수와의 상관분석을 위하여 SPSS (Statistical Package for the Social Sciences) 프로그램(ver.18, IBM Inc.) 을 사용하였다.
성능/효과
2002년부터 2019년까지의 시계열 마스크 맵에 의한 평균 논벼 재배면적은 449,978 ha로 산출되었으며, 통계자료상의 평균 논 경지면적은 991,422 ha로 나타나 MCD12Q1 PFT 기반의 마스크 맵은 전체 논 경지면적의 약 45%에 해당하는 것으로 나타났다(Table 2). 그러나 시계열 마스크 맵에 의한 면적은 MCD12Q1 공간해상도의 한계로 인하여 25 ha 이하의 소규모 논이 제외된 면적이다.
결과를 살펴보면, 스마트 팜 맵의 논에 해당하는 영역과 MCD12Q1 PFT의 ‘Cereal croplands’ 영역의 분포가 매우 유사한 경향을 나타내었다.
본 연구에서는 MCD12Q1 기반의 연차별 시계열 마스크 맵을 작성하고 이를 이용하여 추출한 식생지수와 단수와의 상관분석을 통해 시계열 마스크 맵이 논벼 식생지수와 단수와의 관계에 미치는 영향을 검토하였다. 그 결과, 시계열 마스크 맵에 의한 논벼 재배면적이 실제 논벼 재배면적의 변동 특성을 잘 반영하고 있는 것으로 나타났으며, 마스크 맵에 따른 논벼 NDVI의 연차별 변화도 NDVITM이 NDVIFM 보다 적게 나타나 매년 논벼 재배현황의 변화를 잘 반영하고 있는 것으로 판단되었다. 또한, 논벼(정곡 및 백미) 단수와의 상관계수도 NDVIFM 보다 NDVITM이 높게 나타나 시계열 마스크 맵이 논벼 고유의 식생지수 추출 및 작황 모니터링, 단수 추정의 정확도를 높이는 방법임을 확인하였다.
그림을 살펴보면 마스크 맵에 따른 NDVI의 연차별 증감 추이는 유사하게 나타났으나 연도별 최고치는 시계열 마스크 맵(TM)에 의한 NDVI (이하 NDVITM)가 고정형 마스크 맵(FM)에 의한 NDVI (이하 NDVIFM)와 비교하여 상대적으로 낮은 것으로 나타난 반면에 연도별 최저치는 NDVIFM가 낮게 나타났다.
5는 시계열 마스크 맵에 의한 논벼 재배면적과 통계자료상의 경지면적의 2002년부터 2019년까지의 시계열 변화이다. 두 자료 간 정량적 오차는 발생하였으나, 전체적으로 감소되는 추세는 유사하게 나타나 MCD12Q1 PFT 기반의 시계열 마스크 맵이 실제 논벼 재배면적의 변동 특성을 잘 반영하고 있는 것으로 나타났다.
MCD12Q1 PFT는 공간해상도가 500 m 로써 실제 농경지를 정밀하게 표현하기에는 한계가 있고 논 주위의 기타 토지피복도 포함된다. 따라서 스마트 팜 맵 또는 통계자료 상의 면적과 비교하여 과대추정 된 것으로 나타났으며, 논벼 고유의 식생지수 추출을 위한 마스크 맵으로 바로 활용하기에는 부적합한 것으로 판단된다.
05 유의수준 하에서 유의한 것으로 나타났으며, 단수와의 상관계수는 정곡과 백미 모두 NDVIFM 보다 NDVITM이 높게 나타나 시계열 마스크 맵의 효과를 확인할 수 있었다. 따라서 위성영상 기반의 식생지수 추출 시 시계열 마스크 맵을 적용할 경우 논벼 고유의 식생지수 추출이 가능하며, 작황 모니터링 및 단수 추정의 정확도를 높일 수 있는 방법으로 판단된다.
그 결과, 시계열 마스크 맵에 의한 논벼 재배면적이 실제 논벼 재배면적의 변동 특성을 잘 반영하고 있는 것으로 나타났으며, 마스크 맵에 따른 논벼 NDVI의 연차별 변화도 NDVITM이 NDVIFM 보다 적게 나타나 매년 논벼 재배현황의 변화를 잘 반영하고 있는 것으로 판단되었다. 또한, 논벼(정곡 및 백미) 단수와의 상관계수도 NDVIFM 보다 NDVITM이 높게 나타나 시계열 마스크 맵이 논벼 고유의 식생지수 추출 및 작황 모니터링, 단수 추정의 정확도를 높이는 방법임을 확인하였다.
그림을 살펴보면 마스크 맵에 따른 NDVI의 연차별 증감 추이는 유사하게 나타났으나 연도별 최고치는 시계열 마스크 맵(TM)에 의한 NDVI (이하 NDVITM)가 고정형 마스크 맵(FM)에 의한 NDVI (이하 NDVIFM)와 비교하여 상대적으로 낮은 것으로 나타난 반면에 연도별 최저치는 NDVIFM가 낮게 나타났다. 또한, 전 기간 동안의 NDVITM 범위는 0.351~0.757으로 나타났으나 NDVIFM 범위는 0342~0.785으로 나타나 NDVI의 변동폭이 NDVITM 보다 NDVIFM가 큰 것으로 나타났다. 이는 마스크 맵에 따른 NDVI 추출 기준에 의한 결과로서 고정형 마스크 맵을 적용할 경우 매년 변화하는 논벼 재배현황을 반영하지 못하여 NDVI 추출시 일부 산림, 나지(휴경지) 등의 기타 토지피복이 포함지만, 시계열 마스크 맵의 경우 매년 논벼 재배현황의 변화를 반영하여 기타 토지피복의 NDVI를 최대한 배제하여 추출할 수 있다.
와 단수(정곡 및 백미)간의 상관분석 결과를 나타내고 있다. 모든 상관계수는 0.05 유의수준 하에서 유의한 것으로 나타났으며, 단수와의 상관계수는 정곡과 백미 모두 NDVIFM 보다 NDVITM이 높게 나타나 시계열 마스크 맵의 효과를 확인할 수 있었다. 따라서 위성영상 기반의 식생지수 추출 시 시계열 마스크 맵을 적용할 경우 논벼 고유의 식생지수 추출이 가능하며, 작황 모니터링 및 단수 추정의 정확도를 높일 수 있는 방법으로 판단된다.
후속연구
포장 단위에서 실측한 경우 현장에서 측정한 식생지수와 논벼 단수는 양의 상관성이 높게 나타남에도 불구하고 위성영상에서 추출한 식생지수는 단수와 상관관계가 낮게 나오거나 음의 상관계수를 갖는 경우가 발생하는데 이는 논벼 마스크 맵에 실제 논이 아닌 지역이 포함될 수 있다는 점과 기상조건(구름, 수증기, 연무 등)으로 인하여 내재된 오차 등이 주요 원인으로 생각된다. 따라서 정밀한 농업 관측을 위해서는 시계열 마스크 맵의 적용과 함께 위성영상의 공간해상도 확보, 전처리 및 보정 기술 개발 등의 기술적 방법의 개선이 필요하며, NDVI 외 다양한 식생지수에 대한 검토 및 작물별, 위성영상별 작황 추정에 최적화된 식생지수의 개발 등이 필요할 것으로 생각된다. 또한, 보다 정밀한 시계열 마스크 맵 작성을 위하여 토지피복 및 작물구분을 위한 고급 분류기법의 적용, 농경지 변화탐지 기술의 개발 및 마스크 맵이 단수 추정 정확도에 미치는 영향에 대한 정량화 연구 등 다양한 노력이 필요할 것으로 사료된다.
MCD12Q1 PFT의 ‘Cereal croplands’ 영역에 포함된 기타 토지피복을 제거함으로서 이를 이용하면 실제 논벼 고유의 식생지수 추출이 가능하다. 또한 논벼 재배현황의 연차 간 변이를 반영할 수 있어 보다 정밀한 작황 모니터링이 가능할 것으로 판단된다.
따라서 정밀한 농업 관측을 위해서는 시계열 마스크 맵의 적용과 함께 위성영상의 공간해상도 확보, 전처리 및 보정 기술 개발 등의 기술적 방법의 개선이 필요하며, NDVI 외 다양한 식생지수에 대한 검토 및 작물별, 위성영상별 작황 추정에 최적화된 식생지수의 개발 등이 필요할 것으로 생각된다. 또한, 보다 정밀한 시계열 마스크 맵 작성을 위하여 토지피복 및 작물구분을 위한 고급 분류기법의 적용, 농경지 변화탐지 기술의 개발 및 마스크 맵이 단수 추정 정확도에 미치는 영향에 대한 정량화 연구 등 다양한 노력이 필요할 것으로 사료된다.
참고문헌 (11)
Allen, J.D., 1990. A Look at the Remote Sensing Applications Program of the National Agricultural Statistics Service, Journal of Official Statistics, 6(4): 393-409.
Becker-Reshef, I., E. Vermote, M. Lineman, and C. Justice, 2010. A generalized regression-based model for forecasting winter wheat yields in Kansas and Ukraine using MODIS data, Remote Sensing of Environment, 114(6): 1312-1323.
Ferencz, Cs., P. Bognar, J. Lichtenberger, D. Hamar, G. Tarcsai, G. Timar, G. Molnar, SZ. Pasztor, P. Steinbach, B. Szekely, O.E. Ferencz, and I. FerenczArkos, 2004. Crop Yield Estimation by Satellite Remote Sensing, International Journal of Remote Sensing, 25(20): 4113-4149.
Hong, S.Y., J.T. Lee, S.K. Rim, and J.S. Shin, 1997. Radiometric estimates of grain yields related to crop aboveground net production (ANP) in paddy rice, Proc. of 1997 International Geoscience and Remote Sensing Symposium, Singapore, Aug. 3-8, pp. 1793-1795.
Hong, S.Y., E.Y. Choe, G.Y. Kim, S.K. Kang, Y.H. Kim, and Y.S. Zhang, 2009. A study on Estimating Rice Yield of North Korea using MODIS NDVI, Proc. of 2009 Korea Remote Sensing Symposium, Seoul, Korea, Mar. 25, pp. 116-120.
Hong, S.Y., J.N. Hur, J.B. Ahn, J.M. Lee, B.K. Min, C.K. Lee, Y.H. Kim, K.D. Lee, S.H. Kim, G.Y. Kim, and K.M. Shim, 2012. Estimating Rice Yield Using MODIS NDVI and Meteorological Data in Korea, Korean Journal of Remote Sensing, 28(5): 509-520 (in Korean with English abstract).
KOrean Statistical Information Service (KOSIS), 2020. http://www.kosis.kr/, Accessed on Aug. 7, 2020.
Na, S.I., J.H. Park, and J.K. Park, 2012. Development of Korean Paddy Rice Yield Prediction Model (KRPM) using Meteorological Element and MODIS NDVI, Journal of the Korean Society of Agricultural Engineers, 54(3): 141-148 (in Korean with English abstract).
Na, S.I., S.Y. Hong, Y.H. Kim, K.D. Lee, and S.Y. Jang, 2013. Prediction of Rice Yield in Korea using Paddy Rice NPP index: Application of MODIS data and CASA model, Korean Journal of Remote Sensing, 29(5): 461-476 (in Korean with English abstract).
Ren, J.Q., Z.X. Chen, Q.B. Zhou, and H.J. Tang, 2008. Regional Yield Estimation for Winter Wheat with MODIS-NDVI Data in Shandong, China, International Journal of Applied Earth Observation and Geoinformation, 10(4): 403-413.
Rojas, O., 2007. Operational Maize Yield Model Development and Validation Based on Remote Sensing and Agro-meteorological Data in Kenya, International Journal of Remote Sensing, 28(17): 3775-3793.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.