SEARCH · Search NASA

Results for “Gaussian beam”

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.

210 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS2008

Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Data Processing, Sky Maps, and Basic Results

The Wilkinson Microwave Anisotropy Probe (WMAP) is a Medium-Class Explorer (MIDEX) satellite aimed at elucidating cosmology through full-sky observations of the cosmic microwave background (CMB). The WMAP full-sky maps of the temperature and polarization anisotropy in five frequency bands provide our most accurate view to date of conditions in the early universe. The multi-frequency data facilitate the separation of the CMB signal from foreground emission arising both from our Galaxy and from extragalactic sources. The CMB angular power spectrum derived from these maps exhibits a highly coherent acoustic peak structure which makes it possible to extract a wealth of information about the composition and history of the universe. as well as the processes that seeded the fluctuations. WMAP data have played a key role in establishing ACDM as the new standard model of cosmology (Bennett et al. 2003: Spergel et al. 2003; Hinshaw et al. 2007: Spergel et al. 2007): a flat universe dominated by dark energy, supplemented by dark matter and atoms with density fluctuations seeded by a Gaussian, adiabatic, nearly scale invariant process. The basic properties of this universe are determined by five numbers: the density of matter, the density of atoms. the age of the universe (or equivalently, the Hubble constant today), the amplitude of the initial fluctuations, and their scale dependence. By accurately measuring the first few peaks in the angular power spectrum, WMAP data have enabled the following accomplishments: Showing the dark matter must be non-baryonic and interact only weakly with atoms and radiation. The WMAP measurement of the dark matter density puts important constraints on supersymmetric dark matter models and on the properties of other dark matter candidates. With five years of data and a better determination of our beam response, this measurement has been significantly improved. Precise determination of the density of atoms in the universe. The agreement between the atomic density derived from WMAP and the density inferred from the deuterium abundance is an important test of the standard big bang model. Determination of the acoustic scale at redshift z = 1090. Similarly, the recent measurement of baryon acoustic oscillations (BAO) in the galaxy power spectrum (Eisenstein et al. 2005) has determined the acoustic scale at redshift z approx. 0.35. When combined, these standard rulers accurately measure the geometry of the universe and the properties of the dark energy. These data require a nearly flat universe dominated by dark energy consistent with a cosmological constant. Precise determination of the Hubble Constant, in conjunction with BAO observations. Even when allowing curvature (Omega(sub 0) does not equal 1) and a free dark energy equation of state (w does not equal -1), the acoustic data determine the Hubble constant to within 3%. The measured value is in excellent agreement with independent results from the Hubble Key Project (Freedman et al. 2001), providing yet another important consistency test for the standard model. Significant constraint of the basic properties of the primordial fluctuations. The anti-correlation seen in the temperature/polarization (TE) correlation spectrum on 4deg scales implies that the fluctuations are primarily adiabatic and rule out defect models and isocurvature models as the primary source of fluctuations (Peiris et al. 2003).

Weiland, J.L.↗

NASA NTRS1976

Dielectric breakdown induced by picosecond laser pulses

The damage thresholds of transparent optical materials were investigated. Single picosecond pulses at 1.06 microns, 0.53 microns and 0.35 microns were obtained from a mode locked Nd-YAG oscillator-amplifier-frequency multiplier system. The pulses were Gaussian in space and time and permitted the determination of breakdown thresholds with a reproducibility of 15%. It was shown that the breakdown thresholds are characteristic of the bulk material, which included nine alkali halides, five different laser host materials, KDP, quartz, sapphire and calcium fluoride. The extension of the damage data to the ultraviolet is significant, because some indication was obtained that two- and three-photon absorption processes begin to play a role in determining the threshold. Throughout the visible region of the spectrum the threshold is still an increasing function of frequency, indicating that avalanche ionization is the dominant factor in determining the breakdown threshold. This was confirmed by a detailed study of the damage morphology with a high resolution microscope just above the threshold. The influence of self focusing is discussed, and evidence for beam distortion below the power threshold for complete self focusing is presented, confirming the theory of Marburger.

Smith, W. L.↗

NASA NTRS2010

