Extraction of planting areas of main crops based on sparse representation of time-series leaf area index

Pengxin WANG; Lan XUN; Li LI; Lei WANG; Qingling KONG

doi:10.11834/jrs.20197391

Views : 0 下载量: 111 CSCD: 4 更多指标

PDF
Export
Share
Collection
Album

Extraction of planting areas of main crops based on sparse representation of time-series leaf area index
Vol. 23, Issue 5, Pages: 959-970(2019)
Published： 2019-9 ，

Accepted： 27 May 2018
DOI： 10.11834/jrs.20197391

扫描看全文

王鹏新, 荀兰, 李俐, 王蕾, 孔庆玲. 2019. 基于时间序列叶面积指数稀疏表示的作物种植区域提取. 遥感学报, 23(5): 959–970

Wang P X, Xun L, Li L, Wang L and Kong Q L. 2019. Extraction of planting areas of main crops based on sparse representation of time-series leaf area index. Journal of Remote Sensing, 23(5): 959–970
王鹏新, 荀兰, 李俐, 王蕾, 孔庆玲. 2019. 基于时间序列叶面积指数稀疏表示的作物种植区域提取. 遥感学报, 23(5): 959–970 DOI： 10.11834/jrs.20197391.

Wang P X, Xun L, Li L, Wang L and Kong Q L. 2019. Extraction of planting areas of main crops based on sparse representation of time-series leaf area index. Journal of Remote Sensing, 23(5): 959–970 DOI： 10.11834/jrs.20197391.

摘要

以华北平原黄河以北地区为研究区域，以时间序列叶面积指数LAI(Leaf Area Index)傅里叶变换的谐波特征作为不同作物识别的数据源，利用稀疏表示的分类方法识别2007年—2016年冬小麦、春玉米、夏玉米等主要农作物种植区域。首先利用上包络线Savitzky-Golay滤波分别对2007年—2016年的时间序列MODIS LAI曲线进行重构，进而对重构的年时间序列LAI进行傅里叶变换，以0—5级谐波振幅、1—5级谐波相位作为作物识别的依据，基于各类地物的训练样本，通过在线字典学习算法构建稀疏表示方法的判别字典，对每个待测样本利用正交匹配追踪算法求解稀疏系数，从而计算对应于各类地物的重构误差，根据最小重构误差判定待测样本的作物类型，并对作物识别结果的位置精度进行验证。结果表明，2007年—2016年作物识别的总体精度为77.97%，Kappa系数为0.74，表明本文提出的方法可以用于研究区域主要作物种植区域的提取。

Abstract

Crop mapping is an important component of agriculture monitoring. Accurate information on crop area coverage is vital for food security and the agricultural industry

and the demand for timely crop mapping is high. Previous research indicated that remote sensing technology is a practical and feasible method for agricultural crop area extraction. In this study

the north area of the Yellow River in the North China Plain is chosen as the study area

where the main crops are winter wheat

maize

cotton

and soybean. To obtain the distribution information of crops

the yearly four-day composite MODIS time-series Leaf Area Index (LAI) with 500 m spatial resolution is collected. A total of 92 MODIS LAI images obtained yearly from 2007 to 2016 are used to build time-series LAI curves. To avoid the edge effect of the time-series LAI caused by the Savitzky–Golay filter

the last two phases of LAI images in the last year and the first two phases of LAI images in the next year are added to build the time-series LAI in a year. The Savitzky–Golay filter is then applied on the yearly time-series LAI pixel by pixel to minimize effects of anomalous values caused by atmospheric haze

cloud contamination

and so on. Fourier transform method based on reconstructed LAI is further employed to extract the key parameters. The 11 parameters

including the amplitudes of 0–5 terms and the phases of 1–5 terms

are taken as the features for crop identification. The training samples and verification samples of various crops are obtained through ground investigation and Google Earth images. On the basis of the training samples of various crops

online dictionary learning algorithm is applied to construct the dictionary used to identify the crops. With the dictionary

the orthogonal matching pursuit algorithm is further applied on samples under testing to obtain the sparse representation coefficient. Then the crops are identified according to the minimum reconstruction error

which can be calculated by the dictionary and the coefficient. Therefore

the areas planting winter wheat

spring maize

summer maize

cotton

and orchard from 2007 to 2016 are extracted in the study area. Lastly

the accuracy of the identification results is evaluated yearly by a confusion matrix. Results show that the reconstructed time-series LAI curves are smooth and consistent with crop growth and development characteristics. Overall identification accuracy reaches 77.97% with a Kappa coefficient of 0.74 from 2007 to 2016. User accuracies for individual crops are as follows: winter wheat and summer maize

90.60%; spring maize

73.40%; early summer maize

81.80%; cotton

69.40%; and orchard

81.60%. Annual overall accuracies from 2007 to 2016 range between 70.57% and 83.71% and Kappa coefficients range from 0.66 to 0.81. In conclusion

