基于GF-2高分辨率遥感影像的水体提取方法研究

针对GF-2卫星影像数据的特点,利用单波段阈值法、多波段算法、归一化水体指数(NDWI)阈值法、单波段阈值法与阴影水体指数(SWI)相结合的决策树法对刘家峡地区的水体信息进行了提取,同时提出一种单波段阈值法与增强阴影水体指数(ESWI)相结合的决策树分类法,并对以上几种水体提取的效果进行比较分析,发现单波段阈值法与ESWI相结合的决策树分类法能够有效地消除绝大部分阴影的影响,且提取效果较SWI的效果要好,NDWI与多波段算法提取效果次之,单波段阈值法提取效果较差.      

不同指数模型提取ALOS AVNIR-2影像中水体的敏感性和精度分析

以广东惠州附近水域为研究区,利用ALOS AVNIR-2多光谱数据分析水体和其他主要地物在影像上的光谱特征。构建几种常用的波谱指数,分析阈值对各种指数模型提取水体的敏感性及所能达到的最大精度。发现利用这些指数均可以较好地提取图像中的水体信息,提取效果从好到差依次是:归一化差值水体指数法(NDWI)、近红外波段分割法(NIR)、归一化差值植被指数法(NDVI)和比值植被指数法(RVI)。其中NDWI指数模型,在选择合适阈值的情况下,水体提取的总体精度最高可达98%左右,并且提取过程对阈值影响不敏感。     

基于标记分水岭分割的高分辨率影像水体信息提取

提出一种基于标记分水岭分割的高分辨率卫星融合影像提取水体信息的方法。首先采用pansharpening 融合法获得光谱扭曲小的高分辨率卫星融合影像;其次利用标记分水岭算法对高分辨率卫星融合影像进行分割;最后在分割基础上利用水体指数模型提取水体信息。采用QuickBird高分辨率遥感影像数据试验,与eCognition多尺度分割提取水体方法的结果相比,表明本文方法更快速有效,具有实用推广价值。     

基于GF-1遥感数据的城市水体信息提取方法研究

水体信息的高精度提取是水资源监测、调查与管理等研究领域的关键问题。本文选取国产GF-1卫星影像作为数据源,根据影像中典型地物光谱特征采用归一化水体指数法(NDWI)、多波段谱间关系法对研究区进行了水体提取,并通过结果分析提出了基于主成分分析与多尺度分割技术的综合水体信息提取方法。对以上4种方法进行对比分析,验证结果表明:归一化水体指数法易受阴影信息影响;改进的多波段谱间关系法能较为完整地提取水体信息,但受小范围阴影信息影响;主成分分析综合法总体精度较高;多尺度分割提取法能有效地分离水体信息与非水体信息,水体信息提取的效果最佳。     

综合水体指数的创建

在对NDWI和MNDWI等传统水体指数进行分析的基础上,本文提出了一种创建水体指数的新思路:将多光谱图像进行LBV和HSV等复合变换,综合利用变换前的多波段信息和变换后的新波段信息绘制地物光谱曲线,寻找水体的光谱特性,进行水体指数的创建。本文利用ETM+图像HSV变换后的Sat波段以及LBV变换后的V波段,创建了综合水体指数(SWI),并利用SWI进行了遥感影像水体信息自动提取实验。实验表明:利用SWI能准确提取水体信息,有效减少水体提取结果中的误提取像元。     

综合水体指数的创建

在对NDWI和MNDWI等传统水体指数进行分析的基础上,本文提出了一种创建水体指数的新思路:将多光谱图像进行LBV和HSV等复合变换,综合利用变换前的多波段信息和变换后的新波段信息绘制地物光谱曲线,寻找水体的光谱特性,进行水体指数的创建.本文利用ETM+图像HSV变换后的Sat波段以及LBV变换后的V波段,创建了综合水体指数(SWI),并利用SWI进行了遥感影像水体信息自动提取实验.实验表明:利用SWI能准确提取水体信息,有效减少水体提取结果中的误提取像元.     

ASTER数据的煤矿区水体信息提取

为实现多光谱数据对煤矿区水体信息提取,以萍乡市芦新岭及周边煤矿为试验区,在国内外学者对水体指数法提取水体信息和应用的启发下,对试验区不同地物的波谱特征进行分析。本文首先利用数学统计方法建立一种ASTER数据综合归一化差异水体指数模型SWI ASTER_xy,然后利用这一模型对试验区进行水体信息提取。试验表明,该模型能够有效地提取矿区各类水体,还可以运用在ASTER数据（或其他多光谱数据）对不同区域不同地物的提取,在遥感地物定量提取上具有很好的效果和潜力,具有一定的推广价值。     

