Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian process fitting”

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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

A simulation of air pollution model parameter estimation using data from a ground-based LIDAR remote sensor

One way to obtain estimates of the unknown parameters in a pollution dispersion model is to compare the model predictions with remotely sensed air quality data. A ground-based LIDAR sensor provides relative pollution concentration measurements as a function of space and time. The measured sensor data are compared with the dispersion model output through a numerical estimation procedure to yield parameter estimates which best fit the data. This overall process is tested in a computer simulation to study the effects of various measurement strategies. Such a simulation is useful prior to a field measurement exercise to maximize the information content in the collected data. Parametric studies of simulated data matched to a Gaussian plume dispersion model indicate the trade offs available between estimation accuracy and data acquisition strategy.

Kibler, J. F.↗
NASA NTRS2013

Herschel Observations of the Debris Disc Around HIP 92043

Context: Typical debris discs are composed of particles ranging from several micron sized dust grains to kilometer-sized asteroidal bodies, and their infrared emission peaks at wavelengths 60-100 microns. Recent Herschel DUNES (Dust Around Nearby Stars) observations have identified several debris discs around nearby Sun-like stars (F, G and K spectral type) with significant excess emission only at 160 microns. Aims: We observed HIP 92043 (110 Her, HD 173667) at far-infrared and sub-millimetre wavelengths with Herschel PACS (Photodetector Array Camera and Spectrometer) and SPIRE (Spectral and Photometric Imaging Receiver).Identification of the presence of excess emission from HIP 92043 and the origin and physical properties of any excess was undertaken through analysis of its spectral energy distribution (SED) and the PACS images. Methods: The PACS and SPIRE images were produced using the HIPE (Herschel Interactive Processing Environment) photProject map maker routine. Fluxes were measured using aperture photometry. A stellar photosphere model was scaled to optical and near infrared photometry and subtracted from the far-infared and sub-mm fluxes to determine the presence of excess emission. Source radial profiles were fitted using a 2D Gaussian and compared to a PSF (Point Spread Function) model based on Herschel observations of alpha Boo to check for extended emission. Results: Clear excess emission from HIP 92043 was observed at 70 and 100 microns. Marginal excess was observed at 160 and 250 microns.Analysis of the images reveals that the source is extended at 160 microns. A fit to the source SED is inconsistent with a photosphere and single temperature black body. Conclusions: The excess emission from HIP 92043 is consistent with the presence of an unresolved circumstellar debris disc at 70and 100 microns, with low probability of background contamination. The extended 160 micron emission may be interpreted as an additional cold component to the debris disc or as the result of background contamination along the line of sight. The nature of the 160 micron excess cannot be determined absolutely from the available data, but we favour a debris disc interpretation, drawing parallels with previously identified cold disc sources in the DUNES sample.

debris discs↗
NASA NTRS2004

Synoptic Study of an X-Ray Nova in Outburst (ADP 2000)

In the completed work listed above, we present an X--ray spectral and timing analysis of 4131543-47 during its 2002 outburst based on 49 pointed observations obtained using the Rossi X-ray Timing Explorer (RXTE). The outburst reached a peak intensity of 4.2 Crab in the 2-12 keV band and declined by a factor of 32 throughout the month-long observation. A 21.9+/-0.6 mJy radio flare was detected at 1026.75 MHz two days before the X-ray maximum; the radio source was also detected late in the outburst, after the X-ray source entered the hard state. The X-ray light curve exhibits the classic shape of a rapid rise and an exponential decay. The spectrum is soft and dominated by emission from the accretion disk. The continuum is fit with a multicolor disk blackbody (kT-max = 1.04 keV) and a power-law (Gamma -2.7). Midway through the decay phase, a strong low-frequency quasi-periodic oscillation (QPO; nu = 7.3-8.1 Hz) was present for several days. The spectra feature a broad Fe K line that is asymmetric, suggesting that the line is due to relativistic broadening rather than Gomptonization. Relativistic Laor models provide much better fits to the line than non-relativistic Gaussian models, particularly near the beginning and end of our observations. The line fits yield estimates for the inner disk radius that are within 6 GM/c 2; this result and additional evidence indicates that this black hole may have a non-zero angular momentum. In current work, we are in the process of analyzing a huge set of RXTE data (-100 pointed observations) for H1743-322, which was in eruption throughout most of 2003. The source s behavior is remarkable. It displays all of the canonical X-ray states, sometimes in unusual ways. In this short report, we feature a single important result, namely, the detection of a pair of high-frequency QPOs. The bright flaring behavior is of special interest. During the flares the shape of the power spectrum and the presence of low-frequency QPOs indicate that the source is in the steep power-law (SPL or very high) state. It is in this state that the high-frequency QPOs, which are believed to emanate from the very inner accretion disk, were oberved. Specifically, we observe a 240 Hz QPO at a significance of 4.7sigma and a 162 Hz QPO at 3.7sigma. The higher frequency QPO was discovered by Homan et al. in an independent data set; they also found marginal evidence (2.5sigma)

