利用遥感指数时间序列轨迹监测森林扰动

作为陆地生态系统的主体,森林的碳循环与碳蓄积对研究陆地生态系统起着重要作用,但目前森林扰动资料的缺乏在很大程度上影响着森林碳通量的估算精度。利用1986年-2011年共14期的Landsat TM/ ETM+影像,以江西武宁县为例,使用遥感指数时间序列轨迹分析方法,研究了适用于中国南方森林的扰动监测技术,该技术不仅可以识别森林的扰动变化,同时还可以监测植被的恢复信息。精度分析表明该方法得出的扰动产品的Kappa系数达到0.80,总体精度达到89.7%,表明该方法对武宁县森林扰动具有较好的监测能力。森林扰动特征分析表明武宁县森林在20世纪90年代受扰动最为剧烈,并且扰动主要发生在低海拔地区。     

光谱特征扩展的时间序列Landsat数据地表覆盖分类

时间序列遥感影像常用于地表覆盖监测及其变化监测。然而,利用时序遥感数据—尤其是中分辨率遥感数据监测地表覆盖变化,其方法基本是先对多期影像分别进行监督分类然后对比分类结果。由于这种方法需要对每期遥感影像单独选择分类训练样本,而对于历史影像,常常难以获得可靠的样本数据。本文基于遥感数据定量化处理,尝试利用光谱特征扩展方法对时间序列Landsat数据进行分类:首先,结合一种新的大气校正方法和相对辐射归一化方法,对时间序列Landsat数据进行定量化处理,以消除各期影像之间的辐射差异,获得地表反射率数据。然后,论文选择一期易于获得分类训练样本的反射率数据作为"参考影像",并结合样本数据提取不同地表覆盖类型的光谱特征。最后,将"参考影像"中提取的地物光谱特征,扩展到所有时间序列反射率数据进行分类。论文利用青藏高原玛多地区的5景Landsat数据对本文的方法进行了验证,结果显示:基于光谱特征扩展的分类方法,可有效对定量化处理后的Landsat数据进行分类,分类总体精度为88.35%—94.25%,分类结果和传统的单景监督分类结果具有较好的一致性。此外,研究也发现,"参考影像"和待分类图像获取时间的季相差异会影响其分类的精度。     

长时间序列多源遥感数据的森林干扰监测算法研究进展

时空意义明确的森林干扰和恢复信息是评价森林生态系统碳动态的关键因素之一。然而由于诸多的现实困难,多尺度的森林干扰定量化时空信息相对缺乏。Landsat数据具备光谱、时间和空间分辨率上的优势,以及可以免费获取的特点,使其成为主要的长时间序列动态监测的遥感数据源之一,为长时间周期内提供具有合适的空间细节和时间频率的森林干扰信息成为可能。特别是基于Landsat时间序列堆栈(LTSS)的森林干扰自动分析算法的出现,更为森林生态系统的近实时监测提供强有力的工具。本文全面评述了长时间序列遥感数据准备和预处理技术以及国内外基于遥感数据源的多时相森林干扰监测方法,重点分析了基于Landsat的多种指数监测和自动化方法的优缺点,并总结了其与多源数据结合的扩展应用,最后就现有方法与国内外新的数据、技术手段的关联进行了展望,以期为推广中国本土卫星影像应用于森林干扰监测提供理论借鉴。     

多模型耦合下的植被冠层可燃物含水率遥感反演

植被冠层可燃物含水率FMC(Fuel Moisture Content)是评估野火风险及估算火灾蔓延速率的重要指标。以中国西部6个典型研究区为例,基于辐射传输模型,利用Landsat 5 TM及Landsat 8 OLI遥感数据,开展草原、森林冠层FMC定量反演研究。为克服基于物理模型的病态反演问题、FMC自身的弱敏感性问题及西南森林多具复杂的双层冠层结构问题,研究中考虑了模型参数之间的相关特征,使用多波段遥感数据及耦合辐射传输模型等方法。反演结果显示,总体植被冠层FMC反演精度R~2为0.64,RMSE为44.86%,其中草地冠层FMC的反演精度(R~2=0.64,RMSE=47.57%)略低于森林冠层FMC的反演精度(R~2=0.71,RMSE=30.82%)。为进一步论证该反演结果对野火风险评估的有效性,研究中选取并分析了2011年3月2日于云南大理白族自治州剑川县金华镇金和村森林火灾爆发前、爆发时及灾后该区域植被冠层FMC的变化特征。结果显示,火灾爆发时该地区植被冠层FMC明显低于火灾发生前后(约一月时间)植被冠层FMC,证明了本文FMC反演结果对野火风险评估的有效性。     

