基于超分辨率重建的多时相MODIS与Landsat反射率融合方法

提出一种基于超分辨率重建的MODIS与Landsat反射率图像融合方法,以STARFM算法与超分辨率重建为基础,使用观测的MODIS和Landsat地表反射率图像预测给定时刻的Landsat合成反射率图像。该方法利用基于稀疏表示的超分辨率重建方法对MODIS图像进行分辨率增强,实验结果表明这一操作能够增加原MODIS图像的空间细节,有助于提高STARFM算法的预测精度;另一方面,考虑输入两个基时刻图像相差较大时原STARFM算法预测的反射率会存在"时间平滑"的问题,限制每次只使用一个基时刻MODIS和Landsat图像对进行STARFM预测,使用逐图像块选择策略,从由两个基时刻图像分别进行预测得到的两组预测图像中选择最优的预测,同样得到了优于STARFM算法的预测结果。     

在变化检测中Landsat-7异常数据处理方法试验

为了更好地利用Landsat系列卫星影像进行变化检测,本文对以SLC-off方式获取的Landsat-7影像进行修复处理,使其能够参与变化检测。本文试验对比3种常用异常数据修复模型,认为多影像局部自适应回归分析模型修复方法最适合后续进行的变化检测,并总结了Landsat-7异常数据参与变化检测的流程方法。     

A Landsat surface reflectance dataset for North America, 1990-2000

The Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center has processed and released 2100 Landsat Thematic Mapper and Enhanced Thematic Mapper Plus surface reflectance scenes, providing 30-m resolution wall-to-wall reflectance coverage for North America for epochs centered on 1990 and 2000. This dataset can support decadal assessments of environmental and land-cover change, production of reflectance-based biophysical products, and applications that merge reflectance data from multiple sensors [e.g., the Advanced Spaceborne Thermal Emission and Reflection Radiometer, Multiangle Imaging Spectroradiometer, Moderate Resolution Imaging Spectroradiometer (MODIS)]. The raw imagery was obtained from the orthorectified Landsat GeoCover dataset, purchased by NASA from the Earth Satellite Corporation. Through the LEDAPS project, these data were calibrated, converted to top-of-atmosphere reflectance, and then atmospherically corrected using the MODIS/6S methodology. Initial comparisons with ground-based optical thickness measurements and simultaneously acquired MODIS imagery indicate comparable uncertainty in Landsat surface reflectance compared to the standard MODIS reflectance product (the greater of 0.5% absolute reflectance or 5% of the recorded reflectance value). The rapid automated nature of the processing stream also paves the way for routine high-level products from future Landsat sensors.     

Evaluation of a rule-based compositing technique for Landsat-5 TM and Landsat-7 ETM+ images

Image compositing is a multi-objective optimization process. Its goal is to produce a seamless cloud and artefact-free artificial image. This is achieved by aggregating image observations and by replacing poor and cloudy data with good observations from imagery acquired within the timeframe of interest. This compositing process aims to minimise the visual artefacts which could result from different radiometric properties, caused by atmospheric conditions, phenologic patterns and land cover changes. It has the following requirements: (1) image compositing must be cloud free, which requires the detection of clouds and shadows, and (2) the image composite must be seamless, minimizing artefacts and visible across inter image seams. This study proposes a new rule-based compositing technique (RBC) that combines the strengths of several existing methods. A quantitative and qualitative evaluation is made of the RBC technique by comparing it to the maximum NDVI (MaxNDVI), minimum red (MinRed) and maximum ratio (MaxRatio) compositing techniques. A total of 174 Landsat TM and ETM+ images, covering three study sites and three different timeframes for each site, are used in the evaluation. A new set of quantitative/qualitative evaluation techniques for compositing quality measurement was developed and showed that the RBC technique outperformed all other techniques, with MaxRatio, MaxNDVI, and MinRed techniques in order of performance from best to worst.     

Landsat TM遥感影像中厚云和阴影去除