McClintock, Jeffrey↗
NASA NTRS2000

The High Mosaicity Illusion: Revealing the True Physical Characteristics of Macromolecular Crystals

An experimental system and software have been developed for simultaneously measuring the diffraction resolution and mosaic spread of macromolecular crystals. Hundreds of reflection profiles over a wide resolution range were rapidly measured by using a charge coupled device (CCD) area detector in combination with superfine phi slicing data collection. The contributions of the X-ray beam to the reflection widths were minimized by using a highly-parallel, highly-monochromatic synchrotron source. These contributions and Lorentz effects were evaluated and deconvoluted from the recorded data. Data collection and processing is described. From one degree of superfine phi slice data collected on a crystal of manganese superoxide dismutase the mosaicity of 261 reflections were measured. The average mosaicity was 0.0101 degrees (0.0035) at the full-width-at-half-maximum (FWHM) and ranged from 0.0011 degrees to 0.0188 degrees. Each reflection profile was individually fit with two gaussian profiles with the first gaussian contributing 55% and the second contributing 35% of the reflection. On average, the mosaicity of the first gaussian was 0.0054 degrees (0.0015) and the second was 0.0061 degrees (0.0023). The mosaicity of the crystal was anisotropic with fh, f k, and fl values of 0.0068 degrees, 0.0140 degrees and 0.0046 degrees, respectively at the FWHM. The anisotropic mosaicity analysis indicates that the crystal is the most perfect in the I direction which corresponds to the favored growth direction of the crystal.

Bellamy, Henry↗
NASA NTRS1977

A quasi-linear control theory analysis of timesharing skills

The compliance of the human ankle joint is measured by applying 0 to 50 Hz band-limited gaussian random torques to the foot of a seated human subject. These torques rotate the foot in a plantar-dorsal direction about a horizontal axis at a medial moleolus of the ankle. The applied torques and the resulting angular rotation of the foot are measured, digitized and recorded for off-line processing. Using such a best-fit, second-order model, the effective moment of inertia of the ankle joint, the angular viscosity and the stiffness are calculated. The ankle joint stiffness is shown to be a linear function of the level of tonic muscle contraction, increasing at a rate of 20 to 40 Nm/rad/Kg.m. of active torque. In terms of the muscle physiology, the more muscle fibers that are active, the greater the muscle stiffness. Joint viscosity also increases with activation. Joint stiffness is also a linear function of the joint angle, increasing at a rate of about 0.7 to 1.1 Nm/rad/deg from plantar flexion to dorsiflexion rotation.

Agarwal, G. C.↗
NASA NTRS1991

A hydrodynamic approach to cosmology - Texture-seeded cold dark matter and hot dark matter cosmogonies

