Stray Light Correction Research Articles

Low computational demand stray light correction method for hyperspectral imaging spectrometers

Stray light correction in hyperspectral imaging spectrometers has long been restricted by high computational requirements. This paper presents a low computational demand method, based on the matrix operations for spectrometer stray light correction, using an iterative approach to efficiently correct stray light across both spectral and spatial dimensions. The efficacy of this method is demonstrated through its application to simulated and real images, achieving an overall reduction of stray light by over 50%, with significantly reduced computation time and memory usage compared to the method by Zong et al. based on a sparse matrix [Appl. Opt. 45, 1111 (2006)10.1364/AO.45.001111APOPAI2155-3165]. By enabling stray light correction on general-purpose computers, this method enhances affordability and accessibility, promoting broader use and reducing measurement uncertainties in various hyperspectral imaging applications.

Going beyond hardware limitations with advanced stray light calibration for the Metop-3MI space instrument

Space optical instruments play a pivotal role in enhancing our understanding of the universe and our planet, and are crucial in addressing the urgent challenges posed by climate change. In this context, stray light has emerged as a primary performance limitation. Originating from ghost reflections or scattering, it obscures essential details and introduces false information into images. With the demand for increasingly high-performing instruments, mitigation through hardware optimization is becoming insufficient. We are entering an era where future instruments require a stray light correction algorithm to meet user specifications, necessitating extensive on-ground calibration. This paper examines the Metop-3MI Earth observation instrument, which, with wide field of view, broad spectral range, and multi-polarization capabilities, epitomizes the challenges of stray light calibration and correction. A custom calibration apparatus was constructed to evaluate the complex stray light dependence on field-of-view, wavelength, and polarization. Data were processed, and stray light kernels database was derived, which then fed into a specially developed correction algorithm. Applied to the image of an extended scene, it effectively reduces stray light by a remarkable factor of 91. This achievement sets a new standard for low-stray-light instruments and provides a comprehensive case study for future missions.

Spectral properties of bright deposits in permanently shadowed craters on Ceres

Context. Bright deposits in permanently shadowed craters on Ceres are thought to harbor water ice. However, the evidence for water ice presented thus far is indirect. Aims. We aim to directly detect the spectral characteristics of water ice in bright deposits present in permanently shadowed regions (PSRs) in polar craters on Ceres. Methods. We analyzed narrowband images of four of the largest shadowed bright deposits acquired by the Dawn Framing Camera to reconstruct their reflectance spectra, carefully considering issues such as in-field stray light correction and image compression artifacts. Results. The sunlit portion of a polar deposit known to harbor water ice has a negative (blue) spectral slope of −58 ± 12% µm−1 relative to the background in the visible wavelength range. We find that the PSR bright deposits have similarly blue spectral slopes, consistent with a water ice composition. Based on the brightness and spectral properties, we argue that the ice is likely present as particles of high purity. Other components such as phyllosilicates may be mixed in with the ice. Salts are an unlikely brightening agent given their association with cryovolcanic processes, of which we find no trace. Conclusions. Our spectral analysis strengthens the case for the presence of water ice in permanently shadowed craters on Ceres.

Out-of-field stray light correction in optical instruments: the case of Metop-3MI

Out-of-field stray light correction in optical instruments: the case of Metop-3MI

A physically based correction for stray light in Brewer spectrophotometer data analysis

