Atmospheric correction of ocean color sensors: computing atmospheric diffuse transmittance

Using the reciprocal equation derived by Yang and Gordon [Appl. Opt. 36, 7887-7897 (1997)] for atmospheric diffuse transmittance of the ocean-atmosphere system, I examined the accuracy of an analytical equation proposed by Gordon et al. [Appl. Opt. 22, 20-36 (1983)] in computing the atmospheric diffuse transmittance for wavelengths from 412 to 865 nm for both a pure Rayleigh and a two-layer Rayleigh-aerosol atmosphere overlying a flat Fresnel-reflecting ocean surface. It was found that for viewing angles up to approximately 40 degrees , the analytical formula produces errors usually between 2% and 3% for nonabsorbing and weakly absorbing aerosols and for aerosol optical thicknesses tau(a) 相似文献   

Role of oceanic air bubbles in atmospheric correction of ocean color imagery

Ocean color is the radiance that emanates from the ocean because of scattering by chlorophyll pigments and particles of organic and inorganic origin. Air bubbles in the ocean also scatter light and thus contribute to the water-leaving radiance. This additional water-leaving radiance that is due to oceanic air bubbles could violate the black pixel assumption at near-infrared wavelengths and be attributed to chlorophyll in the visible. Hence, the accuracy of the atmospheric correction required for the retrieval of ocean color from satellite measurements is impaired. A comprehensive radiative transfer code for the coupled atmosphere--ocean system is employed to assess the effect of oceanic air bubbles on atmospheric correction of ocean color imagery. This effect is found to depend on the wavelength-dependent optical properties of oceanic air bubbles as well as atmospheric aerosols.     

Atmospheric correction of ocean color imagery: use of the junge power-law aerosol size distribution with variable refractive index to handle aerosol absorption

When strongly absorbing aerosols are present in the atmosphere, the usual two-step procedure of processing ocean color data-(1) atmospheric correction to provide the water-leaving reflectance (rho(w)), followed by (2) relating rho(w) to the water constituents-fails and simultaneous estimation of the ocean and aerosol optical properties is necessary. We explore the efficacy of using a simple model of the aerosol-a Junge power-law size distribution consisting of homogeneous spheres with arbitrary refractive index-in a nonlinear optimization procedure for estimating the relevant oceanic and atmospheric parameters for case 1 waters. Using simulated test data generated from more realistic aerosol size distributions (sums of log-normally distributed components with different compositions), we show that the ocean's pigment concentration (C) can be retrieved with good accuracy in the presence of weakly or strongly absorbing aerosols. However, because of significant differences in the scattering phase functions for the test and power-law distributions, large error is possible in the estimate of the aerosol optical thickness. The positive result for C suggests that the detailed shape of the aerosol-scattering phase function is not relevant to the atmospheric correction of ocean color sensors. The relevant parameters are the aerosol single-scattering albedo and the spectral variation of the aerosol optical depth. We argue that the assumption of aerosol sphericity should not restrict the validity of the algorithm and suggest an avenue for including colored aerosols, e.g., wind-blown dust, in the procedure. A significant advantage of the new approach is that realistic multicomponent aerosol models are not required for the retrieval of C.     

Atmospheric correction algorithm for hyperspectral remote sensing of ocean color from space

Existing atmospheric correction algorithms for multichannel remote sensing of ocean color from space were designed for retrieving water-leaving radiances in the visible over clear deep ocean areas and cannot easily be modified for retrievals over turbid coastal waters. We have developed an atmospheric correction algorithm for hyperspectral remote sensing of ocean color with the near-future Coastal Ocean Imaging Spectrometer. The algorithm uses lookup tables generated with a vector radiative transfer code. Aerosol parameters are determined by a spectrum-matching technique that uses channels located at wavelengths longer than 0.86 mum. The aerosol information is extracted back to the visible based on aerosol models during the retrieval of water-leaving radiances. Quite reasonable water-leaving radiances have been obtained when our algorithm was applied to process hyperspectral imaging data acquired with an airborne imaging spectrometer.     

Atmospheric correction of ocean color sensors: analysis of the effects of residual instrument polarization sensitivity