combining the harmonic characteristics of the time-series LAI with the sparse representation can effectively identify the areas for planting different crops. The approach developed in this study is feasible for extracting information on main crop distribution in the study area.

关键词

遥感叶面积指数稀疏表示Savitzky-Golay滤波华北平原农作物识别

Keywords

remote sensingLeaf Area Indexsparse representationSavitzky-Golay filterNorth China Plaincropsidentification

references

Arvor D, Jonathan M, Meirelles M S P, Dubreuil V and Durieux L. 2011. Classification of MODIS EVI time series for crop mapping in the state of Mato Grosso, Brazil. International Journal of Remote Sensing, 32(22): 7847–7871

Belgiu M and Csillik O. 2018. Sentinel-2 cropland mapping using pixel-based and object-based time-weighted dynamic time warping analysis. Remote Sensing of Environment, 204: 509–523

Candes E J and Wakin M B. 2008. An introduction to compressive sampling. IEEE Signal Processing Magazine, 25(2): 21–30

Chen J, Jönsson P, Tamura M, Gu Z H, Matsushita B and Eklundh L. 2004. A simple method for reconstructing a high-quality NDVI time-series data set based on the savitzky–golay filter. Remote Sensing of Environment, 91(3/4): 332–344

Chen Y, Nasrabadi N M and Tran T D. 2011. Hyperspectral image classification using dictionary-based sparse representation. IEEE Transactions on Geoscience and Remote Sensing, 49(10): 3973–3985

Conrad C, Colditz R, Dech S, Klein D and Vlek P L G. 2011. Temporal segmentation of MODIS time series for improving crop classification in Central Asian irrigation systems. International Journal of Remote Sensing, 32(23): 8763–8778

Efron B, Hastie T, Johnstone I and Tibshirani R. 2004. Least angle regression. The Annals of Statistics, 32(2): 407–499

方墨人, 田庆久. 2004. 基于MODIS的LAI时间序列谱的地物分类方法研究. 国土资源遥感, 16(3): 5–7, 22

Fang M R and Tian Q J. 2004. The study of land type classification technology based on time-series LAI of MODIS. Remote Sensing for Land and Resources, 16(3): 5–7, 22

Jia K, Liang S L, Wei X Q, Yao Y J, Su Y R, Jiang B and Wang X X. 2014. Land cover classification of Landsat data with phenological features extracted from time series MODIS NDVI data. Remote Sensing, 6(11): 11518–11532

梁守真, 邢前国, 施平, 周迪. 2011. 山东省典型地表覆被NDVI时间序列谐波分析. 生态学杂志, 30(1): 59–65

Liang S Z, Xing Q G, Shi P and Zhou D. 2011. Harmonic analysis on NDVI time series of typical land covers in Shandong Province. Chinese Journal of Ecology, 30(1): 59–65

刘帅, 朱亚杰, 薛磊. 2015. 一种结合稀疏表示和纹理分块的遥感影像超分辨率方法. 武汉大学学报(信息科学版), 40(5): 578–582

Liu S, Zhu Y J and Xue L. 2015. Remote sensing image super-resolution method using sparse representation and classified texture patches. Geomatics and Information Science of Wuhan University, 40(5): 578–582

Mairal J, Bach F, Ponce J and Sapiro G. 2010. Online learning for matrix factorization and sparse coding. The Journal of Machine Learning Research, 11: 19–60

Mallat S G and Zhang Z F. 1993. Matching pursuits with time-frequency dictionaries. IEEE Transactions on Signal Processing, 41(12): 3397–3415

Massey R, Sankey T T, Congalton R G, Yadav K, Thenkabail P S, Ozdogan M and Meador A J S. 2017. MODIS phenology-derived, multi-year distribution of conterminous U.S. crop types. Remote Sensing of Environment, 198: 490–503

莫兴国, 刘苏峡, 林忠辉, 邱建秀. 2011. 华北平原蒸散和GPP格局及其对气候波动的响应. 地理学报, 66(5): 589–598

Mo X G, Liu S X, Lin Z H and Qiu J X. 2011. Patterns of evapotranspiration and GPP and their responses to climate variations over the North China Plain. Acta Geographica Sinica, 66(5): 589–598

Pan Y Z, Li L, Zhang J S, Liang S L, Zhu X F and Sulla-Menashe D. 2012. Winter wheat area estimation from MODIS-EVI time series data using the Crop Proportion Phenology Index. Remote Sensing of Environment, 119: 232–242

Savitzky A and Golay M J E. 1964. Smoothing and differentiation of data by simplified least squares procedures. Analytical Chemistry, 36(8): 1627–1639

