MODIS卫星数据地表反照率反演的简化模式

以内蒙西部地区的MODIS遥感图像数据和地表野外同步观测的光谱数据为例,在野外数据量较少且有定标数据的条件下反演地表反照率。使用6S大气1辐射传输模型进行大气校正,并通过MODTRAN4.0模型获取各波段地表入射光通量和窄波段的地表反照率;在窄波段反照率与宽波段反照率之间存在线性关系的前提下,以各波段的入射光通量占总入射通量的比例作为反演参数,实现窄波段到宽波段的反演。反演结果证明此方法简便可行。     

NOAA AVHRR land surface albedo algorithm development

The primary objective of this research is to develop a surface albedo model for the National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR). The primary test site is the Konza prairie, Kansas (U.S.A.), used by the International Satellite Land Surface Climatology Project (ISLSCP) in the First ISLSCP Field Experiment (FIFE). In this research, high spectral resolution field spectrometer data was analyzed to simulate AVHRR wavebands and to derive surface albedos. Development of a surface albedo algorithm was completed by analysing a combination of satellite, field spectrometer, and ancillary data. Estimated albedos from the field spectrometer data were compared to reference albedos derived using pyranometer data. Variations from surface anisotropy of reflected solar radiation were found to be the most significant albedo-related error. Additional error or sensitivity came from estimation of a shortwave mid-IR reflectance (1.3-4.0 mu m) using the AVHRR red and near-IR bands. Errors caused by the use of AVHRR spectral reflectance to estimate both a total visible (0.4-0.7 mu m) and near-IR (0.7-1.3 mu m) reflectance were small. The solar spectral integration, using the derived ultraviolet, visible, near-IR and SW mid-IR reflectivities, was not sensitive to many clear-sky changes in atmospheric properties and illumination conditions.     

Study of infrared reflection characteristics of aerial target using MODIS data on GPU

To study of the infrared signature of an aerial target, it is required to precisely model the background radiation. Simple empirical models or standard atmospheric models in LOWTRAN/MODTRAN were used in earlier studies. To further precisely model the thermal radiation of earth’s surface and atmospheric radiance/transmittance, the atmospheric profile, land surface temperature, and emissivity, the sea surface temperature retrieved from moderate-resolution imaging spectroradiometer, and the sea surface emissivity model developed by Wu and Smith are utilized in this study. Meanwhile, considering that the reflection of background radiation incident from different directions in each spectral wavelength can be calculated in parallel, implementations using open multi-processing and compute unified device architecture on a multi-core CPU and many-core graphics processing unit (GPU) are presented and speedups of 9\(\times\) and 258\(\times\) are obtained on a platform with dual Xeon E5-2652 CPU and an NVIDIA Tesla K80 GPU card, respectively.     

Mapping sediment-laden sea ice in the Arctic using AVHRR remote-sensing data: Atmospheric correction and determination of reflectances as a function of ice type and sediment load

Exploiting the fact that the spectral characteristics of light backscattered from sediment-laden ice differ substantially from those of clean ice and that sediment tends to accumulate at the ice surface during the first melt season, remote-sensing techniques provide a valuable tool for mapping the extent of particle-laden ice in the Arctic basin and assessing its particulate loading. This study considers two fundamental problems that still need to be addressed in order to make full use of satellite observations for this type of assessment: (i) the effects of the atmosphere on surface reflectances derived from radiances measured by the satellite sensor need to be quantified and ultimately corrected for, and (ii) the spectral reflectance of the ice surface as a function of particle loading and sub-pixel distribution needs to be determined in order to derive quantitative estimates from the at-sensor satellite signal. Here, spectral albedos have been computed for different ice surfaces of variable sediment load with a radiative transfer model for sea ice coupled with an optical model for particulates included in sea ice. In a second step, the role of the atmosphere in modulating the surface reflectance signal is assessed with the aid of an atmospheric radiative transfer model applied to a “standard” Arctic atmosphere and surface boundary conditions as prescribed by the sea ice radiative transfer model. A series of sensitivity studies helps assess differences between top-of-the-atmosphere and true surface reflectance and has been utilized to derive a look-up table for atmospheric correction of Advanced Very High Resolution Radiometer (AVHRR) data over sediment-laden sea ice surfaces. In particular, the effects of solar elevation, viewing geometry, and atmospheric properties are considered. The atmospheric corrections are necessary for certain geometries and surface types. Large discrepancies between raw and corrected data are particularly evident in the derived coverage of clean ice and ice with small sediment loading.     