We provide an analysis of the influence of instrument polarization sensitivity on the radiance measured by spaceborne ocean color sensors. Simulated examples demonstrate the influence of polarization sensitivity on the retrieval of the water-leaving reflectance rho(w). A simple method for partially correcting for polarization sensitivity-replacing the linear polarization properties of the top-of-atmosphere reflectance with those from a Rayleigh-scattering atmosphere-is provided and its efficacy is evaluated. It is shown that this scheme improves rho(w) retrievals as long as the polarization sensitivity of the instrument does not vary strongly from band to band. Of course, a complete polarization-sensitivity characterization of the ocean color sensor is required to implement the correction.     

Effects of stratospheric aerosols and thin cirrus clouds on the atmospheric correction of ocean color imagery: simulations

Using simulations, we determine the influence of stratospheric aerosol and thin cirrus clouds on the performance of the proposed atmospheric correction algorithm for the moderate resolution imaging spectroradiometer (MODIS) data over the oceans. Further, we investigate the possibility of using the radiance exiting the top of the atmosphere in the 1.38-microm water vapor absorption band to remove their effects prior to application of the algorithm. The computations suggest that for moderate optical thicknesses in the stratosphere, i.e., tau(s) < or approximately 0.15, the stratospheric aerosol-cirrus cloud contamination does not seriously degrade the MODIS except for the combination of large (approximately 60 degrees) solar zenith angles and large (approximately 45 degrees) viewing angles, for which multiple-scattering effects can be expected to be particularly severe. The performance of a hierarchy of stratospheric aerosol/cirrus cloud removal procedures for employing the 1.38-microm water vapor absorption band to correct for stratospheric aerosol/cirrus clouds, ranging from simply subtracting the reflectance at 1.38 microm from that in the visible bands, to assuming that their optical properties are known and carrying out multiple-scattering computations of their effect by the use of the 1.38-microm reflectance-derived concentration, are studied for stratospheric aerosol optical thicknesses at 865 nm as large as 0.15 and for cirrus cloud optical thicknesses at 865 nm as large as 1.0. Typically, those procedures requiring the most knowledge concerning the aerosol optical properties (and also the most complex) performed the best; however, for tau(s) < or approximately 0.15, their performance is usually not significantly better than that found by applying the simplest correction procedure. A semiempirical algorithm is presented that permits accurate correction for thin cirrus clouds with tau(s) as large as unity when an accurate estimate of the cirrus cloud scattering phase function is provided, and as large as 0.5 when a coarse approximation to the phase function is used. Given estimates of the stratospheric aerosol optical properties, the implementation of the algorithm by using a set of lookup tables appears to be straightforward.     

Moderate-resolution imaging spectroradiometer ocean color polarization correction

The polarization correction for the Moderate-Resolution Imaging Spectroradiometer (MODIS) instruments on the Terra and Aqua satellites is described. The focus is on the prelaunch polarization characterization and on the derivation of polarization correction coefficients for the processing of ocean color data. The effect of the polarization correction is demonstrated. The radiances at the top of the atmosphere need to be corrected by as much as 3.2% in the 412 nm band. The effect on the water-leaving radiances can exceed 50%. The polarization correction produces good agreement of the MODIS Aqua water-leaving radiance time series with data from another, independent satellite-based ocean color sensor, the Sea-Viewing Wide Field-of-View Sensor (SeaWiFS).     

Sensor-independent approach to the vicarious calibration of satellite ocean color radiometry

The retrieval of ocean color radiometry from space-based sensors requires on-orbit vicarious calibration to achieve the level of accuracy desired for quantitative oceanographic applications. The approach developed by the NASA Ocean Biology Processing Group (OBPG) adjusts the integrated instrument and atmospheric correction system to retrieve normalized water-leaving radiances that are in agreement with ground truth measurements. The method is independent of the satellite sensor or the source of the ground truth data, but it is specific to the atmospheric correction algorithm. The OBPG vicarious calibration approach is described in detail, and results are presented for the operational calibration of SeaWiFS using data from the Marine Optical Buoy (MOBY) and observations of clear-water sites in the South Pacific and southern Indian Ocean. It is shown that the vicarious calibration allows SeaWiFS to reproduce the MOBY radiances and achieve good agreement with radiometric and chlorophyll a measurements from independent in situ sources. We also find that the derived vicarious gains show no significant temporal or geometric dependencies, and that the mission-average calibration reaches stability after approximately 20-40 high-quality calibration samples. Finally, we demonstrate that the performance of the vicariously calibrated retrieval system is relatively insensitive to the assumptions inherent in our approach.     

