基于ASTER数据的地物光谱比辐射率的反演研究

比辐射率光谱表征着一个物体的内部物理和化学特征,是定量化遥感的一个关键参数.本文利用ASTER数据的第10至14波段,根据其数据的特点,基于温度比辐射率分离算法的思想,提出了将比辐射率标准化法(Normalize Emissivity Method,NEM)、经验公式、比值法(Ratio Method)这几个模块结合起来,在迭代的基础上计算出比辐射率的新算法.本文简要分析了模型误差的主要来源,并且把本文的算法与简化的ASTER的TES算法进行了比较.分析表明本文的算法是可行、有效的,而且该算法简单,易于实现,在能够保证精度的情况下运算速度快;同时也说明ASTER遥感数据用于反演地物的比辐射率可以得到比较理想的结果,具有良好的应用前景.     

Radiometric normalization and image mosaic generation of ASTER thermal infrared data: An application to extensive sand sheets and dune fields

Data from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) have a significant advantage over previous datasets because of the combination of high spatial resolution (15-90 m) and enhanced multispectral capabilities, particularly in the thermal infrared (TIR) atmospheric window (8-12 μm) of the Earth where common silicate minerals are more easily identified. However, the 60 km swath width of ASTER can limit the effectiveness of accurately tracing large-scale features, such as eolian sediment transport pathways, over long distances. The primary goal of this paper is to describe a method for generating a seamless and radiometrically accurate ASTER TIR mosaic of atmospherically corrected radiance and from that, extract surface emissivity for arid lands, specifically, sand seas. The Gran Desierto in northern Sonora, Mexico was used as a test location for the radiometric normalization technique because of past remote sensing studies of the region, its compositional diversity, and its size. A linear approach was taken to transform adjacent image swaths into a direct linear relationship between image acquisition dates. Pseudo-invariant features (PIFs) were selected using a threshold of correlation between radiance values, and change-pixels were excluded from the linear regression used to determine correction factors. The degree of spectral correlation between overlapping pixels is directly related to the amount of surface change over time; therefore, the gain and offsets between scenes were based only on regions of high spectral correlation. The result was a series of radiometrically normalized radiance-at-surface images that were combined with a minimum of image edge seams present. These edges were subsequently blended to create the final mosaic. The advantages of this approach for TIR radiance (as opposed to emissivity) data include the ability to: (1) analyze data acquired on different dates (with potentially very different surface temperatures) as one seamless compositional dataset; (2) perform decorrelation stretches (DCS) on the entire dataset in order to identify and discriminate compositional units; and (3) separate brightness temperature from surface emissivity for quantitative compositional analysis of the surface, reducing seam-line error in the emissivity mosaic. The approach presented here is valid for any ASTER-related study of large geographic regions where numerous images spanning different temporal and atmospheric conditions are encountered.     

Validation of GOES and MODIS active fire detection products using ASTER and ETM+ data

In this study we implemented a comprehensive analysis to validate the MODIS and GOES satellite active fire detection products (MOD14 and WFABBA, respectively) and characterize their major sources of omission and commission errors which have important implications for a large community of fire data users. Our analyses were primarily based on the use of 30 m resolution ASTER and ETM+ imagery as our validation data. We found that at the 50% true positive detection probability mark, WFABBA requires four times more active fire area than is necessary for MOD14 to achieve the same probability of detection, despite the 16× factor separating the nominal spatial resolutions of the two products. Approximately 75% and 95% of all fires sampled were omitted by the MOD14 and WFABBA instantaneous products, respectively; whereas an omission error of 38% was obtained for WFABBA when considering the 30-minute interval of the GOES data. Commission errors for MOD14 and WFABBA were found to be similar and highly dependent on the vegetation conditions of the areas imaged, with the larger commission errors (approximately 35%) estimated over regions of active deforestation. Nonetheless, the vast majority (> 80%) of the commission errors were indeed associated with recent burning activity where scars could be visually confirmed in the higher resolution data. Differences in thermal dynamics of vegetated and non-vegetated areas were found to produce a reduction of approximately 50% in the commission errors estimated towards the hours of maximum fire activity (i.e., early-afternoon hours) which coincided with the MODIS/Aqua overpass. Lastly, we demonstrate the potential use of temporal metrics applied to the mid-infrared bands of MODIS and GOES data to reduce the commission errors found with the validation analyses.     

