SEARCH · Search NASA

Results for “retrieval algorithm”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Mechanisms for Diurnal Variability of Global Tropical Rainfall Observed from TRMM

The behavior and various controls of diurnal variability in tropical-subtropical rainfall are investigated using Tropical Rainfall Measuring Mission (TRMM) precipitation measurements retrieved from the three level-2 TRMM standard profile algorithms for the 1998 annual cycle. Results show that diurnal variability characteristics of precipitation are consistent for all three algorithms, providing assurance that TRMM retrievals are producing consistent estimates of rainfall variability. As anticipated, most ocean areas exhibit more rainfall at night, while over most land areas, rainfall peaks during daytime; however, important exceptions are noted. The dominant feature of the oceanic diurnal cycle is a rainfall maximum in late-evening-early-morning (LE-EM) hours, while over land the dominant maximum occurs in the mid- to late afternoon (MLA). In conjunction with these maxima are pronounced seasonal variations of the diurnal amplitudes. Amplitude analysis shows that the diurnal pattern and its seasonal evolution are closely related to the rainfall accumulation pattern and its seasonal evolution. In addition, the horizontal distribution of diurnal variability indicates that for oceanic rainfall, there is a secondary MLA maximum coexisting with the LE-EM maximum at latitudes dominated by large-scale convergence and deep convection. Analogously, there is a preponderancy for an LE EM maximum over land coexisting with the stronger MLA maximum, although it is not evident that this secondary continental feature is closely associated with the large-scale circulation. Neither of the secondary maxima exhibit phase behavior that can be considered semidiurnal in nature. Diurnal rainfall variability over the ocean associated with large-scale convection is clearly an integral component of the general circulation. Phase analysis reveals differences in regional and seasonal features of the diurnal cycle, indicating that underlying forcing mechanisms differ from place to place. This is underscored by the appearance of secondary ocean maxima in the presence of large-scale convection, along with other important features. Among these, there are clear-cut differences between the diurnal variability of seasonal rainfall over the mid-Pacific and Indian Ocean Basins. The mid-Pacific exhibits double maxima in spring and winter but only LE-EM maxima in summer and autumn, while the Indian Ocean exhibits double maxima in spring and summer and only an LE-EM maximum in autumn and winter. There are also evident daytime maxima within the major large-scale marine stratocumulus regions off the west coasts of continents. The study concludes with a discussion concerning how the observational evidence either supports or repudiates possible forcing mechanisms that have been suggested to explain diurnal rainfall variability.

Yang, Song↗

DOE OSTI2018

The Complexity of Bit Retrieval

Bit retrieval is the problem of reconstructing a periodic binary sequence from its periodic autocorrelation, with applications in cryptography and x-ray crystallography. After defining the problem, with and without noise, we describe and compare various algorithms for solving it. A geometrical constraint satisfaction algorithm, relaxed-reflect-reflect, is currently the best algorithm for noisy bit retrieval.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2016

TRMM Data Improvement as Part of the GPM Data Processing

NASA has a long standing commitment to the improvement of its mission datasets. Indeed, data reprocessing is always built into the plans, schedule and budget for the mission data processing system. However, in addition to these ongoing mission reprocessing, NASA also supports a final reprocessing of all the data for a mission upon its completion (known as Phase F). TRMM Phase F started with the end of the TRMM mission in June of 2015. This last reprocessing has two overall goals: improvement of the TRMM mission data products; incorporation of the 17+ years of TRMM data into the ongoing NASA/JAXA GPM data processing. The first goal guarantees that the latest algorithms used for precipitation retrievals will also be used in reprocessing the TRMM data. The second goal ensures that as GPM algorithms are improved, the entire TRMM data will always be reprocessed with each GPM reprocessing. In essence TRMM becomes another of the GPM constellation satellites. This paper will concentrate on presenting the improvements to TMI level 1 data including calibration, geolocation, and emissive antenna corrections. It will describe the format changes that will occur how the TMI level 1C product will be intercalibrated using GMI as the reference calibration. It will also provide an overview of changes in the precipitation radar products as well as the combined TMIPR product.

GPM↗

NASA NTRS2004

Hexagonal Pixels and Indexing Scheme for Binary Images

