Advanced Characterization of Spatial Aspects of Image System Blur
Presentation of SDRD research at an ASU seminar.
SEARCH · Engineering Papers
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.
Presentation of SDRD research at an ASU seminar.
Computer program, which predicts effects of optical aberrations on transmitters and receivers, includes effects of Gaussian pupil functions and utilizes algorithm that permits specification of number and location of output points for computed spread function results.
Polarized fluorescence microscopy is a valuable tool for measuring molecular orientations in biological samples, but techniques for recovering three-dimensional orientations and positions of fluorescent ensembles are limited. We report a polarized dual-view light-sheet system for determining the diffraction-limited three-dimensional distribution of the orientations and positions of ensembles of fluorescent dipoles that label biological structures. We share a set of visualization, histogram, and profiling tools for interpreting these positions and orientations. We model the distributions based on the polarization-dependent efficiency of excitation and detection of emitted fluorescence, using coarse-grained representations we call orientation distribution functions (ODFs). We apply ODFs to create physics-informed models of image formation with spatio-angular point-spread and transfer functions. We use theory and experiment to conclude that light-sheet tilting is a necessary part of our design for recovering all three-dimensional orientations. We use our system to extend known two-dimensional results to three dimensions in FM1-43-labeled giant unilamellar vesicles, fast-scarlet-labeled cellulose in xylem cells, and phalloidin-labeled actin in U2OS cells. Additionally, we observe phalloidin-labeled actin in mouse fibroblasts grown on grids of labeled nanowires and identify correlations between local actin alignment and global cell-scale orientation, indicating cellular coordination across length scales.
Scanning slit across point or line spread function in determination of radial irradiance distribution
Abstract Background To facilitate in silico studies that investigate digital mammography (DM) and breast tomosynthesis (DBT), models replicating the variety in imaging performance of the DM and DBT systems, observed across manufacturers are needed. Purpose The main purpose of this work is to develop generic physics models for direct and indirect detector technology used in commercially available systems, with the goal of making them available open source to manufacturers to further tweak and develop the exact in silico replicas of their systems. Methods We recently reported on an in silico version of the SIEMENS Mammomat Inspiration DM/DBT system using an open‐source GPU‐accelerated Monte Carlo x‐ray imaging simulation code (MC‐GPU). We build on the previous version of the MC‐GPU codes to mimic the imaging performances of two other Food and Drug Administration (FDA)‐approved DM/DBT systems, such as Hologic Selenia Dimensions (HSD) and the General Electric Senographe Pristina (GSP) systems. In this work, we developed a hybrid technique to model the optical spread and signal crosstalk observed in the GSP and HSD systems. MC simulations are used to track each x‐ray photon till its first interaction within the x‐ray detector. On the other hand, the signal spread in the x‐ray detectors is modeled using previously developed analytical equations. This approach allows us to preserve the modeling accuracy offered by MC methods in the patient body, while speeding up secondary carrier transport (either electron–hole pairs or optical photons) using analytical equations in the detector. The analytical optical spread model for the indirect detector includes the depth‐dependent spread and collection of optical photons and relies on a pre‐computed set of point response functions that describe the optical spread as a function of depth. To understand the capabilities of the computational x‐ray detector models, we compared image quality metrics like modulation transfer function (MTF), normalized noise power spectrum (NNPS), and detective quantum efficiency (DQE), simulated with our models against measured data. Please note that the purpose of these comparisons with measured data would be to gauge if the model developed as part of this work could replicate commercially used direct and indirect technology in general and not to achieve perfect fits with measured data. Results We found that the simulated image quality metrics such as MTF, NNPS, and DQE were in reasonable agreement with experimental data. To demonstrate the imaging performance of the three DM/DBT systems, we integrated the detector models with the VICTRE pipeline and simulated DM images of a fatty breast model containing a spiculated mass and a calcium oxalate cluster. In general, we found that the images generated using the indirect model appeared more blurred with a different noise texture and contrast as compared to the systems with direct detectors. Conclusions We have presented computational models of three commercially available FDA‐approved DM/DBT systems, which implement both direct and indirect detector technology. The updated versions of the MC‐GPU codes that can be used to replicate three systems are available in open source format through GitHub.
The performance of optical systems is limited by imperfect optical surfaces that degrade the images below the level set by wave theoretical limits. The central irradiance functions are derived for slit and circular apertures with five distributions of wavefront errors and for a range of maximal wavefront deviations. For practical frequency of occurrence distributions of wavefront deviations, the point spread and the image energy integral functions are determined. Practical performances of optical systems are derived and performance limits discussed.
The spatial structure of the depth of rainfall from a stationary storm event is investigated by using point process techniques. Cells are assumed to be stationary and to be distributed in space either independently according to a Poisson process, or with clustering according to a Neyman-Scott scheme. Total storm rainfall at the centre of each cell is a random variable and rainfall is distributed around the centre in a way specified by a spread function that may incorporate random parameters. The mean, variance and covariance structure of the precipitation depth at a point are obtained for different spread functions. For exponentially distributed centre depth and a spread function having quadratically exponential decay, the total storm depth at any point in the field is shown to have a gamma distribution. The probability of zero rainfall at a point is investigated, as is the stochastic variability of model parameters from storm to storm. Data from the Upper Rio Guaire basin in Venezuela are used in illustration.
The circumgalactic medium (CGM) traced by O vi λλ1031,1037 doublet absorption has been found to concentrate along the projected major and minor axes of the host galaxies. This suggests that O vi traces accreting and outflowing gas, respectively, which are key components of the baryon cycle of galaxies. We investigate this further by examining the kinematics of 29 O vi absorbers associated with galaxies at 0.13<z{sub gal}<0.66 as a function of galaxy color, inclination, and azimuthal angle. Each galaxy was imaged with the Hubble Space Telescope (HST), and the absorption was detected in COS/HST spectra of nearby (D<200 kpc) background quasars. We use the pixel-velocity two-point correlation function to characterize the velocity spread of the absorbers, which is a method used previously for a sample of Mg ii absorber–galaxy pairs. The absorption velocity spread for O vi is more extended than Mg ii, which suggests that the two ions trace differing components of the CGM. Again, in contrast to Mg ii, the O vi absorption velocity spreads are similar regardless of galaxy color, inclination, and azimuthal angle. This indicates that the kinematics of the high-ionization gas is not strongly influenced by the current star formation activity in the galaxy. The kinematic homogeneity of O vi absorption and its tendency to be observed mainly along the projected galaxy major and minor axes is likely due to varying ionization conditions and gas densities about the galaxy. Gas in intermediate azimuthal angles may be ionized out of the O vi phase, possibly resulting in an azimuthal angle dependence of the distribution of gas in higher ionization states.
The equivalent widths of Mg II absorption in the circumgalactic medium (CGM) trace the global star formation rate up to z < 6, are larger for star-forming galaxies than passively evolving galaxies, and decrease with increasing distance from the galaxy. We delve further into the physics involved by investigating gas kinematics and cloud column density distributions as a function of galaxy color, redshift, and projected distance from the galaxy (normalized by galaxy virial radius, D/R{sub vir}). For 39 isolated galaxies at 0.3 < z{sub gal} < 1.0, we have detected Mg II absorption in high-resolution (Δv ≃ 6.6 km s{sup −1}) spectra of background quasars within a projected distance of 7 < D < 190 kpc. We characterize the absorption velocity spread using pixel-velocity two-point correlation functions. Velocity dispersions and cloud column densities for blue galaxies do not differ with redshift nor with D/R{sub vir}. This suggests that outflows continually replenish the CGM of blue galaxies with high velocity dispersion, large column density gas out to large distances. Conversely, absorption hosted by red galaxies evolves with redshift where the velocity dispersions (column densities) are smaller (larger) at z{sub gal} < 0.656. After taking into account larger possible velocities in more massive galaxies, we find that there is no difference in the velocity dispersions or column densities for absorption hosted by red galaxies with D/R{sub vir}. Thus, a lack of outflows in red galaxies causes the CGM to become more quiescent over time, with lower velocity dispersions and larger column densities toward lower z{sub gal}. The quenching of star formation appears to affect the CGM out to D/R{sub vir} = 0.75.
The Off-plane Grating Rocket Experiment (OGRE) is a soft X-ray grating spectrometer to be flown on a suborbital rocket. The payload is designed to obtain the highest-resolution soft X-ray spectrum of Capella to date with a resolution goal of R(lambda/delta lambda) > 2000 at select wavelengths in its 10 - 55 A bandpass of interest. The optical design of the spectrometer realizes a maximum resolution of R 5000 with all components performing optimally, in perfect alignment, and not considering in-flight pointing errors; however, performance errors, misalignments, and in-flight pointing errors work to degrade this performance. A comprehensive line spread function (LSF) error budget has been constructed for the spectrometer to identify errors contributing to the LSF, to determine how each affects the LSF, and to inform performance requirements and alignment tolerances for the spectrometer. In this document, the comprehensive LSF error budget for the OGRE spectrometer is presented and its implications are discussed.
The X-Ray Spectrometer (XRS) will be launched in February 2005 as part of the Astro-E2 mission. It will provide unprecedented throughput and resolving powers particularly at high energies. In this presentation we will describe the in-flight calibration program. The energy scale of the XRS is a complex, non-linear function of the noise and power conditions on the array. It will be calibrated empirically using the bright point sources, Capella and GX301-2. Ground calibration of the line spread function show it to be almost perfectly Gaussian. The in-flight calibration is designed to verify this using the energy scale targets. The effective area curve of the XRS contains discreet edge structure from the mirrors, the optical blocking filters, and the microcalorimeter HgTe absorbers. The effective area Calibration program will simultaneously measure these absorption edges and the global effective area properties using the relatively featureless sources 3C273 and Mrk421. Additional monitoring of any ice build-up on the filters will be conducted using observations of the supernova remnants N132D and E0102.
We present a study comparing O vi λλ1031, 1037 doublet absorption found toward group galaxy environments to that of isolated galaxies. The O vi absorption in the circumgalactic medium (CGM) of isolated galaxies has been studied previously by the “Multiphase Galaxy Halos” survey, where the kinematics and absorption properties of the CGM have been investigated. We extend these studies to group environments. We define a galaxy group as having two or more galaxies with a line-of-sight velocity difference of no more than 1000 km s{sup −1} and located within 350 kpc (projected) of a background quasar sightline. We identified a total of six galaxy groups associated with O vi absorption W {sub r} > 0.06 Å that have a median redshift of 〈z{sub gal}〉=0.1669 and a median impact parameter of 〈D〉=134.1 kpc. An additional 12 non-absorbing groups were identified with a median redshift of 〈z{sub gal}〉=0.2690 and a median impact parameter of 〈D〉=274.0 kpc. We find the average equivalent width to be smaller for group galaxies than for isolated galaxies (3σ). However, the covering fractions are consistent with both samples. We used the pixel-velocity two-point correlation function method and find that the velocity spread of O vi in the CGM of group galaxies is significantly narrower than that of isolated galaxies (10σ). We suggest that the warm/hot CGM does not exist as a superposition of halos; instead, the virial temperature of the halo is hot enough for O vi to be further ionized. The remaining O vi likely exists at the interface between hot diffuse gas and cooler regions of the CGM.
Laser-light spatial-domain scanning function for deconvolution of blurred photographs using point-spread and holographic Fourier transform division filter
STIS long-slit medium resolution spectra reduced in CALSTIS extended-source mode with narrow extraction heights (GWIDTH=3 pixels) show photometric uncertainties of +/- 3% relative to point-source extractions. These uncertainties are introduced through interpolation in the spectral image rectification processing stage, and are correlated with the number of pixel crossings the spectral profile core encounters in the spatial direction. The line-spread-function may be determined as a function of pixel crossing- position from calibration data sub-sampled in the spatial direction. This line spread function will be applied to science data to perform optimal extractions and point- source de-blending. Wavelength and breathing effects will be studied. Viability of the method to de-convolve extended source 'blobs' will be investigated.
Purpose: Scatter-to-primary ratio (SPR) is an important dosimetric quantity that describes the contribution from the scatter photons in an external photon beam. The purpose of this study is to develop an improved analytical formula to describe SPR as a function of circular field size (r) and depth (d) using Monte Carlo (MC) simulation. Methods: MC simulation was performed for Mohan photon spectra (Co-60, 4, 6, 10, 15, 23 MV) using EGSNRC code. Point-spread scatter dose kernels in water are generated. The scatter-to-primary ratio (SPR) is also calculated using MC simulation as a function of field size for circular field size with radius r and depth d. The doses from forward scatter and backscatter photons are calculated using a convolution of the point-spread scatter dose kernel and by accounting for scatter photons contributing to dose before (z' d) reaching the depth of interest, d, where z' is the location of scatter photons, respectively. The depth dependence of the ratio of the forward scatter and backscatter doses is determined as a function of depth and field size. Results: We are able to improve the existing 3-parameter (a, w, d0) empirical formula for SPR by introducing depth dependence for one of the parameter d0, which becomes 0 for deeper depths. The depth dependence of d0 can be directly calculated as a ratio of backscatter-to-forward scatter doses for otherwise the same field and depth. With the improved empirical formula, we can fit SPR for all megavoltage photon beams to within 2%. Existing 3-parameter formula cannot fit SPR data for Co-60 to better than 3.1%. Conclusion: An improved empirical formula is developed to fit SPR for all megavoltage photon energies to within 2%.
Field equations for spreading resistance of variable capacitance diode and point contact variable resistance diode derived in terms of frequency, material characteristics and physical parameters
Purpose: To characterize the modulation transfer function (MTF) of the pre-clinical (phase II) head scanner developed for proton computed tomography (pCT) by the pCT collaboration. To evaluate the spatial resolution achievable by this system. Methods: Our phase II proton CT scanner prototype consists of two silicon telescopes that track individual protons upstream and downstream from a phantom, and a 5-stage scintillation detector that measures a combination of the residual energy and range of the proton. Residual energy is converted to water equivalent path length (WEPL) of the protons in the scanned object. The set of WEPL values and associated paths of protons passing through the object over a 360° angular scan is processed by an iterative parallelizable reconstruction algorithm that runs on GP-GPU hardware. A custom edge phantom composed of water-equivalent polymer and tissue-equivalent material inserts was constructed. The phantom was first simulated in Geant4 and then built to perform experimental beam tests with 200 MeV protons at the Northwestern Medicine Chicago Proton Center. The oversampling method was used to construct radial and azimuthal edge spread functions and modulation transfer functions. The spatial resolution was defined by the 10% point of the modulation transfer function in units of lp/cm. Results: The spatial resolution of the image was found to be strongly correlated with the radial position of the insert but independent of the relative stopping power of the insert. The spatial resolution varies between roughly 4 and 6 lp/cm in both the the radial and azimuthal directions depending on the radial displacement of the edge. Conclusion: The amount of image degradation due to our detector system is small compared with the effects of multiple Coulomb scattering, pixelation of the image and the reconstruction algorithm. Improvements in reconstruction will be made in order to achieve the theoretical limits of spatial resolution.
In contrast to the ordinary contact line problem, virtually no information is available on the similar problem associated with a molten material spreading on a solid which is below the melt's fusion point. The latter is a more complex problem which heat transfer and solidification take place simultaneously with spreading, and requires answers not only for the hot melt's advance speed over the cold solid as a function of contact angle, but also for how one is to predict the point of the molten contact line's arrest by freezing. This issues are of importance in evolving methods of materials processing. The purpose of our work is to develop, based on both experiments and theory, an understanding of the dynamic processes that occur when a molten droplet touches a subcooled solid, spreads partly over it by capillary action, and freezes. We seek answers to the following basic questions. First, what is the relationship between the melt's contact line speed and the apparent (dynamic) contact angle? Secondly, at what point will the contact line modon be arrested by freezing? The talk will describe three components of our work: (1) deposition experiments with small molten droplets; (2) investigation of the dynamics of the molten contact line by means of a novel forced spreading method; and (3) an attempt to provide a theoretical framework for answering the basic questions posed above.