Abstract. Brewer ozone spectrophotometers have become an integral part of the global ground-based ozone monitoring network collecting data since the early 1980s. The double-monochromator Brewer version (MkIII) was introduced in 1992. With the Brewer hardware being so robust, both single- and double-monochromator instruments are still in use. The main difference between the single Brewers and the double Brewers is the much lower stray light in the double instrument. Laser scans estimate the rejection level of the single Brewers to be 10−4.5, while the doubles improve this to 10−8, virtually eliminating the effects of stray light. For a typical single-monochromator Brewer, stray light leads to an underestimation of ozone of approximately 1 % at 1000 DU ozone slant column density (SCD) and can exceed 5 % at 2000 DU, while underestimation of sulfur dioxide reaches 30 DU when no sulfur dioxide is present. This is because even a small additional stray light contribution at shorter wavelengths significantly reduces the calculated SCD at large values. An algorithm for stray light correction based on the physics of the instrument response to stray light (PHYCS) has been developed. The simple assumption is that count rates measured at any wavelength have a contribution from stray light from longer, and thus brighter, wavelengths because of the ozone cross-section gradient leading to a rapid change in intensity as a function of wavelength. Using the longest measured wavelength (320 nm) as a proxy for the overall brightness provides an estimate of this contribution. The sole parameter, on the order of 0.2 % to 0.6 %, that describes the percentage of light at the longest wavelength to be subtracted from all channels is determined by comparing ozone calculations from single- and double-monochromator Brewers making measurements side-by-side. Removing this additional count rate from the signal mathematically before deriving ozone corrects for the extra photons scattering within the instrument that produce the stray light effect. Analyzing historical data from co-located single- and double-monochromator Brewers provides an estimate of how the stray light contribution changes over time in an instrument. The corrected count rates of the measured wavelengths can also be used to improve other calculations: the sulfur dioxide column and the aerosol optical depth, the effective temperature of the ozone layer, or any other products. A multi-platform implementation of PHYCS, rmstray, to correct the count rates for stray light and save the corrected values in a new B-file for use with any existing Brewer data analysis software is available to the global Brewer user community at https://doi.org/10.5281/zenodo.8097038 (Savastiouk and Diémoz, 2023).

Characterization of the East—West Spatial Uniformity for GOES-16/17 ABI Bands Using the Moon

The Advanced Baseline Imager (ABI) is the primary instrument onboard the NOAA Geostationary Operational Environmental Satellite-R Series (GOES-R) satellites, providing continuous weather imagery over the vast area in the Western Hemisphere. It is imperative to ensure consistent calibration accuracy within the instrument’s field of regard (FOR). This paper characterized the spatial uniformity in the east–west (EW) direction for the six ABI visible and near-infrared (VNIR) bands of the first two GOES-R satellites, GOES-16 (G16) and GOES-17 (G17), using a special collection of lunar chasing images during their post-launch testing and post-launch product testing (PLT/PLPT) periods. The EW response versus scan-angle (RVS) is examined with the normalized lunar irradiance ratios at varying scan angles combined from multiple lunar-chasing events. The impacts of straylight from the Earth were found in some of the B01–B03 lunar images. The straylight, including those scattered into the spacelook scenes near the polar regions and those leaked into space near the Moon, can cause RVS variation up to 1% for B01 and to a lesser magnitude for the other two bands. Straylight correction algorithms are applied for the accurate ABI lunar image irradiance calculation. After the corrections, the RVS variation is reduced to less than 0.3% for all the VNIR bands of both G16/17 in full-disk (FD) images. Results of this study also confirm that the Global Space-based Inter-Calibration System (GSICS) Implementation of the ROLO (GIRO) model has high relative accuracy for the ABI VNIR bands when the lunar images are collected within a relatively short time. The method described in this paper can be applied to validate the EW spatial uniformity for imagers on other geostationary satellites, including the recently launched GOES-18 and the future GOES-U satellites.

The Decay Process of an α-configuration Sunspot

The decay of sunspot plays a key role in magnetic flux transportation in solar active regions (ARs). To better understand the physical mechanism of the entire decay process of a sunspot, an α-configuration sunspot in AR NOAA 12411 was studied. Based on the continuum intensity images and vector magnetic field data with stray light correction from Solar Dynamics Observatory/Helioseismic and Magnetic Imager, the area, vector magnetic field and magnetic flux in the umbra and penumbra are calculated with time, respectively. Our main results are as follows: (1) The decay curves of the sunspot area in its umbra, penumbra, and whole sunspot take the appearance of Gaussian profiles. The area decay rates of the umbra, penumbra and whole sunspot are −1.56 MSH day−1, −12.61 MSH day−1 and −14.04 MSH day−1, respectively; (2) With the decay of the sunspot, the total magnetic field strength and the vertical component of the penumbra increase, and the magnetic field of the penumbra becomes more vertical. Meanwhile, the total magnetic field strength and vertical magnetic field strength for the umbra decrease, and the inclination angle changes slightly with an average value of about 20°; (3) The magnetic flux decay curves of the sunspot in its umbra, penumbra, and whole sunspot exhibit quadratic patterns, their magnetic flux decay rates of the umbra, penumbra and whole sunspot are −9.84 × 1019 Mx day−1, −1.59 × 1020 Mx day−1 and −2.60 × 1020 Mx day−1, respectively. The observation suggests that the penumbra may be transformed into the umbra, resulting in the increase of the average vertical magnetic field strength and the reduction of the inclination angle in the penumbra during the decay of the sunspot.

