高分二号遥感影像提取冬小麦空间分布

精细的农作物空间分布数据对于资源、环境、生态、气候变化和粮食安全问题均具有重要的意义,卷积神经网络已经成为从遥感影像中提取农作物空间分布数据的主要方法,但提取结果中的种植区域边缘往往比较粗糙。本文以高分二号遥感影像为数据源,选择冬小麦为提取目标,利用RefineNet模型和最大后验概率模型构建冬小麦遥感提取模型WWRSE(Winter Wheat Remote Sensing Extraction),获取精细的冬小麦空间分布数据。WWRSE模型利用RefineNet网络提取像素的语义特征,使用改进的SoftMax模型生成像素的类别概率向量;以类别概率向量的最大分量与次大分量的差值作为置信度,根据置信度将类别概率向量分为可信和不可信两组,可信组直接使用最大分量对应的类别标签作为相应像素的分类结果;结合最大后验概率模型确定不可信组像素的分类结果。利用随机梯度法对WWRSE模型进行训练。选择SegNet、DeepLab、RefineNet作为对比模型进行实验,WWRSE提取结果的精度为92.9%,比SegNet提高了13.8%,比DeepLab提高了10.9%,比RefineNet提高了8.6%。实验结果表明WWRSE模型在提取冬小麦空间分布数据方面具有一定的优势。WWRSE模型提取的结果能够为大范围冬小麦种植面积统计提供依据。     

FY-3250m分辨率数据的华北平原冬小麦提取

冬小麦是中国最主要的粮食作物之一,利用遥感技术提取冬小麦种植区是遥感应用研究的一个重要方向。2008年以来发射的系列风云三号(FY-3)卫星均携带着中分辨率光谱成像仪(MERSI),该传感器有5个250 m分辨率的波段,波段范围包括可见光、近红外和热红外,观测数据包含丰富的地表信息,为大范围冬小麦种植区提取提供了新的数据源。首先,选取生长季前期多幅高质量的MERSI数据,采用分层提取的方法,对于不同的层次选用与待提取类别最为敏感的特征波段来构建相应的决策树,从而将每一幅影像中冬小麦种植区提取出来,然后,将多幅数据融合为一幅生长季内的冬小麦种植区图。最后,使用野外实地调查的数据进行精度验证,面积提取精度为90.8%。结果表明,在春季返青后,即可做出当季冬小麦种植分布图,为农情监测提供及时的信息支撑。     

基于高分六号WFV数据的冬小麦种植面积提取

冬小麦是我国重要的粮食作物之一,准确获取冬小麦种植面积具有重要的现实意义。为探究高分六号卫星影像进行冬小麦遥感监测的可行性和精确性,本文选取甘肃省崆峒区为研究区,运用红边波段+监督分类中的支持向量机法,提取了2019年崆峒区冬小麦种植面积,并利用混淆矩阵对分类结果进行精度验证。结果表明:提取崆峒区冬小麦种植面积为15045 hm 2,与实际种植面积相比,误差率为1.02%;该模型能有效地提取崆峒区冬小麦,总体精度为98.88%,Kappa系数为0.97;红边波段能有效地提取干扰地物,提取精度比直接使用监督分类高7.88个百分点;GF6影像在提取冬小麦种植面积上具有明显优势。     

统计数据总量约束下全局优化阈值的冬小麦分布制图

