陡峭山区影像的半经验地形校正

本文提出一种新的半经验地形校正模型SCEDIL(Simple topographic Correction using Estimation of Diffuse Light),该模型通过结合DEM与光学影像数据寻找局部区域内完全光照和阴影的水平像元,并以光照、阴影水平像元的平均反射率值估算局部区域散射辐射比,提高了陡峭山区影像的地形校正精度。以高分一号卫星和Landsat ETM+影像为例,从目视判读和定量分析两个方面,比较分析该算法与传统半经验地形校正算法(C、SCS+C)的校正结果。结果表明:(1)对较为平坦的地形,SCEDIL和C、SCS+C校正都有较好的目视结果;对地面起伏较大的陡峭地形,C、SCS+C校正后,原阴影区域易呈现破碎化特征,SCEDIL校正后,原阴影区域过渡较为平滑。(2)SCEDIL校正后,各波段反射率的均值和标准差优于C、SCS+C校正,SCEDIL校正后,影像总分类精度与同类地物光谱信息均一性均优于C和SCS+C校正。SCEDIL半经验地形校正方法能有效地去除影像中的地形干扰,尤其对陡峭地形的校正效果,优于常规地形校正模型。     

Assessment of different topographic correction methods in ALOS AVNIR-2 data over a forest area

Abstract

Because the removal of topographic effects is one the most important pre-processing steps when extracting information from satellite images in digital Earth applications, the problem of differential terrain illumination on satellite imagery has been investigated for at least 20 years. As there is no superior topographic correction method applicable to all areas and all images, a comparison of topographic normalization methods in different regions and images is necessary. In this study, common topographic correction methods were applied on an ALOS AVNIR-2 image of a rugged forest area, and the results were evaluated through different criteria. The results show that the simple correction methods [Cosine, Sun-Canopy-sensor (SCS), and Minnaert correction] are inefficient in exceptionally rough forests. Among the improved correction methods (SCS+C, modified Minnaert, and pixel-based Minnaert), the best result was achieved using a pixel-based Minnaert approach in which a separate correction factor in various slope angles is used. Thus, this method should be considered for topographic correction, especially in forests with severe topography.     

地形抹平与半经验模型的Landsat TM/ETM+地形校正方法

地形校正可以减小地形起伏对地物光谱的影响,提高计算机分类在山区的精度。设计了针对全球土地覆盖分类的Landsat TM/ETM+数据地形校正方法 SCOS(Smoothed COS余弦),首先对地形的坡度角进行抹平处理,很大程度上削弱了地表非朗伯性对地形校正的影响,然后利用简单有效的余弦校正去除地形效应。该方法与其他常用地形校正算法的对比分析是通过对全球不同区域、不同地表覆盖的有代表性的6景Landsat TM/ETM+数据的试验,采用统计分析与目视判读的方式,从过度校正和类内均一性两个方面进行的。结果表明,该方法在目视效果和统计结果上优于常规方法,并且更加简单有效,无需复杂的大气参数及传感器参数,满足全球地表覆盖分类对地形校正的需求。     

A coupled atmospheric and topographic correction algorithm for remotely sensed satellite imagery over mountainous terrain

Radiometric correction is an important issue in the quantitative remote-sensing community. By integrating dark object subtraction (DOS)-based atmospheric correction with physics-based topographic correction, a coupled land surface reflectance retrieval algorithm (coupled atmospheric and topographic correction algorithm, named the CAT algorithm) for rugged mountainous regions is proposed. Terra MODIS-derived atmospheric characterization data (including aerosol optical depth, integrated precipitable water, surface pressure, and ozone concentration) are employed as inputs for the proposed algorithm. A physics-based path radiance estimation model is proposed and embedded in the CAT algorithm, and band-specific per-pixel path radiance values are calculated. After the CAT algorithm was performed, the correlation between reflectance and terrain was dramatically reduced, with correlation coefficients nearly equal zero, especially for the near infrared and short-wave infrared bands, meanwhile the image information content increased over 20%. To provide a comparison with previous studies, two commonly used methods in the literature (DOS + Cosine and DOS + C) were employed. The results of the comparison show that the proposed algorithm performed better in both atmospheric and topographic corrections without empirical regression.     

The effect of atmospheric and topographic correction methods on land cover classification accuracy