A scheme for resampling binaryimage data from a rectangular grid to a regular hexagonal grid and an associated tree-structured pixel-indexing scheme keyed to the level of resolution have been devised. This scheme could be utilized in conjunction with appropriate image-data-processing algorithms to enable automated retrieval and/or recognition of images. For some purposes, this scheme is superior to a prior scheme that relies on rectangular pixels: one example of such a purpose is recognition of fingerprints, which can be approximated more closely by use of line segments along hexagonal axes than by line segments along rectangular axes. This scheme could also be combined with algorithms for query-image-based retrieval of images via the Internet. A binary image on a rectangular grid is generated by raster scanning or by sampling on a stationary grid of rectangular pixels. In either case, each pixel (each cell in the rectangular grid) is denoted as either bright or dark, depending on whether the light level in the pixel is above or below a prescribed threshold. The binary data on such an image are stored in a matrix form that lends itself readily to searches of line segments aligned with either or both of the perpendicular coordinate axes. The first step in resampling onto a regular hexagonal grid is to make the resolution of the hexagonal grid fine enough to capture all the binaryimage detail from the rectangular grid. In practice, this amounts to choosing a hexagonal-cell width equal to or less than a third of the rectangular- cell width. Once the data have been resampled onto the hexagonal grid, the image can readily be checked for line segments aligned with the hexagonal coordinate axes, which typically lie at angles of 30deg, 90deg, and 150deg with respect to say, the horizontal rectangular coordinate axis. Optionally, one can then rotate the rectangular image by 90deg, then again sample onto the hexagonal grid and check for line segments at angles of 0deg, 60deg, and 120deg to the original horizontal coordinate axis. The net result is that one has checked for line segments at angular intervals of 30deg. For even finer angular resolution, one could, for example, then rotate the rectangular-grid image +/-45deg before sampling to perform checking for line segments at angular intervals of 15deg.

Johnson, Gordon G.↗

DOE OSTI2021

Imaging pyrometry for most color cameras using a triple pass filter

A simple combination of the Planck blackbody emission law, optical filters, and digital image processing is demonstrated to enable most commercial color cameras (still and video) to be used as an imaging pyrometer for flames and explosions. The hardware and data processing described take advantage of the color filter array (CFA) that is deposited on the surface of the light sensor array present in most digital color cameras. In this work, a triple-pass optical filter incorporated into the camera lens allows light in three 10-nm wide bandpass regions to reach the CFA/light sensor array. These bandpass regions are centered over the maxima in the blue, green, and red transmission regions of the CFA, minimizing the spectral overlap of these regions normally present. A computer algorithm is used to retrieve the blue, green, and red image matrices from camera memory and correct for remaining spectral overlap. A second algorithm calibrates the corrected intensities to a gray body emitter of known temperature, producing a color intensity correction factor for the camera/filter system. The Wien approximation to the Planck blackbody emission law is used to construct temperature images from the three color (blue, green, red) matrices. A short pass filter set eliminates light of wavelengths longer than 750 nm, providing reasonable accuracy (±10%) for temperatures between 1200 and 6000 K. The effectiveness of this system is demonstrated by measuring the temperature of several systems for which the temperature is known.

47 OTHER INSTRUMENTATION↗

NASA NTRS2023

Integrating NASA Aqua AIRS in a Real-Time NUCAPS Science-to-Applications System to Support Severe Weather Forecasting

In recent years, National Oceanic and Atmospheric Administration (NOAA) Unique Combined Atmospheric Processing System (NUCAPS) hyperspectral infrared satellite sounding retrievals derived from Joint Polar Satellite System (JPSS) polar-orbiting satellites have been documented as observations that add value to weather forecasting applications. NUCAPS is currently the operational algorithm delivering JPSS satellite sounding retrievals to the NOAA National Weather Service (NWS) and is based on the heritage Atmospheric Infrared Sounder (AIRS) Science Team algorithm for processing vertical temperature, moisture, and trace gas retrievals. For the Special Collection on “Twenty Years of Observations from AIRS,” we highlight the methodology we implemented to develop a prototype science-to-applications system to enable real-time processing of AIRS satellite sounding retrievals through the NUCAPS algorithm (i.e., NUCAPS-Aqua) to support weather forecasting applications. The addition of NUCAPS-Aqua to experimental real-time pathways alongside operational JPSS NUCAPS soundings, facilitated assessment of NUCAPS-Aqua at the 2022 Hazardous Weather Testbed (HWT) Spring Experiment. Development of NUCAPS-Aqua described in this technical report includes preservation of microwave observations and calculation of a-priori regression coefficients. Additionally, the real-time processing and challenges with implementing a science-to-applications system are discussed. Two illustrative pre-convective forecasting examples analyzed by NWS forecasters during the 2022 HWT Spring Experiment are highlighted to demonstrate the benefit of NUCAPS-Aqua as (a) special afternoon soundings, (b) an additional observation to assess temporal trends using multiple satellites, and (c) a complement to observational and model analysis.

