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.     

基于ICESat-2的陕西省SRTM和ASTER GDEM数据对比分析研究

SRTM1 DEM和ASTER GDEM V3作为全球公开使用的最新DEM数据,为地形研究、GIS应用、资源管理、环境监测等领域提供重要的数据支持。本研究旨在分析它们在不同地形和地表覆盖物区域的误差分布情况,探究其在地形表达中的优势、劣势。选取陕西省特有地形,结合地形剖面和曲面误差方法,首次利用ICESat-2分析比较了SRTM1 DEM和ASTER GDEM V3数据的垂直精度在卫星轨道、地形地貌、土地覆盖中的误差分布。通过研究发现：STRM1 DEM和ASTER GDEM V3的总体平均误差分别为0.056 m和5.301 m,受轨道因子的影响,SRTM1 DEM垂直精度高于ASTER GDEM V3。坡度对数据的精度影响占比最大,林地的数据精度显示最低,地貌中平原地区影响最小,但误差随地形起伏不断增加。SRTM1 DEM误差曲面跨度小于ASTER GDEM V3,对于地形表达比ASTER GDEM V3更加精确。分析结果进一步验证了两种DEM数据测绘原理不同导致精度差异。     

MODIS spectral signals at a flux tower site: Relationships with high-resolution data, and CO2 flux and light use efficiency measurements

In this study, spectral indices were calculated from single date HyMap (3 m; 126 bands), Hyperion (30 m; 242 bands), ASTER (15/30 m; 9 bands), and a time series of MODIS nadir BRDF-adjusted reflectance (NBAR; 1 km, 7 bands) for a study area surrounding the Tumbarumba flux tower site in eastern Australia. The study involved: a) the calculation of a range of physiologically-based vegetation indices from ASTER, HyMap, Hyperion and MOD43B NBAR imagery over the flux tower site; b) comparison across scales between HyMap, Hyperion and MODIS for the normalized difference water index (NDWI) and the Red-Green ratio; c) analysis of relationships between tower-based flux and light use efficiency (LUE) measurements and seasonal and climatic constraints on growth; and d) examination of relationships between fluxes, LUE and time series of NDVI, NDWI and Red-Green ratio. Strong seasonal patterns of variation were observed in NDWI and Red/Green ratio from MODIS NBAR. Correlations between fine (3 and 30 m) and coarse (1 km) scale indices for a small region around the flux tower site were moderately good for Red/Green ratio, but poor for NDWI. Hymap NDWI values for the understorey canopy were much lower than values for the tree canopy. MODIS NDWI was negatively correlated with CO₂ fluxes during warm and cool seasons. The correlation indicated that surface reflectance, affected by a spectrally bright grassland understorey canopy, was decoupled from growth of trees with access to deep soil moisture. The application of physiologically-based indices at earth observation scale requires careful attention to applicability of band configurations, contribution of vegetation components to reflectance signals, mechanistic relationships between biochemical processes and spectral indexes, and incorporation of ancillary information into any analysis.     

Temperature and emissivity separation from ASTER data for low spectral contrast surfaces

The performance of Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) thermal infrared (TIR) data product algorithms was evaluated for low spectral contrast surfaces (such as vegetation and water) in a test site close to Valencia, Spain. Concurrent ground measurements of surface temperature, emissivity, and atmospheric radiosonde profiles were collected at the test site, which is a thermally homogeneous area of rice crops with nearly full vegetation cover in summer. Using the ground data and the local radiosonde profiles, at-sensor radiances were simulated for the ASTER TIR channels and compared with L1B data (calibrated at-sensor radiances) showing discrepancies up to 3% in radiance for channel 10 at 8.3 μm (equivalently, 2.5 °C in temperature or 7% in emissivity), whereas channel 13 (10.7 μm) yielded a closer agreement (maximum difference of 0.5% in radiance or 0.4 °C in temperature). We also tested the ASTER standard products of land surface temperature (LST) and spectral emissivity generated with the Temperature-Emissivity Separation (TES) algorithm with standard atmospheric correction from both global data assimilation system profiles and climatology profiles. These products showed anomalous emissivity spectra with lower emissivity values and larger spectral contrast (or maximum-minimum emissivity difference, MMD) than expected, and as a result, overestimated LSTs. In this work, a scene-based procedure is proposed to obtain more accurate MMD estimates for low spectral contrast materials (vegetation and water) and therefore a better retrieval of LST and emissivity with the TES algorithm. The method uses various gray-bodies or near gray-bodies with known emissivities and assumes that the calibration and atmospheric correction performed with local radiosonde data are accurate for channel 13. Taking the channel 13 temperature (atmospherically and emissivity corrected) as the true LST, the radiances for the other channels were simulated and used to derive linear relationships between ASTER digital numbers and at-ground radiances for each channel. The TES algorithm was applied to the adjusted radiances and the resulting products showed a closer agreement with the ground measurements (differences lower than 1% in channel 13 emissivities and within ± 0.3 °C in temperature for rice and sea pixels).     

提取黔西南红粘土的ASTER波段比值技术研究