Five-Year Wilkinson Microwave Anisotropy Probe Observations: Data Processing, Sky Maps, and Basic Results

We present new full-sky temperature and polarization maps in five frequency bands from 23 to 94 GHz, based on data from the first five years of the Wilkinson Microwave Anisotropy Probe (WMAP) sky survey. The new maps are consistent with previous maps and are more sensitive. The five-year maps incorporate several improvements in data processing made possible by the additional years of data and by a more complete analysis of the instrument calibration and in-flight beam response. We present several new tests for systematic errors in the polarization data and conclude that W-band polarization data is not yet suitable for cosmological studies, but we suggest directions for further study. We do find that Ka-band data is suitable for use; in conjunction with the additional years of data, the addition of Ka band to the previously used Q- and V-band channels significantly reduces the uncertainty in the optical depth parameter, tau. Further scientific results from the five-year data analysis are presented in six companion papers and are summarized in Section 7 of this paper. With the five-year WMAP data, we detect no convincing deviations from the minimal six-parameter ACDM model: a flat universe dominated by a cosmological constant, with adiabatic and nearly scale-invariant Gaussian fluctuations. Using WMAP data combined with measurements of Type Ia supernovae and Baryon Acoustic Oscillations in the galaxy distribution, we find (68% CL uncertainties): OMEGA(sub b)h(sup 2) = 0.02267(sup +0.00058)(sub -0.00059), OMEGA(sub c)h(sup 2) = 0.1131 plus or minus 0.0034, OMEGA(sub logical and) = 0.726 plus or minus 0.015, ns = .960 plus or minus 0.013, tau = 0.84 plus or minus 0.016, and DELTA(sup 2)(sub R) = (22.445 plus or minus 0.096) x 10(exp -9) at k = 0.002 Mpc(exp -1). From these we derive sigma(sub 8) = 0.812 plus or minus 0.026, H(sub 0) = 70.5 plus or minus 1.3 kilometers per second Mpc(exp -1), OMEGA(sub b) = 0.0456 plus or minus 0.0015, OMEGA(sub c) = .228 plus or minus 0.013, OMEGA(sub m)h(sup 2) = 0.1358(sup +0.0037)(sub -0.0036), z reion = 10.9 plus or minus 1.4, and t(sub 0) = 13.72 plus or minus 0.12 Gyr. The new limit on the tensor-to-scalar ration is r less than 0.22 (95% CL), while the evidence for a running spectral index is insignificant, dn(sub s)/d ln k = -0.028 plus or minus 0.020 (68% CL). We obtain tight, simultaneous limits on the (constant) dark energy equation of state and the spatial curvature of the universe: -0.14 less than 1 + w less than 0.12 (95% CL) and -0.0179 less than OMEGA(sub k) less than 0.0081 (95% CL). The number of relativistic degrees of freedom, expressed in units of the effective number of neutrino species, is found to be N(sub eff) = 4.4 plus or minus 1.5 (69% CL), consistent with the standard value of 3.04. Models with N(sub eff) = 0 are disfavored at greater than 99% confidence. Finally, new limits on physically motivated primordial non-Gaussianity parameters are -9 less than f(sup local)(sub NL) less than 111 (95% CL) and -151 less than f(sup equal)(sub NL) less than 253 (95% CL) for the local and equilateral models, respectively.

Hinshaw, G.↗

NASA NTRS1990

NGC 1058: Gas motions in an extended, quiescent spiral disk