remote sensing↗

NASA NTRSDate not supplied

Development of a Global Reference Surface Reflectance and BRDF Datasets from Geostationary Satellite Observations and AERONET Measurements

Surface reflectances and their dependency on illumination-view geometries (i.e., BRDF) are the foundation of many high-level satellite products for land and water monitoring. Yet it is difficult to evaluate the quality of satellite-based surface reflectances with ground-based measurements due to the spatial scale differences. In order to fill the gap, here we develop a reference dataset of surface reflectance and BRDF at the global AERONET sites with data streams from operational geostationary sensors including Himawari 8/9 AHI, GK-2A AMI, and GOES 16/17 ABI. Taking the top-of-atmosphere (TOA) reflectance and the site measured atmospheric aerosol optical depth (AOD) as the main inputs, we apply the GeoNEX-AC algorithm to performance accurate atmospheric correction and derive 10-minute surface reflectance and daily Ross-Thick-Li-Sparse (RTLS) BRDF parameters at AERONET sites where coincident AOD measurements and TOA observations are available from 2016 (for Himawari) or 2018 (for GOES) onwards. The algorithm ensures that the retrieved surface BRDF parameters, along with the site-measured AOD, allow the atmospheric radiative transfer model, SHARM, accurately simulate the observed TOA reflectance at diurnal and longer time scales. They are our best estimates of the surface optical properties and thus can serve as the “reference” to evaluate the performance of operational atmospheric correction algorithms (where AOD is assumed unknown and needs to be retrieved). The reference BRDF also allow us to evaluate the spectral band ratios between the SWIR (e.g., 2200 nm) and the visible (e.g., 650 nm) regions, which are commonly used in operational atmospheric correction algorithms. Finally, we demonstrate that the reference dataset can be used to develop potential data synergies between different GEO satellites as well as GEO-LEO sensors.

Weile Wang↗

NASA NTRS2005

Converting Soil Moisture Observations to Effective Values for Improved Validation of Remotely Sensed Soil Moisture

We compare soil moisture retrieved with an inverse algorithm with observations of mean moisture in the 0-6 cm soil layer. A significant discrepancy is noted between the retrieved and observed moisture. Using emitting depth functions as weighting functions to convert the observed mean moisture to observed effective moisture removes nearly one-half of the discrepancy noted. This result has important implications in remote sensing validation studies.

Laymon, Charles A.↗

NASA NTRSDate not supplied

NASA GPM Status and Future Activities

The joint U.S.-Japan Global Precipitation Measurement (GPM) mission is approaching a decade of operations, and continues to pursue research, dataset production, and outreach related to precipitation. One key activity over the last year was the release of an improved “Version 07” of all GPM precipitation and latent heating products. This talk summarizes key improvements to the GPM products for which NASA has lead responsibility and provides some examples of the changes between Versions 06 and 07 in algorithm performance. One important operational change that affected Version 07 is that the scanning strategy for the Ka-band radar channel changed in May 2018; all products that depend on Ka were revised to accommodate this change. For example, in Version 07 the Goddard Profiling (GPROF) algorithm has implemented improvements in regions where orographic enhancement and suppression take place and where the surface is snowy/icy, and again covers radiometers reaching back to 1987. The Combined Radar Radiometer Algorithm (CORRA) now incorporates modified drop-size distribution constraints that substantially reduce bias. Revisions to the Convective-Stratiform Heating (CSH) algorithm employ new radiative transfer retrievals as well as accounting for terrain in the vertical coordinates. Each algorithm was adjusted to ensure continuity for each product across the boundary in 2014 between the predecessor Tropical Rainfall Measuring Mission (TRMM) and the GPM Core Observatory. The U.S. Science Team’s Integrated Multi-satellitE Retrievals for GPM (IMERG) was upgraded to account for distortions in the probability density function of regional precipitation rates due to weighted averaging in the Kalman filter used for “morphing” the passive microwave data. The talk will conclude by considering major issues that require continued attention, including the use of machine learning algorithms, the operational challenge of swarms of “small”, perhaps short-lived satellites, and estimates of the remaining lifespan of the Core Observatory.

Global Precipitation Measuremen↗

NASA NTRS2022