Mapping of vegetation in mountain areas based on remote sensing is obstructed by atmospheric and topographic distortions. A variety of atmospheric and topographic correction methods has been proposed to minimize atmospheric and topographic effects and should in principle lead to a better land cover classification. Only a limited number of atmospheric and topographic combinations has been tested and the effect on class accuracy and on different illumination conditions is not yet researched extensively. The purpose of this study was to evaluate the effect of coupled correction methods on land cover classification accuracy. Therefore, all combinations of three atmospheric (no atmospheric correction, dark object subtraction and correction based on transmittance functions) and five topographic corrections (no topographic correction, band ratioing, cosine correction, pixel-based Minnaert and pixel-based C-correction) were applied on two acquisitions (2009 and 2010) of a Landsat image in the Romanian Carpathian mountains. The accuracies of the fifteen resulting land cover maps were evaluated statistically based on two validation sets: a random validation set and a validation subset containing pixels present in the difference area between the uncorrected classification and one of the fourteen corrected classifications. New insights into the differences in classification accuracy were obtained. First, results showed that all corrected images resulted in higher overall classification accuracies than the uncorrected images. The highest accuracy for the full validation set was achieved after combination of an atmospheric correction based on transmittance functions and a pixel-based Minnaert topographic correction. Secondly, class accuracies of especially the coniferous and mixed forest classes were enhanced after correction. There was only a minor improvement for the other land cover classes (broadleaved forest, bare soil, grass and water). This was explained by the position of different land cover types in the landscape. Finally, coupled correction methods showed most efficient on weakly illuminated slopes. After correction, accuracies in the low illumination zone (cos β ≤ 0.65) were improved more than in the moderate and high illumination zones. Considering all results, best overall classification results were achieved after combination of the transmittance function correction with pixel-based Minnaert or pixel-based C-topographic correction. Furthermore, results of this bi-temporal study indicated that the topographic component had a higher influence on classification accuracy than the atmospheric component and that it is worthwhile to invest in both atmospheric and topographic corrections in a multi-temporal study.     

Evaluation of different topographic correction methods for Landsat imagery

The recent free availability of Landsat historical data provides new potentials for land-cover change studies. Multi-temporal studies require a previous radiometric and geometric homogenization of input images, to better identify true changes. Topographic normalization is one of the key steps to create consistent and radiometricly stable multi-temporal time series, since terrain shadows change throughout time. This paper aims to evaluate different methods for topographic correction of Landsat TM-ETM+ data. They were assessed for 15 ETM+ images taken under different illumination conditions, using two criteria: (a) reduction of the standard deviation (SD) for different land-covers and (b) increase in temporal stability of a time series for individual pixels. We observed that results improve when land-cover classes where processed independently when applying the more advanced correction algorithms such as the C-correction and the Minnaert correction. Best results were obtaining for the C-correction and the empiric–statistic correction. Decreases of the SD for bare soil pixels were larger than 100% for the C-correction and the empiric–statistic correction method compared to the other correction methods in the visible spectrum and larger than 50% in the IR region. In almost all tests the empiric–statistic method provided better results than the C-correction. When analyzing the multi-temporal stability, pixels under bad illumination conditions (northern orientation) improved after correction, while a deterioration was observed for pixels under good illumination conditions (southern orientation). Taken this observation into account, a simple but robust method for topographic correction of Landsat imagery is proposed.     

一种顾及空间相关性遥感影像辐射度的地形校正算法

地形校正的目的是消除太阳光照对不规则地面地物辐射值的影响。这种影响会使相似植被类型地物的辐射值发生很大的变化。因此,在地形复杂的地区,地形校正是影像预处理的一个重要步骤。传统的基于单像素的地形校正方法,虽然减小了辐射值的变化,但在太阳入射角低的地区常常出现校正过度的情况。针对这种误差进行分析,提出一种考虑了空间相关性的校正算法,并且利用鄂西地区的Landsat7卫星影像进行的试验证明,该算法优于传统的地形校正模型。     

山地叶面积指数反演理论、方法与研究进展

叶面积指数LAI（Leaf Area Index）是表征叶片疏密程度和冠层结构特征的重要植被参数,在气候变化、作物生长模型以及碳、水循环研究中发挥着重要作用。遥感是获取区域及全球尺度LAI的一个重要手段,当前LAI产品主要基于遥感数据反演得到,但是多数LAI产品算法并未考虑地形特征的影响,导致山地LAI遥感反演精度不确定性大。提高山地LAI遥感反演精度亟需考虑地形因子对冠层反射率的影响,其中山地冠层反射率模型和遥感数据地形校正是提升山地LAI遥感反演精度的关键。本文围绕山地LAI遥感反演理论与方法,综合分析了国内外山地冠层反射率模型和地形校正模型的研究进展,总结了目前山地LAI遥感反演存在的问题,并讨论了未来研究的发展趋势。     

