SEARCH · Engineering Papers

Results for “Dynamic range”

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,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS1998

Resiliency of the Multiscale Retinex Image Enhancement Algorithm

The multiscale retinex with color restoration (MSRCR) continues to prove itself in extensive testing to be very versatile automatic image enhancement algorithm that simultaneously provides dynamic range compression, color constancy, and color rendition, However, issues remain with regard to the resiliency of the MSRCR to different image sources and arbitrary image manipulations which may have been applied prior to retinex processing. In this paper we define these areas of concern, provide experimental results, and, examine the effects of commonly occurring image manipulation on retinex performance. In virtually all cases the MSRCR is highly resilient to the effects of both the image source variations and commonly encountered prior image-processing. Significant artifacts are primarily observed for the case of selective color channel clipping in large dark zones in a image. These issues are of concerning the processing of digital image archives and other applications where there is neither control over the image acquisition process, nor knowledge about any processing done on th data beforehand.

Rahman, Zia-Ur↗

NASA NTRS2000

Retinex Preprocessing for Improved Multi-Spectral Image Classification

The goal of multi-image classification is to identify and label "similar regions" within a scene. The ability to correctly classify a remotely sensed multi-image of a scene is affected by the ability of the classification process to adequately compensate for the effects of atmospheric variations and sensor anomalies. Better classification may be obtained if the multi-image is preprocessed before classification, so as to reduce the adverse effects of image formation. In this paper, we discuss the overall impact on multi-spectral image classification when the retinex image enhancement algorithm is used to preprocess multi-spectral images. The retinex is a multi-purpose image enhancement algorithm that performs dynamic range compression, reduces the dependence on lighting conditions, and generally enhances apparent spatial resolution. The retinex has been successfully applied to the enhancement of many different types of grayscale and color images. We show in this paper that retinex preprocessing improves the spatial structure of multi-spectral images and thus provides better within-class variations than would otherwise be obtained without the preprocessing. For a series of multi-spectral images obtained with diffuse and direct lighting, we show that without retinex preprocessing the class spectral signatures vary substantially with the lighting conditions. Whereas multi-dimensional clustering without preprocessing produced one-class homogeneous regions, the classification on the preprocessed images produced multi-class non-homogeneous regions. This lack of homogeneity is explained by the interaction between different agronomic treatments applied to the regions: the preprocessed images are closer to ground truth. The principle advantage that the retinex offers is that for different lighting conditions classifications derived from the retinex preprocessed images look remarkably "similar", and thus more consistent, whereas classifications derived from the original images, without preprocessing, are much less similar.

Thompson, B.↗

NASA NTRS1999

Observations of Land Surface Variability Using Passive Microwave Sensing

Understanding the global variability of land surface wetness (soil moisture), skin temperature, and related surface fluxes of heat and moisture is key to assessing the importance of the land surface in influencing climate. The feasibility of producing model estimates of these quantities is being studied as part of the International Satellite Land Surface Climatology Project (ISLSCP) Global Soil Wetness Project (GSWP). In the GSWP approach, meteorological observations and analyses are used to drive global circulation models. Satellite measurements can provide independent estimates of key land surface parameters that are needed for initializing and validating the climate models and for monitoring long-term change. Satellite observations of the land surface can also be assimilated into soil models to estimate moisture in the root zone. In our research, passive microwave satellite data recorded during 1978-1987 from the Nimbus-7 Scanning Multichannel Microwave Radiometer (SMMR) are being used to examine spatial and temporal trends in surface soil moisture, vegetation, and temperature. These data include observations at C and X bands (6.6 and 10.7 GHz), which are not available on the current Special Sensor Microwave/Imager (SSM/I) and are precursors to data that will become available from the Advanced Microwave Scanning Radiometer (AMSR) on Advanced Earth Observing Satellite (ADEOS-II) and Earth Observing System (EOS) PM1 in the year 2000. A chart shows a time-series of SMMR-derived surface temperature, T-e and surface soil moisture M, retrieved on a 0.5 deg x 0.5 deg grid and further averaged over a 4 deg x 10 deg study region in the African Sahel. Also shown are National Center for Environmental Prediction (NCEP) model outputs of surface temperature, T-sfc, and soil wetness, Soil-w. The variables have been scaled to have similar dynamic ranges on the plots. The NCEP data from the NCEP Reanalysis Project are monthly averages on a 2.5 deg x 2.5 deg grid averaged over the 4 deg x 10 deg study area. Comparisons of SMMR retrievals with forecast model output show the potential of the satellite data for validating model output and monitoring long-term trends. Continuing work will extend these results to other regions to validate the retrievals more quantitatively. In preparation for the launch of AMSR, field experiments are planned in collaboration with the Global Energy and Water Cycle Experiment (GEWEX) Coordinated Enhanced Observing Period (CEOP) experiments to evaluate the satellite-derived soil moisture measurements and to demonstrate their usefulness for land surface hydrology and climate. Additional information is contained in the original.