Sensitivity of solar-induced fluorescence to spectral stray light in high resolution imaging spectroscopy

Remote sensing of solar-induced chlorophyll fluorescence (SIF) advances our ability to monitor gross primary productivity. Because SIF is usually <5% of recorded canopy radiance, sources of measurement error must be reduced for accurate SIF retrieval. Here we quantify the impact of spectral stray light on SIF retrieval using a high-resolution imaging spectrometer. We first derived the FLD retrieval solution in the presence of stray light under the simplifying assumptions that SIF, reflectance, and stray light are constant across the region of the spectrum used in the retrieval. We show that stray light contributions from the canopy and reference spectra could cancel out if obtained from the same instrument. In practice however, the simplifying assumptions are unlikely to hold because SIF and reflectance vary across regions of the spectrum commonly exploited by the FLD retrieval, and stray light can also vary across the spectrum, as we demonstrate with an empirical example. To quantify the impact of spectral stray light on SIF retrievals in more realistic scenarios, we performed a sensitivity analysis. We used a stray light signal distribution function and measurements of spectral radiance and SIF to generate spectra with known quantities of SIF and stray light over four orders of magnitude of stray light. We then quantified the bias in retrieved SIF using four SIF retrieval approaches—standard FLD, 3FLD, a spectral fitting method (SFM), and a method based on singular value decomposition (SVD)—applied to the red and/or far-red spectral domains. We found that spectral stray light can cause either a positive or negative bias in SIF estimates. In the highest stray light scenario, with spectral stray light one order of magnitude smaller than the radiance signal, the stray light error ranged from < ± 1% of SIF (for O2-A band retrievals with FLD, 3FLD, and SFM methods, and for far-red retrievals in the 748 nm to 756 nm domain with SVD-based and 3FLD methods applied to Fraunhofer lines) to >30% of SIF (O2-B band retrieval with SFM, and 3FLD applied to Fraunhofer lines in the red domain). For context, a bias of 27% of SIF is a comparable magnitude to seasonal variation in SIF observed in the Amazon basin, and to SIF differences among plant functional types. However, when spectral stray light was three orders of magnitude smaller than the radiance signal, mean bias was < ± 2% for all retrieval methods. Our analysis informs when spectral stray light corrections must be performed for SIF remote sensing.

Calibration of the NOMAD-UVIS data

The Ultraviolet and VIsible Spectrometer (UVIS), covering the 200–650 ​nm range, is one of three spectrometers that comprise the NOMAD instrument on the ExoMars 2016 Trace Gas Orbiter (TGO). UVIS can operate in solar occultation, nadir and limb viewing mode and was designed to monitor ozone and aerosols in the Martian atmosphere. Here, we describe the calibration procedure to convert the UVIS raw data into a calibrated data product ready for scientific exploitation. The calibration includes the CCD offset and dark current subtraction, the wavelength assignment, the noise identification and removal, the smearing removal, and the radiance or transmittance conversion. A straylight correction, critical for some parts of the UVIS spectral range, is also applied during the data reduction process, which is described in more detail in two companion papers [Mason et al., 2022; Depiesse et al., In prep] corresponding to two different and independent methods giving consistent results. The solar occultation observations are converted into transmittance and are therefore self-calibrating, while nadir and limb measurements require an absolute radiometric calibration. A comparison with coincident nadir MRO/MARCI measurements is provided as a final validation and generally shows a ±10% agreement on the radiances measured by both instruments. • Description of the calibration of UVIS measurements to produce level 1 data from raw data • CCD non-linearity correction, offset removal, dark current removal, Bad pixel removal, Smearing removal. • A straylight correction is also applied but described in companion papers. • Conversions: Pixel-wavelength assignment, Radiance (nadir) or Transmittance (occultation) conversion. • Radiometric comparison with coincident MRO/MARCI measurements, serving as validation for the radiance calibrated data.

Spectrophotometry of the lunar regolith using the Chang’E-3 Panoramic Camera (PCAM)