Landsat系列卫星地表反射率产品研究进展

Landsat系列卫星为地球资源环境动态监测提供了长达40余年的中高分辨率卫星影像,近年来USGS等实施的Landsat数据共享计划使Landsat系列卫星数据成为应用最广泛的中高分辨率卫星数据。长期以来,Landsat系列卫星数据以灰度值(DN值)的形式提供给用户,受传感器波段响应、具体成像条件(太阳高度角、大气散射和吸收等)差异等的影响,不同传感器、不同时间、不同地点的卫星影像DN值不具有可比性,给Landsat系列卫星数据的定量应用造成了障碍。近年来国内外多家研究机构陆续推出了Landsat系列卫星地表反射率产品,以增强Landsat长时间系列卫星数据的可比性。本文对USGS、马里兰大学、WELD团队和中国科学院等多家机构推出的Landsat地表反射率产品进行了介绍,并对未来的研究方向进行了展望。     

塞罕坝地区高空间分辨率叶面积指数时序估算与变化检测

叶面积指数LAI (Leaf Area Index)是调节植被冠层生理过程的最重要的生物物理变量之一,高空间分辨率时间序列LAI对于植被生长检测、地表过程模拟与区域和全球变化研究至关重要,但是由于数据缺失和反演方法限制,目前还没有时空连续的高分辨率LAI数据产品。本研究提出了一种生成时间连续的高空间分辨率LAI数据的算法,首先对MODIS LAI产品滤波平滑,生成时间序列LAI的上包络曲线,根据上包络曲线提供的变化信息构建LAI动态模型。然后利用地面实测的LAI数据与Landsat反射率数据构建LAI反演的BP (Back Propagation)神经网络模型。将反演得到的高分辨率LAI数据作为LAI观测数据,利用集合卡尔曼滤波EnKF (Ensemble Kalman Filter)方法实时更新动态模型,生成时间连续的30 m空间分辨率LAI数据集。基于该算法生成了塞罕坝地区2000年—2018年长时间序列LAI数据集,利用Prophet深度学习模型进行模拟和预测,根据预测和原始LAI差异,利用支持向量机SVM (Support Vector Machine)方法检测植被干扰状况。结果表明:EnKF算法能够生成时空连续的高空间分辨率LAI数据,估算结果与地面测量值一致性较高,R~2为0.9498,RMSE为0.1577,在区域尺度上与Landsat LAI参考值较为吻合,R~2高于0.87,RMSE低于0.61。Prophet与SVM模型检测到研究区2009年,2010年,2013年,2014年,2015年植被受干扰较为严重,主要由于年降水量偏少和林区作业砍伐造成,检测结果与当地降水量与砍伐数据吻合。本文提出的算法可用于大范围高时空LAI数据反演和植被变化检测,对塞罕坝乃至全国林区规划管理具有重要的参考价值。     

集成Landsat OLI和机载LiDAR条带数据的亚热带森林生物量制图

以亚热带天然次生林为研究对象,借助一个条带的少量LiDAR点云数据和覆盖整个研究区的免费Landsat OLI多光谱数据,并结合地面实测数据,探索森林生物量低成本高精度制图方法。首先,提取了OLI和LiDAR特征变量,并与地上和地下生物量进行相关分析以筛选变量;然后,借助LiDAR数据覆盖区的样地和条带LiDAR数据构建"LiDAR生物量模型";再从LiDAR反演生物量的结果中进行采样,结合OLI特征变量构建"LiDAR-OLI模型";最后,与单独使用OLI多光谱数据建立的"OLI估算模型"结果进行比较,分析精度并验证新方法的效果。结果表明,"LiDAR-OLI模型"对地上和地下生物量的模型拟合效果较好且均优于"OLI模型",且其交叉验证的精度也较高并优于"OLI模型",从而证明了新方法的可靠性及有效性。本研究为主、被动遥感技术在中小尺度上协同反演森林参数提供了实验基础,也为基于全覆盖免费OLI多光谱数据及条带LiDAR数据的低成本森林生物量制图探索了技术路线。     