An algorithm for the retrieval of 30-m snow-free albedo from Landsat surface reflectance and MODIS BRDF

We present a new methodology to generate 30-m resolution land surface albedo using Landsat surface reflectance and anisotropy information from concurrent MODIS 500-m observations. Albedo information at fine spatial resolution is particularly useful for quantifying climate impacts associated with land use change and ecosystem disturbance. The derived white-sky and black-sky spectral albedos may be used to estimate actual spectral albedos by taking into account the proportion of direct and diffuse solar radiation arriving at the ground. A further spectral-to-broadband conversion based on extensive radiative transfer simulations is applied to produce the broadband albedos at visible, near infrared, and shortwave regimes. The accuracy of this approach has been evaluated using 270 Landsat scenes covering six field stations supported by the SURFace RADiation Budget Network (SURFRAD) and Atmospheric Radiation Measurement Southern Great Plains (ARM/SGP) network. Comparison with field measurements shows that Landsat 30-m snow-free shortwave albedos from all seasons generally achieve an absolute accuracy of ±0.02-0.05 for these validation sites during available clear days in 2003-2005, with a root mean square error less than 0.03 and a bias less than 0.02. This level of accuracy has been regarded as sufficient for driving global and regional climate models. The Landsat-based retrievals have also been compared to the operational 16-day MODIS albedo produced every 8-days from MODIS on Terra and Aqua (MCD43A). The Landsat albedo provides more detailed landscape texture, and achieves better agreement (correlation and dynamic range) with in-situ data at the validation stations, particularly when the stations include a heterogeneous mix of surface covers.     

Leaf level detection of solar induced chlorophyll fluorescence by means of a subnanometer resolution spectroradiometer

A leaf level investigation on the spectral signature of Phaseolus vulgaris was undertaken by using a very high spectral resolution spectroradiometer featuring full width at half maximum of 0.06 nm and spectral range of 635.5-802.5 nm. High spectral resolution allows detection of leaf reflected and emitted radiance fields in two narrow absorption bands at 687 and 760 nm, respectively, where solar irradiance is strongly reduced owing to molecular oxygen absorption of the terrestrial atmosphere. The flux emitted due to chlorophyll fluorescence was measured using the Fraunhofer line depth principle by spectrally modelling the signal, capitalizing on the high resolution of the spectroradiometer devices. An experiment was conducted on two potted bean plants. One was maintained in good health for use as a reference while the other was treated with a photosystem II inhibitor. Collected spectra show that the fluorescence emission produces a pair of characteristic peaks superimposed on the typical leaf-specific reflectance curve. The magnitude of the fluorescence signal of the herbicide-treated leaf was four times greater than that of the control plant, thus indicating damage to the photosynthetic apparatus of the plant.     

Spectral calibration and atmospheric correction of ultra-fine spectral and spatial resolution remote sensing data. Application to CASI-1500 data

Imaging spectrometers operating in the solar spectrum measure the upwelling reflected solar radiation, and are an important tool in the bio/geochemical characterization of the Earth system. Surface reflectance is usually the starting point for the retrieval of biophysical parameters from remote measurements. Reliable radiometric and spectral calibrations and accurate atmospheric correction are mandatory in the interpretation of the surface reflectance. A complete surface reflectance retrieval scheme specifically designed for ultra-fine spectral resolution (bandwidth from 10 to 2 nm) and spatial resolution (pixel size less than 10 m) imaging spectrometers is presented in this work. The assessment of the spectral calibration is coupled to the removal of the atmospheric distortion so that maps of surface reflectance are derived, as well as columnar water vapor (CWV) maps, estimations of aerosol optical thickness (AOT) and updated sensor gain coefficients and spectral calibration. Radiative transfer calculations are performed by an optimized version of the MODTRAN4 code, which is run before processing each image. The method is tested against Compact Airborne Spectrographic Imager (CASI) 1500 images acquired during the ESA SENtinel-2 and FLuorescence EXperiment (SEN2FLEX) campaigns held in the Barrax (La Mancha, Spain) study site during June and July 2005. A peak-to-peak spectral shift variation of up to 2.3 nm is detected in CASI. Concerning atmospheric products and surface reflectance retrievals, an extensive validation is performed using ground-based measurements. A good correlation between ground measurements and CASI-derived AOT and CWV is found, with a Pearson correlation coefficient r² up to 0.71 and 0.74, respectively. The subsequent surface reflectance retrievals also hold a good correspondence with ground-based measurements.     