Atmospheric correction of SeaWiFS imagery for turbid coastal and inland waters: comment

Ruddick et al. [Appl. Opt. 39, 897 (2000)] extended the standard SeaWiFS atmospheric-correction algorithm for use over turbid coastal and inland waters. However, Ruddick's method is based on the assumption of a spatially homogeneous constant ratio for the water-leaving reflectances normalized by the sun-sea atmospheric transmittance at 765 and 865 nm. Such first-order b(b)/a model-based assumption can result in an inaccuracy for highly turbid water. Using the first- and second-order b(b)/(a + b(b)) models as well as the second-order b(b)/a model (which, more realistically, do not assume spatial homogeneity ratio), we suggest using the modified assumption, R(8)(-1) = alpha0 R(7)(-1) + (l1Q)(-1) (1 - alpha0), instead of Ruddick's assumption, in SeaWiFS atmospheric-correction algorithms for highly turbid waters.     

Atmospheric correction of SeaWiFS imagery for turbid coastal and inland waters

The standard SeaWiFS atmospheric correction algorithm, designed for open ocean water, has been extended for use over turbid coastal and inland waters. Failure of the standard algorithm over turbid waters can be attributed to invalid assumptions of zero water-leaving radiance for the near-infrared bands at 765 and 865 nm. In the present study these assumptions are replaced by the assumptions of spatial homogeneity of the 765:865-nm ratios for aerosol reflectance and for water-leaving reflectance. These two ratios are imposed as calibration parameters after inspection of the Rayleigh-corrected reflectance scatterplot. The performance of the new algorithm is demonstrated for imagery of Belgian coastal waters and yields physically realistic water-leaving radiance spectra. A preliminary comparison with in situ radiance spectra for the Dutch Lake Markermeer shows significant improvement over the standard atmospheric correction algorithm. An analysis is made of the sensitivity of results to the choice of calibration parameters, and perspectives for application of the method to other sensors are briefly discussed.     

Sources and assumptions for the vicarious calibration of ocean color satellite observations

Spaceborne ocean color sensors require vicarious calibration to sea-truth data to achieve accurate water-leaving radiance retrievals. The assumed requirements of an in situ data set necessary to achieve accurate vicarious calibration were set forth in a series of papers and reports developed nearly a decade ago, which were embodied in the development and site location of the Marine Optical BuoY (MOBY). Since that time, NASA has successfully used data collected by MOBY as the sole source of sea-truth data for vicarious calibration of the Sea-viewing Wide field-of-view Sensor (SeaWiFS) and Moderate Resolution Imaging Spectroradiometer instruments. In this paper, we make use of the 10-year, global time series of SeaWiFS measurements to test the sensitivity of vicarious calibration to the assumptions inherent in the in situ requirements (e.g., very low chlorophyll waters, hyperspectral measurements). Our study utilized field measurements from a variety of sources with sufficient diversity in data collection methods and geophysical variability to challenge those in situ restrictions. We found that some requirements could be relaxed without compromising the ability to vicariously calibrate to the level required for accurate water-leaving radiance retrievals from satellite-based sensors.     

Dynamic range and sensitivity requirements of satellite ocean color sensors: learning from the past

Sensor design and mission planning for satellite ocean color measurements requires careful consideration of the signal dynamic range and sensitivity (specifically here signal-to-noise ratio or SNR) so that small changes of ocean properties (e.g., surface chlorophyll-a concentrations or Chl) can be quantified while most measurements are not saturated. Past and current sensors used different signal levels, formats, and conventions to specify these critical parameters, making it difficult to make cross-sensor comparisons or to establish standards for future sensor design. The goal of this study is to quantify these parameters under uniform conditions for widely used past and current sensors in order to provide a reference for the design of future ocean color radiometers. Using measurements from the Moderate Resolution Imaging Spectroradiometer onboard the Aqua satellite (MODISA) under various solar zenith angles (SZAs), typical (L_typical) and maximum (L_max) at-sensor radiances from the visible to the shortwave IR were determined. The L_typical values at an SZA of 45° were used as constraints to calculate SNRs of 10 multiband sensors at the same L_typical radiance input and 2 hyperspectral sensors at a similar radiance input. The calculations were based on clear-water scenes with an objective method of selecting pixels with minimal cross-pixel variations to assure target homogeneity. Among the widely used ocean color sensors that have routine global coverage, MODISA ocean bands (1?km) showed 2-4 times higher SNRs than the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) (1?km) and comparable SNRs to the Medium Resolution Imaging Spectrometer (MERIS)-RR (reduced resolution, 1.2?km), leading to different levels of precision in the retrieved Chl data product. MERIS-FR (full resolution, 300?m) showed SNRs lower than MODISA and MERIS-RR with the gain in spatial resolution. SNRs of all MODISA ocean bands and SeaWiFS bands (except the SeaWiFS near-IR bands) exceeded those from prelaunch sensor specifications after adjusting the input radiance to L_typical. The tabulated L_typical, L_max, and SNRs of the various multiband and hyperspectral sensors under the same or similar radiance input provide references to compare sensor performance in product precision and to help design future missions such as the Geostationary Coastal and Air Pollution Events (GEO-CAPE) mission and the Pre-Aerosol-Clouds-Ecosystems (PACE) mission currently being planned by the U.S. National Aeronautics and Space Administration (NASA).     