Global GNSS-RO Electron Density in the Lower Ionosphere

Lack of instrument sensitivity to low electron density (Ne) concentration makes it difficult to measure sharp Ne vertical gradients (four orders of magnitude over 30 km) in the D/E-region. A robust algorithm is developed to retrieve global D/E-region Ne from the high-rate GNSS radio occultation (RO) data, to improve spatio-temporal coverage using recent SmallSat/CubeSat constellations.The new algorithm removes F-region contributions in the RO excess phase profile by fitting a linear function to the data below the D-region.The new GNSS-RO observations reveal many interesting features in the diurnal, seasonal, solar-cycle, and magnetic-field-dependent variations in the Ne morphology. While the D/E-region Ne is a function of solar zenith angle (c), it exhibits strong latitudinal variations for the same c with a distribution asymmetric about noon. In addition, large longitudinal variations are observed along the same magnetic field pitch angle. The summer-midlatitude Ne and sporadic-E(Es) show a distribution similar to each other. The distribution of auroral electron precipitation correlates better with the pitch angle from the magnetosphere than from one at100 km. Finally, anew TEC retrieval technique is developed for the high-rate RO data with a top reaching at least 120 km. For better characterization of the E-to-F transition in Ne and more accurate TEC retrievals, it is recommended to have all GNSS-RO acquisition routinely up to 220 km.

GNSS radio occultation↗

NASA NTRS2014

Simultaneous Retrieval of Aerosol Properties and Clear-Sky Direct Radiative Effect Over the Global Ocean from MODIS

A unified satellite algorithm is presented to simultaneously retrieve aerosol properties (aerosol optical depth; AOD and aerosol type) and clear-sky shortwave direct radiative effect (hereafter, DREA) over ocean. The algorithm is applied to Moderate Resolution Imaging spectroradiometer (MODIS) observations for a period from 2003 to 2010 to assess the DREA over the global ocean. The simultaneous retrieval utilizes lookup table (LUT) containing both spectral reflectances and solar irradiances calculated using a single radiative transfer model with the same aerosol input data. This study finds that aerosols cool the top-of-atmosphere (TOA) and bottom-of-atmosphere (BOA) by 5.2 +/- 0.5 W/sq m and 8.3 W/sq m, respectively, and correspondingly warm the atmosphere (hereafter, ATM) by 3.1 W/sq m. These quantities, solely based on the MODIS observations, are consistent with those of previous studies incorporating chemical transport model simulations and satellite observations. However, the DREAs at BOA and ATM are expected to be less accurate compared to that of TOA due to low sensitivity in retrieving aerosol type information, which is related with the atmospheric heating by aerosols, particularly in low AOD conditions; consequently, the uncertainties could not be quantified. Despite the issue in the aerosol type information, the present method allows us to confine the DREA attributed only to fine-mode dominant aerosols, which are expected to be mostly anthropogenic origin, in the range from 1.1 W/sq m to 1.3 W/sq m at TOA. Improvements in size-resolved AOD and SSA retrievals from current and upcoming satellite instruments are suggested to better assess the DREA, particularly at BOA and ATM, where aerosol absorptivity induces substantial uncertainty.

MODIS↗

NASA NTRS2014

Retrieval of Aerosol Optical Properties under Thin Cirrus from MODIS

Retrieval of aerosol optical properties using shortwave bands from passive satellite sensors, such as MODIS, is typically limited to cloud-free areas. However, if the clouds are thin enough (i.e. thin cirrus) such that the satellite-observed reflectance contains signals under the cirrus layer, and if the optical properties of this cirrus layer are known, the TOA reflectance can be corrected for the cirrus layer to be used for retrieving aerosol optical properties. To this end, we first correct the TOA reflectances in the aerosol bands (0.47, 0.55, 0.65, 0.86, 1.24, 1.63, and 2.12 micron for ocean algorithm and 0.412, 0.47, and 0.65 micron for deep blue algorithm) for the effects of thin cirrus using 1.38 micron reflectance and conversion factors that convert cirrus reflectance in 1.38 micron band to those in aerosol bands. It was found that the conversion factors can be calculated by using relationships between reflectances in 1.38 micron band and minimum reflectances in the aerosol bands (Gao et al., 2002). Refer to the example in the figure. Then, the cirrus-corrected reflectance can be calculated by subtracting the cirrus reflectance from the TOA reflectance in the optically thin case. A sensitivity study suggested that cloudy-sky TOA reflectances can be calculated with small errors in the form of simple linear addition of cirrus-only reflectances and clear-sky reflectances. In this study, we correct the cirrus signals up to TOA reflectance at 1.38 micron of 0.05 where the simple linear addition is valid without extensive radiative transfer simulations. When each scene passes the set of tests shown in the flowchart, the scene is corrected for cirrus contamination and passed into aerosol retrieval algorithms.