一种基于城市水体指数与分形几何算法的OLI遥感影像水体提取方法

针对严重污染的城市水体与道路、建筑物、阴影等易于混分,以及遥感水体提取结果不连续、存在斑点问题,本文以广州市流溪河与东江水系为研究对象,基于2016年与2017年OLI遥感影像,采用本文新提出的城市水体指数法（CWI）,同时结合分形几何算法,通过设置形状面积等特征,实现城市复杂环境下的水体信息的自动提取。并与单通道算法、改进的归一化差异水体指数（MNDWI）算法、支持向量机法（SVM）与光谱角度法的水体提取结果进行对比分析。结果表明：SVM算法出现大量斑点,其次为MNDWI水体指数算法,光谱角度算法与单通道算法斑点较少,但水体提取结果不连续,部分河道漏分。本文提出的算法能够克服山体阴影、道路、建筑物等影响,实现城市污染水体以及一般水体连续、准确提取。本文的提取结果可为水资源调查、洪水灾害预测评估、水利规划、环境监测等工作提供基础数据支撑。     

基于ASTER数据的归一化差异水体指数的建立及其应用

在分析归一化差异水体指数的基础上,根据ASTER数据不同地物的光谱特征,提出了基于ASTER数据的归一化差异化水体指数NDWIASTER1、4。并且利用这一指数进行水体信息的提取,发现这是一种简单易行结果客观有效的提取水体的方法。     

一种基于城市水体指数与分形几何算法的OLI遥感影像水体提取方法

针对严重污染的城市水体与道路、建筑物、阴影等易于混分,以及遥感水体提取结果不连续、存在斑点问题,本文以广州市流溪河与东江水系为研究对象,基于2016年与2017年OLI遥感影像,采用本文新提出的城市水体指数法(CWI),同时结合分形几何算法,通过设置形状面积等特征,实现城市复杂环境下的水体信息的自动提取。并与单通道算法、改进的归一化差异水体指数(MNDWI)算法、支持向量机法(SVM)与光谱角度法的水体提取结果进行对比分析。结果表明:SVM算法出现大量斑点,其次为MNDWI水体指数算法,光谱角度算法与单通道算法斑点较少,但水体提取结果不连续,部分河道漏分。本文提出的算法能够克服山体阴影、道路、建筑物等影响,实现城市污染水体以及一般水体连续、准确提取。本文的提取结果可为水资源调查、洪水灾害预测评估、水利规划、环境监测等工作提供基础数据支撑。     

Estimate of heavy metals in soil and streams using combined geochemistry and field spectroscopy in Wan-sheng mining area,Chongqing, China

Heavy metals contaminated soils and water will become a major environmental issue in the mining areas. This paper intends to use field hyper-spectra to estimate the heavy metals in the soil and water in Wan-sheng mining area in Chongqing. With analyzing the spectra of soil and water, the spectral features deriving from the spectral of the soils and water can be found to build the models between these features and the contents of Al, Cu and Cr in the soil and water by using the Stepwise Multiple Linear Regression (SMLR). The spectral features of Al are: 480 nm, 500 nm, 565 nm, 610 nm, 680 nm, 750 nm, 1000 nm, 1430 nm, 1755 nm, 1887 nm, 1920 nm, 1950 nm, 2210 nm, 2260 nm; The spectral features of Cu are: 480 nm, 500 nm, 610 nm, 750 nm, 860 nm, 1300 nm, 1430 nm, 1920 nm, 2150 nm, 2260 nm; And the spectral features of Cr are: 480 nm, 500 nm, 610 nm, 715 nm, 750 nm, 860 nm, 1300 nm, 1430 nm, 1755 nm, 1920 nm, 1950 nm. With these features, the best models to estimate the heavy metals in the study area were built according to the maximal R². The R² of the models of estimating Al, Cu and Cr in the soil and water are 0.813, 0.638, 0.604 and 0.742, 0.584, 0.513 respectively. And the gradient maps of these three types of heavy metals’ concentrations can be created by using the Inverse distance weighted (IDW).The gradient maps indicate that the heavy metals in the soil have similar patterns, but in the North-west of the streams in the study area, the contents are of great differences. These results show that it is feasible to predict contaminated heavy metals in the soils and streams due to mining activities by using the rapid and cost-effective field spectroscopy.     