大范围、长时间和高精度农作物空间分布基础农业科学数据的准确获取对资源、环境、生态、气候变化和国家粮食安全等问题研究具有重要现实意义和科学意义。本文针对传统阈值法农作物识别过程中阈值设置存在灵巧性差和自动化程度低等弱点,以中国粮食主产区黄淮海平原内河北省衡水市景县为典型实验区,首次将全局优化算法应用于阈值模型中阈值优化选取,开展了利用全局优化算法改进基于阈值检测的农作物分布制图方法创新研究。以冬小麦为研究对象,国产高分一号(GF-1)为主要遥感数据源,在作物面积统计数据为总量控制参考标准和全局参数优化的复合型混合演化算法SCE-UA (Shuffled Complex Evolution-University of Arizona)支持下,提出利用时序NDVI数据开展阈值模型阈值参数自动优化的冬小麦空间分布制图方法。最终,获得实验区冬小麦阈值模型最优参数,并利用优化后的阈值参数对冬小麦空间分布进行提取。通过地面验证表明,利用本研究所提方法获取的冬小麦识别结果分类精度均达到较高水平。其中冬小麦识别结果总量精度达到了99.99%,证明本研究所提阈值模型参数优化方法冬小麦提取分类结果总量控制效果良好;同时,与传统的阈值法、最大似然和支持向量机等分类方法相比,本研究所提阈值模型参数优化法区域冬小麦作物分类总体精度和Kappa系数分别都有所提高,其中,总体精度分别提高4.55%、2.43%和0.15%,Kappa系数分别提高0.12、0.06和0.01,这体现出SCE-UA全局优化算法对提高阈值模型冬小麦空间分布识别精度具有一定优势。以上研究结果证明了利用本研究所提基于作物面积统计数据总量控制以及SCE-UA全局优化算法支持下阈值模型参数优化作物分布制图方法的有效性和可行性,可获得高精度冬小麦作物空间分布制图结果,这对提高中国冬小麦空间分布制图精度和自动化水平具有一定意义,也可为农作物面积农业统计数据降尺度恢复重建和大范围区域作物空间分布制图研究提供一定技术参考。     

Comparison of three remotely sensed drought indices for assessing the impact of drought on winter wheat yield

ABSTRACT

Agricultural drought threatens food security. Numerous remote-sensing drought indices have been developed, but their different principles, assumptions and physical quantities make it necessary to compare their suitability for drought monitoring over large areas. Here, we analyzed the performance of three typical remote sensing-based drought indices for monitoring agricultural drought in two major agricultural production regions in Shaanxi and Henan provinces, northern China (predominantly rain-fed and irrigated agriculture, respectively): vegetation health index (VHI), temperature vegetation dryness index (TVDI) and drought severity index (DSI). We compared the agreement between these indices and the standardized precipitation index (SPI), soil moisture, winter wheat yield and National Meteorological Drought Monitoring (NMDM) maps. On average, DSI outperformed the other indices, with stronger correlations with SPI and soil moisture. DSI also corresponded better with soil moisture and NMDM maps. The jointing and grain-filling stages of winter wheat are more sensitive to water stress, indicating that winter wheat required more water during these stages. Moreover, the correlations between the drought indices and SPI, soil moisture, and winter wheat yield were generally stronger in Shaanxi province than in Henan province, suggesting that remote-sensing drought indices provide more accurate predictions of the impacts of drought in predominantly rain-fed agricultural areas.     

基于图优化的紧耦合双目视觉/惯性/激光雷达SLAM方法

基于单一传感器的同时定位与地图构建技术已经逐渐不能满足移动机器人、无人机及自动驾驶车辆等智能移动载体日益复杂的应用场景。为了进一步提升移动载体在复杂环境下的定位与建图性能,基于多传感器融合的SLAM技术成为目前研究的热点内容。本文提出了一种基于图优化的紧耦合双目视觉/惯性/激光雷达SLAM方法(S-VIL SLAM),该方法在视觉惯性系统中引入激光雷达原始观测,基于滑动窗口实现了IMU量测、视觉特征及激光点云特征的多源数据联合非线性优化。利用视觉与激光雷达的互补特性设计了视觉增强的激光雷达闭环优化算法,进一步提升了多源融合SLAM系统的全局定位与建图精度。为了验证本文算法的性能,利用自主搭建的集成多传感器的硬件采集平台在室外场景下进行了车载试验。试验结果表明,本文提出的紧耦合双目视觉/惯性/激光雷达里程计相比于紧耦合双目视觉惯性里程计和激光雷达里程计定位定姿性能显著提升,视觉增强的激光雷达闭环优化算法能够在大尺度场景下有效探测出轨迹中的闭环信息,并实现高精度的全局位姿图优化,经过闭环优化的点云地图具有良好的分辨率和全局一致性。     

Envisat ASAR的区域森林-非森林制图

Envisat卫星ASAR传感器的双极化数据对区域森林监测十分有效。通过分别采用SRTM DEM和Landsat TM图像对地形起伏区域和平坦区域的SAR图像进行地理编码,发展了一种SAR图像自动预处理方法。基于冬季单时相ASAR数据的HH(水平发射,水平接收)、HV(水平发射,垂直接收)极化比值和HV极化图像,提出了一种面向对象的森林-非森林分类方法。将之应用于中国东北森林/非森林制图,分类总体精度、森林用户精度和生产者精度分别为83.7%,85.6%和75.7%。结果表明,本文提出的方法十分适合区域森林-非森林制图的业务化运行。     