Njoku, Eni G.↗

NASA NTRS2000

Earth Observing-1 Advanced Land Imager: Radiometric Response Calibration

The Advanced Land Imager (ALI) is one of three instruments to be flown on the first Earth Observing mission (EO-1) under NASA's New Millennium Program (NMP). ALI contains a number of innovative features, including a wide field of view optical design, compact multispectral focal plane arrays, non-cryogenic HgCdTe detectors for the short wave infrared bands, and silicon carbide optics. This document outlines the techniques adopted during ground calibration of the radiometric response of the Advanced Land Imager. Results from system level measurements of the instrument response, signal-to-noise ratio, saturation radiance, and dynamic range for all detectors of every spectral band are also presented.

Mendenhall, J. A.↗

NASA NTRS2001

Sources, Sinks, and the Distribution of OH in the Lower Stratosphere

Extensive measurement campaigns by the NASA ER-2 research aircraft have obtained a nearly pole-to-pole database of the species that control HOx (OH + HO2) chemistry. The wide dynamic range of these in situ measurements provides an opportunity to demonstrate empirically the mechanisms that control the HOx system. Measurements in the lower stratosphere show a remarkably tight correlation of OH concentration with the solar zenith angle (SZA). This correlation is nearly invariant over latitudes ranging from 70 S to 90 N and all seasons. An analysis of the production and loss of HOx in terms of the rate determining steps of reaction sequences developed by Johnston and Podolske and Johnston and Kinnison is used to clarify the behavior of the system and to directly test our understanding of the system with observations. Calculations using in situ measurements show that the production rate of HOx is proportional to O3 and ultraviolet radiation flux. The loss rate is proportional to the concentration and the partitioning of NOy (reactive nitrogen) and the concentration of HO2. In the absence of heterogeneous reactions, the partitioning of NOy is controlled by O3 and NOx and the concentration of HO2 is controlled by NOy and O3, so that the removal rate of OH is buffered against changes in the correlation of O3 and NOy. The heterogeneous conversion of NO2 to HNO3 is not an important net source of HO, because production and removal sequences are nearly balanced. Changes in NOy partitioning resulting from heterogeneous chemistry have a large effect on the loss rates of HOx, but little or no impact on the measured abundance of OH. The enhanced loss rates at high NO2/HNO3 are offset in the data set examined here by enhanced production rates resulting from increased photolysis rates resulting from the decreased O3 column above the ER-2.

Hanisco, T. F.↗

NASA NTRS2000

Detection of Objects Hidden in Highly Scattering Media Using Time-Gated Imaging Methods

Non-intrusive and non-invasive optical imaging techniques has generated great interest among researchers for their potential applications to biological study, device characterization, surface defect detection, and jet fuel dynamics. Non-linear optical parametric amplification gate (NLOPG) has been used to detect back-scattered images of objects hidden in diluted Intralipid solutions. To directly detect objects hidden in highly scattering media, the diffusive component of light needs to be sorted out from early arrived ballistic and snake photons. In an optical imaging system, images are collected in transmission or back-scattered geometry. The early arrival photons in the transmission approach, always carry the direct information of the hidden object embedded in the turbid medium. In the back-scattered approach, the result is not so forth coming. In the presence of a scattering host, the first arrival photons in back-scattered approach will be directly photons from the host material. In the presentation, NLOPG was applied to acquire time resolved back-scattered images under the phase matching condition. A time-gated amplified signal was obtained through this NLOPG process. The system's gain was approximately 100 times. The time-gate was achieved through phase matching condition where only coherent photons retain their phase. As a result, the diffusive photons, which were the primary contributor to the background, were removed. With a large dynamic range and high resolution, time-gated early light imaging has the potential for improving rocket/aircraft design by determining jets shape and particle sizes. Refinements to these techniques may enable drop size measurements in the highly scattering, optically dense region of multi-element rocket injectors. These types of measurements should greatly enhance the design of stable, and higher performing rocket engines.