Researchers investigate in detail the motion of gas in the galaxy NGC 1058 using the very large array (VLA) to map the emission in the 21-cm line. This galaxy is so nearly face-on that the contribution to the line width due to the variation of the rotational velocity across the D-array beam is small compared with the random z-motion of the gas. Researchers confirm results of earlier studies (Lewis 1987, A. and A. Suppl., 63, 515; van der Kruit and Shostak 1984, A. and A., 134, 258) of the galaxy's total neutral hydrogen (HI) and kinematics, including the fact that the rotation curve drops faster than Keplerian at the outer edge of the disk, which is interpreted as a fortuitous twist of the plane of rotation in the outer disk. However, their very high velocity resolution (2.58 km s(exp -1) after Hanning smoothing) coupled with good spatial resolution, allows researchers to measure more accurately the line width, and even to some extent its shape, throughout the disk. One of the most interesting results of this study is the remarkable constancy of the line width in the outer disk. From radius 90 to 210 seconds the Gaussian velocity dispersion (sigma sub nu) of the 21-cm line has a mean value of 5.7 km s(exp -1) (after correcting for the spectral resolution) with a dispersion of less than 0.9 km s(exp -1) (after correcting for the spectral resolution) with a dispersion of less than 0.9 km s(exp -1). Translating this directly into a kinetic temperature (Doppler temperature): T sub Dopp equals 121K (sigma sub mu exp 2/(km s(exp -1) (exp 2) gives 4000 K, with a dispersion of less than 1500 K over the outer disk. This constancy is observed even when comparing the spiral arms versus inter-arm regions, which in the radius range from 100 to 150 seconds the surface density modulates (defined as the ratio N sub peak -N sub trough/N sub peak + N sub trough) from 0.5 to 0.25 in the range 150 to 200 seconds.

Hanson, Margaret Murray↗

NASA NTRS1993

Measurements of the anisotropy of the cosmic background radiation at 0.5 deg scale near the star Mu Pegasi

Results are presented from the third flight of the MAX experiment, an attitude-controlled balloon-borne millimeter-wave telescope with a 0.5 deg beam, a 1 deg chop, and a three-channel bolometric photometer. Several hours of high-quality data were obtained during a flight on 1991 June 5, including long integrations to search for CBR anisotropy, two separate measurements of dust in the Galactic plane, a brief scan of the Coma Cluster to search for the Sunyaev-Zel'dovich (SZ) effect, and a number of important systematic tests. Data from one of the long CBR integrations, carried out in a region of sky near the star Mu Pegasi, are presented. The primary structure in the data is shown to be emission from Galactic dust via its spectrum and correlation with the IRAS 100/micron map. Several approaches are used to fit this dust component and remove it from the data. An upper limit to CBR anisotropy of deltaT/T less than 2.5 x 10 exp -5 is obtained for a Gaussian autocorrelation function with coherence angle omega(c) = 25'. This limit is significantly higher than the measurement sensitivity of deltaT/T about 1 x 10 exp -5 due to the presence of residual structure in the data after removal of the dust component.

Meinhold, P.↗

NASA NTRS2000

Surface Change Detection Using Large Footprint Laser Altimetry

Laser altimeters provide a precise and accurate method for mapping topography at fine horizontal and vertical scales. A laser altimeter provides range by measuring the roundtrip flight time of a short pulse of laser light from the laser altimeter instrument to the target surface. The range is then combined with laser beam pointing knowledge and absolute position knowledge to provide an absolute measurement of the surface topography. Newer generations of laser altimeters measure the range by recording the shape and time of the outgoing and received laser pulses. The shape of the return pulse can also provide unique information about the vertical structure of material such as vegetation within each laser footprint. Distortion of the return pulse is caused by the time-distributed reflections adding together and representing the vertical distribution of surfaces within the footprint. Larger footprints (10 - 100m in diameter) can support numerous target surfaces and thus provide the potential for producing complex return pulses. Interpreting the return pulse from laser altimeters has evolved from simple timing between thresholds, range-walk corrections, constant-fraction discriminators, and multi-stop time interval units to actual recording of the time varying return pulse intensity - the return waveform. Interpreting the waveform can be as simple as digitally thresholding the return pulse, calculating a centroid, to fitting one or more gaussian pulse-shapes to the signal. What we present here is a new technique for using the raw recorded return pulse as a raw observation to detect centimeter-level vertical topographic change using large footprint airborne and spaceborne laser altimetry. We use the correlation of waveforms from coincident footprints as an indication of the similarity in structure of the waveforms from epoch to epoch, and assume that low correlation is an indicator of vertical structure or elevation change. Thus, using vertically and horizontally geolocated waveforms as raw observables (i.e., waveforms tied to a common reference ellipsoid), we assess whether epoch-to-epoch vertical ground motion results in a decrease in the correlation of coincident waveforms over time, and whether this can be used to quantify the magnitude of the deformation. Results of computer models and an example over an area of eroded beachfront will be presented.