提出了一种新的利用多时相Landsat TM影像数据进行的厚云及其阴影去除的方法。该方法通过分析厚云及其阴影的光谱特征, 设计了厚云和云阴影识别模型。该算法的实现是采用图像配准技术、非监督分类、像元替换等运算, 计算出厚云和云阴影区域的TM影像替换数据, 进而得到消除或者减少云影响的TM遥感影像。试验结果表明本文提出的厚云及其阴影去除方法效果很好, 能消除或者弱化云对TM影像数据的影响。     

A Cloud Removal Algorithm to Generate Cloud and Cloud Shadow Free Images Using Information Cloning

One of the main problems of optical remote sensing is clouds and cloud shadows caused by specific atmospheric conditions during data acquisition. These features limit the usage of acquired images and increase the difficulty in data analysis, such as normalized difference vegetation index values, misclassification, and atmospheric correction. Accurate detection and reliable cloning of cloud and cloud shadow features in satellite images are very useful processes for optical remote sensing applications. In this study, an automated cloud removal algorithm to generate cloud and cloud shadow free images from multitemporal Landsat-8 images is introduced. Cloud and cloud shadow areas are classified by using process-based rule set developed by using spectral and spatial features after applying simple linear iterative clustering superpixel segmentation algorithm to the image to find cloud pixel groups easily and correctly. Segmentation-based cloud detection method gives better results than pixel-based for detection of cloud and cloud shadow patches. After detection of clouds and cloud shadows, cloud-free images are created by cloning cloudless regions from multitemporal dataset. Spectral and structural consistency are preserved by considering spectral features and seasonal effects while cloning process. Statistical similarity tests are applied to find best cloud-free image to use for cloning process. Cloning results are tested with the structural similarity index metric to evaluate the performance of cloning algorithm.     

Cloud removal for remotely sensed images by similar pixel replacement guided with a spatio-temporal MRF model

Cloud cover is generally present in remotely sensed images, which limits the potential of the images for ground information extraction. Therefore, removing the clouds and recovering the ground information for the cloud-contaminated images is often necessary in many applications. In this paper, an effective method based on similar pixel replacement is developed to solve this task. A missing pixel is filled using an appropriate similar pixel within the remaining region of the target image. A multitemporal image is used as the guidance to locate the similar pixels. A pixel-offset based spatio-temporal Markov random fields (MRF) global function is built to find the most suitable similar pixel. The proposed method was tested on MODIS and Landsat images and their land surface temperature products, and the experiments verify that the proposed method can achieve highly accurate results and is effective at dealing with the obvious atmospheric and seasonal differences between multitemporal images.     

融合多光谱影像的高光谱影像厚云去除方法

云遮挡对高光谱影像的应用造成了不可忽视的影响。现有云去除方法通常利用时域近邻的同源影像提供辅助信息。然而,高光谱影像(如GF-5和EO-1高光谱影像)较低的时间分辨率导致同源辅助影像中可能存在较大的地物覆盖变化。时间分辨率更高的多光谱影像(如Landsat 8 OLI影像)能提供时间上更接近于高光谱云影像的辅助信息,从而减少地物覆被变化带来的影响。为应对高光谱和多光谱波段之间差异较大的问题,本文基于空谱随机森林(spatial-spectral-based random forest,SSRF)方法,提出一种利用多光谱影像(Landsat 8 OLI影像)对高光谱影像进行厚云去除的方法,将其简记为SSRF_M。SSRF_M较强的非线性拟合能力使其能够综合利用多光谱影像所有波段的有效数据对各个高光谱波段进行重建。本文使用GF-5和EO-1高光谱影像进行模拟云去除试验,视觉和定量评价结果均表明,与利用时间间隔更长的同源辅助影像的方法相比,本文方法能获得更高精度的云下信息重建结果。     

Selection and quality assessment of Landsat data for the North American forest dynamics forest history maps of the US

