基于FastICA算法和MODIS数据的水稻面积提取  

Extracting the paddy rice area from MODIS imagery by FastICA algorithm

作　　者：耿利宁[1,2] 景元书[1,2] 杨沈斌[1] 浩宇 

机构地区：[1]南京信息工程大学江苏省农业气象重点实验室,江苏南京210044 [2]南京信息工程大学应用气象学院,江苏南京210044

出　　处：《大气科学学报》2015年第6期819-826,共8页Transactions of Atmospheric Sciences

基　　金：公益性行业(气象)科研专项(GYHY201306036;GYHY201306035;GYHY20090622);"十二五"农村领域国家科技计划课题(2011BAD32B01);江苏省"青蓝工程"资助项目

摘　　要：以苏、皖、赣三省为研究区域,采用Fast ICA算法从MODIS数据中提取2010年水稻种植面积,并验证该算法在混合像元分解中的有效性。在对2010年46景8 d合成地表反射率产品数据进行预处理的基础上,结合MODIS土地利用产品和平滑滤波算法,构建耕地类型像元的ILSW和INDV时相变化曲线。依据ILSW和INDV曲线在水稻移栽期前后的变化规律,并根据由各地区水稻INDV时相曲线计算得到水稻相似性指数,从MODIS影像中提取水稻像元。采用Fast ICA算法对潜在水稻像元水稻生长期内的INDV时相曲线进行分解,计算每个像元的水稻丰度,绘制水稻丰度图,获取研究区各省水稻分布和种植面积。利用统计年鉴数据和样方资料对Fast ICA算法提取的水稻面积进行了验证。结果显示:采用水稻相似性曲线有利于提高稻田识别效率,所获取的水稻分布与实际情况吻合;Fast ICA算法能够分解不同地区水稻INDV时相曲线;与统计资料比较,江苏、安徽、江西三省水稻面积的提取精度分别为86.4%、87.9%、51.5%。江西水稻面积提取误差主要出现在地形起伏较大的山区。This paper is aimed to extract the paddy rice area of Jiangsu,Anhui and Jiangxi provinces in2010 by using multi-temporal MODIS imagery,and to investigate the efficiency of Fast Independent Component Analysis（ Fast ICA） algorithm in dealing with mixed pixels. In order to map the paddy rice area,the time series of Land Surface Water Content Index（ ILSW） and Normalized Difference Vegetation Index（ INDV） are first produced from MODIS data. Then they are smoothed in order to reduce contamination by clouds. Meanwhile,the pixels with non-crop classes are removed from the index-images using the Land Cover Type data,which is one of the public-distributed MODIS products. The mapping accuracy of the paddy rice area is always affected by the spatial variation of rice planting schedules and rice species. To solve this issue,Rice Similarity Index（ RSI） is calculated from INDVand used to establish standard rice growing curves for different parts of the study area according to the field survey. Combining the standardrice growing curves and the time series of index-images,a rice mapping algorithm is developed and performed to extract the rice planting area. The results show that the obtained rice distribution map is well consistent with the actual situation. However,the retrieved rice planting area is overestimated in most part of the study area,which can be primarily ascribed to the mixed pixels. In this paper,the Fast ICA algorithm is adopted to decompose the mixed pixels using time series of INDV. Before applying the Fast ICA algorithm,the images of times series of INDVare masked by the obtained paddy rice map. Menwhile,several parameters of the Fast ICA are optimized by comparisons. As a result,a map of paddy rice abundance fraction is retrieved for each province by Fast ICA. According to the analysis and comparison with the statistic data,the average accuracies of the obtained rice area are 86. 4%,87. 9% and 51. 5% in Jiangsu,Anhui and Jiangxi provinces,respectively. The larger area error in Jiangxi province

关 键 词：水稻面积提取 相似性曲线 混合像元 丰度 

分 类 号：S511[农业科学—作物学] S127