Galland, Pierre A.↗

NASA NTRS2001

Demodulation System for Fiber Optic Bragg Grating Dynamic Pressure Sensing

Fiber optic Bragg gratings have been used for years to measure quasi-static phenomena. In aircraft engine applications there is a need to measure dynamic signals such as variable pressures. In order to monitor these pressures a detection system with broad dynamic range is needed. This paper describes an interferometric demodulator that was developed and optimized for this particular application. The signal to noise ratio was maximized through temporal coherence analysis. The demodulator was incorporated in a laboratory system that simulates conditions to be measured. Several pressure sensor configurations incorporating a fiber optic Bragg grating were also explored. The results of the experiments are reported in this paper.

Lekki, John D.↗

NASA NTRS2001

Optimizing Orbit-Instrument Configuration for Global Precipitation Mission (GPM) Satellite Fleet

Following the scientific success of the Tropical Rainfall Measuring Mission (TRMM) spearheaded by a group of NASA and NASDA scientists, their external scientific collaborators, and additional investigators within the European Union's TRMM Research Program (EUROTRMM), there has been substantial progress towards the development of a new internationally organized, global scale, and satellite-based precipitation measuring mission. The highlights of this newly developing mission are a greatly expanded scope of measuring capability and a more diversified set of science objectives. The mission is called the Global Precipitation Mission (GPM). Notionally, GPM will be a constellation-type mission involving a fleet of nine satellites. In this fleet, one member is referred to as the "core" spacecraft flown in an approximately 70 degree inclined non-sun-synchronous orbit, somewhat similar to TRMM in that it carries both a multi-channel polarized passive microwave radiometer (PMW) and a radar system, but in this case it will be a dual frequency Ku-Ka band radar system enabling explicit measurements of microphysical DSD properties. The remainder of fleet members are eight orbit-synchronized, sun-synchronous "constellation" spacecraft each carrying some type of multi-channel PMW radiometer, enabling no worse than 3-hour diurnal sampling over the entire globe. In this configuration the "core" spacecraft serves as a high quality reference platform for training and calibrating the PMW rain retrieval algorithms used with the "constellation" radiometers. Within NASA, GPM has advanced to the pre-formulation phase which has enabled the initiation of a set of science and technology studies which will help lead to the final mission design some time in the 2003 period. This presentation first provides an overview of the notional GPM program and mission design, including its organizational and programmatic concepts, scientific agenda, expected instrument package, and basic flight architecture. Following this introduction, we focus specifically on the last topic, that being an analysis which leads to an optimal flight architecture dictated in part by science requirements but constrained by allowable orbital mechanics, instrument scan patterns, and antenna aperture properties. Because the optimal architecture involves an interplay between orbit mechanics and instrument specifications, it is important to recognize that in attempting to serve various scientific themes, the final optimal architecture will represent a compromise concerning dynamic range, spatial resolution, sampling interval, pointing, beam coincidence, and measurement uncertainty. Moreover, cost becomes a major factor in seeking the optimal architecture through the pathways of antenna and instrument scan designs, as well as propulsion requirements associated with the orbit heights of various "constellation" members. Although the results presented at the IGARSS-2001 meeting will likely not be the fully refined flight architecture specifications, they are expected to be nearly complete.

Smith, Eric A.↗

NASA NTRS1994

Nonlinear Optical Image Processing with Bacteriorhodopsin Films

The transmission properties of some bacteriorhodopsin film spatial light modulators are uniquely suited to allow nonlinear optical image processing operations to be applied to images with multiplicative noise characteristics. A logarithmic amplitude transmission feature of the film permits the conversion of multiplicative noise to additive noise, which may then be linearly filtered out in the Fourier plane of the transformed image. The bacteriorhodopsin film displays the logarithmic amplitude response for write beam intensities spanning a dynamic range greater than 2.0 orders of magnitude. We present experimental results demonstrating the principle and capability for several different image and noise situations, including deterministic noise and speckle. Using the bacteriorhodopsin film, we successfully filter out image noise from the transformed image that cannot be removed from the original image.