Albedo of vegetated surfaces: its variability with differing irradiances

The albedo of four vegetated surfaces was investigated to derive its variability with differing distributions of the irradiance. The results are based on measured values of the spectral biconical reflectance factor, which are combined with calculated spectral irradiances for low and high atmospheric turbidity. The solar zenith angle is varied from 0° to 80°. The derived spectral albedos are then integrated with respect to wavelength in order to achieve the albedo. It is found that the variability of the albedo with respect to the atmospheric turbidity is less than 0.01 in nearly all cases. The variability of the albedo with respect to the solar elevation angle, however, is larger than 0.02 in many cases. For solar elevation angles from 20° to 60°, the variability of the albedo of the four surfaces can be represented by a mean curve which fits the individual variabilities with an accuracy of 0.015.     

Normalized difference haze index: a new spectral index for monitoring urban air pollution

Haze is an undesirable meteorological and environmental phenomenon that can cause enormous harm to the environment, people's lives and health, and economic activities. This study focuses on Nanjing, Yangzhou and Suzhou in the lower reaches of the Yangtze River valley, China, which have suffered from the adverse effects of hazy weather in recent years. The spectral influence of haze on surface features was determined through analysis of the spectral variations of surface covers between hazy and haze-free days. On the basis of the established relationship, a new index called the normalized difference haze index (NDHI) was derived using moderate resolution imaging spectroradiometer (MODIS) data from winter 2008–2009. Correlation analysis of the derived NDHI with in situ observed PM₁₀ (particulate matter with diameter <10 μm) data reveals that NDHI over water bodies has a coefficient of 0.74, 0.57 and 0.67 with PM₁₀ for Nanjing, Yangzhou and Suzhou, respectively. It is concluded that NDHI is a reliable indicator of air pollution. It can be used as a new method of effectively monitoring air pollution from remotely sensed data.     

Impacts of spectral band difference effects on radiometric cross-calibration between satellite sensors in the solar-reflective spectral domain

In order for quantitative applications to make full use of the ever-increasing number of Earth observation satellite systems, data from the various imaging sensors involved must be on a consistent radiometric scale. This paper reports on an investigation of radiometric calibration errors due to differences in spectral response functions between satellite sensors when attempting cross-calibration based on near-simultaneous imaging of common ground targets in analogous spectral bands, a commonly used post-launch calibration methodology. Twenty Earth observation imaging sensors (including coarser and higher spatial resolution sensors) were considered, using the Landsat solar reflective spectral domain as a framework. Scene content was simulated using spectra for four ground target types (Railroad Valley Playa, snow, sand and rangeland), together with various combinations of atmospheric states and illumination geometries. Results were obtained as a function of ground target type, satellite sensor comparison, spectral region, and scene content. Overall, if spectral band difference effects (SBDEs) are not taken into account, the Railroad Valley Playa site is a “good” ground target for cross calibration between most but not all satellite sensors in most but not all spectral regions investigated. “Good” is defined as SBDEs within ± 3%. The other three ground target types considered (snow, sand and rangeland) proved to be more sensitive to uncorrected SBDEs than the RVPN site overall. The spectral characteristics of the scene content (solar irradiance, surface reflectance and atmosphere) are examined in detail to clarify why spectral difference effects arise and why they can be significant when comparing different imaging sensor systems. Atmospheric gas absorption features are identified as being the main source of spectral variability in most spectral regions. The paper concludes with recommendations on spectral data and tools that would facilitate cross-calibration between multiple satellite sensors.     

Seasonal landscape and spectral vegetation index dynamics in the Brazilian Cerrado: An analysis within the Large-Scale Biosphere-Atmosphere Experiment in Amazônia (LBA)

