利用Sentinel-2影像超分辨率重建的红树林冠层氮含量反演

氮素是植被整个生命周期的必要元素,红树林冠层氮素含量（CNC）遥感估算对红树林健康监测具有重要意义。以广东湛江高桥红树林保护区为研究区,本文旨在基于Sentinel-2影像超分辨率重建技术进行红树林CNC估算和空间制图。研究首先基于三次卷积重采样、Sen2Res和SupReMe算法实现Sentinel-2影像从20 m分辨率到10 m的重建;然后以重建后的影像和原始20 m影像为数据源构建40个相关植被指数,采用递归特征消除法（SVM-RFE）确定CNC估算的最优变量组合,进而构建CNC反演的核岭回归（KRR）模型;最后选取最优模型实现CNC制图。研究结果表明：基于Sen2Res和SupReMe超分辨率算法的重建影像不仅与原始影像具有很高的光谱一致性,且明显提高了影像的清晰度和空间细节。红树林CNC反演波段主要集中在红（B4）、红边（B5）、近红外波段（B8a）以及短波红外波段（B11和B12）,与“红边波段”相关的植被指数（RSSI和TCARIre1/OSAVI）也是红树林CNC反演的有效变量。基于3种方法重建后10 m的影像构建的模型反演精度（R²val>0.579）均优于原始20 m的影像（R²val=0.504）;基于Sen2Res算法重建影像构建的反演模型拟合精度（R²val=0.630,RMSE_val=5.133,RE_val=0.179）与基于三次卷积重采样重建影像的模型拟合精度（R²val=0.640,RMSE_val=5.064,RE_val=0.179）基本相当,前者模型验证精度（R²cv=0.497,RMSE_cv=5.985,RE_cv=0.214）较高且模型变量选择数量最为合理。综合重建影像光谱细节及模型精度,基于Sen2Res算法重建的Sentinel-2影像在红树林CNC估算中具有良好的应用潜力,能为区域尺度红树林冠层健康状况的精细监测提供有效的方法借鉴和数据支撑。     

Sentinel-2A与Landsat 8O LI逐像元辐射归一化方法研究

考虑不同传感器光谱响应函数差异及不同地物类型反射率光谱的差异,提出了一种逐像元辐射归一化方法,并以2017年7月17日内蒙古达里诺尔湖地区准同步过境的Sentinel-2A及Landsat 8数据为例,对两类数据可见-近红外波段（VNIR）地表反射率结果进行归一化。首先采用Sen2cor方法及NASA官方提供大气校正算法,分别对Sentinel-2A及Landsat 8 OLI影像进行大气校正并重采样到同一空间分辨率;然后基于光谱库计算匹配因子并构建图像与光谱库之间的匹配转换模型,实现像元尺度上从Sentinel-2影像到Landsat 8影像地表反射率相似波段之间的转换。结果表明,经逐像元归一化的影像相比原始影像及经HLS光谱归一化的影像,与Landsat 8 VNIR波段的相关性明显提高,辐射一致性增强。该转换模型为多源中高分辨率遥感图像高精度辐射归一化提供了新思路。     

Sentinel-2数据的冬小麦地上干生物量估算及评价

作物生物量快速精确的监测对于农业资源的合理利用与农田的精准管理具有重要意义。近年来,遥感技术因其独特的优势已被广泛用于作物生物量的估算中。本文主要针对不同宽波段植被指数在冬小麦生物量(文中的生物量均是指地上干生物量)估算方面的表现进行探索。首先利用欧洲空间局最新的Sentinel-2A卫星数据提取出17种常见的植被指数,之后分别构建其与相应时期内采集的冬小麦地上生物量间的最优估算模型,通过分析两者间的相关性与敏感性,获取适宜进行生物量估算的指数。最后,绘制了研究区的生物量空间分布图。结果表明,所选的植被指数均与生物量显著相关。其中,红边叶绿素指数(CI_(re))与生物量的估算精度最高(决定性系数R~2为0.83;均方根误差RMSE为180.29 g·m~(–2))。虽然相关性较高,但部分指数,如归一化差值植被指数(NDVI)等在生物量较高时会出现饱和现象,从而导致生物量的低估。而加入红边波段的指数不仅能够延缓指数的饱和趋势,而且能够提高反演精度。此外,通过敏感性分析发现,归一化差值指数和比值指数分别在作物生长的早期和中后期对生物量的变化保持较高的敏感性。由于红边比值指数(SR_(re))和MERIS叶绿素敏感指数(MTCI)在冬小麦全生长季内一直对生物量的变化保持高灵敏性,二者是生物量估算中最为稳定的指数。     