Hydrodynamical simulations of galaxy formation in a texture-seeded cosmology are presented, with attention given to Omega = 1 galaxies dominated by both hot dark matter (HDM) and cold dark matter (CDM). The simulations include both gravitational and hydrodynamical physics with a detailed treatment of collisional and radiative thermal processes, and use a cooling criterion to estimate galaxy formation. Background radiation fields and Zel'dovich-Sunyaev fluctuations are explicitly computed. The derived galaxy mass function is well fitted by the observed Schechter luminosity function for a baryonic M/L of 3 and total M/L of 60 in galaxies. In both HDM and CDM texture scenarios, the 'galaxies' and 'clusters' are significantly more strongly correlated than the dark matter due to physical bias processes. The slope of the correlation function in both cases is consistent with observations. In contrast to Gaussian models, peaks in the dark matter density distributrion are less correlated than average.

Cen, R. Y.↗
NASA NTRS2022

Modeling the East-West Asymmetry of Energetic Particle Fluence in Large Solar Energetic Particle Events Using the iPATH Model

It has been noted that in large solar energetic particle (SEP) events, the peak intensities show an East-West asymmetry with respect to the source flare locations. Using the 2D improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, we investigate the origin of this longitudinal trend. We consider multiple cases with different solar wind speeds and eruption speeds of the coronal mass ejections (CMEs) and fit the longitudinal distributions of time-averaged fluence by symmetric/asymmetric Gaussian functions with three time intervals of 8, 24 and 48 hr after the flare onset time respectively. The simulation results are compared with a statistical study of three-spacecraft events. We suggest that the East-West asymmetry of SEP fluence and peak intensity can be primarily caused the combined effect of an extended shock acceleration process and the evolution of magnetic field connection to the shock front. Our simulations show that the solar wind speed and the CME speed are important factors determining the East-West fluence asymmetry.

Zheyi Ding↗
NASA NTRS2021

The Asymmetric Inner Disk of the Herbig Ae Star HD 163296 in the Eyes of VLTI/MATISSE: Evidence for a Vortex?

Context.A complex environment exists in the inner few astronomical units of planet-forming disks. High-angular-resolution observa-tions play a key role in our understanding of the disk structure and the dynamical processes at work.Aims.In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-massstar HD 163296 from early VLTI/MATISSE observations taken in theL- andN-bands. We put special emphasis on the detection ofpotential disk asymmetries.Methods.We use simple geometric models to fit the interferometric visibilities and closure phases. Our models include a smoothedring, a flat disk with an inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetriesas well. We also perform numerical hydrodynamical simulations of the inner edge of the disk.Results.Our modeling reveals a significant brightness asymmetry in theL-band disk emission. The brightness maximum of the asym-metry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio inthe azimuthal variation is3.5±0.2. Comparing our result on the location of the asymmetry with other interferometric measurements,we confirm that the morphology of ther<0.3au disk region is time-variable. We propose that this asymmetric structure, located in ornear the inner rim of the dusty disk, orbits the star. To find the physical origin of the asymmetry, we tested a hypothesis where a vortexis created by Rossby wave instability, and we find that a unique large-scale vortex may be compatible with our data. The half-lightradius of theL-band-emitting region is0.33±0.01au, the inclination is52◦+5◦−7◦, and the position angle is143◦±3◦. Our models predictthat a non-negligible fraction of theL-band disk emission originates inside the dust sublimation radius forμm-sized grains. Refractorygrains or large (&10μm-sized) grains could be the origin of this emission.N-band observations may also support a lack of smallsilicate grains in the innermost disk (r.0.6au), in agreement with our findings fromL-band data.

J Varga↗
NASA NTRS1989

Laser transit anemometer software development program

Algorithms were developed for the extraction of two components of mean velocity, standard deviation, and the associated correlation coefficient from laser transit anemometry (LTA) data ensembles. The solution method is based on an assumed two-dimensional Gaussian probability density function (PDF) model of the flow field under investigation. The procedure consists of transforming the data ensembles from the data acquisition domain (consisting of time and angle information) to the velocity space domain (consisting of velocity component information). The mean velocity results are obtained from the data ensemble centroid. Through a least squares fitting of the transformed data to an ellipse representing the intersection of a plane with the PDF, the standard deviations and correlation coefficient are obtained. A data set simulation method is presented to test the data reduction process. Results of using the simulation system with a limited test matrix of input values is also given.