The Brazilian Cerrado biome comprises a vertically structured mosaic of grassland, shrubland, and woodland physiognomies with distinct phenology patterns. In this study, we investigated the utility of spectral vegetation indices in differentiating these physiognomies and in monitoring their seasonal dynamics. We obtained high spectral resolution reflectances, during the 2000 wet and dry seasons, over the major Cerrado types at Brasilia National Park (BNP) using the light aircraft-based Modland Quick Airborne Looks (MQUALS) package, consisting of a spectroradiometer and digital camera. Site-intensive biophysical and canopy structural measurements were made simultaneously at each of the Cerrado types including Cerrado grassland, shrub Cerrado, wooded Cerrado, Cerrado woodland, and gallery forest. We analyzed the spectral reflectance signatures, their first derivative analogs, and convolved spectral vegetation indices (VI) over all the Cerrado physiognomies. The high spectral resolution data were convolved to the MODIS, AVHRR, and ETM⁺ bandpasses and converted to the normalized difference vegetation index (NDVI) and the enhanced vegetation index (EVI) to simulate their respective sensors. Dry and wet season comparisons of the measured biophysical attributes were made with the reflectance and VI data for the different Cerrado physiognomies. We found that three major domains of Cerrado could be distinguished with the dry and wet season spectral signatures and vegetation indices. The EVI showed a higher sensitivity to seasonality than the NDVI; however, both indices displayed seasonal variations that were approximately one-half that found with the measured landscape green cover dynamics. Inter-sensor comparisons of seasonal dynamics, based on spectral bandpass properties, revealed the ETM⁺-simulated VIs had the best seasonal discrimination capability, followed by MODIS and AVHRR. Differences between sensor bandpass-derived VI values, however, varied with Cerrado type and between dry and wet seasons, indicating the need for inter-sensor VI translation equations for effective multi-sensor applications.     

Application of tilt modulation of the scatter cross-sections to distinguish between different non-Gaussian rough surfaces: unified full wave approach

The tilt modulations of the like- and cross-polarized cross-sections for arbitrarily oriented resolution cells are determined using the unified full-wave approach. A broad family of non-Gaussian rough surfaces characterized by the gamma surface height probability density functions of order K are considered. Furthermore a Pearson-Moskowitz surface height spectral density function is assumed for the sea surface. The surface height autocorrelation function is also assumed to be non-Gaussian. An arbitrarily oriented mean plane associated with the resolution cell is characterized by tilt angles in and perpendicular to a fixed reference plane of incidence. The ‘tilt modulation’ of the scattering cross-sections is determined as functions of the wavelength of the incident field 𝛌₀ and the backscatter angle 0ⁱ ₀. Each resolution cell represents the real (or synthetic) radar footprint. The size of the resolution cell, orientation, and the statistical characteristics of the non-Gaussian surface determine the radar cross-section. Radar remote sensing options are discussed to distinguish between random rough surfaces characterized by different probability density functions.     

Comparison between in situ and MODIS-derived spectral reflectances of snow and sea ice in the Amundsen Sea,Antarctica

The spectral albedo and directional reflectance of snow and sea ice were measured on sea ice of various types, including nilas, grey ice, pancake ice, multi-year pack ice, and land-fast ice in the Ross, Amundsen and Bellingshausen seas during a summer cruise in February through March 2000. Measurements were made using a spectroradiometer that has 512 channels in the visible and near-infrared (VNIR) region in which 16 of the 36 bands of the Moderate Resolution Imaging Spectroradiometer (MODIS) are covered. Directional reflectance is also retrieved from the MODIS radiometrically calibrated data (Level 1B) concurrently acquired from the first National Aeronautics and Space Administration (NASA) Earth Observing System (EOS) satellite, Terra. The locations of the ground ice stations are identified accurately on the MODIS images, and the spectral albedo and directional reflectance values at the 16 VNIR MODIS bands are extracted for those pixel locations. MODIS-derived reflectance is then corrected for the intervening atmosphere whose parameters are retrieved from the MODIS atmospheric profiles product (MOD07_L2) for the same granule. The corresponding spectral albedo and directional reflectance with the same viewing geometry as MODIS are derived from our ground-based spectroradiometer measurements. Because the footprint of the ground spectroradiometer is much smaller than the pixel sizes of MODIS images, the averaged spectral reflectance and albedo in the vicinity of each ice station are simulated for the corresponding MODIS pixel from the ground spectral measurements by weighting over different surface types (various ice types and open water). An accurate determination of ice concentration is important in deriving ground reflectance of a simulated pixel from in situ measurements. The best agreement between the in situ and MODIS measurements was found when the ground had 10/10 ice concentration (discrepancy range 0.2–11.69%, average 4.8%) or was oneice-type dominant (discrepancy range 0.8–16.9%, average 6.2%). The more homogeneous the ground surface and the less variable the ground topography, the more comparable between the in situ and satellite-derived reflectance is expected.     