Step-wise Land-class Elimination Approach for extracting mixed-type built-up areas of Kolkata megacity

The extraction of urban built-up areas is an important aspect of urban planning and understanding the complex drivers and biophysical mechanism of urban climate processes. However, built-up area extraction using Landsat data is a challenging task due to spatio-temporal dynamics and spatially intermixed nature of Land Use and Land Cover (LULC) in the cities of the developing countries, particularly in tropics. In the light of advantages and drawbacks of the Normalized Difference Built-up Index (NDBI) and Built-up Area Extraction Method (BAEM), a new and simple method i.e. Step-wise Land-class Elimination Approach (SLEA) is proposed for rapid and accurate mapping of urban built-up areas without depending exclusively on the band specific normalized indices, in order to pursue a more generalized approach. It combines the use of a single band layer, Normalized Difference Vegetation Index (NDVI) image and another binary image obtained through Logit model. Based on the spectral designation of the satellite image in use, a particular band is chosen for identification of water pixels. The Double-window Flexible Pace Search (DFPS) approach is employed for finding the optimum threshold value that segments the selected band image into water and non-water categories. The water pixels are then eliminated from the original image. The vegetation pixels are similarly identified using the NDVI image and eliminated. The residual pixels left after elimination of water and vegetation categories belong either to the built-up areas or to bare land categories. Logit model is used for separation of the built-up areas from bare lands. The effectiveness of this method was tested through the mapping of built-up areas of the Kolkata Metropolitan Area (KMA), India from Thematic Mapper (TM) images of 2000, 2005 and 2010, and Operational Land Imager (OLI) image of 2015. Results of the proposed SLEA were 95.33% accurate on the whole, while those derived by the NDBI and BAEM approaches returned an overall accuracy of 83.67% and 89.33%, respectively. Comparisons of the results obtained using this method with those obtained from NDBI and BAEM approaches demonstrate that the proposed approach is quite reliable. The SLEA generates new patterns of evidence and hypotheses for built-up areas extraction research, providing an integral link with statistical science and encouraging trans-disciplinary collaborations to build robust knowledge and problem solving capacity in urban areas. It also brings landscape architecture, urban and regional planning, landscape and ecological engineering, and other practice-oriented fields to bear together in processes for identifying problems and analyzing, synthesizng, and evaluating desirable alternatives for urban change. This method produced very accurate results in a more efficient manner compared to the earlier built-up area extraction approaches for the landscape and urban planning.     

Use of field reflectance data for crop mapping using airborne hyperspectral image

Recent developments in hyperspectral remote sensing technologies enable acquisition of image with high spectral resolution, which is typical to the laboratory or in situ reflectance measurements. There has been an increasing interest in the utilization of in situ reference reflectance spectra for rapid and repeated mapping of various surface features. Here we examined the prospect of classifying airborne hyperspectral image using field reflectance spectra as the training data for crop mapping. Canopy level field reflectance measurements of some important agricultural crops, i.e. alfalfa, winter barley, winter rape, winter rye, and winter wheat collected during four consecutive growing seasons are used for the classification of a HyMAP image acquired for a separate location by (1) mixture tuned matched filtering (MTMF), (2) spectral feature fitting (SFF), and (3) spectral angle mapper (SAM) methods. In order to answer a general research question “what is the prospect of using independent reference reflectance spectra for image classification”, while focussing on the crop classification, the results indicate distinct aspects. On the one hand, field reflectance spectra of winter rape and alfalfa demonstrate excellent crop discrimination and spectral matching with the image across the growing seasons. On the other hand, significant spectral confusion detected among the winter barley, winter rye, and winter wheat rule out the possibility of existence of a meaningful spectral matching between field reflectance spectra and image. While supporting the current notion of “non-existence of characteristic reflectance spectral signatures for vegetation”, results indicate that there exist some crops whose spectral signatures are similar to characteristic spectral signatures with possibility of using them in image classification.     

Understanding the Unique Spectral Signature of Winter Rape