Aims. We present a multiband photometric study of the lunar mare surface based on in situ measurements using the Panoramic Camera on board the Chang'E-3 mission. We aim to obtain photometric properties and to investigate the differences in measuring sites in various conditions. Methods. Data were collected with a vast range of phase angles from 0.5° to 160°, probing 0–5° phases to constrain the opposition effect (OE) and phases above 120° to constrain the phase function. Stray light and spectral corrections were conducted to calibrate the regolith’s bidirectional reflectance distribution function. We employed the Hapke photometric model to describe the regolith’s photometric properties. Results. Phase functions and Hapke parameters at three measuring sites were retrieved. The mare regolith at the measuring sites had a single-scattering albedo (SSA) of 0.17-0.21 for the green channel and exhibited backscattering and a strong OE. The coherent backscatter opposition effect (CBOE) represented the majority of OE, with a proportion of around 80% for pristine regolith. The phase ratio curves were arch-shaped, and the arch positions of the curves were associated with grain size.

New Reprocessing towards Life-Time Quality-Consistent Suomi NPP OMPS Nadir Sensor Data Records (SDR): Calibration Improvements and Impact Assessments on Long-Term Quality Stability of OMPS SDR Data Sets

The Nadir Mapper (NM) and Nadir Profiler (NP) within the Ozone Mapping and Profiler Suites (OMPS) are ultraviolet spectrometers to measure Earth radiance and Solar irradiance spectra from 300–380 nm and 250–310 nm, respectively. The OMPS NM and NP instruments flying on the Suomi-NPP (SNPP) satellite have provided over ten years of operational Sensor Data Records (SDRs) data sets to support a variety of OMPS Environmental Data Record (EDR) applications. However, the discrepancies of quality remain in the operational OMPS SDR data prior to 28 June 2021 due to changes in calibration algorithms associated with the calibration coefficient look-up tables (LUTs) during this period. In this study, we present results for the newly (v2) reprocessed SNPP OMPS NM and NP SDR data prior to 30 June 2021, which uses consistent calibration tables with improved accuracy. Compared with a previous (v1) reprocessing, this new reprocessing includes the improvements associated with the following updated tables or error correction: an updated stray light correction table for the NM, an off-nadir geolocation error correction for the NM, an artificial offset error correction in the NM dark processing code, and biweekly solar wavelength LUTs for the NP. This study further analyzes the impact of each improvement on the quality of the OMPS SDR data by taking advantage of the existing OMPS SDR calibration/validation studies. Finally, this study compares the v2 reprocessed OMPS data sets with the operational and the v1 reprocessed data sets. The results demonstrate that the new reprocessing significantly improves the accuracy and consistency of the life-time SNPP OMPS NM and NP SDR data sets. It also advances the uniformity of the data over the dichroic range from 300 to 310 nm between the NM and NP. The normalized radiance differences at the same wavelength between the NM and NP observations are reduced from 0.001 order (v1 reprocessing) or 0.01 order (operational processing) to 0.001 order or smaller. The v2 reprocessed data are archived in the NOAA CLASS data center with the same format as the operational data.

Stray Light Correction Algorithm for High Performance Optical Instruments: The Case of Metop-3MI

Stray light is a critical aspect for high performance optical instruments. When stray light control by design is insufficient to reach the performance requirement, correction by post-processing must be considered. This situation is encountered, for example, in the case of the Earth observation instrument 3MI, whose stray light properties are complex due to the presence of many ghosts distributed on the detector array. We implement an iterative correction method and discuss its convergence properties. Spatial and field binning can be employed to reduce the computation time but at the cost of a decreased performance. Interpolation of the stray light properties is required to achieve high performance correction. For that, two methods are proposed and tested. The first interpolate the stray light in the field domain while the second applies a scaling operation based on a local symmetry assumption. Ultimately, the scaling method is selected and a stray light reduction by a factor of 58 is obtained at 2σ (129 at 1σ) for an extended scene illumination.

A simple method of UV stray light correction for field spectrometers in Ground Validation Sites

Field spectrometers are widely applied in ground validation sites for remote sensing. However, the control of stray light is one of the most important factors to ensure accurate measurements, since the spectral distribution of target source differs significantly from the laboratory calibration source. Here, a simple UV stray light correction method for continuously distributed wide-spectrum light sources was established by using a set of bandpass filters. This analysis enabled a simple, low-cost method for correcting for the observed straylight.

