Threshold Based Regional Cropland Classification from Synthetic Aperture Radar using the Coefficient of Variation Method
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
The TRMM Precipitation Radar is well suited to statistical methods in that the measurements over any given region are sparsely sampled in time. Moreover, the instantaneous rain rate estimates are often of limited accuracy at high rain rates because of attenuation effects and at light rain rates because of receiver sensitivity. For the estimation of the time-averaged rain characteristics over an area both errors are relevant. By enlarging the space-time region over which the data are collected, the sampling error can be reduced. However. the bias and distortion of the estimated rain distribution generally will remain if estimates at the high and low rain rates are not corrected. In this paper we use the TRMM PR data to investigate the behavior of 2 statistical methods the purpose of which is to estimate the rain rate over large space-time domains. Examination of large-scale rain characteristics provides a useful starting point. The high correlation between the mean and standard deviation of rain rate implies that the conditional distribution of this quantity can be approximated by a one-parameter distribution. This property is used to explore the behavior of the area-time-integral (ATI) methods where fractional area above a threshold is related to the mean rain rate. In the usual application of the ATI method a correlation is established between these quantities. However, if a particular form of the rain rate distribution is assumed and if the ratio of the mean to standard deviation is known, then not only the mean but the full distribution can be extracted from a measurement of fractional area above a threshold. The second method is an extension of this idea where the distribution is estimated from data over a range of rain rates chosen in an intermediate range where the effects of attenuation and poor sensitivity can be neglected. The advantage of estimating the distribution itself rather than the mean value is that it yields the fraction of rain contributed by the light and heavy rain rates. This is useful in estimating the fraction of rainfall contributed by the rain rates that go undetected by the radar. The results at high rain rates provide a cross-check on the usual attenuation correction methods that are applied at the highest resolution of the instrument.
The Antarctic Ice Sheet is surrounded by large floating ice shelves that spread under their own weight into the ocean. Ice shelf rigidity depends on ice temperature and fabrics, and is influenced by ice flow and the delicate balance between bottom and surface accumulation. Here, we use an inverse control method to infer the rigidity of the Ronne Ice Shelf that best matches observations of ice velocity from satellite radar interferometry. Ice rigidity, or flow law parameter B, is shown to vary between 300 and 900 kPa a(sup 1/3). Ice is softer along the side margins due to frictional heating, and harder along the outflow of large glaciers, which advect cold continental ice. Melting at the bottom surface of the ice shelf increases its rigidity, while freezing decreases it. Accurate numerical modelling of ice shelf flow must account for this spatial variability in mechanical characteristics.
Purpose: Respiratory gating is commonly used to reduce motion artifacts in positron emission tomography (PET). Clinically established methods for respiratory gating in PET require contact to the patient or a direct optical line between the sensor and the patient’s torso and time consuming preparation. In this work, a contactless method for capturing a respiratory signal during PET is presented based on continuous-wave radar. Methods: The proposed method relies on the principle of emitting an electromagnetic wave and detecting the phase shift of the reflected wave, modulated due to the respiratory movement of the patient’s torso. A 24 GHz carrier frequency was chosen allowing wave propagation through plastic and clothing with high reflections at the skin surface. A detector module and signal processing algorithms were developed to extract a quantitative respiratory signal. The sensor was validated using a high precision linear table. During volunteer measurements and [{sup 18}F] FDG PET scans, the radar sensor was positioned inside the scanner bore of a PET/computed tomography scanner. As reference, pressure belt (one volunteer), depth camera-based (two volunteers, two patients), and PET data-driven (six patients) signals were acquired simultaneously and the signal correlation was quantified. Results: The developed system demonstrated a high measurement accuracy for movement detection within the submillimeter range. With the proposed method, small displacements of 25 μm could be detected, not considerably influenced by clothing or blankets. From the patient studies, the extracted respiratory radar signals revealed high correlation (Pearson correlation coefficient) to those derived from the external pressure belt and depth camera signals (r = 0.69–0.99) and moderate correlation to those of the internal data-driven signals (r = 0.53–0.70). In some cases, a cardiac signal could be visualized, due to the representation of the mechanical heart motion on the skin. Conclusions: Accurate respiratory signals were obtained successfully by the proposed method with high spatial and temporal resolution. By working without contact and passing through clothing and blankets, this approach minimizes preparation time and increases the convenience of the patient during the scan.
A representative soil surface roughness parameterization needed for the retrieval of soil moisture from active microwave satellite observation is difficult to obtain through either in-situ measurements or remote sensing-based inversion techniques. Typically, for the retrieval of soil moisture, temporal variations in surface roughness are assumed to be negligible. Although previous investigations have suggested that this assumption might be reasonable for natural vegetation covers (Moran et al. 2002, Thoma et al. 2006), insitu measurements over plowed agricultural fields (Callens et al. 2006) have shown that the soil surface roughness can change considerably over time. This paper reports on the temporal stability of surface roughness effects on radar observations and soil moisture retrieved from these radar observations collected once a week during a corn growth cycle (May 10th - October 2002). The data set employed was collected during the Optimizing Production Inputs for Economic and Environmental Enhancement (OPE3) field campaign covering this 2002 corn growth cycle and consists of dual-polarized (HH and VV) L-band (1.6 GHz) acquired at view angles of 15, 35, and 55 degrees. Cross-polarized L baud radar data were also collected as part of this experiment, but are not used in the analysis reported on here. After accounting for vegetation effects on radar observations, time-invariant optimum roughness parameters were determined using the Integral Equation Method (IEM) and radar observations acquired over bare soil and cropped conditions (the complete radar data set includes entire corn growth cycle). The optimum roughness parameters, soil moisture retrieval uncertainty, temporal distribution of retrieval errors and its relationship with the weather conditions (e.g. rainfall and wind speed) have been analyzed. It is shown that over the corn growth cycle, temporal roughness variations due to weathering by rain are responsible for almost 50% of soil moisture retrieval uncertainty depending on the sensing configuration. The effects of surface roughness variations are found to be smallest for observations acquired at a view angle of 55 degrees and HH polarization. A possible explanation for this result is that at 55 degrees and HH polarization the effect of vertical surface height changes on the observed radar response are limited because the microwaves travel parallel to the incident plane and as a result will not interact directly with vertically oriented soil structures.
Weather radars, recording information about precipitation around the globe, will soon be significantly upgraded. Most of today s weather radars transmit and receive microwave energy with horizontal orientation only, but upgraded systems have the capability to send and receive both horizontally and vertically oriented waves. These enhanced "dual-polarimetric" (DP) radars peer into precipitation and provide information on the size, shape, phase (liquid / frozen), and concentration of the falling particles (termed hydrometeors). This information is valuable for improved rain rate estimates, and for providing data on the release and absorption of heat in the atmosphere from condensation and evaporation (phase changes). The heating profiles in the atmosphere influence global circulation, and are a vital component in studies of Earth s changing climate. However, to provide the most accurate interpretation of radar data, the radar must be properly calibrated and data must be quality controlled (cleaned) to remove non-precipitation artifacts; both of which are challenging tasks for today s weather radar. The DP capability maximizes performance of these procedures using properties of the observed precipitation. In a notable paper published in 2005, scientists from the Cooperative Institute for Mesoscale Meteorological Studies (CIMMS) at the University of Oklahoma developed a method to calibrate radars using statistically averaged DP measurements within light rain. An additional publication by one of the same scientists at the National Severe Storms Laboratory (NSSL) in Norman, Oklahoma introduced several techniques to perform quality control of radar data using DP measurements. Following their lead, the Topical Rainfall Measuring Mission (TRMM) Satellite Validation Office at NASA s Goddard Space Flight Center has fine-tuned these methods for specific application to the weather radar at Kwajalein Island in the Republic of the Marshall Islands, approximately 2100 miles southwest of Hawaii and 1400 miles east of Guam in the tropical North Pacific Ocean. This tropical oceanic location is important because the majority of rain, and therefore the majority of atmospheric heating, occurs in the tropics where limited ground-based radar data are available.
Determination of the antenna pattern is important for a spaceborne Synthetic Aperture Radar such as Shuttle Imaging Radar-B (SIR-B). For SIR-B the antenna was so large that apart from one section, no complete pattern could be measured on the ground. Attempts were made to measure the pattern while the shuttle was in space by using ground receivers and active radar calibrators. The method used and described is a supplement to these measurements. The vertical pattern of an antenna can be extracted from radar signals returned from regions whose scattering coefficients versus incidence angle characteristics are suitably flat and uniform. The method used shows that the main vertical lobe of the SIR-B antenna is slightly wider than previously reported (6.9 deg at 3 dB points versus 6.2 deg used in radiometric corrections).
Using a two-dimensional linear least-squares method applied to Doppler radar data, we test the viability of determining F-factor remotely. The ultimate application of such an algorithm will be supplying real-time F-factor maps, derived from ground-based Doppler radars to air traffic control personnel and pilots. Data from NASA deployments to the MIT/Lincoln Lab TDWR testbed radar in Orlando in 1991 and 1992 along with NASA deployments to the NCAR TDWR testbed radar in Denver are examined. Preliminary analyses show that the two-dimensional method correlates reasonably well with in situ measurements. Several effects, independent of the method used, act to reduce the correlation to less than one. These include time differences between radar and aircraft data, vertical misalignment between the aircraft and the radar beam, different spatial resolution scales between aircraft and radar data, inhomogeneous radar beam filling, noise in radar data that eludes filtering, and phase lag between time and space due to low pass filtering of the aircraft data. In the final assessment, it appears that a shear-based F-factor algorithm is preferable to the currently implemented TDWR algorithms which lack any local shear estimates.
A method for analyzing multiple Doppler radar data is briefly described. A feature of this method is a new automated objective unfolding scheme which unfolds Doppler velocities for all the radars simultanceously rather than one radar at a time. The method has been used to study a region of widespread precipitation trailing a squall line (anvil), using tripple Doppler radar data to calculate horizontal winds, divergence, and vertical air velocity over an area 120 km by 120 km by 9 km deep. The method yields vertical profiles of divergence which verify the extended VAD method results of Matejka and Srivastava (1981) for the same anvil, and provides a more detailed picture of the variation of winds, divergence, and vertical air velocity in the anvil behind the squall line.
Advanced Air Mobility (AAM) aircraft require precision approach and landing systems (PALS) in several types of environments such as urban, suburban, and rural. It is difficult to implement current state-of-the-art methods approved for automated approach and landing for AAM operations. However, existing technology and systems that use vision, IR, radar, and GPS methods provide baseline perception and sensing requirements for AAM aircraft approach and landing. This paper focuses on vision-based PAL and computer vision feature correspondence methods to demonstrate a baseline navigation system while adhering to the Federal Aviation Administration requirements and regulations. The coplanar algorithm determines pose estimation, which feeds into an Extended Kalman filter that combines IMU with vision to create a sensor fusion navigation solution for GPS-denied environments. The state estimate leads to glideslope and localizer error computations, which will be pertinent for designing and deriving guidance laws and control laws for AAM PALS. The IMU and vision navigation solution provides promising simulation results for AAM PALS. This paper builds on previous work by incorporating high fidelity simulations with computer graphics rendering to demonstrate a distributed sensor network to track an AAM aircraft during approach and landing to compare with the aircraft's onboard navigation solution.
A computerized method developed by Rosenfeld (1987) is used to track the three-dimensional structure of the rainshafts of about 3000 summer afternoon convective rain cells in the semiarid region of central South Africa. By assuming a fixed relationship between the radar reflectivity factor and the rain intensity, time- and area-integrated rain volumes of each individual rainshaft are obtained at the cloudbase level and at a lower level. These results are used to obtain the cumulative fractional evaporation of the falling rain. This permits the cumulative evaporation of the falling rain to be obtained in a quantitative way as a function of the fall distance from the cloud base level and the rain intensity at the cloud base.
Bistatic-radar data reduction and Explorer 35 reflectivity and bandwidth tables with polarization spectra
The connection of the radar characteristics of cloud systems with the genetics of their formation, the microphysical structure, and the spatial extent of the clouds is examined. The relation of the probability of detecting clouds at various distances by their microphysical structure and the capabilities of the equipment are also considered. A method of discriminating cloud systems is proposed on the basis of obtained relations.
High-resolution multi-frequency radars when networked together, form a robust and seamless monitoring system for high-valued items. The multi-look-direction radar network investigated by LLNL and DSI, operates at 77 GHz and 5.8GHz, and is capable of detection at centimeter resolution movement of objects of interest. This dual frequency radar (unlike optical methods) can operate under various environments with conditions such a smoke, dust, and other natural or manmade obscurations. In addition, the multiband network radar is based on low-cost FCC approved EM specifications and can be scaled with the dimensions of the room and objects under consideration by user-defined inputs. Furthermore, the underlying radar signal processing based on software-defined radar engine is adaptive to its operational environment, where operators can input the size of the facility and object of interest and the distance between the object and radar. Finally, this multi-frequency radar employs a host of change detection and machine learning algorithms to reach performance levels in terms of high probability of detection and low false alarm rates while the monitoring state of the object can be reported remotely.
The absolute accuracy of instrumentation radar systems, refined methods of calibrating these systems, and the techniques employed in processing the associated data. A world-wide network of C-Band instrumentation radars augmented by lasers and other tracking instrumentation systems were used. The NASA WFC AN/FPQ-6 instrumentation radar and the AN/FPS-16 instrumentation radar also located at NASA WFC were the primary instruments used in the accuracy and calibration evaluations. The results achieved at WFC were then disseminated to other Ranges where they were verified, augmented and used as part of routine operations.
We present continuous estimates of snow and firn density, layer depth and accumulation from a multi-channel, multi-offset, ground-penetrating radar traverse. Our method uses the electromagnetic velocity, estimated from waveform travel-times measured at common-midpoints between sources and receivers. Previously, common-midpoint radar experiments on ice sheets have been limited to point observations. We completed radar velocity analysis in the upper ~2 m to estimate the surface and average snow density of the Greenland Ice Sheet. We parameterized the Herron and Langway (1980) firn density and age model using the radar-derived snow density, radar-derived surface mass balance (2015–2017) and reanalysis-derived temperature data. We applied structure-oriented filtering to the radar image along constant age horizons and increased the depth at which horizons could be reliably interpreted. We reconstructed the historical instantaneous surface mass balance, which we averaged into annual and multidecadal products along a 78 km traverse for the period 1984–2017. We found good agreement between our physically constrained parameterization and a firn core collected from the dry snow accumulation zone, and gained insights into the spatial correlation of surface snow density.
The motivation for this research is to develop a precipitation classification and rain rate estimation method using cloud radar-only measurements for Atmospheric Radiation Measurement (ARM) long-term cloud observation analysis, which are crucial and unique for studying cloud lifecycle and precipitation features under different weather and climate regimes. Based on simultaneous and collocated observations of the Ka-band ARM zenith radar (KAZR), two precipitation radars (NCAR S-PolKa and Texas A&M University SMART-R), and surface precipitation during the DYNAMO/AMIE field campaign, a new cloud radar-only based precipitation classification and rain rate estimation method has been developed and evaluated. The resulting precipitation classification is equivalent to those collocated SMART-R and S-PolKa observations. Both cloud and precipitation radars detected about 5% precipitation occurrence during this period. The convective (stratiform) precipitation fraction is about 18% (82%). The 2-day collocated disdrometer observations show an increased number concentration of large raindrops in convective rain compared to dominant concentration of small raindrops in stratiform rain. The composite distributions of KAZR reflectivity and Doppler velocity also show two distinct structures for convective and stratiform rain. These indicate that the method produces physically consistent results for two types of rain. The cloud radar-only rainfall estimation is developed based on the gradient of accumulative radar reflectivity below 1 km, near-surface Ze, and collocated surface rainfall (R) measurement. The parameterization is compared with the Z-R exponential relation. The relative difference between estimated and surface measured rainfall rate shows that the two-parameter relation can improve rainfall estimation.
In writing the integral equations for the median mass diameter and number concentration, or comparable parameters of the raindrop size distribution, it is apparent that the forms of the equations for dual-polarization and dual-wavelength radar data are identical when attenuation effects are included. The differential backscattering and extinction coefficients appear in both sets of equations: for the dual-polarization equations, the differences are taken with respect to polarization at a fixed frequency while for the dual-wavelength equations, the differences are taken with respect to frequency at a fixed polarization. An alternative to the integral equation formulation is that based on the k-Z (attenuation coefficient-radar reflectivity factor) parameterization. This-technique was originally developed for attenuating single-wavelength radars, a variation of which has been applied to the TRMM Precipitation Radar data (PR). Extensions of this method have also been applied to dual-polarization data. In fact, it is not difficult to show that nearly identical equations are applicable as well to dualwavelength radar data. In this case, the equations for median mass diameter and number concentration take the form of coupled, but non-integral equations. Differences between this and the integral equation formulation are a consequence of the different ways in which attenuation correction is performed under the two formulations. For both techniques, the equations can be solved either forward from the radar outward or backward from the final range gate toward the radar. Although the forward-going solutions tend to be unstable as the attenuation out to the range of interest becomes large in some sense, an independent estimate of path attenuation is not required. This is analogous to the case of an attenuating single-wavelength radar where the forward solution to the Hitschfeld-Bordan equation becomes unstable as the attenuation increases. To circumvent this problem, the equations can be expressed in the form of a final-value problem so that the recursion begins at the far range gate and proceeds inward towards the radar. Solving the problem in this way traditionally requires estimates of path attenuation to the final gate: in the case of orthogonal linear polarizations, the attenuations at horizontal and vertical polarizations (same frequency) are required while in the dual-wavelength case, attenuations at the two frequencies (same polarization) are required.