Abbiss, John B.↗
NASA NTRS2018

Learning Model Structural Uncertainty with Gaussian Processes

The advent of commercially available quantum computers has marked the beginning of quantum computing as a reality. Both quantum gate and annealing computers have been released by major computer hardware companies. In this work, the D-Wave 2XTM quantum annealing computer housed at the NASA Advanced Systems computational facility is investigated to accelerate Machine Learning (ML) for image registration. NASA collects large amounts of images over the globe remotely using space-based monitoring. Images of a fixed areas of the land surface are taken over time. Due to the orbit of the sensors, the viewing angles deviate slightly, and it is necessary to align or register the images precisely to create image time series over the land surface. Unaligned images can lead to substantial analysis errors. These time-series are then used in modeling Earth Systems models such as hydrological, weather, and carbon monitoring models. In this work, we consider the Moderate Resolution Image Spectrometer (MODIS) data collected by the NASA's terra satellite. Artificial Neural Networks (ANNs) is a natural fit for ML modelling of images. Several successes have been reported using machine learning related to image processing. We investigate the use of ML to register MODIS images. ANNs are investigated in combination with a Restricted Boltzmann Machines (RBM) as an auto-encoder. We will present results showing the accuracy and efficiency of this approach.The D-Wave 2XTM quantum annealer samples the ground-state wave-function of a spin-Ising systems with quadratic interactions between qubits and a Chimera connectivity. The system sits in a ~15 mK thermal bath. One can think of the system as being placed in the ground state initially and subject to thermal excitations governed by Boltzmann statistics. If this is assumed true, one can use the statistics from the D-Wave 2XTM to train RBMs. Generating statistics for training Boltzmann machines is an NP-hard problem and constitutes the largest compute cost. We investigate the use of the D-Wave 2XTM to accelerate the training of the RBMs in our ANNs and report on the results.

Kouatchou, Jules↗
NASA NTRS1999

Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data

Estimates of solar normal mode frequencies from helioseismic observations can be improved by using Multitaper Spectral Analysis (MTSA) to estimate spectra from the time series, then using wavelet denoising of the log spectra. MTSA leads to a power spectrum estimate with reduced variance and better leakage properties than the conventional periodogram. Under the assumption of stationarity and mild regularity conditions, the log multitaper spectrum has a statistical distribution that is approximately Gaussian, so wavelet denoising is asymptotically an optimal method to reduce the noise in the estimated spectra. We find that a single m-upsilon spectrum benefits greatly from MTSA followed by wavelet denoising, and that wavelet denoising by itself can be used to improve m-averaged spectra. We compare estimates using two different 5-taper estimates (Stepian and sine tapers) and the periodogram estimate, for GONG time series at selected angular degrees l. We compare those three spectra with and without wavelet-denoising, both visually, and in terms of the mode parameters estimated from the pre-processed spectra using the GONG peak-fitting algorithm. The two multitaper estimates give equivalent results. The number of modes fitted well by the GONG algorithm is 20% to 60% larger (depending on l and the temporal frequency) when applied to the multitaper estimates than when applied to the periodogram. The estimated mode parameters (frequency, amplitude and width) are comparable for the three power spectrum estimates, except for modes with very small mode widths (a few frequency bins), where the multitaper spectra broadened the modest compared with the periodogram. We tested the influence of the number of tapers used and found that narrow modes at low n values are broadened to the extent that they can no longer be fit if the number of tapers is too large. For helioseismic time series of this length and temporal resolution, the optimal number of tapers is less than 10.

Komm, R. W.↗
NASA NTRS1991

Preliminary results from the COBE differential microwave radiometers - Large angular scale isotropy of the cosmic microwave background