Downie, John D.↗

NASA NTRS1994

Contrast Gain Control Model Fits Masking Data

We studied the fit of a contrast gain control model to data of Foley (JOSA 1994), consisting of thresholds for a Gabor patch masked by gratings of various orientations, or by compounds of two orientations. Our general model includes models of Foley and Teo & Heeger (IEEE 1994). Our specific model used a bank of Gabor filters with octave bandwidths at 8 orientations. Excitatory and inhibitory nonlinearities were power functions with exponents of 2.4 and 2. Inhibitory pooling was broad in orientation, but narrow in spatial frequency and space. Minkowski pooling used an exponent of 4. All of the data for observer KMF were well fit by the model. We have developed a contrast gain control model that fits masking data. Unlike Foley's, our model accepts images as inputs. Unlike Teo & Heeger's, our model did not require multiple channels for different dynamic ranges.

Watson, Andrew B.↗

NASA NTRS2001

Effects of the Ionosphere on Passive Microwave Remote Sensing of Ocean Salinity from Space

Among the remote sensing applications currently being considered from space is the measurement of sea surface salinity. The salinity of the open ocean is important for understanding ocean circulation and for modeling energy exchange with the atmosphere. Passive microwave remote sensors operating near 1.4 GHz (L-band) could provide data needed to fill the gap in current coverage and to complement in situ arrays being planned to provide subsurface profiles in the future. However, the dynamic range of the salinity signal in the open ocean is relatively small and propagation effects along the path from surface to sensor must be taken into account. In particular, Faraday rotation and even attenuation/emission in the ionosphere can be important sources of error. The purpose or this work is to estimate the magnitude of these effects in the context of a future remote sensing system in space to measure salinity in L-band. Data will be presented as a function of time location and solar activity using IRI-95 to model the ionosphere. The ionosphere presents two potential sources of error for the measurement of salinity: Rotation of the polarization vector (Faraday rotation) and attenuation/emission. Estimates of the effect of these two phenomena on passive remote sensing over the oceans at L-band (1.4 GHz) are presented.

LeVine, D. M.↗

NASA NTRS2001

A Multi-Band Analytical Algorithm for Deriving Absorption and Backscattering Coefficients from Remote-Sensing Reflectance of Optically Deep Waters

A multi-band analytical (MBA) algorithm is developed to retrieve absorption and backscattering coefficients for optically deep waters, which can be applied to data from past and current satellite sensors, as well as data from hyperspectral sensors. This MBA algorithm applies a remote-sensing reflectance model derived from the Radiative Transfer Equation, and values of absorption and backscattering coefficients are analytically calculated from values of remote-sensing reflectance. There are only limited empirical relationships involved in the algorithm, which implies that this MBA algorithm could be applied to a wide dynamic range of waters. Applying the algorithm to a simulated non-"Case 1" data set, which has no relation to the development of the algorithm, the percentage error for the total absorption coefficient at 440 nm a (sub 440) is approximately 12% for a range of 0.012 - 2.1 per meter (approximately 6% for a (sub 440) less than approximately 0.3 per meter), while a traditional band-ratio approach returns a percentage error of approximately 30%. Applying it to a field data set ranging from 0.025 to 2.0 per meter, the result for a (sub 440) is very close to that using a full spectrum optimization technique (9.6% difference). Compared to the optimization approach, the MBA algorithm cuts the computation time dramatically with only a small sacrifice in accuracy, making it suitable for processing large data sets such as satellite images. Significant improvements over empirical algorithms have also been achieved in retrieving the optical properties of optically deep waters.

Lee, Zhong-Ping↗

NASA NTRS1996

Determination of Small Animal Long Bone Properties Using Densitometry