Determination of the factors governing soil erodibility using hyperspectral visible and near-infrared reflectance spectroscopy

Soil erodibility, which is difficult to estimate and upscaling, was determined in this study using multiple spectral models of soil properties (soil organic matter (SOM), water-stable aggregates (WSA) > 0.25 mm, the geometric mean radius (D_g)). Herein, the soil erodibility indicators were calculated, and soil properties were quantitatively analyzed based on laboratory simulation experiments involving two selected contrasting soils. In addition, continuous wavelet transformation was applied to the reflectance spectra (350–2500 nm) of 65 soil samples from the study area. To build the relationship, the soil properties that control erodibility were identified prior to the spectral analysis. In this study, the SOM, D_g and WSA >0.25 mm were selected to represent the most significant soil properties controlling erodibility and describe the erodibility indicator based on a logarithmic regression model as a function of SOM or WSA > 0.25 mm. Five, six and three wavelet features were observed to calibrate the estimated soil properties model, and the best performance was obtained with a combination feature regression model for SOM (R² = 0.86, p < 0.01), D_g (R² = 0.79, p < 0.01) and WSA >0.25 mm (R² = 0.61, p < 0.01), respectively. One part of the wavelet features captured amplitude variations in the broad shape of the reflectance spectra, and another part captured variations in the shape and depth of the soil dry substances. The wavelet features for the validated dataset used to predict the SOM, WSA >0.25 mm and D_g were not significantly different compared with the calibrated dataset. The synthesized spectral models of soil properties, and the formation of a new equation for soil erodibility transformed from the spectral models of soil properties are presented in this study. These results show that a spectral analytical approach can be applied to complex datasets and provide new insights into emerging dynamic variation with erodibility estimation.     

Determination of the crop row orientations from Formosat-2 multi-temporal and panchromatic images

This paper presents a technique developed for the retrieval of the orientation of crop rows, over anthropic lands dedicated to agriculture in order to further improve estimate of crop production and soil erosion management. Five crop types are considered: wheat, barley, rapeseed, sunflower, corn and hemp. The study is part of the multi-sensor crop-monitoring experiment, conducted in 2010 throughout the agricultural season (MCM’10) over an area located in southwestern France, near Toulouse. The proposed methodology is based on the use of satellite images acquired by Formosat-2, at high spatial resolution in panchromatic and multispectral modes (with spatial resolution of 2 and 8 m, respectively). Orientations are derived and evaluated for each image and for each plot, using directional spatial filters (45° and 135°) and mathematical morphology algorithms. “Single-date” and “multi-temporal” approaches are considered. The single-date analyses confirm the good performances of the proposed method, but emphasize the limitation of the approach for estimating the crop row orientation over the whole landscape with only one date. The multi-date analyses allow (1) determining the most suitable agricultural period for the detection of the row orientations, and (2) extending the estimation to the entire footprint of the study area. For the winter crops (wheat, barley and rapeseed), best results are obtained with images acquired just after harvest, when surfaces are covered by stubbles or during the period of deep tillage (0.27 > R² > 0.99 and 7.15° > RMSE > 43.02°). For the summer crops (sunflower, corn and hemp), results are strongly crop and date dependents (0 > R² > 0.96, 10.22° > RMSE > 80°), with a well-marked impact of flowering, irrigation equipment and/or maximum crop development. Last, the extent of the method to the whole studied zone allows mapping 90% of the crop row orientations (more than 45,000 ha) with an error inferior to 40°, associated to a confidence index ranging from 1 to 5 for each agricultural plot.     

Exploring hyperspectral bands and estimation indices for leaf nitrogen accumulation in wheat