Using the NASA Earth Exchange platform, the North American Forest Dynamics (NAFD) project mapped forest history wall-to-wall, annually for the contiguous US (1986–2010) using the Vegetation Change Tracker algorithm. As with any effort to identify real changes in remotely sensed time-series, data gaps, shifts in seasonality, misregistration, inconsistent radiometry and cloud contamination can be sources of error. We discuss the NAFD image selection and processing stream (NISPS) that was designed to minimize these sources of error. The NISPS image quality assessments highlighted issues with the Landsat archive and metadata including inadequate georegistration, unreliability of the pre-2009 L5 cloud cover assessments algorithm, missing growing-season imagery and paucity of clear views. Assessment maps of Landsat 5–7 image quantities and qualities are presented that offer novel perspectives on the growing-season archive considered for this study. Over 150,000+ Landsat images were considered for the NAFD project. Optimally, one high quality cloud-free image in each year or a total of 12,152 images would be used. However, to accommodate data gaps and cloud/shadow contamination 23,338 images were needed. In 220 specific path-row image years no acceptable images were found resulting in data gaps in the annual national map products.     

Evaluation of directional normalization methods for Landsat TM/ETM+ over primary Amazonian lowland forests

Biodiversity mapping in extensive tropical forest areas poses a major challenge for the interpretation of Landsat images, because floristically clearly distinct forest types may show little difference in reflectance. In such cases, the effects of the bidirectional reflection distribution function (BRDF) can be sufficiently strong to cause erroneous image interpretation and classification. Since the opening of the Landsat archive in 2008, several BRDF normalization methods for Landsat have been developed. The simplest of these consist of an empirical view angle normalization, whereas more complex approaches apply the semi-empirical Ross–Li BRDF model and the MODIS MCD43-series of products to normalize directional Landsat reflectance to standard view and solar angles. Here we quantify the effect of surface anisotropy on Landsat TM/ETM+ images over old-growth Amazonian forests, and evaluate five angular normalization approaches. Even for the narrow swath of the Landsat sensors, we observed directional effects in all spectral bands. Those normalization methods that are based on removing the surface reflectance gradient as observed in each image were adequate to normalize TM/ETM+ imagery to nadir viewing, but were less suitable for multitemporal analysis when the solar vector varied strongly among images. Approaches based on the MODIS BRDF model parameters successfully reduced directional effects in the visible bands, but removed only half of the systematic errors in the infrared bands. The best results were obtained when the semi-empirical BRDF model was calibrated using pairs of Landsat observation. This method produces a single set of BRDF parameters, which can then be used to operationally normalize Landsat TM/ETM+ imagery over Amazonian forests to nadir viewing and a standard solar configuration.     

Global operational land imager Landsat-8 reflectance-based active fire detection algorithm

A global operational land imager (GOLI) Landsat-8 daytime active fire detection algorithm is presented. It utilizes established contextual active fire detection approaches but takes advantage of the significant increase in fire reflectance in Landsat-8 band 7 (2.20?μm) relative to band 4 (0.66?μm). The detection thresholds are fixed and based on a statistical examination of 39 million non-burning Landsat-8 pixels. Multi-temporal tests based on band 7 reflectance and relative changes in normalized difference vegetation index in the previous six months are used to reduce commissions errors. The probabilities of active fire detection for the GOLI and two recent Landsat-8 active fire detection algorithms are simulated to provide insights into their performance with respect to the fire size and temperature. The algorithms are applied to 11 Landsat-8 images that encompass a range of burning conditions and environments. Commission and omission errors are assessed by visual interpretation of detected active fire locations and by examination of the Landsat-8 images and higher spatial resolution Google Earth imagery. The GOLI algorithm has lower omission and comparable commission errors than the recent Landsat-8 active fire detection algorithms. The GOLI algorithm has demonstrable potential for global application and is suitable for implementation with other Landsat-like reflective wavelength sensors.     

基于支持向量机的遥感影像厚云及云阴影去除

本文提出了一种基于支持向量机的遥感影像厚云及云阴影去除方法。首先利用支持向量机的学习性能检测影像中的云层,并利用太阳角度信息,判定云阴影区域,得到云层和云阴影的二值图。再对影像进行支持向量值轮廓波变换,利用云层和云阴影二值图生成的选择矩阵,对变换系数进行多层镶嵌,完成云层及云阴影的初去除。对影像镶嵌未能去除的云层及云阴影,通过统计学补偿的方法进行修复。最后重构图像并进行中值滤波实现厚云及云阴影去除。仿真实验表明,该方法能更好地再现云层覆盖区域的地物信息,去云后的图像具有更好的光滑度和清晰度。     

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波段的相关性明显提高,辐射一致性增强。该转换模型为多源中高分辨率遥感图像高精度辐射归一化提供了新思路。     