Assessment of bone structural property changes due to loading regimens or pharmacological treatment typically requires destructive mechanical testing and sectioning. Our group has accurately and non-destructively estimated three dimensional cross-sectional areal properties (principal moments of inertia, Imax and Imin, and principal angle, Theta) of human cadaver long bones from pixel-by-pixel analysis of three non-coplanar densitometry scans. Because the scanner beam width is on the order of typical small animal diapbyseal diameters, applying this technique to high-resolution scans of rat long bones necessitates additional processing to minimize errors induced by beam smearing, such as dependence on sample orientation and overestimation of Imax and Imin. We hypothesized that these errors are correctable by digital image processing of the raw scan data. In all cases, four scans, using only the low energy data (Hologic QDR-1000W, small animal mode), are averaged to increase image signal-to-noise ratio. Raw scans are additionally processed by interpolation, deconvolution by a filter derived from scanner beam characteristics, and masking using a variable threshold based on image dynamic range. To assess accuracy, we scanned an aluminum step phantom at 12 orientations over a range of 180 deg about the longitudinal axis, in 15 deg increments. The phantom dimensions (2.5, 3.1, 3.8 mm x 4.4 mm; Imin/Imax: 0.33-0.74) were comparable to the dimensions of a rat femur which was also scanned. Cross-sectional properties were determined at 0.25 mm increments along the length of the phantom and femur. The table shows average error (+/- SD) from theory of Imax, Imin, and Theta) over the 12 orientations, calculated from raw and fully processed phantom images, as well as standard deviations about the mean for the femur scans. Processing of phantom scans increased agreement with theory, indicating improved accuracy. Smaller standard deviations with processing indicate increased precision and repeatability. Standard deviations for the femur are consistent with those of the phantom. We conclude that in conjunction with digital image enhancement, densitometry scans are suitable for non-destructive determination of areal properties of small animal bones of comparable size to our phantom, allowing prediction of Imax and Imin within 2.5% and Theta within a fraction of a degree. This method represents a considerable extension of current methods of analyzing bone tissue distribution in small animal bones.

Breit, Gregory A.↗

NASA NTRS2002

Microgravity Compatible Reagentless Instrumentation for Detection of Dissolved Organic Acids and Alcohols in Potable Water

The Organic Acid and Alcohol Monitor (OAAM) program has resulted in the successful development of a computer controlled prototype analyzer capable of accurately determining aqueous organic acids and primary alcohol concentrations over a large dynamic range with high sensitivity. Formic, acetic, and propionic acid were accurately determined at concentrations as low as 5 to 10 micrograms/L in under 20 minutes, or as high as 10 to 20 mg/L in under 30 minutes. Methanol, ethanol, and propanol were determined at concentrations as low as 20 to 100 micrograms/L, or as high as 10 mg/L in under 30 minutes. Importantly for space based application, the OAAM requires no reagents or hazardous chemicals to perform these analyses needing only power, water, and CO2 free purge gas. The OAAM utilized two membrane processes to segregate organic acids from interfering ions. The organic acid concentration was then determined based upon the conductiometric signal. Separation of individual organic acids was accomplished using a chromatographic column. Alcohols are determined in a similar manner after conversion to organic acids by sequential biocatalytic and catalytic oxidation steps. The OAAM was designed to allow the early diagnosis of under performing or failing sub-systems within the Water Recovery System (WRS) baselined for the International Space Station (ISS). To achieve this goal, several new technologies were developed over the course of the OAAM program.

Akse, James R.↗

NASA NTRS1998

Effect of Free Stream Turbulence on the Flow-Induced Background Noise of In-Flow Microphones