Stray-Light Correction Methodology for the Precision Solar Spectroradiometer

A stray-light correction methodology for the Precision Solar Spectroradiometer (PSR) is presented. The correction is based on laboratory-measured line spread functions also taking into account the radiation from the 2nd and 3rd grating orders. The efficiency of the correction is validated on solar and lamp measurement data. The results are compared to those obtained with a PSR equipped with an order-sorting filter and with a Precision Filter Radiometer.

SNPP VIIRS Day Night Band: Ten Years of On-Orbit Calibration and Performance

Aboard the polar-orbiting SNPP satellite, the VIIRS instrument has been in operation since launch in October 2011. It is a visible and infrared radiometer with a unique panchromatic channel capability designated as a day-night band (DNB). This channel covers wavelengths from 0.5 to 0.9 µm and is designed with a near-constant spatial resolution for Earth observations 24 h a day. The DNB operates at 3 gain stages (low, middle, and high) to cover a large dynamic range. An onboard solar diffuser (SD) is used for calibration in the low gain stage, and to enable the derivation of gain ratios between the different stages. In this paper, we present the SNPP VIIRS DNB calibration performed by the NASA VIIRS characterization support team (VCST). The DNB calibration algorithms are described to generate the calibration coefficient look up tables (LUTs) for the latest NASA Level 1B Collection 2 products. We provide an evaluation of DNB on-orbit calibration performance. This activity supports the NASA Earth science community by delivering consistent VIIRS sensor data products via the Land Science Investigator-led Processing Systems, including the SD degradation applied for DNB calibrations in detector gain and gain ratio trending. The DNB stray light contamination and its correction are highlighted. Performance validations are presented using comparisons to the calibration methods employed by NOAA’s operational Interface Data Processing Segment. Further work on stray light corrections is also discussed.

Correction of atmospheric stray light in restored slit spectra

Context. A long-standing issue in solar ground-based observations has been the contamination of data due to stray light, which is particularly relevant in inversions of spectropolarimetric data. Aims. We aim to build on a statistical method of correcting stray-light contamination due to residual high-order aberrations and apply it to ground-based slit spectra. Methods. The observations were obtained at the Swedish Solar Telescope, and restored using the multi-frame blind deconvolution restoration procedure. Using the statistical properties of seeing, we created artificially degraded synthetic images generated from magneto-hydrodynamic simulations. We then compared the synthetic data with the observations to derive estimates of the amount of the residual stray light in the observations. In the final step, the slit spectra were deconvolved with a stray-light point spread function to remove the residual stray light from the observations. Results. The RMS granulation contrasts of the deconvolved spectra were found to increase to approximately 12.5%, from 9%. Spectral lines, on average, were found to become deeper in the granules and shallower in the inter-granular lanes, indicating systematic changes to gradients in temperature. The deconvolution was also found to increase the redshifts and blueshifts of spectral lines, suggesting that the velocities of granulation in the solar photosphere are higher than had previously been observed.

Raw EPIC Data Calibration

Earth Polychromatic Imaging Camera (EPIC) raw level-0 (L0) data in one channel is a 12-bit 2,048 × 2,048 pixels image array plus auxiliary data such as telemetry, temperature, etc. The EPIC L1a processor applies a series of correction steps on the L0 data to convert them into corrected count rates (level-1a or L1a data): Dark correction, Enhanced pixel detection, Read wave correction, Latency correction, Non-linearity correction, Temperature correction, Conversion to count rates, Flat fielding, and Stray light correction. L1a images should have all instrumental effects removed and only need to be multiplied by one single number for each wavelength to convert counts to radiances, which are the basis for all higher-level EPIC products, such as ozone and sulfur dioxide total column amounts, vegetation index, cloud, aerosol, ocean surface, and vegetation properties, etc. This paper gives an overview of the mathematics and the pre-launch and on-orbit calibration behind each correction step.

Mission-Long Recalibrated Science Quality Suomi NPP VIIRS Radiometric Dataset Using Advanced Algorithms for Time Series Studies