Downscaling of Land Surface Temperature Maps in the Texas High Plains with the TsHARP Method

High-resolution evapotranspiration (ET) maps can assist demand-based irrigation management. Development of high-resolution daily ET maps requires high-resolution land surface temperature (LST) images. Earth-observing satellite sensors such as the Landsat 5 Thematic Mapper (TM) and MODerate resolution Imaging Spectroradiometer (MODIS) provide thermal images that are coarser than simultaneously acquired visible and near-infrared images. In this study, we evaluated the TsHARP downscaling technique for its capability to downscale coarser LST images using finer resolution normalized difference vegetation index (NDVI) data. The TsHARP technique was implemented to downscale seven coarser scale (240, 360, 480, 600, 720, 840, and 960 m) synthetic images to a 120 m LST image. The TsHARP was also evaluated for downscaling a coarser 960 m LST image to 240 m to mimic MODIS datasets. Comparison between observed 120 m LST images and 120 m LST images downscaled from coarser 240, 360, 480, 600, 720, 840, and 960 m images yielded root mean square errors of 1.0, 1.3, 1.5, 1.6, 1.7, 1.8, and 1.9°C, respectively. This indicates that the TsHARP method can be used for downscaling coarser (960 m) MODIS-based LST images using finer Landsat (120 m) or MODIS (240 m)-derived NDVI images. However, the TsSHARP method should be evaluated further with real datasets before using it for an operational ET remote sensing program for irrigation scheduling purposes.     

An experimental comparison between KELM and CART for crop classification using Landsat-8 OLI data

The operational land imager (OLI) is the latest instrument in the Landsat series of satellite imagery, which officially began normal operations on 30 May 2013. The OLI includes two bands that are not on the thematic mapper series of sensors aboard Landsat-5 and 7; a cirrus band and a coastal/aerosol band. This paper compares the classification and regression tree and the kernel-based extreme learning machine (KELM) for mapping crops in Hokkaido, Japan, using OLI data, except the cirrus band and the pan band. The OLI data acquired on 8 July 2013 was used for crop classification of beans, beets, grassland, maize, potatoes and winter wheat. The KELM algorithm performed better in this study and achieved overall accuracies of 90.1%. According to the Jeffries–Matusita (J–M) distances, the short wavelength infrared band provides the greater contribution (the highest value was observed for band 6 in OLI data).     

噪声调节主成分变换下的遥感图像薄云去除

针对在光学遥感图像中的地球表面反射率由于云覆盖所造成的衰减问题,本文提出了一种基于噪声调整主成分变换(NAPCT)模型的薄云去除方法.首先,通过噪声估计构建NAPCT除云模型,实现云图像的转换;然后,利用NAPCT第一成分(NAPC1)提供的云分布信息,通过合并云掩膜,对NAPC1上的浑浊像素进行校正;最后,将逆变换应...     

Comparison of eMODIS and MOD/MYD13A2 NDVI products during 2012–2014 spring green-up periods in Alaska and northwest Canada