Remote sensing of structural complexity indices for habitat and species distribution modeling

Spatial distribution models are increasingly used in ecological studies, but are limited by the poor accuracy of remote sensing (RS) for mapping microhabitat (< 0.1 ha) features. Mapping accuracy can be improved by combining advanced RS image-processing techniques with microhabitat data expressed as a structural complexity index (SCI). To test this idea, we used principal components analysis (PCA) and an additive SCI method developed for forest ecology (calculated by re-scaling and summing representative structural variables) to summarize 13 microhabitat-scale (0.04 ha) vegetation structure attributes describing the rare mountain bongo antelope's (Tragelaphus eurycerus isaaci) habitat in Kenya's Aberdare mountains. Microhabitat data were collected in 127 plots: 37 related to bongo habitat use, 90 from 1 km-spaced grid points representing overall habitat availability and bongo non-presence. We then assessed each SCI's effectiveness for discerning microhabitat variability and bongo habitat selection, using Wilcoxon Rank Sum tests for differences in mean SCI scores among plots divided into 4 vegetation classes, and the Area Under the Curve (AUC) of receiver operating characteristics from logistic regressions. We also examined the accuracy of predicted SCI scores resulting from regression models based on variables derived from a) ASTER imagery processed with spectral mixture and texture analysis, b) an SRTM DEM and c) rainfall data, using the 90 grid plots for model training and the bongo plots as an independent test dataset. Of the five SCIs derived, two performed best: the PCA-derived Canopy Structure Index (CSI) and an additive index summarizing 8 structural variables (AI₈). CSI and AI₈ showed significant differences between 5 of 6 vegetation class pairs, strong abilities to distinguish bongo-selected from available habitat (AUCs = 0.71 (CSI); 0.70 (AI₈)), and predicted scores 60-110% more accurate than reported by other studies using RS to quantify individual microhabitat structural attributes (CSI model R² = 0.51, RMSE = 0.19 (training) and 0.21 (test); AI₈ model R² = 0.46, RMSE = 0.17 (training) and 0.19 (test)). Repeating the Wilcoxon tests and logistic regressions with RS-predicted SCI values showed that AI₈ most effectively preserved the patterns found with the observed SCIs. These results demonstrate that SCIs effectively characterize microhabitat structure and selection, and boost microhabitat mapping accuracy when combined with enhanced RS image-processing techniques. This approach can improve distribution models and broaden their applicability, makes RS more relevant to applied ecology, and shows that processing field data to be more compatible with RS can improve RS-based habitat mapping accuracy.     

Atmospheric correction of optical imagery from MODIS and Reanalysis atmospheric products