Driven by significant technological developments in the hyperspectral imaging, material mapping using reference spectra has received renewed interest of the remote sensing community. The applicability of reference spectral signatures in image classification depends mainly on the material type and its spectral signature behaviour. Identification and spectral characterization of materials which exhibit unique spectral behaviour is the first step in this approach. Consequently there have been active researches for the identification of surface materials which exhibit unique spectral signatures. The uniqueness of reflectance signature of winter rape relative to its co-occurring crop species was reported in this study. Reflectance spectral libraries constructed from field spectral reflectance measurements collected over five agricultural crops (alfalfa, winter barley, winter rape, winter rye, and winter wheat) during four subsequent growing seasons were classified by the linear discriminant analysis (LDA). Further, the reference field spectral database was used for the spectral feature fitting and classification of a historical HyMAP airborne hyperspectral imagery acquired at a separate site, by spectral library search. Results indicate the existence of a meaningful spectral matching between image and field spectra for winter rape and demonstrate the potential for transferring spectral library for hyperspectral image classification. The observed consistency in the discrimination of winter rape demonstrates experimentally the fundamental principle of remote sensing which suggests the theoretical existence of unique spectral signatures for materials which can be incorporated as reference spectral signatures for hyperspectral image classification.     

Assessing the performance of two unsupervised dimensionality reduction techniques on hyperspectral APEX data for high resolution urban land-cover mapping

Despite the high richness of information content provided by airborne hyperspectral data, detailed urban land-cover mapping is still a challenging task. An important topic in hyperspectral remote sensing is the issue of high dimensionality, which is commonly addressed by dimensionality reduction techniques. While many studies focus on methodological developments in data reduction, less attention is paid to the assessment of the proposed methods in detailed urban hyperspectral land-cover mapping, using state-of-the-art image classification approaches. In this study we evaluate the potential of two unsupervised data reduction techniques, the Autoassociative Neural Network (AANN) and the BandClust method – the first a transformation based approach, the second a feature-selection based approach – for mapping of urban land cover at a high level of thematic detail, using an APEX 288-band hyperspectral dataset. Both methods were tested in combination with four state-of-the-art machine learning classifiers: Random Forest (RF), AdaBoost (ADB), the multiple layer perceptron (MLP), and support vector machines (SVM). When used in combination with a strong learner (MLP, SVM) BandClust produces classification accuracies similar to or higher than obtained with the full dataset, demonstrating the method’s capability of preserving critical spectral information, required for the classifier to successfully distinguish between the 22 urban land-cover classes defined in this study. In the AANN data reduction process, on the other hand, important spectral information seems to be compromised or lost, resulting in lower accuracies for three of the four classifiers tested. Detailed analysis of accuracies at class level confirms the superiority of the SVM/Bandclust combination for accurate urban land-cover mapping using a reduced hyperspectral dataset. This study also demonstrates the potential of the new APEX sensor data for detailed mapping of land cover in spatially and spectrally complex urban areas.     

超像元尺度下的多光谱图像超分辨率洪水淹没制图

超分辨率制图SRM (Super-resolution Mapping)技术可以有效地处理遥感图像中的混合像元,获得准确的地物类别分布信息。目前,SRM技术已经成功地应用于多光谱图像洪水淹没定位中,称为超分辨率洪水淹没制图SRFIM (Super-resolution Flood Inundation Mapping)。然而,现有的SRFIM方法往往基于像元尺度空间相关性,这种空间相关性考虑设定的矩形窗内的像元之间的空间关系,但实际情况下淹没区域与非淹没区域的形状是不规则的,因此这种像元尺度空间相关性不够准确,影响最终的洪水淹没制图精度。为了解决这一问题,提出了超像元尺度空间相关性下的多光谱图像超分辨率洪水淹没制图SSSC-SRFIM (Super-resolution Flood Inundation Mapping for Multispectral Image Based on Super-pixel Scale Spatial Correlation)。在SSSC-SRFIM中,首先利用双立方插值改善原始粗糙多光谱图像,获得改善后的图像,并利用光谱解混方法对改善后的图像进行光谱解混,获得具有每个亚像元属于淹没类别概率值的丰度图像;然后利用主成分分析法提取改善后图像的第一主成分,并利用基于多分辨率的图像分割算法分割第一主成分,获得不规则形状的超像元;再者将丰度图像与超像元进行整合计算,并引入随机游走算法计算各个超像元之间的空间相关性;最后,依据超像元空间相关性,利用基于类别单元的类别方法将淹没区域或非淹没区域标签分配给每个亚像元中,得到最终的洪水淹没制图结果。利用两个Landsat 8 OLI多光谱图像对该方法进行了评价。结果表明,与传统的SRFIM方法相比,本文提出的SSSC-SRFIM方法具有更好的效果。     