Preliminary but precise micowave maps are presented of the sky, and thus of the early universe, derived as the first results from the Differential Microwave Radiometers experiment aboard COBE. The dipole anisotropy attributed to the motion of the solar system with respect to the CMB reference frame shows strongly in all six sky maps and is consistent with a Doppler-shifted thermal spectrum. The best-fitted dipole has amplitude 3.3 + or - 0.2 mK in the direction (alpha, delta) = 11.2 h + or - 0.2 h, -7 deg + or - 2 deg (J2000) or (l,b) = 265 deg + or - 2 deg, 48 deg + or - 2 deg. There is no clear evidence in the maps for any other large angular-scale feature. Limits on Delta T/T0 of 3 x 10 to the -5th (T0 = 2.735 K), 4 x 10 to the -5th, and 4 x 10 to the -5th are found for the rms quadrupole amplitude, monochromatic fluctuations, and Gaussian fluctuations, respectively. These measurements place the most severe constraints to date on many potential physical processes in the early universe.

Smoot, G. F.↗
NASA NTRS2015

Intensity Conserving Spectral Fitting

The detailed shapes of spectral line profiles provide valuable information about the emitting plasma, especially when the plasma contains an unresolved mixture of velocities, temperatures, and densities. As a result of finite spectral resolution, the intensity measured by a spectrometer is the average intensity across a wavelength bin of non-zero size. It is assigned to the wavelength position at the center of the bin. However, the actual intensity at that discrete position will be different if the profile is curved, as it invariably is. Standard fitting routines (spline, Gaussian, etc.) do not account for this difference, and this can result in significant errors when making sensitive measurements. Detection of asymmetries in solar coronal emission lines is one example. Removal of line blends is another. We have developed an iterative procedure that corrects for this effect. It can be used with any fitting function, but we employ a cubic spline in a new analysis routine called Intensity Conserving Spline Interpolation (ICSI). As the name implies, it conserves the observed intensity within each wavelength bin, which ordinary fits do not. Given the rapid convergence, speed of computation, and ease of use, we suggest that ICSI be made a standard component of the processing pipeline for spectroscopic data.

Klimchuk, J. A.↗
NASA NTRS2004

Characterizing Oscillatory Bursts in Single-Trial EEG Data

Oscillatory bursts in numerous bands ranging from low (theta) to high frequencies (e.g., gamma) undoubtedly play an important role in cortical dynamics. Largely because of the inadequacy of existing analytic techniques. however, oscillatory bursts and their role in cortical processing remains poorly understood. To study oscillatory bursts effectively one must be able to isolate them and characterize them in the single trial. We describe a series of straightforward analysis techniques that produce useful indices of burst characteristics. First, stimulus-evoked responses are estimated using Differentially Variable Component Analysis (dVCA), and are subtracted from the single-trial. The single-trial characteristics of the evoked responses are stored to identify possible correlations with burst activity. Time-frequency (T-F), or wavelet, analyses are then applied to the single trial residuals. While T-F plots have been used in recent studies to identify and isolate bursts, we go further by fitting each burst in the T-F plot with a two-dimensional Gaussian. This provides a set of burst characteristics, such as, center time. burst duration, center frequency. frequency dispersion. and amplitude, all of which contribute to the accurate characterization of the individual burst. The burst phase can also be estimated. Burst characteristics can be quantified with several standard techniques (e.g.. histogramming and clustering), as well as Bayesian techniques (e.g., blocking) to allow a more parametric description analysis of the characteristics of oscillatory bursts, and the relationships of specific parameters to cortical excitability and stimulus integration.

Knuth, K. H.↗
NASA NTRS1997

Quantitative Simulations of MST Visual Receptive Field Properties Using a Template Model of Heading Estimation