Aerosol polarization effects on atmospheric correction and aerosol retrievals in ocean color remote sensing

The current ocean color data processing system for the Sea-viewing Wide Field-of-View Sensor (SeaWiFS) and the moderate resolution imaging spectroradiometer (MODIS) uses the Rayleigh lookup tables that were generated using the vector radiative transfer theory with inclusion of the polarization effects. The polarization effects, however, are not accounted for in the aerosol lookup tables for the ocean color data processing. I describe a study of the aerosol polarization effects on the atmospheric correction and aerosol retrieval algorithms in the ocean color remote sensing. Using an efficient method for the multiple vector radiative transfer computations, aerosol lookup tables that include polarization effects are generated. Simulations have been carried out to evaluate the aerosol polarization effects on the derived ocean color and aerosol products for all possible solar-sensor geometries and the various aerosol optical properties. Furthermore, the new aerosol lookup tables have been implemented in the SeaWiFS data processing system and extensively tested and evaluated with SeaWiFS regional and global measurements. Results show that in open oceans (maritime environment), the aerosol polarization effects on the ocean color and aerosol products are usually negligible, while there are some noticeable effects on the derived products in the coastal regions with nonmaritime aerosols.     

Analytic model of ocean color

Ocean color is determined by spectral variations in reflectance at the sea surface. In the analytic model presented here, reflectance at the sea surface is estimated with the quasi-single-scattering approximation that ignores transspectral processes. The analytic solutions we obtained are valid for a vertically homogeneous water column. The solution provides a theoretical expression for the dimensionless, quasi-stable parameter (r), with a value of ~0.33, that appears in many models in which reflectance at the sea surface is expressed as a function of absorption coefficient (a) and backscattering coefficient (b(b)). In the solution this parameter is represented as a function of the mean cosines for downwelling and upwelling irradiances and as the ratio of the upward-scattering coefficient to the backscattering coefficient. Implementation of the model is discussed for two cases: (1) that in which molecular scattering is the main source of upwelling light, and (2) that in which particle scattering is responsible for all the upwelled light. Computations for the two cases are compared with Monte Carlo simulations, which accounts for processes not considered in the analytic model (multiple scattering, and consequent depth-dependent changes in apparent optical properties). The Monte Carlo models show variations in reflectance with the zenith angle of the incident light. The analytic model can be used to reproduce these variations fairly well for the case of molecular scattering. For the particle-scattering case also, the analytic and Monte Carlo models show similar variations in r with zenith angle. However, the analytic model (as implemented here) appears to underestimate r when the value of the backscattering coefficient b(b) increases relative to the absorption coefficient a. The errors also vary with the zenith angle of the incident light field, with the maximum underestimate being approximately 0.06 (equivalent to relative errors from 12 to 17%) for the range of b(b)/a studied here. One implication of this result is that the model could also be used to obtain approximate solutions for the Q factor, defined for a given look angle as the ratio of the upwelling irradiance at the surface to the upwelling radiance at the surface at that angle. This is a quantity that is important in remote-sensing applications of ocean-color models. An advantage of the model discussed here is that its implementation requires inputs that are in principle accessible only in a remote-sensing context.     

Scene-based nonuniformity correction and enhancement: pixel statistics and subpixel motion