In this paper we analyze the differences obtained in the atmospheric correction of optical imagery covering bands located in the Visible and Near Infra-Red (VNIR), Short-Wave Infra-Red (SWIR) and Themal-Infrared (TIR) spectral regions when atmospheric profiles extracted from different sources are used. In particular, three sensors were used, Compact High Resolution Imaging Spectrometer (CHRIS), Advanced Spaceborne Thermal Emission and Reflection radiometer (ASTER) and Landsat5 Thematic Mapper (TM), whereas four atmospheric profiles sources were considered: i) local soundings launched near the sensor overpass time, ii) Moderate Resolution Radiometer (MODIS) atmospheric profiles product (MOD07), iii) Atmospheric Correction Parameter Calculator (ACPC) generated by the National Center for Environmental Prediction (NCEP) and iv) Modified Atmospheric Profiles from Reanalysis Information (MAPRI), which includes data from NCEP and National Center of Atmospheric Research (NCAR) Reanalysis project but interpolated to 34 atmospheric levels and resampled to 0.5° × 0.5°. MODIS aerosol product (MOD04) was also used to extract Aerosol Optical Thickness (AOT) values at 550 nm. Analysis was performed for three test dates (12th July 2003, 18th July 2004 and 13th July 2005) over an agricultural area in Spain. Results showed that air temperature vertical profiles were similar for the four sources, whereas dew point temperature profiles showed significant differences at some particular levels. Atmospheric profiles were used as input to MODTRAN4 radiative transfer code in order to compute atmospheric parameters involved in atmospheric correction, with the aim of retrieving surface reflectances in the case of VNIR and SWIR regions, and Land Surface Temperature (LST) in the case of the TIR region. For the VNIR and SWIR region, significant differences depending on the atmospheric profile used were not found, particularly in the Visible region in which the AOT content is the main parameter involved in the atmospheric correction. In the case of TIR, differences depending on the atmospheric profile used were appreciable, since in this case the main parameter involved in the atmospheric correction is the water vapor content, which depends on the vertical profile. In terms of LST retrieval from ASTER data (2004 test case), all profiles provided satisfactory results compared to the ones obtained when using a local sounding, with errors of 0.3 K for ACPC and MAPRI cases and 0.7 K for MOD07. When retrieving LST from TM data (2005 test case), errors for MOD07 and MAPRI were 0.6 and 0.9 K respectively, whereas ACPC provided an error of 2 K. The results presented in this paper show that the different atmospheric profile sources are useful for accurate atmospheric correction when local soundings are not available. In particular, MOD07 product provides atmospheric information at the highest spatial resolution, 5 km, although its use is limited from 2000 to present, whereas MAPRI provides historical information from 1970 to present, but at lower spatial resolution.     

ASTER Level 1A数据产品文件解析及应用

深入解析了ASTER数据的结构,研究如何正确读取数据中卫星的位置、速度、时间、姿态角、姿态变化率等与影像定位有关的数据,并研究了这些数据的变化规律。根据ASTER数据的特点,给出了该数据在辐射校正及影像定位方面的应用方法,使得影像能更好的在这些方面得到应用。     

ETM+和ASTER数据在遥感信息提取中的对比研究

遥感蚀变信息提取是找矿的一个重要技术手段。本文选择位于秘鲁南部阿雷基帕(AREQUIPA)省境内的萨卡纳(CERCANA)和伊卡(ICA)省境内的Moarcona铁矿区作为本文的两个研究区,从分析地物光谱出发,利用ETM+和ASTER卫星影像数据,通过主成份分析法和比值分析法分别对两个研究区进行泥化蚀变信息提取和铁染蚀变信息提取,并对两者的提取结果进行对比分析。最后结果表明,相较于ETM+数据,ASTER数据在矿化蚀变信息的提取方面具有更大的优势。     

Effects of soil composition and mineralogy on remote sensing of crop residue cover