Blair, J. Bryan↗

NASA NTRS2017

Looking at A 0535+26 at Low Luminosities with NuSTAR

We report on two NuSTAR observations of the high-mass X-ray binary A 0535+26 taken toward the end of its normal 2015 outburst at very low 3-50 keV luminosities of approximately 1.4 times x 10 (sup 36) ergs per second and approximately 5 times x 10 (sup 35) ergs per second, which are complemented by nine Swift observations. The data clearly confirm indications seen in earlier data that the source's spectral shape softens as it becomes fainter. The smooth exponential rollover at high energies seen in the first observation evolves to a much more abrupt steepening of the spectrum at 20-30 keV. The continuum evolution can be nicely described with emission from a magnetized accretion column, modeled using the compmag model modified by an additional Gaussian emission component for the fainter observation. Between the two observations, the optical depth changes from 0.75 plus or minus 0.04 to 0.56 plus 0.01 (sup) minus 0.04 (sub), the electron temperature remains constant, and there is an indication that the column decreases in radius. Since the energy-resolved pulse profiles remain virtually unchanged in shape between the two observations, the emission properties of the accretion column reflect the same accretion regime. This conclusion is also confirmed by our result that the energy of the cyclotron resonant scattering feature (CRSF) at approximately 45 keV is independent of the luminosity, implying that the magnetic field in the region in which the observed radiation is produced is the same in both observations. Finally, we also constrain the evolution of the continuum parameters with the rotational phase of the neutron star. The width of the CRSF could only be constrained for the brighter observation. Based on Monte Carlo simulations of CRSF formation in single accretion columns, its pulse phase dependence supports a simplified fan beam emission pattern. The evolution of the CRSF width is very similar to that of the CRSF depth, which is, however, in disagreement with expectations.

Ballhausen, Ralf↗

NASA NTRS2013

Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results