A shadow- eliminated vegetation index (SEVI) for removal of self and cast shadow effects on vegetation in rugged terrains

ABSTRACT

The effect of terrain shadow, including the self and cast shadows, is one of the main obstacles for accurate retrieval of vegetation parameters by remote sensing in rugged terrains. A shadow- eliminated vegetation index (SEVI) was developed, which was computed from only red and near-infrared top-of-atmosphere reflectance without other heterogeneous data and topographic correction. After introduction of the conceptual model and feature analysis of conventional wavebands, the SEVI was constructed by ratio vegetation index (RVI), shadow vegetation index (SVI) and adjustment factor (f (Δ)). Then three methods were used to validate the SEVI accuracy in elimination of terrain shadow effects, including relative error analysis, correlation analysis between the cosine of solar incidence angle (cosi) and vegetation indices, and comparison analysis between SEVI and conventional vegetation indices with topographic correction. The validation results based on 532 samples showed that the SEVI relative errors for self and cast shadows were 4.32% and 1.51% respectively. The coefficient of determination between cosi and SEVI was only 0.032 and the coefficient of variation (std/mean) for SEVI was 12.59%. The results indicate that the proposed SEVI effectively eliminated the effect of terrain shadows and achieved similar or better results than conventional vegetation indices with topographic correction.     

Issues in the application of Digital Surface Model data to correct the terrain illumination effects in Landsat images

The accuracy of topographic correction of Landsat data based on a Digital Surface Model (DSM) depends on the quality, scale and spatial resolution of the DSM data used and the co-registration between the DSM and the satellite image. A physics-based bidirectional reflectance distribution function (BRDF) and atmospheric correction model in conjunction with a 1-second DSM was used to conduct the analysis in this paper. The results show that for the examples used from Australia, the 1-second DSM, can provide an effective product for this task. However, it was found that some remaining artefacts in the DSM data, originally due to radar shadow, can still cause significant local errors in the correction. Where they occur, false shadows and over-corrected surface reflectance factors can be observed. More generally, accurate co-registration between satellite images and DSM data was found to be critical for effective correction. Mis-registration by one or two pixels could lead to large errors of retrieved surface reflectance factors in gully and ridge areas. Using low-resolution DSM data in conjunction with high-resolution satellite images will also fail to correct significant terrain components where they occur at the finer scales of the satellite images. DSM resolution appropriate to the resolution of satellite image and the roughness of the terrain is needed for effective results, and the rougher the terrain, the more critical will be the accurate registration.     

ALI遥感数据的岩矿反射率反演方法研究

岩矿反射率反演是遥感蚀变信息提取中最基础、最关键的一个过程。本介绍了卫星EO-1携带的先进陆地成像仪(ALI)为及其定标方式,分析了两种重要的岩矿反射率反演方法:表观反射率的计算和地表真实反射率的反演。通过分析表明:表观反射率实现关键是参数的获取,地表真实反射率反演关键是模型的选择和参数设置。就反演效果来说,基于大气辐射传输的FLAASH模型优于公式计算的表观反射率,这对应用遥感数据进行岩矿蚀变信息提取的研究人员有一定的参考价值。     

An improved ANUDEM method combining topographic correction and DEM interpolation

Void filling and anomaly replacement in conjunction with auxiliary sources of data have been widely used to improve the quality of existing problematic high-resolution digital elevation models. However, the traditional interpolation methods used for this purpose have always failed to eliminate the discrepancies between different data-sets. In this paper, an improved ANUDEM method is presented for DEM interpolation, which incorporates the idea of topographic correction using high correlation of topological structure between contour lines (CLs) from multi-scale digital elevation models (DEM). Firstly, the terrain topological structure is extracted from the CLs of a low-resolution DEM. The topographic surface correction is then undertaken based on the extracted structure, which recovers the topographic information of the sharp depressions and eminences to fit the high-resolution representation. Finally, the breaklines of the terrain surface are distilled and integrated into the denser DEM generation. The experiments undertaken confirmed the superiority of the proposed method over the other DEM interpolation methods. It is shown that the proposed method can provide results with a higher accuracy, as well as a better visual quality.     