利用Sentinel-2A多光谱成像仪与Landsat8陆地成像仪影像进行普陀山岛植被分类效果比较

卫星遥感技术可用于海岛资源调查。Sentinel-2A与Landsat 8两颗卫星都可免费提供空间分辨率较高的多光谱遥感影像,在海岛调查中的应用潜力较大。本文以浙江舟山普陀山岛为例开展了针对这两种影像在海岛植被分类中的应用效果的研究,分别利用Sentinel-2A多光谱成像仪（MSI）和Landsat 8陆地成像仪（OLI）影像基于最大似然法分类获得了该岛阔叶林、针阔混交林、针叶林、灌丛、草丛等植被及其他地物的分布情况,并进行了精度检验,结果表明MSI的总体分类精度略高于OLI。     

Spectral unmixing based spatiotemporal downscaling fusion approach

Time-series remote sensing data are important in monitoring land surface dynamics. Due to technical limitations, satellite sensors have a trade-off between temporal, spatial and spectral resolutions when acquiring remote sensing images. In order to obtain remote sensing images with high spatial resolution and high temporal frequency, spatiotemporal fusion methods have been developed. In this paper, we propose a Linear Spectral Unmixing-based Spatiotemporal Data Fusion Model (LSUSDFM) for spatial and temporal data fusion. In this model, the endmember abundance of the low-resolution image pixel is calculated based on that of the high-resolution image by the spectral mixture analysis. The endmember spectrum signals of low-resolution images are then calculated continuously within an optimized moving window. Subsequently, the fused image is reconstructed according to the endmember spectrum and its corresponding abundance map. A simulated dataset and real satellite images are used to test the fusion model, and the fusion results are compared with a current spectral unmixing based downscaling fusion model (SUDFM). Our experimental work shows that, compared to the SUDFM, the proposed LSUSDFM can achieve better quality and accuracy of fused images, especially in effectively eliminating the “plaque” phenomenon in the results by the SUDFM. The LSUSDFM has great potential in generating images with both high spatial resolution and high temporal frequency, as well as increasing the number of spectral bands of the high spatial resolution data.     

Visualizing the ill-posedness of the inversion of a canopy radiative transfer model: A case study for Sentinel-2

Monitoring biophysical and biochemical vegetation variables in space and time is key to understand the earth system. Operational approaches using remote sensing imagery rely on the inversion of radiative transfer models, which describe the interactions between light and vegetation canopies. The inversion required to estimate vegetation variables is, however, an ill-posed problem because of variable compensation effects that can cause different combinations of soil and canopy variables to yield extremely similar spectral responses. In this contribution, we present a novel approach to visualise the ill-posed problem using self-organizing maps (SOM), which are a type of unsupervised neural network. The approach is demonstrated with simulations for Sentinel-2 data (13 bands) made with the Soil-Leaf-Canopy (SLC) radiative transfer model. A look-up table of 100,000 entries was built by randomly sampling 14 SLC model input variables between their minimum and maximum allowed values while using both a dark and a bright soil. The Sentinel-2 spectral simulations were used to train a SOM of 200 × 125 neurons. The training projected similar spectral signatures onto either the same, or contiguous, neuron(s). Tracing back the inputs that generated each spectral signature, we created a 200 × 125 map for each of the SLC variables. The lack of spatial patterns and the variability in these maps indicate ill-posed situations, where similar spectral signatures correspond to different canopy variables. For Sentinel-2, our results showed that leaf area index, crown cover and leaf chlorophyll, water and brown pigment content are less confused in the inversion than variables with noisier maps like fraction of brown canopy area, leaf dry matter content and the PROSPECT mesophyll parameter. This study supports both educational and on-going research activities on inversion algorithms and might be useful to evaluate the uncertainties of retrieved canopy biophysical and biochemical state variables.     