We present the final nine-year maps and basic results from the Wilkinson Microwave Anisotropy Probe (WMAP) mission. The full nine-year analysis of the time-ordered data provides updated characterizations and calibrations of the experiment. We also provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate cosmic microwave background (CMB) anisotropy from foreground emission, and both types of signals are analyzed in detail.We provide new point source catalogs as well as new diffuse and point source foreground masks. An updated template-removal process is used for cosmological analysis; new foreground fits are performed, and new foreground reduced are presented.We nowimplement an optimal C(exp -1)1 weighting to compute the temperature angular power spectrum. The WMAP mission has resulted in a highly constrained Lambda-CDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that big bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (N(sub eff) = 3.84 +/- 0.40). The model fit also implies that the age of the universe is (sub 0) = 13.772 +/- 0.059 Gyr, and the fit Hubble constant is H(sub 0) = 69.32 +/- 0.80 km/s/ Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity, reported earlier by the WMAP team, now has high statistical significance (n(sub s) = 0.9608+/-0.0080); and the universe is close to flat/Euclidean (Omega = −0.0027+0.0039/−0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter Lambda-CDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further.With no significant anomalies and an adequate goodness of fit, the inflationary flat Lambda-CDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

dark matter↗

NASA NTRS2010

Mineralogy and Thermal Properties of V-Type Asteroid 956 Elisa: Evidence for Diogenitic Material from the Spitzer IRS (5-35 Micrometers) Spectrum

We present the thermal infrared (5-35 micrometer) spectrum of 956 Elisa as measured by the Spitzer Infrared Spectrograph ("IRS"; Houck,1.R. et .11. [20041. Astrophys, 1. SuppL 154, 18-24) together with new ground-based lightcurve data and near-IR spectra. From the visible lightcurve photometry, we determine a rotation period of 16.494 +/- 0.001 h, identify the rotational phase of the Spitzer observations, and estimate the visible absolute magnitude (Hv) at that rotational phase to be 12.58 +/- 0.04. From radiometric analysis of the thermal flux spectrum, we find that at the time of observation 956 Elisa had a projected radius of 5.3 +/- 0.4 km with a visible albedo pv = 0.142+/- 0.022, significantly lower than that of the prototype V-type asteroid, 4 Vesta. (This corresponds to a radius of 5.2 +/- 0.4 km at lightcurve mean.) Analysis with the standard thermal model (STM) results in a sub-solar temperature of 292.3 +/- 2.8 K and beaming parameter eta = 1.16 +/- 0.05. Thermophysical modeling places a lower limit of 20 J m(exp -2)K(exp -1)s(exp -1/2) on the thermal inertia of the asteroid's surface layer (if the surface is very smooth) but more likely values fall between 30 and 150 J m(exp -2)K(exp -1)s(exp -1/2) depending on the sense of rotation. The emissivity spectrum, calculated by dividing the measured thermal flux spectrum by the modeled thermal continuum, exhibits mineralogically interpretable spectral features within the 9-12 micrometer reststrahlen band, the 15-16.5 micrometer Si-O-Si stretching region, and the 16-25 micrometer reststrahlen region that are consistent with pyroxene of diogenitic composition: extant diogenitic pyroxenes fall within the narrow compositional range W0(sub 2+/-1)En(sub 74+/-2)Fs(sub 24+/-1). Spectral deconvolution of the 9-12 micrometer reststrahlen features indicates that up to approximately 20% olivine may also be present, suggesting an olivine-diogenite-like mineralogy. The mid-IR spectrum is inconsistent with non-cumulate eucrite as the major component on the surface of 956 Elisa, although cumulate eucrite material may be present at abundances lower than that of the diogenite component. Analysis of new near-IR spectra of 956 Elisa with the Modified Gaussian Model (MGM; Sunshine, J,M., Pieters, C.M., Pratt, S.F.[1990]. J. Geophys. Res, 95 (May), 6955-6966) results in two pyroxene compositions: 75% magnesian low-Ca pyroxene and 25% high-Ca pyroxene. High-Ca pyroxene is not evident in the mid-IR data, but may belong to a component that is underrepresented in the mid-IR spectrum either because of its spatial distribution on the asteroid or because of its particle size. High-Ca pyroxenes that occur as exsolution lamellae may also be more evident spectrally in the NIR than in the mid-IR. In any case, we find that the mid-IR spectrum of 956 Elisa is dominated by emission from material of diogenite- like composition, which has very rarely been observed among asteroids.

Lim, Lucy F.↗

NASA NTRS1988

Approximation theory for LQG (Linear-Quadratic-Gaussian) optimal control of flexible structures

An approximation theory is presented for the LQG (Linear-Quadratic-Gaussian) optimal control problem for flexible structures whose distributed models have bounded input and output operators. The main purpose of the theory is to guide the design of finite dimensional compensators that approximate closely the optimal compensator. The optimal LQG problem separates into an optimal linear-quadratic regulator problem and an optimal state estimation problem. The solution of the former problem lies in the solution to an infinite dimensional Riccati operator equation. The approximation scheme approximates the infinite dimensional LQG problem with a sequence of finite dimensional LQG problems defined for a sequence of finite dimensional, usually finite element or modal, approximations of the distributed model of the structure. Two Riccati matrix equations determine the solution to each approximating problem. The finite dimensional equations for numerical approximation are developed, including formulas for converting matrix control and estimator gains to their functional representation to allow comparison of gains based on different orders of approximation. Convergence of the approximating control and estimator gains and of the corresponding finite dimensional compensators is studied. Also, convergence and stability of the closed-loop systems produced with the finite dimensional compensators are discussed. The convergence theory is based on the convergence of the solutions of the finite dimensional Riccati equations to the solutions of the infinite dimensional Riccati equations. A numerical example with a flexible beam, a rotating rigid body, and a lumped mass is given.

Gibson, J. S.↗

NASA NTRS2019

Star Tracker Accuracy Improvement and Optimization for Attitude Measurement in Three-Axis

High precision attitude measurement systems obviate the need for the beacon from the receiver making it possible for the spacecraft to beam a laser communications signal to a ground station without the ground station advertising its location. The research presented targets new detection and estimation methods to improve the accuracy in locating stars on a focal plane detector, and an understanding of the effects of changes in the optics design parameters and aberration, including defocus, on the navigation solution itself. This understanding can lead to an optimization of the attitude solution with respect to those optics realm parameter changes. The methodology discussed includes the development of a model of a current star tracker system. Using this model, multiple algorithms are implemented, including a multi-hypothesis method (MHT), to detect and estimate the position of the stars on the focal plane detector. It will be shown that using the MHT for detection and estimation, a greater accuracy can be found for each star estimation from more traditional detection and estimation algorithms. The approach then uses the model to develop statistics of the star tracker and the attitude estimation outputs to understand the accuracy, or variance, of the system's attitude solution. This solution is repeated for a range of defocus aberration, and a lower limit to the variance of the attitude solution is shown. A Cramer Rao lower bound solution is derived for the star tracker system and the results are compared to the Monte Carlo analysis from the model and shown to correlate very well. The approach uses a star image not as a Gaussian spot on the focal plane as done in previous work, and use of an image that includes the effects of aberrations of the optic system, and the effects of under-sampling and noise from the focal plane detector as well. Analysis includes exploring a star tracker's accuracy improvement through the combination of focus error and under-sampling effects alone, possibly contradicting conventional wisdom and approaches.

Lichter, Michael J.↗

NASA NTRS2004

Determining Sizes of Particles in a Flow from DPIV Data

A proposed method of measuring the size of particles entrained in a flow of a liquid or gas would involve utilization of data from digital particle-image velocimetry (DPIV) of the flow. That is to say, with proper design and operation of a DPIV system, the DPIV data could be processed according to the proposed method to obtain particle sizes in addition to particle velocities. As an additional benefit, one could then compute the mass flux of the entrained particles from the particle sizes and velocities. As in DPIV as practiced heretofore, a pulsed laser beam would be formed into a thin sheet to illuminate a plane of interest in a flow field and the illuminated plane would be observed by means of a charge-coupled device (CCD) camera aimed along a line perpendicular to the illuminated plane. Unlike in DPIV as practiced heretofore, care would be taken to polarize the laser beam so that its electric field would lie in the illuminated plane, for the reason explained in the next paragraph. The proposed method applies, more specifically, to transparent or semitransparent spherical particles that have an index of refraction different from that of the fluid in which they are entrained. The method is based on the established Mie theory, which describes the scattering of light by diffraction, refraction, and specular reflection of light by such particles. In the case of a particle illuminated by polarized light and observed in the arrangement described in the preceding paragraph, the Mie theory shows that the image of the particle on the focal plane of the CCD camera includes two glare spots: one attributable to light reflected toward the camera and one attributable to light refracted toward the camera. The distance between the glare spots is a known function of the size of the particle, the indices of refraction of the particle material, and design parameters of the camera optics. Hence, the size of a particle can be determined from the distance between the glare spots. The proposed method would be implemented in an algorithm that would automatically identify, and measure the distance between, the glare spots for each particle for which a suitable image has been captured in a DPIV image frame. The algorithm (see figure) would begin with thresholding of data from the entire image frame to reduce noise, thereby facilitating discrimination of particle images from the background and aiding in the separation of overlapping particles. It is important not to pick a threshold level so high that the light intensity between a given pair of glare spots does not fall below the threshold value, leaving the glare spots disconnected. The image would then be scanned in a sequence of rows and columns of pixels to identify groups of adjacent pixels that contain nonzero brightnesses and that are surrounded by pixels of zero brightness. Each such group would be assumed to constitute the image of one particle. Each such group would be further analyzed to determine whether the image was saturated; saturated particle images must be rejected because the locations of glare spots in saturated images cannot accurately be determined. Within each unsaturated particle image, the centroids (deemed to be the locations) of the glare spots would be determined by means of gradients of brightness distributions and three-point horizontal and three-point vertical Gaussian estimates based on the brightness values of the brightest pixels and the pixels adjacent to them. If the brightness of a given particle image contained only one peak, then it would be assumed that a second glare spot did not exist and that image would be rejected.

Wernet, M. P.↗