Estimation of atmospheric particulate matter based on MODIS haze optimized transformation

This research is an attempt to simulate the relationship between haze optimized transformation (HOT) and aerosol optical thickness (AOT), and explore the influence of typical ground covers on this relationship using the 6S atmospheric radiative transfer model for the Chinese city of Nanjing. The HOT data were derived from moderate resolution imaging spectroradiometer (MODIS) satellite images recorded in the winter and spring seasons of December 2007–May 2009. They were analysed in conjunction with ground observed atmospheric particulate matter (PM) data so as to establish their quantitative relationship. Such a relationship may open a new avenue for remotely estimating atmospheric PM based on HOT. The results obtained indicate that HOT is related positively to AOT. This relationship is most accurately depicted by a second-order polynomial equation. Although built-up areas, waterbodies, and vegetation have differing HOT values, all of them bear a close and consistent correlation with AOT. HOT of built-up areas, waterbodies, and vegetative surfaces derived from MODIS images is also positively correlated with PM₁₀ (PM with diameter <10 μm), which was measured near the surface. The second-order polynomial equation has a coefficient of determination (R²) value of 0.375 (built-up), 0.344 (water), and 0.362 (vegetation) and a root mean squared error (RMSE) of 0.0258, 0.0264, and 0.0261, respectively. The closeness in R² value and RMSE for different ground covers suggests that correlation is marginally affected by the ground cover. It is thus concluded that HOT can be used as a reliable alternative for estimating PM₁₀ from MODIS data.     

Lidar and spectroradiometer measurements of atmospheric aerosol optical characteristics over an urban area in Sofia,Bulgaria

A series of campaigns involving a systematic investigation of the atmosphere over an urban area of Sofia city were carried out. A European Aerosol Research Lidar Network (EARLINET) scanning aerosol lidar, a spectroradiometer, a standard sun photometer and a ground meteorological station were used in the observations. Multiple aerosol layers of variable thickness (200–600 m) were observed systematically in the planetary boundary layer (PBL) over the study area and the experimental data were compared with theoretical data. A study of the optical characteristics of the atmospheric aerosol, including the extinction coefficient, aerosol optical depth (AOD) and Angstrom parameters α and β, was performed and their variations followed during the convective boundary layer (CBL) formation. Values of the AOD obtained using the different instruments during simultaneous measurements were compared. Preliminary results show that the AOD values recorded by the sun photometer and those calculated on the basis of the spectroradiometer data are higher than those retrieved from the lidar data. Determination of the atmospheric optical depth and extinction coefficient using a ground-based spectral instrument is a relatively simple and inexpensive method of monitoring the total aerosol content in the atmosphere as well as the air quality over the region.     

High resolution planetary albedos—values and variability

The spectral variability of high-resolution planetary albedo is shown to be a useful tool in surface-type analysis. The possibility of using comparison between planetary albedos in selected wavelength regions for cirrus cloud identification has been investigated. The results of an analysis of digital data streams from the polar orbiter satellite, TIROS-N, in three wavelength regions. Channel 1 (0.55–0.90 μm), Channel 2 (0.725–1.10 μm), and Channel 3 (3.55–3.93 μm) are described. These preliminary findings suggest that a spectral relationship exists for high-resolution planetary albedos in clear sky situations which permits surface type/cover discrimination. However the variability of planetary albedo in regions of cirrus cloud cover appears to be a strong function of the underlying surface or cloud type particularly when the cirrus cloud is tenous and/or it overlies highly reflective surfaces or clouds.     

Validation of MISR land surface broadband albedo