红粘土型金矿是黔西南的一个重要的金矿类型,我们以往的研究实践表明:Aleks Kalinowski等[1]研制的红粘土型金矿ASTER波段比值4/5不适合黔西南红粘土的信息提取,因此,需要研究适合该地区的ASTER波段比值。ASTER波段比值组合2/1、4/3和(5+7)/6彩色合成后被成功地应用到研究区的红粘土提取中。比值的选择是建立在对矿区红粘土光谱曲线特征的分析和ASTER图像预处理的基础上,该比值组合能够有效提取出ASTER图像上的矿区红粘土信息,并可以扩展应用到周边地区,寻找到更多含金红粘土的分布区,这充分显示了ASTER波段比值彩色合成技术在红粘土提取中的作用。     

Object-based crop identification using multiple vegetation indices, textural features and crop phenology

Crop identification on specific parcels and the assessment of soil management practices are important for agro-ecological studies, greenhouse gas modeling, and agrarian policy development. Traditional pixel-based analysis of remotely sensed data results in inaccurate identification of some crops due to pixel heterogeneity, mixed pixels, spectral similarity, and crop pattern variability. These problems can be overcome using object-based image analysis (OBIA) techniques, which incorporate new spectral, textural and hierarchical features after segmentation of imagery. We combined OBIA and decision tree (DT) algorithms to develop a methodology, named Object-based Crop Identification and Mapping (OCIM), for a multi-seasonal assessment of a large number of crop types and field status.In our approach, we explored several vegetation indices (VIs) and textural features derived from visible, near-infrared and short-wave infrared (SWIR) bands of ASTER satellite scenes collected during three distinct growing-season periods (mid-spring, early-summer and late-summer). OCIM was developed for 13 major crops cultivated in the agricultural area of Yolo County in California, USA. The model design was built in four different scenarios (combinations of three or two periods) by using two independent training and validation datasets and the best DTs resulted in an error rate of 9% for the three-period model and between 12 and 16% for the two-period models. Next, the selected DT was used for the thematic classification of the entire cropland area and mapping was then evaluated applying the confusion matrix method to the independent testing dataset that reported 79% overall accuracy. OCIM detected intra-class variations in most crops attributed to variability from local crop calendars, tree-orchard structures and land management operations. Spectral variables (based on VIs) contributed around 90% to the models, although textural variables were necessary to discriminate between most of the permanent crop-fields (orchards, vineyard, alfalfa and meadow). Features extracted from late-summer imagery contributed around 60% in classification model development, whereas mid-spring and early-summer imagery contributed around 30 and 10%, respectively. The Normalized Difference Vegetation Index (NDVI) was used to identify the main groups of crops based on the presence and vigor of green vegetation within the fields, contributing around 50% to the models. In addition, other VIs based on SWIR bands were also crucial to crop identification because of their potential to detect field properties like moisture, vegetation vigor, non-photosynthetic vegetation and bare soil. The OCIM method was built using interpretable rules based on physical properties of the crops studied and it was successful for object-based feature selection and crop identification.     

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数据中提取蚀变异常

选择朱拉扎嘎金矿为试验区,ASTER数据为数据源,利用光谱匹配技术进行蚀变异常提取实验。重点阐述光谱匹配技术(SM)应用于ASTER数据提取蚀变异常工作中的两项关键技术:内部平均相对反射率法(IARR)大气校正和典型蚀变矿物波谱对应于ASTER数据波段重采样。最终通过光谱匹配技术提取了4组蚀变异常,分别为褐铁矿化蚀变异常、高岭土化蚀变异常、碳酸盐化蚀变组合异常和含OH-蚀变组合异常。     

Spectral characterization and ASTER-based lithological mapping of an ophiolite complex: A case study from Neyriz ophiolite, SW Iran

The Neyriz ophiolite occurs along the Zagros suture zone in SW Iran, and is part of a 3000-km obduction belt thrusting over the edge of the Arabian continent during the late Cretaceous. This complex typically consists of altered dunites and peridotites, layered and massive gabbros, sheeted dykes and pillow lavas, and a thick sequence of radiolarites. Reflectance and emittance spectra of Neyriz ophiolite rock samples were measured in the laboratory and their spectra were used as endmembers in a spectral feature fitting (SFF) algorithm. Laboratory spectral reflectance measurements of field samples showed that in the visible through shortwave infrared (VNIR-SWIR) wavelength region the ultramafic and gabbroic rocks are characterized by ferrous-iron and Fe, MgOH spectral features, and the pillow lavas and radiolarites are characterized by spectral features of ferric-iron and AlOH. The laboratory spectral emittance spectra also revealed a wide wavelength range of SiO spectral features for the ophiolite rock units. After continuum removal of the spectra, the SFF classification method was applied to the VNIR + SWIR 9-band stack, and to the 11-band data set of SWIR and TIR data sets of the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) sensor, using field spectra as training sets for evaluating the potential of these data sets in discriminating ophiolite rock units. Output results were compared with the geological map of the area and field observations, and were assessed by the use of confusion matrices. The assessment showed, in terms of kappa coefficient, that the SFF classification method with continuum removal applied to the SWIR data achieved excellent results, which were distinctively better than those obtained using VNIR + SWIR data and TIR data alone.     

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

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