The management of crop residues (non-photosynthetic vegetation) in agricultural fields influences soil erosion and soil carbon sequestration. Remote sensing methods can efficiently assess crop residue cover and related tillage intensity over many fields in a region. Although the reflectance spectra of soils and crop residues are often similar in the visible, near infrared, and the lower part of the shortwave infrared (400-1900 nm) wavelength region, specific diagnostic chemical absorption features are evident in the upper shortwave infrared (1900-2500 nm) region. Two reflectance band height indices used for estimating residue cover are the Cellulose Absorption Index (CAI) and the Lignin-Cellulose Absorption (LCA) index, both of which use reflectances in the upper shortwave infrared (SWIR). Soil mineralogy and composition will affect soil spectral properties and may limit the usefulness of these spectral indices in certain areas. Our objectives were to (1) identify minerals and soil components with absorption features in the 2000 nm to 2400 nm wavelength region that would affect CAI and LCA and (2) assess their potential impact on remote sensing estimates of crop residue cover. Most common soil minerals had CAI values ≤ 0.5, whereas crop residues were always > 0.5, allowing for good contrast between soils and residues. However, a number of common soil minerals had LCA values > 0.5, and, in some cases, the mineral LCA values were greater than those of the crop residues, which could limit the effectiveness of LCA for residue cover estimation. The LCA of some dry residues and live corn canopies were similar in value, unlike CAI. Thus, the Normalized Difference Vegetation Index (NDVI) or similar method should be used to separate out green vegetation pixels. Mineral groups, such as garnets and chlorites, often have wide ranges of CAI and LCA values, and thus, mineralogical analyses often do not identify individual mineral species required for precise CAI estimation. However, these methods are still useful for identifying mineral soils requiring additional scrutiny. Future advanced multi- and hyperspectral remote sensing platforms should include CAI bands to allow for crop residue cover estimation.     

A sequential model for disaggregating near-surface soil moisture observations using multi-resolution thermal sensors

A sequential model is developed to disaggregate microwave-derived soil moisture from 40 km to 4 km resolution using MODIS (Moderate Imaging Spectroradiometer) data and subsequently from 4 km to 500 m resolution using ASTER (Advanced Scanning Thermal Emission and Reflection Radiometer) data. The 1 km resolution airborne data collected during the three-week National Airborne Field Experiment 2006 (NAFE'06) are used to simulate the 40 km pixels, and a thermal-based disaggregation algorithm is applied using 1 km resolution MODIS and 100 m resolution ASTER data. The downscaled soil moisture data are subsequently evaluated using a combination of airborne and in situ soil moisture measurements. A key step in the procedure is to identify an optimal downscaling resolution in terms of disaggregation accuracy and sub-pixel soil moisture variability. Very consistent optimal downscaling resolutions are obtained for MODIS aboard Terra, MODIS aboard Aqua and ASTER, which are 4 to 5 times the thermal sensor resolution. The root mean square error between the 500 m resolution sequentially disaggregated and ground-measured soil moisture is 0.062 vol./vol. with a bias of − 0.045 vol./vol. and values ranging from 0.08 to 0.40 vol./vol.     

Evaluation of ASTER and MODIS land surface temperature and emissivity products using long-term surface longwave radiation observations at SURFRAD sites

Land surface temperature (LST) and emissivity are key parameters in estimating the land surface radiation budget, a major controlling factor of global climate and environmental change. In this study, Terra Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) and Aqua MODerate resolution Imaging Spectroradiometer (MODIS) Collection 5 LST and emissivity products are evaluated using long-term ground-based longwave radiation observations collected at six Surface Radiation Budget Network (SURFRAD) sites from 2000 to 2007. LSTs at a spatial resolution of 90 m from 197 ASTER images during 2000-2007 are directly compared to ground observations at the six SURFRAD sites. For nighttime data, ASTER LST has an average bias of 0.1 °C and the average bias is 0.3 °C during daytime. Aqua MODIS LST at 1 km resolution during nighttime retrieved from a split-window algorithm is evaluated from 2002 to 2007. MODIS LST has an average bias of − 0.2 °C. LST heterogeneity (defined as the Standard Deviation, STD, of ASTER LSTs in 1 × 1 km² region, 11 × 11 pixel in total) and instrument calibration error of pyrgeometer are key factors impacting the ASTER and MODIS LST evaluation using ground-based radiation measurements. The heterogeneity of nighttime ASTER LST is 1.2 °C, which accounts for 71% of the STD of the comparison, while the heterogeneity of the daytime LST is 2.4 °C, which accounts for 60% of the STD. Collection 5 broadband emissivity is 0.01 larger than that of MODIS Collection 4 products and ASTER emissivity. It is essential to filter out the abnormal low values of ASTER daily emissivity data in summer time before its application.