MODIS↗

NASA NTRS1996

Ku-Band Ocean Backscatter Functions for Surface Wind Retrieval

For remote sensing of ocean surface wind, Ku-band scatterometer technology has been demonstrated and will be used to monitor global ocean wind fields from spaceborne sensors such as the NASA Scatterometer (NSCAT). Algorithms for wind velocity retrieval from scatterometer data utilize absolute radar returns and relative azimuth modulations of backscatter. The NUSCAT-SWADE data base is used.

ocean↗

NASA NTRS2012

Estimating Sea Surface Salinity and Wind Using Combined Passive and Active L-Band Microwave Observations

Several L-band microwave radiometer and radar missions have been, or will be, operating in space for land and ocean observations. These include the NASA Aquarius mission and the Soil Moisture Active Passive (SMAP) mission, both of which use combined passive/ active L-band instruments. Aquarius s passive/active L-band microwave sensor has been designed to map the salinity field at the surface of the ocean from space. SMAP s primary objectives are for soil moisture and freeze/thaw detection, but it will operate continuously over the ocean, and hence will have significant potential for ocean surface research. In this innovation, an algorithm has been developed to retrieve simultaneously ocean surface salinity and wind from combined passive/active L-band microwave observations of sea surfaces. The algorithm takes advantage of the differing response of brightness temperatures and radar backscatter to salinity, wind speed, and direction, thus minimizing the least squares error (LSE) measure, which signifies the difference between measurements and model functions of brightness temperatures and radar backscatter. The algorithm uses the conjugate gradient method to search for the local minima of the LSE. Three LSE measures with different measurement combinations have been tested. The first LSE measure uses passive microwave data only with retrieval errors reaching 1 to 2 psu (practical salinity units) for salinity, and 1 to 2 m/s for wind speed. The second LSE measure uses both passive and active microwave data for vertical and horizontal polarizations. The addition of active microwave data significantly improves the retrieval accuracy by about a factor of five. To mitigate the impact of Faraday rotation on satellite observations, the third LSE measure uses measurement combinations invariant under the Faraday rotation. For Aquarius, the expected RMS SSS (sea surface salinity) error will be less than about 0.2 psu for low winds, and increases to 0.3 psu at 25 m/s wind speed for warm waters (25 C). To achieve the required 0.2 psu accuracy, the impact of sea surface roughness (e.g. wind-generated ripples) on the observed brightness temperature has to be corrected to better than one tenth of a degree Kelvin. With this algorithm, the accuracy of retrieved wind speed will be high, varying from a few tenths to 0.6 m/s. The expected direction accuracy is also excellent (less than 10 ) for mid to high winds, but degrades for lower speeds (less than 7 m/s).

Yueh, Simon H.↗

NASA NTRS2000

Retrieved Vertical Profiles of Latent Heat Release Using TRMM Rainfall Products

This paper represents the first attempt to use TRMM rainfall information to estimate the four dimensional latent heating structure over the global tropics for February 1998. The mean latent heating profiles over six oceanic regions (TOGA COARE IFA, Central Pacific, S. Pacific Convergence Zone, East Pacific, Indian Ocean and Atlantic Ocean) and three continental regions (S. America, Central Africa and Australia) are estimated and studied. The heating profiles obtained from the results of diagnostic budget studies over a broad range of geographic locations are used to provide comparisons and indirect validation for the heating algorithm estimated heating profiles. Three different latent heating algorithms, the Goddard Convective-Stratiform (CSH) heating, the Goddard Profiling (GPROF) heating, and the Hydrometeor heating (HH) are used and their results are intercompared. The horizontal distribution or patterns of latent heat release from the three different heating retrieval methods are quite similar. They all can identify the areas of major convective activity (i.e., a well defined ITCZ in the Pacific, a distinct SPCZ) in the global tropics. The magnitude of their estimated latent heating release is also not in bad agreement with each other and with those determined from diagnostic budget studies. However, the major difference among these three heating retrieval algorithms is the altitude of the maximum heating level. The CSH algorithm estimated heating profiles only show one maximum heating level, and the level varies between convective activity from various geographic locations. These features are in good agreement with diagnostic budget studies. By contrast, two maximum heating levels were found using the GPROF heating and HH algorithms. The latent heating profiles estimated from all three methods can not show cooling between active convective events. We also examined the impact of different TMI (Multi-channel Passive Microwave Sensor) and PR (Precipitation Radar) rainfall information on latent heating structures.