Effect of topographic correction on forest change detection using spectral trend analysis of Landsat pixel-based composites

Pixel-based image compositing enables production of large-area surface reflectance images that are largely devoid of clouds, cloud shadows, or haze. Change detection with spectral trend analysis uses a dense time series of images, such as pixel-based composites, to quantify the year, amount, and magnitude of landscape changes. Topographically-related shadows found in mountainous terrain may confound trend-based forest change detection approaches. In this study, we evaluate the impact of topographic correction on trend-based forest change detection outcomes by comparing the amount and location of changes identified on an image composite with and without a topographic correction. Moreover, we evaluated two different approaches to topographic correction that are relevant to pixel-based image composites: the first corrects each pixel according to the day of year (DOY) the pixel was acquired, whilst the second corrects all pixels to a single reference date (August 1st), which was also the target date for generating the pixel-based image composite. Our results indicate that a greater area of change is detected when no topographic correction is applied to the image composite, however, the difference in change area detected between no correction and either the DOY or the August 1st correction is minor and less than 1% (0.54–0.85%). The spatial correspondence of these different approaches is 96.2% for the DOY correction and 97.7% for the August 1st correction. The largest differences between the correction processes occur in valleys (0.71–1.14%), upper slopes (0.71–1.09%), and ridges (0.73–1.09%). While additional tests under different conditions and in other environments are encouraged, our results indicate that topographic correction may not be justified in change detection routines computing spectral trends from pixel-based composites.     

基于DEM的ETM+图像辐射校正及汉江流域反照率的计算

对Gilabert[1]大气校正方法进行改进,提出了一种简单易行的辐射校正方法。该方法只需从ETM⁺图像上获取两个不同海拔高程的暗像元,再依据DEM进行大气校正,反演出各个高程上的大气及辐射参数,并在此基础上进行地形校正,进而获取LandsatETM⁺1～5和7波段的地表反射率,通过对各窄波段地表反射率的线性组合获取反照率。结果表明,该方法可行,精度较高。     

SCS+C地形辐射校正模型的应用分析研究

在对有森林覆盖的山区影像进行地形辐射校正时，基于太阳-冠层-传感器（SCS）几何关系的校正模型优于基于太阳-地形-传感器（STS）几何关系的模型。SCS校正模型解释了树木不依赖于地形、观测角和光照入射角而具有向地性生长的本质特性，但在某些地形区域，SCS与余弦校正同样存在过度校正的问题。为了解决这个问题，研究者在SCS校正模型中引入C校正系数来解释散射辐射项，提出了SCS+C校正模型。以北京密云Landsat 5影像为数据源，通过目视判别、直方图、定量的统计参数和地物光谱曲线对比等方法，对SCS+C校正模型与传统的余弦校正、C校正和SCS校正模型进行了对比。结果表明，4种方法均能在很大程度上消除地形阴影，更好地反映阴影区域的细节信息； 从总体的光谱特性保真程度来说，余弦和SCS校正都因过度校正问题表现较差，SCS+C校正最好，C校正次之。     

Automatic sub-pixel co-registration of Landsat-8 Operational Land Imager and Sentinel-2A Multi-Spectral Instrument images using phase correlation and machine learning based mapping

ABSTRACT

This study investigates misregistration issues between Landsat-8/ Operational Land Imager and Sentinel-2A/ Multi-Spectral Instrument at 30?m resolution, and between multi-temporal Sentinel-2A images at 10?m resolution using a phase-correlation approach and multiple transformation functions. Co-registration of 45 Landsat-8 to Sentinel-2A pairs and 37 Sentinel-2A to Sentinel-2A pairs were analyzed. Phase correlation proved to be a robust approach that allowed us to identify hundreds and thousands of control points on images acquired more than 100 days apart. Overall, misregistration of up to 1.6 pixels at 30?m resolution between Landsat-8 and Sentinel-2A images, and 1.2 pixels and 2.8 pixels at 10?m resolution between multi-temporal Sentinel-2A images from the same and different orbits, respectively, were observed. The non-linear random forest regression used for constructing the mapping function showed best results in terms of root mean square error (RMSE), yielding an average RMSE error of 0.07?±?0.02 pixels at 30?m resolution, and 0.09?±?0.05 and 0.15?±?0.06 pixels at 10?m resolution for the same and adjacent Sentinel-2A orbits, respectively, for multiple tiles and multiple conditions. A simpler 1st order polynomial function (affine transformation) yielded RMSE of 0.08?±?0.02 pixels at 30?m resolution and 0.12?±?0.06 (same Sentinel-2A orbits) and 0.20?±?0.09 (adjacent orbits) pixels at 10?m resolution.     