基于AutoCAD的遥感影像处理方法及其在电力工程中的应用

卫星遥感和低空遥感等多平台多传感器遥感技术的快速发展,使得多层次遥感影像数据越来越丰富,越来越廉价,使其在电力工程的前期选址选线、可行性研究、施工图等设计过程中的使用越来越广泛。针对大多数电力工程主要采用AutoCAD平台而丢失了遥感影像的地理信息(图像的地理坐标,分辨率等),本文提出了一种新的方法,实现AutoCAD平台对遥感影像带地理参考插入,并以插件的形式实现与系统平台的对接,实现了将遥感影像引入到现有基于AuotoCAD电力工程设计平台之中,简化制图流程,提高设计水平,达到一个平台即能完成所有操作的目的,并通过实际工程应用验证了方法的可靠与高效。     

MODIS与TM冬小麦种植面积遥感测量一致性研究——小区域实验研究

针对大范围冬小麦种植面积遥感测量业务化运行中所存在的不同空间分辨率遥感影像的获取能力、空间分辨率与测量精度之间相互制约的现实问题,选择典型实验区,采取全覆盖的多时相低分辨率MODIS数据与中高分辨率TM样区数据相结合的方法,以支持向量机（SVM）为主要技术手段,通过选取不同比例的TM样本量,对MODIS进行混合像元分解,并对MODIS冬小麦测量结果与TM测量结果进行一致性分析,进而提出了一套可用于大范围冬小麦种植面积遥感测量业务化推广的识别和精度检验方法。研究结果表明：当TM样本量大于5%时,相对于TM识别结果而言,MODIS像元精度可以稳定在82%以上;当TM样本量大于40%时,区域精度可以稳定在97%。因此,从理论上讲,在实际业务化运行中,只要能够获得监测区40%的中高分辨率的影像,采用本研究提出的多尺度冬小麦种植面积测量方法,基本上可以满足业务化运行的精度要求。     

A Method to Determine Appropriate Spatial Resolution for Hard Image Classification

This study analyzed the relationship between the spatial resolution and the hard classification effect based on pixel-based image classification, and then discussed how to determine appropriate spatial resolution. Thematic maps of winter wheat derived from 250 m MODIS image, 19.5 m China-Brazil Earth Resources Satellite (CBERS) image, and 2.44 m QuickBird image were used to examine the classification effect as a case study. It indicated that the “Pareto Boundaries” and the “within-class variability” could be used to determine the coarsest and the highest resolution for hard classification, respectively. The methods proposed in this study should be useful to guide how to select appropriate spatial resolution for land cover mapping.     

Optimization of spectral indices and long-term separability analysis for classification of cereal crops using multi-spectral RapidEye imagery

Crop monitoring using remotely sensed image data provides valuable input for a large variety of applications in environmental and agricultural research. However, method development for discrimination between spectrally highly similar crop species remains a challenge in remote sensing. Calculation of vegetation indices is a frequently applied option to amplify the most distinctive parts of a spectrum. Since no vegetation index exist, that is universally best-performing, a method is presented that finds an index that is optimized for the classification of a specific satellite data set to separate two cereal crop types. The η² (eta-squared) measure of association – presented as novel spectral separability indicator – was used for the evaluation of the numerous tested indices. The approach is first applied on a RapidEye satellite image for the separation of winter wheat and winter barley in a Central German test site. The determined optimized index allows a more accurate classification (97%) than several well-established vegetation indices like NDVI and EVI (<87%). Furthermore, the approach was applied on a RapidEye multi-spectral image time series covering the years 2010–2014. The optimized index for the spectral separation of winter barley and winter wheat for each acquisition date was calculated and its ability to distinct the two classes was assessed. The results indicate that the calculated optimized indices perform better than the standard indices for most seasonal parts of the time series. The red edge spectral region proved to be of high significance for crop classification. Additionally, a time frame of best spectral separability of wheat and barley could be detected in early to mid-summer.     