Tao, W.-K.↗

NASA NTRS2019

An Approach to Retrieve BRDF from Satellite and Airborne Measurements of Surface-Reflected Radiance Based on Decoupling of Atmospheric Radiative Transfer and Surface Reflection

Bi-directional Reflection Distribution Function (BRDF) defines anisotropy of the surface reflection. It is required to specify the boundary condition for radiative transfer (RT) modeling. Measurements of reflected radiance by satellite- and air-borne sensors provide information about anisotropy of surface reflection. Atmospheric correction needs to be performed to derive BRDF from the reflected radiance. Common approach for BRDF retrievals consists of the use of kernel-based BRDF and RT modeling that needs to be done anew at every step of the iterative process. The kernels’ weights are obtained by minimization of the difference between measured and modeled radiance. This study develops a new method of retrieving kernel-based BRDF that requires RT calculations to be done only once. The method employs the exact analytical expression of radiance at any atmospheric level through the solutions of two auxiliary atmosphere-only RT problems and the surface-reflected radiance at the surface level. The latter is related to BRDF and solutions of the auxiliary RT problems by a Fredholm integral equation of the second kind. The approach requires to perform RT calculations one time before the iterations. It can use observations taken at different atmospheric conditions assuming that surface conditions remain unchanged during the time span of observations. The algorithm accurately catches zero weights of the kernels that may be a concern if the number of kernels is greater than 3 in current mainstream approaches. The study presents numerical tests of the BRDF retrieval algorithm for various surface and atmospheric conditions.

Radkevich, Alexander↗

NASA NTRSDate not supplied

TPSAS-NF1676L-29044-DND

In Langley NASA, Clouds and the Earth’s Radiant Energy System (CERES) and Moderate Resolution Imaging Spectroradiometer (MODIS) are merged with Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) on the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) and CloudSat Cloud Profiling Radar (CPR). The CERES merged product (C3M) matches up to three CALIPSO footprints with each MODIS pixel along its ground track. It then assigns the nearest CloudSat footprint to each of those MODIS pixels. The cloud properties from MODIS, retrieved using the CERES algorithms, are included in C3M with the matched CALIPSO and CloudSat products along with radiances from 18 MODIS channels. The dataset is used to validate the CERES retrieved MODIS cloud properties and the computed TOA and surface flux difference using MODIS or CALIOP/CloudSAT retrieved clouds. This information is then used to tune the computed fluxes to match the CERES observed TOA flux. The capability of viewing C3M product on browsers will be available on CERES ordering tool web site in the future, https://ceres.larc.nasa.gov

David R Doelling↗

NASA NTRS1998

Cloud Filtering Using a Bi-Spectral Spatial Coherence Approach

The research in this paper focuses on describing a technique developed for cloud filtering using a bi-spectral approach on GOES-8/9 Imager data. The application was developed for use with infrared retrievals of geophysical parameters in mind, where cloud cover contaminates the derived product. However, numerous potential applications of the technique exist. The technique will be described and a preliminary validation of the algorithm will be presented. Although initially based on the spatial coherence approach from Coakley and Brethereton (1982), it has evolved to utilize a difference image of the 11 and 3.9 micrometer channels on the GOES-8/9 Imager. This image is very similar to that produced by Nelson and Ellrod (1996). During the daytime the technique makes use of the varying solar reflectance in the 3.9 micrometer channel by clouds and land to identify cloudy pixels. While at night, the technique makes use of the varying emissivity of the clouds in the scene to discriminate between clear and cloudy pixels. The technique applies three basic threshold tests to produce the final cloud filtered image: 1) a standard deviation threshold to detect the spatial variance in the scene, 2) a difference threshold between adjacent pixels, and 3) a simple infrared temperature threshold. The first test is applied to the entire image at once, then in a second pass the next two tests are applied. The final infrared temperature threshold is only meant to identify the coldest clouds that might pass the previous tests. The technique performs well during the daytime, while nighttime performance is degraded but is promising. The technique has proven to be robust and shows great promise of meeting its original goal of cloud filtering for use in an infrared retrieval algorithm for use in climate studies.

Guillory, Anthony R.↗