Fusion of multispectral and panchromatic Satellite images using the curvelet transform

A useful technique in various applications of remote sensing involves the fusion of different types of satellite images, namely multispectral (MS) satellite images with a high spectral and low spatial resolution and panchromatic (Pan) satellite image with a low spectral and high spatial resolution. Recent studies show that wavelet-based image fusion provides high-quality spectral content in fused images. However, the results of most wavelet-based methods of image fusion have a spatial resolution that is less than that obtained via the Brovey, intensity-hue-saturation, and principal components analysis methods of image fusion. We introduce an improved method of image fusion which is based on the amelioration de la resolution spatiale par injection de structures (ARSIS) concept using the curvelet transform, because the curvelet transform represents edges better than wavelets. Because edges are fundamental in image representation, enhancing the edges is an effective means of enhancing spatial resolution. Curvelet-based image fusion has been used to merge a Landsat Enhanced Thematic Mapper Plus Pan and MS image. The proposed method simultaneously provides richer information in the spatial and spectral domains.     

Remote estimation of crop and grass chlorophyll and nitrogen content using red-edge bands on Sentinel-2 and -3

Sentinel-2 is planned for launch in 2014 by the European Space Agency and it is equipped with the Multi Spectral Instrument (MSI), which will provide images with high spatial, spectral and temporal resolution. It covers the VNIR/SWIR spectral region in 13 bands and incorporates two new spectral bands in the red-edge region, which can be used to derive vegetation indices using red-edge bands in their formulation. These are particularly suitable for estimating canopy chlorophyll and nitrogen (N) content. This band setting is important for vegetation studies and is very similar to the ones of the Ocean and Land Colour Instrument (OLCI) on the planned Sentinel-3 satellite and the Medium Resolution Imaging Spectrometer (MERIS) on Envisat, which operated from 2002 to early 2012. This paper focuses on the potential of Sentinel-2 and Sentinel-3 in estimating total crop and grass chlorophyll and N content by studying in situ crop variables and spectroradiometer measurements obtained for four different test sites. In particular, the red-edge chlorophyll index (CI_red-edge), the green chlorophyll index (CI_green) and the MERIS terrestrial chlorophyll index (MTCI) were found to be accurate and linear estimators of canopy chlorophyll and N content and the Sentinel-2 and -3 bands are well positioned for deriving these indices. Results confirm the importance of the red-edge bands on particularly Sentinel-2 for agricultural applications, because of the combination with its high spatial resolution of 20 m.     

融合高时空分辨率数据估算植被净初级生产力

植被净初级生产力NPP(Net Primary Production)遥感估算与分析,有赖于高时空分辨率的遥感数据,但目前中高分辨率的遥感数据受卫星回访周期及天气的影响,在中国南方地区难以获取连续时间序列的数据,从而影响了高精度的区域植被净初级生产力的遥感估算。为此,提出一种基于多源遥感数据时空融合技术与CASA模型估算高时空分辨率NPP的方法。首先,利用多源遥感数据,即Landsat8 OLI数据与MODIS13Q1数据,采用遥感数据时空融合方法,获得了时间序列的Landsat8 OLI融合数据;然后,基于Landsat8 OLI时空融合数据,并采用CASA模型,以长株潭城市群核心区为例,进行区域植被NPP的遥感估算。研究结果表明,基于时间序列Landsat融合数据估算的30m分辨率的NPP具有良好的空间细节信息,且估算值与实测值的相关系数达0.825,与实测NPP数据保持了较好的一致性。     