Suomi NPP has been successfully operating since its launch on 28 October 2011. As one of the major payloads, along with microwave and infrared sounders (Advanced Technology Microwave Sounder (ATMS), Cross-track Infrared Sounder (CrIS)), and ozone mapping/profiling (OMPS) instruments, the Visible Infrared Imaging Radiometer Suite (VIIRS) has performed for well beyond its mission design life. Its data have been used for a variety of applications for nearly 30 environmental data products, including global imagery twice daily with 375 and 750 m resolutions, clouds, aerosol, cryosphere, ocean color and sea-surface temperature, a number of land products (vegetation, land-cover, fire and others), and geophysical and social economic studies with nightlights. During the early days of VIIRS operational calibration and data production, there were inconsistencies in both algorithms and calibration inputs, for several reasons. While these inconsistencies have less impact on nowcasting and near real-time applications, they introduce challenges for time series analysis due to calibration artifacts. To address this issue, we developed a comprehensive algorithm, and recalibrated and reprocessed the Suomi NPP VIIRS radiometric data that have been produced since the launch. In the recalibration, we resolved inconsistencies in the processing algorithms, terrain correction, straylight correction, and anomalies in the thermal bands. To improve the stability of the reflective solar bands, we developed a Kalman filtering model to incorporate onboard solar, lunar, desert site, inter-satellite calibration, and a deep convective cloud calibration methodology. We further developed and implemented the Solar Diffuser Surface Roughness Rayleigh Scattering model to account for the sensor responsivity degradation in the near infrared bands. The recalibrated dataset was validated using vicarious sites and alternative methods, and compared with independent processing from other organizations. The recalibrated radiometric dataset (namely, the level 1b or sensor data records) also incorporates a bias correction for the reflective solar bands, which not only addresses known calibration biases, but also allows alternative calibrations to be applied if so desired. The recalibrated data have been proven to be of high quality, with much improved stability (better than 0.3%) and accuracy (by up to 2%). The recalibrated radiance data are now available from 2012 to 2020 for users and will eventually be archived on the NOAA CLASS database.

NOAA-20 Visible Infrared Imaging Radiometer Suite day–night band on-orbit calibration and performance

The NOAA-20 Visible Infrared Imaging Radiometer Suite (VIIRS) instrument has been successfully operating on orbit since November 28, 2017. The day–night band (DNB) onboard NOAA-20 VIIRS is a panchromatic channel covering wavelengths from 0.5 to 0.9 μm, capable of observing the Earth scene in visible/near-infrared spectral range at a spatial resolution of 750 m. The DNB operates at low-, mid-, or high-radiometric gain stages, and it uses an onboard solar diffuser (SD) panel for low-gain stage calibration. The SD observations also provide a means to compute gain ratios between low-to-mid and mid-to-high-gain stages. With their large dynamic range and high sensitivity, the DNB detectors can make observations during both daytime and nighttime. We provide an assessment of the DNB on-orbit performance and behavior in the first two-year mission period and beyond. The calibration methodology used by the VIIRS Characterization Support Team in support of the NASA Earth science community has been described. The trending of on-board calibrators dark-offsets, SD gains, and gain ratios, and signal-to-noise ratio at minimum radiance have been analyzed, especially during key events such as the nadir and cryocooler doors opening. Furthermore, we performed intercomparison studies between Suomi National Polar-orbiting Partnership and NOAA-20 instruments and evaluated DNB radiometric calibration and characterization, including the SD degradation, detector gains, gain ratios, and straylight correction, as well as the calibration comparison between the NOAA Interface Data Processing Segment look-up-tables and our delivery results.

Use of Moon Observations for Characterization of Sentinel-3B Ocean and Land Color Instrument

During the commissioning of the Sentinel-3B satellite, a single lunar observation was performed to assess the possible use of the moon for characterization and validation of onboard instruments. The observation was carried out in stable orientation after a roll maneuver, allowing the moon to be imaged by the Earth view of instruments. Data acquired by the Ocean Land Color Instrument (OLCI) allowed inflight verification of stray-light correction (SLC) performed by the Mission Performance Centre (MPC), and assessment of radiometric behavior of instrument in comparison with lunar irradiance models performed in cooperation between European Space Research and Technology Centre (ESTEC) and European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT). This paper describes the results of those activities along with the proposed update of stray-light correction developed with the use of lunar data.