地表反射率地形校正物理模型与效果评价方法研究进展

地形校正是准确获取地形复杂区遥感反射率的重要步骤,对提高山区地表遥感参数定量化反演精度,扩大遥感产品应用广度具有重要意义。从20世纪80年代开始,国内外学者开始对准确获取山区地表遥感反射率进行研究,建立了多种地形校正模型来减少或消除遥感图像中地形效应影响,减少同种地表类型的反射率差异,并将地形校正模型分为经验模型和物理模型。根据构建物理模型时是否考虑地表非朗伯体特性,将物理模型分为朗伯体假设模型和非朗伯体假设模型,本文分别从朗伯体假设模型和非朗伯体假设校正模型展开叙述。从两类模型构建的理论基础,模型特点,局限性等几方面进行分析和讨论,描述了两类模型的发展历程,系统阐述了朗伯体假设模型和非朗伯体假设模型的适用性和不足,剖析了目前地形校正模型存在的问题与挑战。同时,本文也比较了应用于地形校正的效果评价方法,并展望了地形校正方法和地形校正评价方法的未来主要发展方向。     

Terrain characterization using SRTM data

Earth’s surface possesses relief because the geomorphic processes operate at different rates, and geologic structure plays a dominant role in the evolution of landforms (Thornbury, 1954). The spatial pattern of relief yields the topographic mosaic of a terrain and is normally extracted from the topographical maps which are available at various scales. As cartographic abstractions are scale dependent, topographical maps are rarely good inputs for terrain analysis. Currently, the shuttle radar topography mission (SRTM) provides one of the most complete, highest resolution digital elevation model (DEM) of the Earth. It is an ideal data-set for precise terrain analysis and topographic characterization in terms of the nature of altimetric distribution, relief aspects, patterns of lineaments and surface slope, topographic profiles and their visualisation, correlation between geology and topography, hypsometric attributes and finally, the hierarchy of terrain sub-units. The present paper extracts the above geomorphic features and terrain character of part of the Chotonagpur plateau and the Dulung River basin therein using SRTM data.     

基于LiDAR的农田地形环境三维重建方法

为克服传统农田土地平整测量方法耗时费力的特点,提出采用LiDAR技术对农田地形进行重建的探索性研究。通过HDL-32E型激光雷达等搭建了系统的硬件平台,应用C++语言编写了系统数据的采集程序;在此基础上对激光雷达所采集数据进行了标定,研究了农田地形重建系统中不同坐标系的转换方法;同时基于最小值去噪法设计了更适用于农田地形点云去噪的均值限差去噪法。通过对比在农田起伏较大区域不同坡度范围内RTK与激光雷达所测单元个数,对系统精度进行了评价;最后实现了车载农田地形重建系统的界面显示、应用与精度评估。结果表明,在10°~15°、25°~30°大坡度范围内激光雷达所获农田地形更为丰富,精度更高。该方法重建的农田地形模型点云数据和原始农田地形点云数据投影面积逼近度可达93%,验证了本文研究方法应用于农田地形环境重建的可行性,同时为今后的土地精细平整工作提供了理论参考与依据。     

包络线去除的丘陵地区遥感影像阴影信息重建

中国西南丘陵常态山和喀斯特山交错分布,遥感影像普遍存在山体阴影,分布零散且无规律,基于DEM的地形校正模型(C校正等)虽然算法成熟、易于操作,但在复杂地形区存在误差。引入基于相似像元包络线的阴影校正方法(CR校正),按照阴影提取、包络线去除、相似像元寻找和阴影亮度重建的步骤,采用西南丘陵地区Landsat 8 OLI影像进行验证实验。结果表明:CR校正后,阴影区的视觉特征与邻近非阴影区趋于一致,阴影像元亮度有明显提升;校正后影像主要波段标准差减小,与非阴影区参考光谱的相对均方根误差在2.919%以内,最低仅为0.516%;自动分类精度从43.59%提高到61.57%,CR校正有效提高了有阴影的丘陵地区遥感影像质量。