为了解决MODIS数据中普遍存在的混合像元问题,该文利用2008年和2009年多时相的MODIS13Q1影像,以经过优化的N-FINDR算法进行线性混合像元分解提取冬小麦种植面积,各省的误差均控制在正负4%左右。利用同期多时相的HJ-1星分类数据作为参考值,在试验区域选择14个均匀分布的样区验证混合像元分解结果。结果显示6个样区的相对误差在10%以内,其余8个样区的误差基本在15%左右。该研究可为冬小麦种植面积的监测提供参考。
Winter wheat is one of the main food crops in the north of China. It is significant to monitor winter wheat planting areas for China’s grain policy and economic planning. The MODIS products are outstanding with the characteristics of large area coverage, frequent repeat, and free access to download. It offers a valuable application on long-term and large-area detection of winter wheat. Because of the coarse spatial resolution of MODIS products, the mixed pixels become the common problem existing in MODIS data. Therefore, it is necessary to solve the problem of mixed pixels in crop area extraction with MODIS data, In this study, we chose the Huanghuaihai Plain (including Hebei province, Shandong province, Henan province, Beijing, and Tianjin) as the study area, and used multi-temporal MODIS data in 2008 and 2009 to extract the winter wheat area with an optimized N-FINDR algorithm and linear unmixing method. In a traditional N-FINDR algorithm, all pixels in the image would be traversed to find the pixel group that can form a simplex with the maximum area. The optimized N-FINDR algorithm we used simplifies the procedure by finding the points set that can form a triangle with the maximum area in a two-dimensional plane composed by any two bands first, then the vertex of the triangle are taken as candidate endmembers, and final endmembers are obtained by traversing all the candidate endmembers. In order to find points set in a two-dimensional plane, we used the convex hull property of a polygon with rotating calipers. This optimized algorithm can improve time complexity from O(n3 ) to O(n2 ). Comparing this with national statistical data in 2009, the relative error of the extracted winter wheat planting area was less than 4% for each province. The results showed that the method we used was applicable for winter wheat area extraction on a large scale. In order to further validate the results, we selected 14 sample areas, and multi-temporal HJ-1 data at same period were taken to produce the winter wheat pl