时序遥感影像的时空特征在土地分类中的应用

基于Landsat时间序列数据的土地覆盖检测成为当前研究热点,但基于时间序列数据空间纹理特征的应用及不同时序特征重要性评估较少.基于时间序列Landsat8数据,在时序光谱特征、指数特征和地形特征基础上引入时序纹理特征,利用随机森林算法建立八种分类模型,对北京密云区进行土地覆盖分类并比较其分类精度,进而基于袋外(OOB)误差方法和基尼指数进行特征变量重要性评估.相比加入归一化建筑指数(Normalized Difference Built-up Index,NDBI)或归一化植被指数(Normalized Difference Vegetation Index,NDVI)样本特征,时序纹理特征的加入使总体精度分别提高1.88％和2.12％;最优分类模型中灰度共生矩阵(Gray-level Co-occurrence Matrix,GLCM)熵参数在纹理特征中较为重要,GLCM差异性参数和GLCM相关性参数其次.本文为进一步挖掘影像的时空特征、提高土地覆盖制图精度提供新思路.     

Landsat长时间序列数据格式统一与反射率转换方法实现

介绍了一种长时间序列遥感影像预处理程序,即陆地卫星生态系统干扰自适应处理系统(landsat ecosystem disturbance adaptive processing system,LEDAPS)。该程序通过使用MODTRAN太阳能输出模型,校正太阳方位、日地距离、TM或ETM+带通以及太阳辐照度,将定标影像转换为表观(top-of-atmosphere,TOA)反射率影像,并将通过浓密植被(dark dense vegetation,DDV)算法插值生成的气溶胶光学厚度(aerosol optical thickness,AOT)以及通过相关资料获得的臭氧(O3)浓度、大气压及水汽值等用于6S辐射传输模型,生成地表反射率产品。以LEDAPS可处理的标准数据Landsat7 ETM+和统一格式后的非标准数据Landsat5 TM影像为例,介绍了长时间(1987—2011年)序列数据的选择、格式统一以及算法的实现过程,同时给出了校正后影像效果评价的方法。结果表明,标准数据和非标准数据经过LEDAPS处理后生成的地表反射率产品能有效降低大气中O3、水汽及气溶胶等对影像真实反射率的影响,为土地覆盖变化和干扰因素等的长时间序列监测和生物物理参数的遥感反演提供科学产品,有助于在国内形成处理长时间序列影像数据的准则。     

应用随机辐射传输模型反演云南松林分郁闭度

随机辐射传输模型可用于模拟水平分布不均一森林的辐射传输过程。本文以云南松林分为研究对象，提出一种应用随机辐射传输模型的郁闭度反演方法。该方法以随机辐射传输模型中参数与林分郁闭度的定量关系为基础，提出了针对云南松的冠型等效模型，构建了郁闭度和卫星反射率（GF-1和Landsat 8卫星影像）的查找表，并实施了反演。基于野外实测的30个样地进行了郁闭度数据验证，并和基于NDVI回归模型的反演方法进行对比。结果表明，反演结果能够较准确反映云南松林分郁闭状况（R²=0.8345，RMSE=0.0688），通过冠型修正能够降低反演误差，冠型等效模型是合理的。反演方法机理清晰且适用范围广，研究成果可为大面积森林郁闭度反演提供模型和方法支持。     

改进的VCT长时间序列森林干扰方法

针对传统的VCT(Vegetation Change Tracker)模型进行了改进,通过引入双季节数据,提高了森林干扰提取精度。以湖南省株洲市天元区和渌口区为例,基于Google Earth Engine搭建了VCT模型,构建了长时间序列Landsat 5/8卫星数据堆栈(LTSS),实现了研究区1989—2019年以来的森林干扰历史重建,分析了研究区森林干扰时空分布特征和变化趋势等。实验结果表明:基于双季节数据的VCT算法可有效提取森林干扰信息,干扰提取的总体精度达到81.33%,较传统的VCT方法提高了4.44%;研究区1989—2019年干扰面积总体呈上升趋势,干扰主要体现在城镇扩张、农业开垦和交通网建设3方面。实验结果表明了该算法的有效性,可应用于我国南方森林干扰提取研究。     

North Carolina’s forest disturbance and timber production assessed using time series Landsat observations