Land surface broadband albedo is a critical variable for many scientific applications. Due to the scarcity of spectral albedo measurements of the Earth's surface environments, it is useful to construct broadband albedo from spectral albedo data obtained by multi‐angle satellite observations. The Multi‐angle Imaging SpectroRadiometer (MISR) onboard NASA's Earth Observing System (EOS) Terra satellite provides land surface albedo products from multi‐angular observations; however, the products have not been comprehensively validated. We convert MISR spectral albedos to total shortwave albedos and validate them using ground measurements at different validation sites. For most surface types, a published narrowband to broadband conversion formula was used, but a new conversion formula for snow and ice covered sites is developed in this study where the spectral range of the instrument is different. Several comparisons are made: (1) between MISR directional‐hemispherical reflectance (DHR) or albedo and MODIS (Moderate Resolution Imaging Spectroradiometer) DHR; and (2) between MISR spectral DHR and bi‐hemispherical reflectance (BHR). The results show that: (1) both the value and the temporal trends of the MISR shortwave albedo and the ground measured shortwave albedo are in good agreement, with the exception of the snow and ice sites; (2) the MISR DHR conforms well to MODIS DHR; and (3) the values of MISR DHR and BHR are nearly identical.     

Topographic information of sand dunes as extracted from shading effects using Landsat images

Topographic variations affect the reflectance properties of the Earth's surface and are often removed in remote sensing studies, especially when significant terrain variations exist. In this study, however, we show that shading effects assessed by Landsat can be treated as a signal that stores important topographic information, especially when the spectral characteristics of a surface are homogenous. The coastal transverse dunes of the Ashdod area, and the desert linear dunes of Nizzana (both located in Israel), were selected to investigate the abovementioned idea. The dune heights in these areas are 10 m on average (relative to their surroundings) and have maximum slopes of 33°. An innovative method for extracting slope, aspect, and height data for sand dunes using Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) images was developed, based on the regularity and periodicity of dunes' landscapes. Using two Landsat images representing different sun zenith and azimuth angles, reflectance values of each image were converted to cos(i) values (i=incident angle between the surface normal and the solar beam radiation), applying histogram matching methods. The slope and aspect of each pixel were determined as those that give the best prediction of the observed value of cos(i). Height profiles were then extracted, using simple trigonometric relationships. The accuracies of heights and slopes along selected profile lines were to the order of 1 m and 3°, respectively (at a spatial resolution of 15 m). Best results were obtained when the images included one from the summer and the other from the winter, corresponding to maximum difference in solar zenith and azimuth angles. Errors in heights were attributed to surface heterogeneity (e.g., presence of biogenic soil crusts in the rainy season), geometric correction errors, cast shadows, and Bidirectional Reflectance Distribution Function (BRDF) effects. Comparison to Advanced Thermal Emission and Reflection Radiometer (ASTER) 3D information showed that the proposed method is better in representing the topographic variation of the area than the digital elevation model (DEM) produced by ASTER.     

Narrowband to broadband conversions of land surface albedo: II. Validation

In the first paper of this series, we developed narrowband to broadband albedo conversion formulae for a series of sensors. These formulae were determined based on extensive radiative transfer simulations under different surface and atmospheric conditions. However, it is important to validate the simulation results using independent measurement data. In this paper, the validation results for three broadband albedos (total-shortwave, -visible and -near-IR albedos) using ground measurement of several cover types on five different days at Beltsville, MD are presented. Results show that the conversion formulae in the previous paper are very accurate and the average residual standard errors of the resulting broadband albedos for most sensors are around 0.02, which meets the required accuracy for land surface modeling.     

基于高光谱遥感影像的大气纠正:用AVIRIS数据评价大气纠正模块FLAASH

高光谱遥感影像由于集中了高光谱分辨率和高空间分辨率的优点,在对地观测中具有不可替代的优势。实际应用当中,往往需要从遥感影像获取地物的地表反射率信息,这就要求首先从影像中去除大气的影响,即进行大气纠正及补偿。目前,对遥感影像进行大气纠正的算法有很多,详细介绍了基于遥感影像自身信息的大气纠正模块FLAASH(Fast Line of Sight Atmospheric Analysis f Spectral Hypercubes)所涉及的算法,并利用该模块对AVIRIS航空遥感影像进行了大气纠正, 对不同的结果进行了分析对比,从而对该算法进行了初步的评价。