Accurate monitoring of vegetation dynamics is required to understand the inter-annual variability and long term trends in terrestrial carbon exchange in tundra and boreal ecoregions. In western North America, two Normalized Vegetation Index (NDVI) products based on spectral reflectance data from the Moderate Resolution Imaging Spectroradiometer (MODIS) are available. The MOD/MYD13A2 NDVI product is available as a 16-day composite product in a sinusoidal projection as global hdf tiles. The eMODIS Alaska NDVI product is available as a 7-day composite geotif product in a regional equal area conic projection covering Alaska and the entire Yukon River Basin. These two NDVI products were compared for the 2012–2014 late May–late June spring green-up periods in Alaska and the Yukon Territory. Relative to the MOD/MYD13A2 NDVI product, it is likely that the eMODIS NDVI product contained more cloud-contaminated NDVI values. For example, the MOD/MYD13A2 product flagged substantially fewer pixels as “good quality” in each 16-day composite period compared to the corresponding MODIS Alaska NDVI product from a 7-day composite period. During the spring green-up period, when field-based NDVI increases, the eMODIS NDVI product averaged 43 percent of pixels that declined by at least 0.05 NDVI between 2 composite periods, consistent with cloud-contamination problems, while the MOD/MYD13A2 NDVI averaged only 6 percent of pixels. Based on a cloudy Landsat-8 scene, the eMODIS compositing process selected 23 percent pixels, while the MOD/MYD13A2 compositing process selected less than 0.003 percent pixels. Based on the results, it appears that the MOD/MYD13A2 NDVI product is superior for scientific applications based on NDVI phenology in the tundra and boreal regions of northwestern North America.     

Vicarious Cross-Calibration of the China Environment Satellite Using Nearly Simultaneously Observations of Landsat-7 ETM+ Sensor

To understand the absolute radiometric calibration accuracy of the HJ-A CCD-1 sensors, image from these sensors were compared to nearly simultaneously image from Landsat-7 ETM+ sensors. Although the HJ-A CCD-1 sensor has almost the same wavelength of each central band and band width as Landsat-7 ETM+ sensor, there is slightly difference in spectral response function (SRF). The impacts of SRF difference effects would produce ~2 % uncertainty in predicting reflectance of HJ-A CCD-1 sensor using Landsat-7 ETM+ sensor. The reflectance observed by satellite at top-of-atmosphere generally depends on its’ geometric conditions. The results reveal that the impacts of geometrical conditions would impact on the vicarious cross-calibration accuracy, which should be removed. The performances of cross-calibration are calibrated and validated by four image pairs collected from Yellow River Delta, China, and Qingdao City, China, at four independent times. The results indicate that the HJ-A CCD-1 sensors can be cross calibrated to the Landsat-7 ETM+ sensors to within an accuracy of 3.99 % (denoted by Relative Root Mean Square Error) of each other in all bands except band 4, which has a 6.33 % difference.     

基于光谱分析的MODIS云检测算法研究

总结了云、雪、植被、沙漠、水体等目标的光谱特性以及云和不同目标间的光谱差异，结合MODIS影像的波段范围，提出了一种对不同下垫面通用的多光谱云检测算法，试验证明了算法的有效性。     

Comparison of the Landsat Surface Reflectance Climate Data Record (CDR) and manually atmospherically corrected data in a semi-arid European study area

This study contributes to the quality assessment of atmospherically corrected Landsat surface reflectance data that are routinely generated by the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS). This dataset, named Landsat Surface Reflectance Climate Data Record (Landsat CDR), is available at global scale and offers unprecedented opportunities to land monitoring and management services that require atmospherically corrected Earth observation (EO) data. Our assessment is based on the comparison of the Landsat CDR data against a set of Landsat and DEIMOS-1 images processed to a high degree of accuracy using an industry-standard atmospheric correction algorithm (ATCOR-2). The software package has been used for many years and its correction procedures can be considered consolidated and well-established. The dataset of Landsat and DEIMOS-1 images was acquired over a semi-arid agricultural area located in Lower Austria and was independently corrected by using a manual fine-tuning of ATCOR-2 parameters to reach the highest possible accuracy. Results show a very good correspondence of the surface reflectance in each of the six reflective spectral channels as well as for the NDVI (Normalized Difference Vegetation Index). An additional comparison against a NDVI time series from MODIS revealed also a good correspondence. Coefficients of determination (R²) between the two multi-year and multi-seasonal Landsat/DEIMOS datasets range between 0.91 (blue band) and 0.98 (nIR, SWIR-1 and SWIR-2). The results obtained for our semi-arid test site in Austria confirm previous findings and suggest that automatic atmospheric procedures, such as the one implemented by LEDAPS are accurate enough to be used in land monitoring services that require consistent multi-temporal surface reflectance data.