利用Sentinel-2A影像的面向对象特征湿地决策树分类方法

苏州市湿地众多、类型多样化、周围环境复杂,使用传统的遥感分类方法很难得到精度较高的湿地分类结果。研究了面向对象特征的湿地决策树分类方法,以苏州市澄湖地区为研究区域,使用欧空局的Sentinel-2A影像,先将研究区域分为湿地水体、植被和非植被3大类型,再分别构建鱼塘、河流、湖泊、农田和裸地等面向对象特征,据此实现湿地遥感分类。研究结果表明,该方法能够有效利用遥感影像提供的光谱特征、几何特征和纹理特征等多种丰富信息,产生较高的分类精度,总体分类精度可达80.67%,Kappa系数为77.80%。与传统的基于中低分辨率遥感影像的分类方法相比,该方法可以有效提取湿地不同地物对象的几何结构和纹理等特征,在提高湿地分类精度的同时实现对大面积湿地的快速动态监测。     

混沌蜂群优化的NSST域多光谱与全色图像融合

为有效融合多光谱图像的光谱信息和全色图像的空间细节信息,提出了一种基于混沌蜂群优化和改进脉冲耦合神经网络(PCNN)的非下采样Shearlet变换(NSST)域图像融合方法。首先对多光谱图像进行Intensity-HueSaturation(IHS)变换,全色图像的直方图按照多光谱图像亮度分量的直方图进行匹配;然后分别对多光谱图像的亮度分量和新全色图像进行NSST变换,对低频分量使用改进加权融合算法进行融合,以互信息作为适应度函数,利用混沌蜂群算法找到最优加权系数。对高频分量采用改进脉冲耦合神经网络(PCNN)方法进行融合,再经NSST逆变换和IHS逆变换得到融合图像。本文方法在主观视觉效果和信息熵、光谱扭曲度等客观定量评价指标上优于基于IHS变换、基于非下采样Contourlet变换(NSCT)和非负矩阵分解(NMF)、基于NSCT和PCNN等5种融合方法。本文方法在提升图像空间分辨率的同时,有效地保留了光谱信息。     

遥感影像时—空融合的“点”—“线”—“面”质量评价

遥感影像时—空融合可集成多源数据高空间分辨率和高时间分辨率互补优势,生成时间连续的高空间分辨率影像,在遥感影像的动态监测与时序分析等方面具有重要应用价值.然而,现有多数研究往往基于单一数据产品对时—空融合算法进行评价,而在实际生产应用中,需要验证算法在多种遥感产品数据的融合表现;此外,目前研究大多基于“单点时刻”进行评...     

利用最优估计理论进行多光谱与全色影像融合

应用小波变换分解遥感影像,利用遥感影像自身的先验信息——空间分辨率确定高频域融合过程中的权值,使用最小二乘估计与小波重构完成影像融合。实验结果表明,相对于参考的其他融合方法,此方法在注入全色影像空间细节和保持多光谱影像的光谱信息方面性能更佳。     

Multitemporal fusion of Landsat/TM and ENVISAT/MERIS for crop monitoring

Monitoring Earth dynamics using current and future satellites is one of the most important objectives of the remote sensing community. The exploitation of image time series from sensors with different characteristics provides new opportunities to increase the knowledge about environmental changes and to support many operational applications. This paper presents an image fusion approach based on multiresolution and multisensor regularized spatial unmixing. The approach yields a composite image with the spatial resolution of the high spatial resolution image while retaining the spectral and temporal characteristics of the medium spatial resolution image. The approach is tested using images from Landsat/TM and ENVISAT/MERIS instruments, but is general enough to be applied to other sensor pairs. The potential of the proposed spatial unmixing approach is illustrated in an agricultural monitoring application where Landsat temporal profiles from images acquired over Albacete, Spain, in 2004 and 2009 are complemented with MERIS fused images. The resulting spatial resolution from Landsat allows monitoring small and medium size crops at the required scale while the fine spectral and temporal resolution from MERIS allow a more accurate determination of the crop type and phenology as well as capturing rapidly varying land-cover changes.     