Hyperspectral sensing can provide an effective means for fast and non-destructive estimation of leaf nitrogen (N) status in crop plants. The objectives of this study were to design a new method to extract hyperspectral spectrum information, to explore sensitive spectral bands, suitable bandwidth and best vegetation indices based on precise analysis of ground-based hyperspectral information, and to develop regression models for estimating leaf N accumulation per unit soil area (LNA, g N m⁻²) in winter wheat (Triticum aestivum L.). Three field experiments were conducted with different N rates and cultivar types in three consecutive growing seasons, and time-course measurements were taken on canopy hyperspectral reflectance and LNA under the various treatments. Then, normalized difference spectral indices (NDSI) and ratio spectral indices (RSI) based on the original spectrum and the first derivative spectrum were constructed within the range of 350–2500 nm, and their relationships with LNA were quantified. The results showed that both LNA and canopy hyperspectral reflectance in wheat changed with varied N rates, with consistent patterns across different cultivars and seasons. The sensitive spectral bands for LNA existed mainly within visible and near infrared regions. The best spectral indices for estimating LNA in wheat were found to be NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆), and the regression models based on the above four spectral indices were formulated as Y = 26.34x^1.887, Y = 5.095x − 6.040, Y = 0.609 e^3.008x and Y = 0.388x^1.260, respectively, with R² greater than 0.81. Furthermore, expanding the bandwidth of NDSI (R₈₆₀, R₇₂₀) and RSI (R₉₉₀, R₇₂₀) from 1 nm to 100 nm at 1 nm interval produced the LNA monitoring models with similar performance within about 33 nm and 23 nm bandwidth, respectively, over which the statistical parameters of the models became less stable. From testing of the derived equations, the model for LNA estimation on NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆) gave R² over 0.79 with more satisfactory performance than previously reported models and physical models in wheat. It can be concluded that the present hyperspectral parameters of NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆) can be reliably used for estimating LNA in winter wheat.     

Nonlinear regional economic dynamics: continuous-time specification,estimation and stability analysis

This paper presents an innovative approach to the study of regional economic dynamics within a nonlinear continuous-time econometric framework—a generalized specification of the Lotka–Volterra system of equations. This specification, which accounts for interdependent behavior of three industrial sectors and spillover effects of activities in neighboring regions, is employed in an analysis of five Italian regions between 1980 and 2003. For these regions, we report estimation results, characterize the varying systems dynamics, analyze the models’ local and global stability properties, and determine via sensitivity analyses which structural features appear to exert the greatest influence on these properties.

Study on Recovering the Earth＇s Potential Field Based on GOCE Gradiometry

Given the second radial derivative Vrr（P） ｜δs of the Earth＇s gravitational potential V（P） on the surface δS corresponding to the satellite altitude, by using the fictitious compress recovery method, a fictitious regular harmonic field rrVrr（P）^＊ and a fictitious second radial gradient field V：（P） in the domain outside an inner sphere Ki can be determined, which coincides with the real field V（P） in the domain outside the Earth. Vrr^＊（P）could be further expressed as a uniformly convergent expansion series in the domain outside the inner sphere, because rrV（P）^＊ could be expressed as a uniformly convergent spherical harmonic expansion series due to its regularity and harmony in that domain. In another aspect, the fictitious field V^＊（P） defined in the domain outside the inner sphere, which coincides with the real field V（P） in the domain outside the Earth, could be also expressed as a spherical harmonic expansion series. Then, the harmonic coefficients contained in the series expressing V^＊（P） can be determined, and consequently the real field V（P） is recovered. Preliminary simulation calculations show that the second radial gradient field Vrr（P） could be recovered based only on the second radial derivative V（P）｜δs given on the satellite boundary. Concerning the final recovery of the potential field V（P） based only on the boundary value Vrr （P）｜δs, the simulation tests are still in process.     

An assessment of chlorophyll-a concentration spatio-temporal variation using Landsat satellite data,in a small tropical reservoir

Traditional approaches to monitoring aquatic systems are often limited by the need for data collection which often is time-consuming, expensive and non-continuous. The aim of the study was to map the spatio-temporal chlorophyll-a concentration changes in Malilangwe Reservoir, Zimbabwe as an indicator of phytoplankton biomass and trophic state when the reservoir was full (year 2000) and at its lowest capacity (year 2011), using readily available Landsat multispectral images. Medium-spatial resolution (30 m) Landsat multispectral Thematic Mapper TM 5 and ETM+ images for May to December 1999–2000 and 2010–2011 were used to derive chlorophyll-a concentrations. In situ measured chlorophyll-a and total suspended solids (TSS) concentrations for 2011 were employed to validate the Landsat chlorophyll-a and TSS estimates. The study results indicate that Landsat-derived chlorophyll-a and TSS estimates were comparable with field measurements. There was a considerable wet vs. dry season differences in total chlorophyll-a concentration, Secchi disc depth, TSS and turbidity within the reservoir. Using Permutational multivariate analyses of variance (PERMANOVA) analysis, there were significant differences (p < 0.0001) for chlorophyll-a concentration among sites, months and years whereas TSS was significant during the study months (p < 0.05). A strong positive significant correlation among both predicted TSS vs. chlorophyll-a and measured vs. predicted chlorophyll-a and TSS concentrations as well as an inverse relationship between reservoir chlorophyll-a concentrations and water level were found (p < 0.001 in all cases). In conclusion, total chlorophyll-a concentration in Malilangwe Reservoir was successfully derived from Landsat remote sensing data suggesting that the Landsat sensor is suitable for real-time monitoring over relatively short timescales and for small reservoirs. Satellite data can allow for surveying of chlorophyll-a concentration in aquatic ecosystems, thus, providing invaluable data in data scarce (limited on site ground measurements) environments.     