We previously developed a template model of primate visual self-motion processing that proposes a specific set of projections from MT-like local motion sensors onto output units to estimate heading and relative depth from optic flow. At the time, we showed that that the model output units have emergent properties similar to those of MSTd neurons, although there was little physiological evidence to test the model more directly. We have now systematically examined the properties of the model using stimulus paradigms used by others in recent single-unit studies of MST: 1) 2-D bell-shaped heading tuning. Most MSTd neurons and model output units show bell-shaped heading tuning. Furthermore, we found that most model output units and the finely-sampled example neuron in the Duffy-Wurtz study are well fit by a 2D gaussian (sigma approx. 35deg, r approx. 0.9). The bandwidth of model and real units can explain why Lappe et al. found apparent sigmoidal tuning using a restricted range of stimuli (+/-40deg). 2) Spiral Tuning and Invariance. Graziano et al. found that many MST neurons appear tuned to a specific combination of rotation and expansion (spiral flow) and that this tuning changes little for approx. 10deg shifts in stimulus placement. Simulations of model output units under the same conditions quantitatively replicate this result. We conclude that a template architecture may underlie MT inputs to MST.

Stone, Leland S.↗
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 NTRS1998

Pursuit Latency for Chromatic Targets

The temporal dynamics of eye movement response to a change in direction of stimulus motion has been used to compare the processing speeds of different types of stimuli (Mulligan, ARVO '97). In this study, the pursuit response to colored targets was measured to test the hypothesis that the slow response of the chromatic system (as measured using traditional temporal sensitivity measures such as contrast sensitivity) results in increased eye movement latencies. Subjects viewed a small (0.4 deg) Gaussian spot which moved downward at a speed of 6.6 deg/sec. At a variable time during the trajectory, the dot's direction of motion changed by 30 degrees, either to the right or left. Subjects were instructed to pursue the spot. Eye movements were measured using a video ophthalmoscope with an angular resolution of approximately 1 arc min and a temporal sampling rate of 60 Hz. Stimuli were modulated in chrominance for a variety of hue directions, combined with a range of small luminance increments and decrements, to insure that some of the stimuli fell in the subjects' equiluminance planes. The smooth portions of the resulting eye movement traces were fit by convolving the stimulus velocity with an exponential having variable onset latency, time constant and amplitude. Smooth eye movements with few saccades were observed for all stimuli. Pursuit responses to stimuli having a significant luminance component are well-fit by exponentials having latencies and time constants on the order of 100 msec. Increases in pursuit response latency on the order of 100-200 msec are observed in response to certain stimuli, which occur in pairs of complementary hues, corresponding to the intersection of the stimulus section with the subjects' equiluminant plane. Smooth eye movements can be made in response to purely chromatic stimuli, but are slower than responses to stimuli with a luminance component.

Mulligan, Jeffrey B.↗
NASA NTRS2009

Hard X-Ray Flare Source Sizes Measured with the Ramaty High Energy Solar Spectroscopic Imager

Ramaty High Energy Solar Spectroscopic Imager (RHESSI) observations of 18 double hard X-ray sources seen at energies above 25 keV are analyzed to determine the spatial extent of the most compact structures evident in each case. The following four image reconstruction algorithms were used: Clean, Pixon, and two routines using visibilities maximum entropy and forward fit (VFF). All have been adapted for this study to optimize their ability to provide reliable estimates of the sizes of the more compact sources. The source fluxes, sizes, and morphologies obtained with each method are cross-correlated and the similarities and disagreements are discussed. The full width at half-maximum (FWHM) of the major axes of the sources with assumed elliptical Gaussian shapes are generally well correlated between the four image reconstruction routines and vary between the RHESSI resolution limit of approximately 2" up to approximately 20" with most below 10". The FWHM of the minor axes are generally at or just above the RHESSI limit and hence should be considered as unresolved in most cases. The orientation angles of the elliptical sources are also well correlated. These results suggest that the elongated sources are generally aligned along a flare ribbon with the minor axis perpendicular to the ribbon. This is verified for the one flare in our list with coincident Transition Region and Coronal Explorer (TRACE) images. There is evidence for significant extra flux in many of the flares in addition to the two identified compact sources, thus rendering the VFF assumption of just two Gaussians inadequate. A more realistic approximation in many cases would be of two line sources with unresolved widths. Recommendations are given for optimizing the RHESSI imaging reconstruction process to ensure that the finest possible details of the source morphology become evident and that reliable estimates can be made of the source dimensions.

Dennis, Brian R.↗