Winter wheat yield forecasting in Ukraine based on Earth observation,meteorological data and biophysical models

Ukraine is one of the most developed agriculture countries and one of the biggest crop producers in the world. Timely and accurate crop yield forecasts for Ukraine at regional level become a key element in providing support to policy makers in food security. In this paper, feasibility and relative efficiency of using moderate resolution satellite data to winter wheat forecasting in Ukraine at oblast level is assessed. Oblast is a sub-national administrative unit that corresponds to the NUTS2 level of the Nomenclature of Territorial Units for Statistics (NUTS) of the European Union. NDVI values were derived from the MODIS sensor at the 250 m spatial resolution. For each oblast NDVI values were averaged for a cropland map (Rainfed croplands class) derived from the ESA GlobCover map, and were used as predictors in the regression models. Using a leave-one-out cross-validation procedure, the best time for making reliable yield forecasts in terms of root mean square error was identified. For most oblasts, NDVI values taken in April–May provided the minimum RMSE value when comparing to the official statistics, thus enabling forecasts 2–3 months prior to harvest. The NDVI-based approach was compared to the following approaches: empirical model based on meteorological observations (with forecasts in April–May that provide minimum RMSE value) and WOFOST crop growth simulation model implemented in the CGMS system (with forecasts in June that provide minimum RMSE value). All three approaches were run to produce winter wheat yield forecasts for independent datasets for 2010 and 2011, i.e. on data that were not used within model calibration process. The most accurate predictions for 2010 were achieved using the CGMS system with the RMSE value of 0.3 t ha⁻¹ in June and 0.4 t ha⁻¹ in April, while performance of three approaches for 2011 was almost the same (0.5–0.6 t ha⁻¹ in April). Both NDVI-based approach and CGMS system overestimated winter wheat yield comparing to official statistics in 2010, and underestimated it in 2011. Therefore, we can conclude that performance of empirical NDVI-based regression model was similar to meteorological and CGMS models when producing winter wheat yield forecasts at oblast level in Ukraine 2–3 months prior to harvest, while providing minimum requirements to input datasets.     

利用Landsat TM和IRS P6 LISS III卫星影像及GIS对耕地和耕种模式的时空动态分析(英文)

The study of the spatial patterns and temporal changes of cropland is important to understand the underlying factors and the functional effects of the agricultural landscape. On the other hand, crop dynamics mapping is essential to know the overall agro-spatial diversity of the area. Therefore, this paper addressed a spatio-temporal analysis of cropland and cropping pattern change in the Bogra district of Bangladesh over the last 16 years (between 1988/89 and 2004/05). In this paper, crop mapping from multi-temporal and multi-sensor satellite images was described. Landsat TM and IRS P6 LISS III satellite images were used with GIS for spatial dynamics of cropland and cropping pattern change analysis. First, seasonal cropland maps were derived from object-based classification of satellite images, then two-date classified image differencing with GIS overlay technique and decision rules were applied. Cropping pattern change was analyzed in a spatial and quantitative way for the 16 years and for this, Integrated Land and Water Information System (ILWIS) and Land Change Modular (LCM) of IDRISI Andes were used. The results showed that in the area, mono crop cultivation was found in summer, but in winter, areas under different crop cultivation had changed dramatically. Change analysis showed that the changes mainly occurred in the north northwest and southwest of the areas, and during the time the highest change area was found under the rice-potato pattern.      

非监督分类的冬小麦种植信息提取模型

为了解决在区域冬小麦种植信息遥感提取过程中监督学习算法存在的需要地面样本数据支持、流程复杂、人为干扰因素多及自动化程度低等问题,本文以非监督分类为核心,结合多尺度技术,提出了一种新的非监督分类冬小麦种植信息提取模型。选取河北省辛集市为典型试验区,以2014年高分一号数据为数据源,对本文提出的模型进行实例验证。试验结果表明：该模型的Kappa系数为0.88,整体精度为94.00%;对于研究区内的冬小麦,在无需训练样本、人为干扰因素少等条件下,该模型具有与监督分类相似的提取精度,能够满足冬小麦种植信息地面遥感监测的需求。     

基于时间序列叶面积指数稀疏表示的作物种植区域提取

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