Wood products provide a relatively long-term carbon storage mechanism. Due to lack of consistent datasets on these products, however, it is difficult to determine their carbon contents. The main purpose of this study was to quantify forest disturbance and timber product output (TPO) using time series Landsat observations for North Carolina. The results revealed that North Carolina had an average forest disturbance rate of 178,000 ha per year from 1985 to 2010. The derived disturbance products were found to be highly correlated with TPO survey data, explaining up to 87% of the total variance of county level industrial roundwood production. State level TPO estimates derived using the Landsat-based disturbance products tracked those derived from ground-based survey data closely. The TPO modeling approach developed in this study complements the ground-based TPO surveys conducted by the US Forest Service. It allows derivation of TPO estimates for the years that did not have TPO survey data, and may be applicable in other regions or countries where at least some ground-based survey data on timber production are available for model development and dense time series Landsat observations exist for developing annual forest disturbance products.     

Forest cover classification using Landsat ETM+ data and time series MODIS NDVI data

Forest cover plays a key role in climate change by influencing the carbon stocks, the hydrological cycle and the energy balance. Forest cover information can be determined from fine-resolution data, such as Landsat Enhanced Thematic Mapper Plus (ETM+). However, forest cover classification with fine-resolution data usually uses only one temporal data because successive data acquirement is difficult. It may achieve mis-classification result without involving vegetation growth information, because different vegetation types may have the similar spectral features in the fine-resolution data. To overcome these issues, a forest cover classification method using Landsat ETM+ data appending with time series Moderate-resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) data was proposed. The objective was to investigate the potential of temporal features extracted from coarse-resolution time series vegetation index data on improving the forest cover classification accuracy using fine-resolution remote sensing data. This method firstly fused Landsat ETM+ NDVI and MODIS NDVI data to obtain time series fine-resolution NDVI data, and then the temporal features were extracted from the fused NDVI data. Finally, temporal features combined with Landsat ETM+ spectral data was used to improve forest cover classification accuracy using supervised classifier. The study in North China region confirmed that time series NDVI features had significant effects on improving forest cover classification accuracy of fine resolution remote sensing data. The NDVI features extracted from time series fused NDVI data could improve the overall classification accuracy approximately 5% from 88.99% to 93.88% compared to only using single Landsat ETM+ data.     

知识引导的稀疏时间序列遥感数据拟合

在多云多雨的地区,光学遥感存在着获取无云数据困难的难题,这会导致时间序列应用中可用数据匮乏。因此,本文面向稀疏时间序列遥感数据,根据噪声造成遥感影像上归一化差分植被指数(NDVI)被低估的事实,提出了一种知识引导的拟合方法。首先,在遥感影像预处理的基础上,利用先验知识和时序差分法对噪声进行识别和剔除;然后,采用高斯二阶模型对原始数据进行拟合;最后,根据拟合残差更新权重,进行迭代拟合,重复上述过程直至获得稳定的结果。本文以Landsat 8 OLI作为数据源,对浙江省杭州地区的森林数据进行拟合,结果表明:在稀疏时间序列数据的情况下,本文方法与MODIS数据拟合结果的相关系数达到0.92,关键时点(如NDVI峰值点等)的时间误差在5 d;相比当前主流方法的0.88与8 d具有更高的精度。     

Normalized algorithm for mapping and dating forest disturbances and regrowth for the United States

Forest disturbances such as harvesting, wildfire and insect infestation are critical ecosystem processes affecting the carbon cycle. Because carbon dynamics are related to time since disturbance, forest stand age that can be used as a surrogate for major clear-cut/fire disturbance information has recently been recognized as an important input to forest carbon cycle models for improving prediction accuracy. In this study, forest disturbances in the USA for the period of ∼1990–2000 were mapped using 400+ pairs of re-sampled Landsat TM/ETM scenes in 500m resolution, which were provided by the Landsat Ecosystem Disturbance Adaptive Processing System project. The detected disturbances were then separated into two five-year age groups, facilitated by Forest Inventory and Analysis (FIA) data, which was used to calculate the area of forest regeneration for each county in the USA.     

Forest cover trends from time series Landsat data for the Australian continent