Mapping forest windthrows using high spatial resolution multispectral satellite images

Wind disturbances represent the main source of damage in European forests, affecting them directly (windthrows) or indirectly due to secondary damages (insect outbreaks and forest fires). The assessment of windthrows damages is very important to establish adequate management plans and remote sensing can be very useful for this purpose. Many types of optical remote sensing data are available with different spectral, spatial and temporal resolutions, and many options are possible for data acquisition, i.e. immediately after the event or after a certain time. The objective of this study is to compare the windthrows mapping capabilities of two multispectral satellite constellations (i.e. Sentinel-2 and PlanetScope) characterized by very different spectral, spatial and temporal resolutions, and to evaluate the impact of the acquisition conditions on the mapping results. The analysed area, with an extent of 732 km², is located in the Trentino-South Tyrol region (Italy) which was affected by the Vaia storm on the 27th-30th of October 2018, causing serious forest damages. The change vector analysis technique was used to detect the windthrows. For each data source, two pairs of images were considered: 1) pre- and post- event images acquired as close as possible to the event; 2) pre- and post- event images acquired at optimal conditions, i.e. at similar phenological state and similar illumination conditions. The results obtained with the two satellite constellations are very similar despite their different resolutions. Data acquired in optimal conditions allowed having the best detection rate (accuracy above 80 %), while data acquired just after the event showed many limitations. Improved spatial resolution (PlanetScope data) allows for a better delineation of the borders of the windthrow areas and of the detection of smaller windthrow patches.     

Estimating mineral content of indigenous browse species using laboratory spectroscopy and sentinel-2 imagery

Trees provide low-cost organic inputs, with the potential to improve livelihoods for rural communities. Understanding foliar nutrients of tree species is crucial for integration of trees into agroecosystems. The study explored nitrogen (N), phosphorus (P), potassium (K) and calcium (Ca) concentrations of nine browse species collected from the bushveld region of South Africa using wet analysis and laboratory spectroscopy in the region 400–2500 nm, along with partial least squares (PLS) regression. We further explore the relationship between canopy reflectance of Sentinel-2 image and foliar N, P, K & Ca. Laboratory spectroscopy was significant for N estimation, while satellite imagery also revealed useful information about the estimation of nitrogen at landscape level. Nitrogen was highly correlated with spectral reflectance (R² = 0.72, p < 0.05) for winter and (R² = 0.88, p < 0.05) for summer, whilst prediction of phosphorus potassium and calcium were considered not accurate enough to be of practical use. Modelling the relationship using Sentinel-2 data showed lower correlations for nitrogen (R² = 0.44, p < 0.05) and the other nutrients when compared to the dried samples. The findings indicate that there is potential to assess and monitor resource quality of indigenous trees using nitrogen as key indicator. This multi-level remote sensing approach has promise for providing rapid plant nutrient analyses at different scales.     

Retrieving soil moisture in rainfed and irrigated fields using Sentinel-2 observations and a modified OPTRAM approach