On the multivariate total least-squares approach to empirical coordinate transformations. Three algorithms

The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (Y−E _Y = (X−E _X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E _Y and E _X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.     

Validation of the SEBS-derived sensible heat for FY3A/VIRR and TERRA/MODIS over an alpine grass region using LAS measurements

In this study, sensible heat (H) calculation using remote sensing data over an alpine grass landscape is conducted from May to September 2010, and the calculation is validated using LAS (large aperture scintillometers) measurements. Data from two remote sensing sensors (FY3A-VIRR and TERRA-MODIS) are analysed. Remote sensing data, combined with the ground meteorological observations (pressure, temperature, wind speed, humidity) are fed into the SEBS (Surface Energy Balance System) model. Then the VIRR-derived sensible heat (VIRR_SEBS_H) and MODIS-derived sensible heat (MODIS_SEBS_H) are compared with the LAS-estimated H, which are obtained at the respective satellite overpass time. Furthermore, the similarities and differences between the VIRR_SEBS_H and MODIS_SEBS_H values are investigated. The results indicate that VIRR data quality is as good as MODIS data for the purpose of H estimation. The root mean square errors (rmse) of the VIRR_SEBS_H and MODIS_SEBS_H values are 45.1098 W/m² (n = 64) and 58.4654 W/m² (n = 71), respectively. The monthly means of the MODIS_SEBS_H are marginally higher than those of VIRR_SEBS_H because the satellite overpass time of the TERRA satellite lags by 25 min to that of the FT3A satellite. Relative evaporation (EFr), which is more time-independent, shows a higher agreement between MODIS and VIRR. Many common features are shared by the VIRR_SEBS_H and the MODIS_SEBS_H, which can be attributed to the SEBS model performance. In May–June, H is over-estimated with more fluctuations and larger rmse, whereas in July–September, H is under-estimated with fewer fluctuations and smaller rmse. Sensitivity analysis shows that potential temperature gradient (delta_T) plays a dominant role in determining the magnitude and fluctuation of H. The largest rmse and over-estimation in H occur in June, which could most likely be attributed to high delta_T, high wind speed, and the complicated thermodynamic state during the transitional period when bare land transforms to dense vegetation cover.     

Plant phenolics and absorption features in vegetation reflectance spectra near 1.66 μm

Past laboratory and field studies have quantified phenolic substances in vegetative matter from reflectance measurements for understanding plant response to herbivores and insect predation. Past remote sensing studies on phenolics have evaluated crop quality and vegetation patterns caused by bedrock geology and associated variations in soil geochemistry. We examined spectra of pure phenolic compounds, common plant biochemical constituents, dry leaves, fresh leaves, and plant canopies for direct evidence of absorption features attributable to plant phenolics. Using spectral feature analysis with continuum removal, we observed that a narrow feature at 1.66 μm is persistent in spectra of manzanita, sumac, red maple, sugar maple, tea, and other species. This feature was consistent with absorption caused by aromatic CH bonds in the chemical structure of phenolic compounds and non-hydroxylated aromatics. Because of overlapping absorption by water, the feature was weaker in fresh leaf and canopy spectra compared to dry leaf measurements. Simple linear regressions of feature depth and feature area with polyphenol concentration in tea resulted in high correlations and low errors (% phenol by dry weight) at the dry leaf (r² = 0.95, RMSE = 1.0%, n = 56), fresh leaf (r² = 0.79, RMSE = 2.1%, n = 56), and canopy (r² = 0.78, RMSE = 1.0%, n = 13) levels of measurement. Spectra of leaves, needles, and canopies of big sagebrush and evergreens exhibited a weak absorption feature centered near 1.63 μm, short ward of the phenolic compounds, possibly consistent with terpenes. This study demonstrates that subtle variation in vegetation spectra in the shortwave infrared can directly indicate biochemical constituents and be used to quantify them. Phenolics are of lesser abundance compared to the major plant constituents but, nonetheless, have important plant functions and ecological significance. Additional research is needed to advance our understanding of the spectral influences of plant phenolics and terpenes relative to dominant leaf biochemistry (water, chlorophyll, protein/nitrogen, cellulose, and lignin).