In perennial and natural vegetation systems, monitoring changes in vegetation over time is of fundamental interest for identifying and quantifying impacts of management and natural processes. Subtle changes in vegetation cover can be identified by calculating the trends of a vegetation density index over time. In this paper, we apply such an index-trends approach, which has been developed and applied to time series Landsat imagery in rangeland and woodland environments, to continental-scale monitoring of disturbances within forested regions of Australia. This paper describes the operational methods used for the generation of National Forest Trend (NFT) information, which is a time-series summary providing visual indication of within-forest vegetation changes (disturbance and recovery) over time at 25 m resolution. This result is based on a national archive of calibrated Landsat TM/ETM+ data from 1989 to 2006 produced for Australia's National Carbon Accounting System (NCAS). The NCAS was designed in 1999 initially to provide consistent fine-scale classifications for monitoring forest cover extent and changes (i.e. land use change) over the Australian continent using time series Landsat imagery. NFT information identifies more subtle changes within forested areas and provides a capacity to identify processes affecting forests which are of primary interest to ecologists and land managers. The NFT product relies on the identification of an appropriate Landsat-based vegetation cover index (defined as a linear combination of spectral image bands) that is sensitive to changes in forest density. The time series of index values at a location, derived from calibrated imagery, represents a consistent surrogate to track density changes. To produce the trends summary information, statistical summaries of the index response over time (such as slope and quadratic curvature) are calculated. These calculated index responses of woody vegetation cover are then displayed as maps where the different colours indicate the approximate timing, direction (decline or increase), magnitude and spatial extent of the changes in vegetation cover. These trend images provide a self-contained and easily interpretable summary of vegetation change at scales that are relevant for natural resource management (NRM) and environmental reporting.     

Update and spatial extension of strategic forest inventories using time series remote sensing and modeling

Up-to-date forest inventory information relating the characteristics of managed and natural forests is fundamental to sustainable forest management and required to inform conservation of biodiversity and assess climate change impacts and mitigation opportunities. Strategic forest inventories are difficult to compile over large areas and are often quickly outdated or spatially incomplete as a function of their long production cycle. As a consequence, automated approaches supported by remotely sensed data are increasingly sought to provide exhaustive spatial coverage for a set of core attributes in a timely fashion. The objective of this study was to demonstrate the integration of current remotely-sensed data products and pre-existing jurisdictional inventory data to map four forest attributes of interest (stand age, dominant species, site index, and stem density) for a 55 Mha study region in British Columbia, Canada. First, via image segmentation, spectrally homogenous objects were derived from Landsat surface-reflectance pixel composites. Second, a suite of Landsat-based predictors (e.g., spectral indices, disturbance history, and forest structure) and ancillary variables (e.g., geographic, topographic, and climatic) were derived for these units and used to develop predictive models of target attributes. For the often difficult classification of dominant species, two modelling approaches were compared: (a) a global Random Forests model calibrated with training samples collected over the entire study area, and (b) an ensemble of local models, each calibrated with spatially constrained local samples. Accuracy assessment based upon independent validation samples revealed that the ensemble of local models was more accurate and efficient for species classification, achieving an overall accuracy of 72% for the species which dominate 80% of the forested areas in the province. Results indicated that site index had the highest agreement between predicted and reference (R² = 0.74, %RMSE = 23.1%), followed by stand age (R² = 0.62, %RMSE = 35.6%), and stem density (R² = 0.33, %RMSE = 65.2%). Inventory attributes mapped at the image-derived unit level captured much finer details than traditional polygon-based inventory, yet can be readily reassembled into these larger units for strategic forest planning purposes. Based upon this work, we conclude that in a multi-source forest monitoring program, spatially localized and detailed characterizations enabled by time series of Landsat observations in conjunction with ancillary data can be used to support strategic inventory activities over large areas.     

兼容多源数据的森林湿地信息提取和比较分析

针对以光谱特征差异为依据,提取森林湿地信息精度低的问题,该文采用兼容多源数据的分类回归树(CART)提取方法,并以大沾河国家森林湿地进行实证研究。基于Landsat8遥感数据、Radarsat-2极化雷达数据和地形辅助数据,采用SPM软件分别构建3种特征变量组合的CART决策树模型,并获取分类规则,最后根据规则对研究区的森林湿地信息进行提取。结果表明:3种特征变量组合中,兼容光谱、纹理、雷达与地形辅助数据的CART决策树的森林湿地信息提取精度最高,用户精度和制图精度分别达到了88.46%和82.14%。研究结果体现了雷达数据与地形辅助数据有助于提取森林湿地信息。