We propose a framework for scene-based nonuniformity correction (NUC) and nonuniformity correction and enhancement (NUCE) that is required for focal-plane array-like sensors to obtain clean and enhanced-quality images. The core of the proposed framework is a novel registration-based nonuniformity correction super-resolution (NUCSR) method that is bootstrapped by statistical scene-based NUC methods. Based on a comprehensive imaging model and an accurate parametric motion estimation, we are able to remove severe/structured nonuniformity and in the presence of subpixel motion to simultaneously improve image resolution. One important feature of our NUCSR method is the adoption of a parametric motion model that allows us to (1) handle many practical scenarios where parametric motions are present and (2) carry out perfect super-resolution in principle by exploring available subpixel motions. Experiments with real data demonstrate the efficiency of the proposed NUCE framework and the effectiveness of the NUCSR method.     

Initial analysis of ocean color data from the ocean color and temperature scanner. I. Imagery analysis

The ocean color and temperature scanner (OCTS) collected global ocean color data from November 1996 to June 1997. Analyses of OCTS imagery indicate three features that impair scientific research uses: (1) band misalignments, (2) image striping, and (3) image noise. These are due to (1) band offsets in the sensor design, (2) detector radiometric response variability, and (3) primarily cloud contamination, respectively. Methods are analyzed to ameliorate the effects of each that facilitate use of OCTS ocean color data for quantitative scientific analyses.     

Correction of satellite imagery over mountainous terrain

A method for the radiometric correction of satellite imagery over mountainous terrain has been developed to remove atmospheric and topographic effects. The algorithm accounts for horizontally varying atmospheric conditions and also includes the height dependence of the atmospheric radiance and transmittance functions to simulate the simplified properties of a three-dimensional atmosphere. A database has been compiled that contains the results of radiative transfer calculations (atmospheric transmittance, path radiance, direct and diffuse solar flux) for a wide range of weather conditions. A digital elevation model is used to obtain information about surface elevation, slope, and orientation. Based on the Lambertian assumption the surface reflectance in rugged terrain is calculated for the specified atmospheric conditions. Regions with extreme illumination geometries sensitive to BRDF effects can be optionally processed separately. The method is restricted to high spatial resolution satellite sensors with a small swath angle such as the Landsat thematic mapper and Systeme pour l'Observation de la Terre high resolution visible, since some simplifying assumptions were made to reduce the required image processing time.     

Estimation of saturated pixel values in digital color imaging

Pixel saturation, in which the incident light at a pixel causes one of the color channels of the camera sensor to respond at its maximum value, can produce undesirable artifacts in digital color images. We present a Bayesian algorithm that estimates what the saturated channel's value would have been in the absence of saturation. The algorithm uses the nonsaturated responses from the other color channels, together with a multivariate normal prior that captures the correlation in response across color channels. The prior may be estimated directly from the image data, since most image pixels are not saturated. Given the prior and the responses of the nonsaturated channels, the algorithm returns the optimal expected mean square estimate for the true response. Extensions of the algorithm to the case in which more than one channel is saturated are also discussed. Both simulations and examples with real images are presented to show that the algorithm is effective.     

Atmospheric correction of ocean-color sensors: effects of the Earth's curvature

We investigate the influence of the curvature of the Earth on a proposed atmospheric-correction scheme for the Sea-Viewing Wide-Field-of-View Sensor (SeaWiFS) by simulating the radiance exiting the top of a spherical-shell atmosphere and inserting the result into the proposed correction algorithm. The error in the derived water-leaving reflectance suggests that the effects of the curvature are negligible for solar zenith angles (θ(0)) ≤ 70°. Furthermore, for θ(0) > 70° the error in atmospheric correction can usually be reduced if the molecular-scattering component of the top of the atmosphere reflectance (ρ(r)) is computed with a spherical-shell atmosphere radiative transfer code. Also, for θ(0) > 70° the error in atmospheric correction in a spherical-shell atmosphere, when ρ(r) is computed with a spherical-shell model, can be predicted reasonably well from computations made with plane-parallel atmosphere radiative transfer codes. This implies that studies aimed at improving atmospheric correction can be made assuming plane-parallel geometry and that the investigator can be confident when θ(0)> 70° that any improvements will still be valid for a spherical-shell atmosphere as long as ρ(r) is computed in spherical-shell geometry. Finally, a scheme for computing ρ(r) in a spherical-shell atmosphere in a relatively simple manner is developed.