Skakun S, Franch B, Vermote E, Roger J C, Becker-Reshef I, Justice C and Kussul N. 2017. Early season large-area winter crop mapping using MODIS NDVI data, growing degree days information and a Gaussian mixture model. Remote Sensing of Environment, 195: 244–258

宋琳, 程咏梅, 赵永强. 2012. 基于稀疏表示模型和自回归模型的高光谱分类. 光学学报, 32(3): 0330003

Song L, Cheng Y M and Zhao Y Q. 2012. Hyper-spectrum classification based on sparse representation model and auto-regressive model. Acta Optica Sinica, 32(3): 0330003

宋相法, 焦李成. 2012. 基于稀疏表示及光谱信息的高光谱遥感图像分类. 电子与信息学报, 34(2): 268–272

Song X F and Jiao L C. 2012. Classification of hyperspectral remote sensing image based on sparse representation and spectral information. Journal of Electronics and Information Technology, 34(2): 268–272

唐鹏钦, 吴文斌, 姚艳敏, 杨鹏. 2011. 基于小波变换的华北平原耕地复种指数提取. 农业工程学报, 27(7): 220–225

Tang P Q, Wu W B, Yao Y M and Yang P. 2011. New method for extracting multiple cropping index of North China Plain based on wavelet transform. Transactions of the CSAE, 27(7): 220–225

唐晓晴, 刘亚洲, 陈骏龙. 2016. 基于稀疏表示的遥感图像分类方法改进. 计算机工程, 42(3): 254–258, 265

Tang X Q, Liu Y Z and Chen J L. 2016. Improvement of remote sensing image classification method based on sparse representation. Computer Engineering, 42(3): 254–258, 265

Tropp J A and Gilbert A C. 2007. Signal recovery from random measurements via orthogonal matching pursuit. IEEE Transactions on Information Theory, 53(12): 4655–4666

王鹏新, 荀兰, 李俐, 解毅, 王蕾. 2017. 基于时间序列叶面积指数傅里叶变换的作物种植区域提取. 农业工程学报, 33(21): 207–215

Wang P X, Xun L, Li L, Xie Y and Wang L. 2017. Extraction of planting areas of main crops based on Fourier transformed characteristics of time series leaf area index products. Transactions of the Chinese Society of Agricultural Engineering, 33(21): 207–215

王学, 李秀彬, 谈明洪, 辛良杰. 2015. 华北平原2001-2011年冬小麦播种面积变化遥感监测. 农业工程学报, 31(8): 190–199

Wang X, Li X B, Tan M H and Xin L J. 2015. Remote sensing monitoring of changes in winter wheat area in North China Plain from 2001 to 2011. Transactions of the Chinese Society of Agricultural Engineering, 31(8): 190–199

Wardlow B D and Egbert S L. 2008. Large-area crop mapping using time-series MODIS 250 m NDVI data: an assessment for the U. S. Central Great Plains. Remote Sensing of Environment, 112(3): 1096–1116

Whelen T and Siqueira P. 2017. Use of time-series L-band UAVSAR data for the classification of agricultural fields in the San Joaquin Valley. Remote Sensing of Environment, 193: 216–224

Wright J, Yang A Y, Ganesh A, Sastry S S and Ma Y. 2009. Robust face recognition via sparse representation. IEEE Transactions on Pattern Analysis and Machine Intelligence, 31(2): 210–227

吴文斌, 杨鹏, 唐华俊, 周清波, Ryosuke S, 张莉. 2009. 两种NDVI时间序列数据拟合方法比较. 农业工程学报, 25(11): 183–188

Wu W B, Yang P, Tang H J, Zhou Q B, Ryosuke S and Zhang L. 2009. Comparison of two fitting methods of NDVI time series datasets. Transactions of the CSAE, 25(11): 183–188

夏琴, 邢帅, 马东洋, 莫德林, 李鹏程, 葛忠孝. 2016. 遥感卫星影像K-SVD稀疏表示去噪. 遥感学报, 20(3): 441–449

Xia Q, Xing S, Ma D Y, Mo D L, Li P C and Ge Z X. 2016. An improved K-SVD-based denoising method for remote sensing satellite images. Journal of Remote Sensing, 20(3): 441–449

尹雯, 李元祥, 周则明, 刘世前. 2013. 基于稀疏表示的遥感图像融合方法. 光学学报, 33(4): 0428003

Yin W, Li Y X, Zhou Z M and Liu S Q. 2013. Remote sensing image fusion based on sparse representation. Acta Optica Sinica, 33(4): 0428003

Alert me when the article has been cited

提交

GF-1 leaf area index product across China based on three-dimensional stochastic radiation transfer model

Spectral unmixing method considering endmember variability of vegetation

Improving the performance of smartphone-derived crop leaf area index

Review of forest Leaf Area Index retrieval over rugged terrain based on remotely sensed data

Estimation of forest leaf area index based on spectrally corrected airborne LiDAR pulse penetration index by intensity of point cloud