Surface soil water content plays an important role in driving the exchange of latent and sensible heat between the atmosphere and land surface through transpiration and evaporation processes, regulating key physiological processes affecting plants growth. Given the high impact of water scarcity on yields, and of irrigated agriculture on the overall withdrawal rate of freshwater, it is important to define models that help to improve water resources management for agricultural purposes, and to optimize rainfed crop yield. Recent advances in satellite-based remote sensing have led to valuable solutions to estimate soil water content based on microwave or optical/thermal-infrared data. This study aims at improving soil water content estimation at high spatial and temporal resolution, by means of the Optical Trapezoid Model (OPTRAM) driven by Copernicus Sentinel-2 data. Two different model variations were considered, based on linear and nonlinear parameters constraints, and validated against in situ soil water content measurements made with time domain reflectometry (TDR) on irrigated maize in central Italy and on rainfed maize and pasture in northern Italy. For the first site the non-linear model shows a better correlation between measured and estimated soil water content values (r = 0.80) compared to the linear model (r = 0.73). In both cases the modeled soil moisture tends to overestimate the measured values at medium to high water content level, while both models underestimate soil moisture at low water content level. Estimated versus measured normalized surface soil water for rainfed pasture plots from nonlinear OPTRAM parametrized based on irrigated maize parameterization (SIM1), and site-specific parametrization for rainfed pasture (SIM2), indicate that both models (SIM1 and SIM2) are comparable for rotational grazing pasture (RMSEsim1 = 0.0581 vs. RMSEsim2 = 0.0485 cm3 cm-3) and the continuous grazing pasture (RMSEsim1 = 0.0485 vs. RMSEsim2 = 0.0602 cm3 cm-3), while for the rainfed maize plots SIM1 shows lower RMSE (average for all plots RMSE = 0.0542 cm3 cm-3) compared to the site-specific calibration model (SIM2 – average for all plots RMSE = 0.0645 cm3 cm-3). Finally, OPTRAM estimations are close to in situ measurement values while Surface Soil Moisture at 1 km (SSM1 km) tends to underestimate the measurements during maize crop growing season. Soil moisture retrieval from high-resolution Sentinel-2 optical images allows water stress conditions to be effectively mapped, supporting decision making in irrigation scheduling and other crop management.     

Comparing methods for mapping canopy chlorophyll content in a mixed mountain forest using Sentinel-2 data

The Sentinel-2 Multi-Spectral Imager (MSI) has three spectral bands centered at 705, 740, and 783 nm wavelengths that exploit the red-edge information useful for quantifying plant biochemical traits. This sensor configuration is expected to improve the prediction accuracy of vegetation chlorophyll content. In this work, we assessed the performance of several statistical and physical-based methods in retrieving canopy chlorophyll content (CCC) from Sentinel-2 in a heterogeneous mixed mountain forest. Amongst the algorithms presented in the literature, 13 different vegetation indices (VIs), a non-parametric statistical approach, and two radiative transfer models (RTM) were used to assess the CCC prediction accuracy. A field campaign was conducted in July 2017 to collect in situ measurements of CCC in Bavarian forest national park, and the cloud-free Sentinel-2 image was acquired on 13 July 2017. The leave-one-out cross-validation technique was used to compare the VIs and the non-parametric approach. Whereas physical-based methods were calibrated using simulated data and validated using the in situ reference dataset. The statistical-based approaches, such as the modified simple ratio (mSR) vegetation index and the partial least square regression (PLSR) outperformed all other techniques. As such the modified simple ratio (mSR3) (665, 865) gave the lowest cross-validated RMSE of 0.21 g/m² (R² = 0.75). The PLSR resulted in the highest R² of 0.78, and slightly higher RMSE =0.22 g/m² than mSR3. The physical-based approach-INFORM inversion using look-up table resulted in an RMSE =0.31 g/m², and R² = 0.67. Although mapping CCC using these methods revealed similar spatial distribution patterns, over and underestimation of low and high CCC values were observed mainly in the statistical approaches. Further validation using in situ data from different terrestrial ecosystems is imperative for both the statistical and physical-based approaches' effectiveness to quantify CCC before selecting the best operational algorithm to map CCC from Sentinel-2 for long-term terrestrial ecosystems monitoring across the globe.     

遥感高时空融合方法的研究进展及应用现状

针对遥感图像的"时空矛盾",评述了当前解决这一问题最主要的方法即遥感时空信息融合的方法,包括基于变化模型的融合、基于重建模型的融合以及基于学习模型的融合。通过分析各个模型的研究现状,指出了每种模型方法的优劣,特别重点介绍了影响较大的自适应时空融合方法的理论以及对其的改进算法。同时本文总结了当前时空融合模型在长时间序列模拟以及大区域数据集生成等方面的实际应用的效果,以及分析了影响时空融合结果的主要因素。最后基于这些问题和影响因素提出了今后时空融合模型发展的目标和方向。