When making noise measurements of sound sources in flow using microphones immersed in an air stream or wind tunnel, the factor limiting the dynamic range of the measurement is, in many cases, the noise caused by the flow over the microphone. To lower this self-noise, and to protect the microphone diaphragm, an aerodynamic microphone forebody is usually mounted on the tip of the omnidirectional microphone. The microphone probe is then pointed into the wind stream. Even with a microphone forebody, however, the self-noise persists, prompting further research in the area of microphone forebody design for flow-induced self-noise reduction. The magnitude and frequency characteristics of in-flow microphone probe self-noise is dependent upon the exterior shape of the probe and on the level of turbulence in the onset flow, among other things. Several recent studies present new designs for microphone forebodies, some showing the forbodies' self-noise characteristics when used in a given facility. However, these self-noise characteristics may change when the probes are used in different facilities. The present paper will present results of an experimental investigation to determine an empirical relationship between flow turbulence and self-noise levels for several microphone forebody shapes as a function of frequency. As a result, the microphone probe self-noise for these probes will be known as a function of freestream turbulence, and knowing the freestream turbulence spectra for a given facility, the probe self-noise can be predicted. Flow-induced microphone self-noise is believed to be related to the freestream. turbulence by three separate mechanisms. The first mechanism is produced by large scale, as compared to the probe size, turbulence which appears to the probe as a variation in the angle of attack of the freestream. flow. This apparent angle of attack variation causes the pressure along the probe surface to fluctuate, and at the location of the sensor orifice this fluctuating surface pressure is sensed by the diaphragm as noise. The second mechanism is caused by the convection of smaller sized turbulence, on the order of the probe cross-section, which passes nearby or strikes the probe giving rise to a fluctuating pressure at the sensor orifice. And, the third mechanism is related to fine scale turbulence through its effects on boundary layer growth and transition to a turbulent boundary layer. The method for relating the probe self-noise to the freestream turbulence will be based on the method of K. J. Young5 from Boeing, who developed the technique and presented flow noise results for a Bruel & Kjaer Type 0385, 1/4 inch (6.35 mm) nose cone. The experimental set-up used in the present experiment is similar to that of Young and is described in the present paper. Finally, flow noise predictions are made using the empirical correlations. These predictions are then compared with actual flow noise measurements made in the National Full-Scale Aerodynamics Complex 40- by 80-Foot Wind Tunnel at NASA Ames Research Center.

Allen, Christopher S.↗

NASA NTRS2002

4D PhaseCam(Trade Mark) Capabilities: Modal Analysis and Multiple-Wavelength Mirror Phasing

The PhaseCam is a dynamic phase shifting interferometer system manufactured by 4D Technology Corporation that is capable of very fast data acquisition. This rapid acquisition extends the capability of conventional interferometry to enable measurement in unstable environments, the generation of phase movies of surface shape and to facilitate modal analysis of structures. The PhaseCam hardware and software have been modified for MSFC to include synchronous modal optical measurement and analysis. These modifications will be discussed and data presented. The dynamic range of a phase shifting measurement is limited by local slope and pixel sampling to lambda/4 wave steps. Two-wavelength techniques can increase the effective measurement wavelength from microns to tens of centimeters and permit the phasing of mirror segments. A two wavelength PhaseCam will be discussed and measurement results presented.

Millerd, James E.↗

NASA NTRS2002

A Large Sparse Aperture Densified Pupil Hypertelescope Concept for Ground Based Detection of Extra-Solar Earth-Like Planets

A concept is presented for a large (10 - 30 meter) sparse aperture hyper telescope to image extrasolar earth-like planets from the ground in the presence of atmospheric seeing. The telescope achieves high dynamic range very close to bright stellar sources with good image quality using pupil densification techniques. Active correction of the perturbed wavefront is simplified by using 36 small flat mirrors arranged in a parabolic steerable array structure, eliminating the need for large delat lines and operating at near-infrared (1 - 3 Micron) wavelengths with flats comparable in size to the seeing cells.

Gezari, D.↗

NASA NTRS2002

Vibration Reduction of Helicopter Blade Using Variable Dampers: A Feasibility Study

In the report, the investigation of controlling helicopter-blade lead-lag vibration is described. Current practice of adding passive damping may be improved to handle large dynamic range of the blade with several peaks of vibration resonance. To minimize extra-large damping forces that may damage the control system of blade, passive dampers should have relatively small damping coefficients, which in turn limit the effectiveness. By providing variable damping, a much larger damping coefficient to suppress the vibration can be realized. If the damping force reaches the maximum allowed threshold, the damper will be automatically switched into the mode with smaller damping coefficient to maintain near-constant damping force. Furthermore, the proposed control system will also have a fail-safe feature to guarantee the basic performation of a typical passive damper. The proposed control strategy to avoid resonant regions in the frequency domain is to generate variable damping force in combination with the supporting stiffness to manipulate the restoring force and conservative energy of the controlled blade system. Two control algorithms are developed and verified by a prototype variable damper, a digital controller and corresponding algorithms. Primary experiments show good potentials for the proposed variable damper: about 66% and 82% reductions in displacement at 1/3 length and the root of the blade respectively.

Lee, George C.↗