Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “point spread functions”

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 145 records · Page 8

An Efficient End-To-End Model for the Kepler Photometer

The Kepler Mission is designed to characterize the frequency of Earth-sized planets in the habitable zones of solar-like stars in the solar galactic neighborhood by observing >100,000 main-sequence stars in a >100 square degree field of view (FOV) and seeking evidence of transiting planets. As part of the system engineering effort, we have developed an End-To-End Model (ETEM) of the photometer to better characterize the expected performance of the instrument and to guide us in making design trades. This model incorporates engineering information such as the point spread function, time histories of pointing offsets, operating temperature, quantization noise, the effects of shutterless readout, and read noise. Astrophysical parameters, such as a realistic distribution of stars vs. magnitude for the chosen FOV, zodiacal light, and cosmic ray events are also included. For a given set of design and operating parameters, ETEM generates pixel time series for all pixels of interest for a single CCD channel of the photometer. These time series are then processed to form light curves for the target stars and the impact of various noise sources on the combined differential photometric precision can be determined. This model is of particular value when investigating the effects of noise sources that cannot be easily subjected to direct analysis, such as residual pointing offsets, thermal drift or cosmic ray effects. This version of ETEM features extremely efficient computation times relative to the previous version while maintaining a high degree of fidelity with respect to the realism of the relevant phenomena.

performance modeling↗

Thick Disks in the Hubble Space Telescope Frontier Fields

Thick disk evolution is studied using edge-on galaxies in two Hubble Space Telescope Frontier Field Parallels. The galaxies were separated into 72 clumpy types and 35 spiral types with bulges. Perpendicular light profiles in F435W, F606W, and F814W ( B , V , and I ) passbands were measured at 1 pixel intervals along the major axes and fitted to sech{sup 2} functions convolved with the instrument line spread function (LSF). The LSF was determined from the average point spread function of ∼20 stars in each passband and field, convolved with a line of uniform brightness to simulate disk blurring. A spread function for a clumpy disk was also used for comparison. The resulting scale heights were found to be proportional to galactic mass, with the average height for a 10{sup 10±0.5} M {sub ⊙} galaxy at z = 2 ± 0.5 equal to 0.63 ± 0.24 kpc. This value is probably the result of a blend between thin and thick disk components that cannot be resolved. Evidence for such two-component structure is present in an inverse correlation between height and midplane surface brightness. Models suggest that the thick disk is observed best between the clumps, and there the average scale height is 1.06 ± 0.43 kpc for the same mass and redshift. A 0.63 ± 0.68 mag V − I color differential with height is also evidence for a mixture of thin and thick components.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-CAMPUS-I-06: Y90 PET/CT for the Instantaneous Determination of Both Target and Non-Target Absorbed Doses Following Hepatic Radioembolization

Purpose The process of converting Yttrium-90 (Y90) PET/CT images into 3D absorbed dose maps will be explained. The simple methods presented will allow the medical physicst to analyze Y90 PET images following radioembolization and determine the absorbed dose to tumor, normal liver parenchyma and other areas of interest, without application of Monte-Carlo radiation transport or dose-point-kernel (DPK) convolution. Methods Absorbed dose can be computed from Y90 PET/CT images based on the premise that radioembolization is a permanent implant with a constant relative activity distribution after infusion. Many Y90 PET/CT publications have used DPK convolution to obtain 3D absorbed dose maps. However, this method requires specialized software limiting clinical utility. The Local Deposition method, an alternative to DPK convolution, can be used to obtain absorbed dose and requires no additional computer processing. Pixel values from regions of interest drawn on Y90 PET/CT images can be converted to absorbed dose (Gy) by multiplication with a scalar constant. Results There is evidence that suggests the Local Deposition method may actually be more accurate than DPK convolution and it has been successfully used in a recent Y90 PET/CT publication. We have analytically compared dose-volume-histograms (DVH) for phantom hot-spheres to determine the difference between the DPK and Local Deposition methods, as a function of PET scanner point-spread-function for Y90. We have found that for PET/CT systems with a FWHM greater than 3.0 mm when imaging Y90, the Local Deposition Method provides a more accurate representation of DVH, regardless of target size than DPK convolution. Conclusion Using the Local Deposition Method, post-radioembolization Y90 PET/CT images can be transformed into 3D absorbed dose maps of the liver. An interventional radiologist or a Medical Physicist can perform this transformation in a clinical setting, allowing for rapid prediction of treatment efficacy by comparison to published tumoricidal thresholds.

60 APPLIED LIFE SCIENCES↗

Bootstrap estimation of the effect of instrument response function uncertainty on the reconstruction of fusion neutron sources

Neutron imagers are important diagnostics for the inertial confinement fusion implosions at the National Ignition Facility. They provide two- and three-dimensional reconstructions of the neutron source shape that are key indicators of the overall performance. To interpret the shape results properly, it is critical to estimate the uncertainty in those reconstructions. There are two main sources of uncertainties: limited neutron statistics, leading to random errors in the reconstructed images, and incomplete knowledge of the instrument response function (the pinhole-dependent point spread function). While the statistical errors dominate the uncertainty for lower yield deuterium-tritium (DT) shots, errors due to the instrument response function dominate the uncertainty for DT yields on the order of 10 16 neutrons or higher. In this work, a bootstrapping method estimates the uncertainty in a reconstructed image due to the incomplete knowledge of the instrument response function. The main reconstruction is created from the fixed collection of pinhole images that are best aligned with the neutron source. Additional reconstructions are then built using subsets of that collection of images. Variations in the shapes of these additional reconstructions originate solely from uncertainties in the instrument response function, allowing us to use them to provide an additional systematic uncertainty estimate.

47 OTHER INSTRUMENTATION↗

Scatter and Blur Corrections for High-Energy X-Ray Radiography

High-energy X-ray radiography is useful as a highly penetrating method for imaging through dense materials. However, the primary modes of interaction of X-rays at these energies involve scattering or the production of secondary high-energy photons, which can interfere with the image. In addition, detector blurring, often resulting from scatter within the detector, can reduce image sharpness. Both of these processes can be mitigated with the use of convolution kernels, with the main challenge being that the proper kernel to use is not known, particularly for the scatter contribution. By radiographing solid slabs of uniform attenuation, we show that point spread functions and material-specific point scatter functions can be determined to significantly reduce the effect of detector blurring and object scatter. Constraining the fits to the slabs and uniform transmission within the slabs is sufficient to recover these functions. A functional form that reproduces the angular distribution of high-energy bremsstrahlung X-rays is presented for recovering point scatter functions. In conclusion, the method is applied to radiographs of objects from bremsstrahlung X-ray sources operating at 4- and 7.5-MV endpoint energies and a significant increase in sharpness is observed.

Blind deconvolution↗

A NEW METHOD FOR FINDING POINT SOURCES IN HIGH-ENERGY NEUTRINO DATA

The IceCube collaboration has reported the first detection of high-energy astrophysical neutrinos, including ∼50 high-energy starting events, but no individual sources have been identified. It is therefore important to develop the most sensitive and efficient possible algorithms to identify the point sources of these neutrinos. The most popular current method works by exploring a dense grid of possible directions to individual sources, and identifying the single direction with the maximum probability of having produced multiple detected neutrinos. This method has numerous strengths, but it is computationally intensive and because it focuses on the single best location for a point source, additional point sources are not included in the evidence. We propose a new maximum likelihood method that uses the angular separations between all pairs of neutrinos in the data. Unlike existing autocorrelation methods for this type of analysis, which also use angular separations between neutrino pairs, our method incorporates information about the point-spread function and can identify individual point sources. We find that if the angular resolution is a few degrees or better, then this approach reduces both false positive and false negative errors compared to the current method, and is also more computationally efficient up to, potentially, hundreds of thousands of detected neutrinos.

79 ASTRONOMY AND ASTROPHYSICS↗

Application of information theory to the design of line-scan imaging systems

Information theory is used to formulate a single figure of merit for assessing the performance of line scan imaging systems as a function of their spatial response (point spread function or modulation transfer function), sensitivity, sampling and quantization intervals, and the statistical properties of a random radiance field. Computational results for the information density and efficiency (i.e., the ratio of information density to data density) are intuitively satisfying and compare well with experimental and theoretical results obtained by earlier investigators concerned with the performance of TV systems.

Huck, F. O.↗

Hubble Space Telescope: Optical telescope assembly handbook. Version 1.0

The Hubble Space Telescope is described along with how its design affects the images produced at the Science Instruments. An overview is presented of the hardware. Details are presented of the focal plane, throughput of the telescope, and the point spread function (image of an unresolved point source). Some detailed simulations are available of this, which might be useful to observers in planning their observations and in reducing their data.

Burrows, Chris↗

Quantifying Uncertainties Through Advanced Theoretical Analysis for Image and Signal Reconstruction, NLV-02-15, Year 2 of 2

Presentation describes project to (1) adapt accelerated matrix-free linear system solvers for applying new statistical approaches to very large–scale problems, (2) characterize the degree to which new statistical methods are superior to classical methods in respect to the sensitivity to “tuning parameters,” and (3) develop a new mathematical theory and associated numerical methods for computing the point spread function of an x-ray imaging system.

97 MATHEMATICS AND COMPUTING↗

NEAR-INFRARED IMAGING OF A z = 6.42 QUASAR HOST GALAXY WITH THE HUBBLE SPACE TELESCOPE WIDE FIELD CAMERA 3

We report on deep near-infrared F125W (J) and F160W (H) Hubble Space Telescope Wide Field Camera 3 images of the z = 6.42 quasar J1148+5251 to attempt to detect rest-frame near-ultraviolet emission from the host galaxy. These observations included contemporaneous observations of a nearby star of similar near-infrared colors to measure temporal variations in the telescope and instrument point-spread function (PSF). We subtract the quasar point source using both this direct PSF and a model PSF. Using direct subtraction, we measure an upper limit for the quasar host galaxy of m{sub J} > 22.8 and m{sub H} > 23.0 AB mag (2 {sigma}). After subtracting our best model PSF, we measure a limiting surface brightness from 0.''3 to 0.''5 radius of {mu}{sub J} > 23.5 and {mu}{sub H} > 23.7 AB mag arcsec{sup -2} (2 {sigma}). We test the ability of the model subtraction method to recover the host galaxy flux by simulating host galaxies with varying integrated magnitude, effective radius, and Sersic index, and conducting the same analysis. These models indicate that the surface brightness limit ({mu}{sub J} > 23.5 AB mag arcsec{sup -2}) corresponds to an integrated upper limit of m{sub J} > 22-23 AB mag, consistent with the direct subtraction method. Combined with existing far-infrared observations, this gives an infrared excess log (IRX) > 1.0 and corresponding ultraviolet spectral slope {beta} > -1.2 {+-} 0.2. These values match those of most local luminous infrared galaxies, but are redder than those of almost all local star-forming galaxies and z {approx_equal} 6 Lyman break galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

CALIBRATION OF THE NuSTAR HIGH-ENERGY FOCUSING X-RAY TELESCOPE

We present the calibration of the Nuclear Spectroscopic Telescope Array (NuSTAR) X-ray satellite. We used the Crab as the primary effective area calibrator and constructed a piece-wise linear spline function to modify the vignetting response. The achieved residuals for all off-axis angles and energies, compared to the assumed spectrum, are typically better than ±2% up to 40 keV and 5%–10% above due to limited counting statistics. An empirical adjustment to the theoretical two-dimensional point-spread function (PSF) was found using several strong point sources, and no increase of the PSF half-power diameter has been observed since the beginning of the mission. We report on the detector gain calibration, good to 60 eV for all grades, and discuss the timing capabilities of the observatory, which has an absolute timing of ±3 ms. Finally, we present cross-calibration results from two campaigns between all the major concurrent X-ray observatories (Chandra, Swift, Suzaku, and XMM-Newton), conducted in 2012 and 2013 on the sources 3C 273 and PKS 2155-304, and show that the differences in measured flux is within ∼10% for all instruments with respect to NuSTAR.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Compound Poisson Generator Approach to Point-source Inference in Astrophysics

Abstract The identification and description of point sources is one of the oldest problems in astronomy, yet even today the correct statistical treatment for point sources remains one of the field’s hardest problems. For dim or crowded sources, likelihood-based inference methods are required to estimate the uncertainty on the characteristics of the source population. In this work, a new parametric likelihood is constructed for this problem using compound Poisson generator (CPG) functionals that incorporate instrumental effects from first principles. We demonstrate that the CPG approach exhibits a number of advantages over non-Poissonian template fitting (NPTF)—an existing method—in a series of test scenarios in the context of X-ray astronomy. These demonstrations show that the effect of the point-spread function, effective area, and choice of point-source spatial distribution cannot, generally, be factorized as they are in NPTF, while the new CPG construction is validated in these scenarios. Separately, an examination of the diffuse-flux emission limit is used to show that most simple choices of priors on the standard parameterization of the population model can result in unexpected biases: when a model comprising both a point-source population and diffuse component is applied to this limit, nearly all observed flux will be assigned to either the population or to the diffuse component. A new parameterization is presented for these priors that properly estimates the uncertainties in this limit. In this choice of priors, CPG correctly identifies that the fraction of flux assigned to the population model cannot be constrained by the data.

79 ASTRONOMY AND ASTROPHYSICS↗

Distinguishing one from many using super-resolution compressive sensing

We present that distinguishing whether a signal corresponds to a single source or a limited number of highly overlapping point spread functions (PSFs) is a ubiquitous problem across all imaging scales, whether detecting receptor-ligand interactions in cells or detecting binary stars. Super-resolution imaging based upon compressed sensing exploits the relative sparseness of the point sources to successfully resolve sources which may be separated by much less than the Rayleigh criterion. However, as a solution to an underdetermined system of linear equations, compressive sensing requires the imposition of constraints which may not always be valid. One typical constraint is that the PSF is known. However, the PSF of the actual optical system may reflect aberrations not present in the theoretical ideal optical system. Even when the optics are well characterized, the actual PSF may reflect factors such as non-uniform emission of the point source (e.g. fluorophore dipole emission). As such, the actual PSF may differ from the PSF used as a constraint. Similarly, multiple different regularization constraints have been suggested including the l 1 -norm, l 0 -norm, and generalized Gaussian Markov random fields (GGMRFs), each of which imposes a different constraint. Other important factors include the signal-to-noise ratio of the point sources and whether the point sources vary in intensity. In this work, we explore how these factors influence super-resolution image recovery robustness, determining the sensitivity and specificity. In conclusion, we determine an approach that is more robust to the types of PSF errors present in actual optical systems.

42 ENGINEERING↗

Performance of Kitt Peak’s Mayall 4-meter telescope during DESI commissioning

In preparation for the Dark Energy Spectroscopic Instrument (DESI), a new top end was installed on the Mayall 4-meter telescope at Kitt Peak National Observatory. The refurbished telescope and the DESI instrument were successfully commissioned on sky between 2019 October and 2020 March. Here we describe the pointing, tracking and imaging performance of the Mayall telescope equipped with its new DESI prime focus corrector, as measured by six guider cameras sampling the outer edge of DESI's focal plane. Analyzing ~500,000 guider images acquired during commissioning, we find a median delivered image FWHM of 1.1 arcseconds (in the r-band at 650 nm), with the distribution extending to a best-case value of ~0.6 arcseconds. The point spread function is well characterized by a Moffat profile with a power-law index of β ~ 3.5 and little dependence of β on FWHM. The shape and size of the PSF delivered by the new corrector at a field angle of 1.57 degrees are very similar to those measured with the old Mayall corrector on axis. Here, we also find that the Mayall achieves excellent pointing accuracy (several arcseconds RMS) and minimal open-loop tracking drift (< 1 milliarcsecond per second), improvements on the telecope's pre-DESI performance. In the future, employing DESI's active focus adjustment capabilities will likely further improve the Mayall/DESI delivered image quality.

47 OTHER INSTRUMENTATION↗

A 189 MHz, 2400 deg{sup 2} POLARIZATION SURVEY WITH THE MURCHISON WIDEFIELD ARRAY 32-ELEMENT PROTOTYPE

We present a Stokes I, Q and U survey at 189 MHz with the Murchison Widefield Array 32 element prototype covering 2400 deg{sup 2}. The survey has a 15.6 arcmin angular resolution and achieves a noise level of 15 mJy beam{sup -1}. We demonstrate a novel interferometric data analysis that involves calibration of drift scan data, integration through the co-addition of warped snapshot images, and deconvolution of the point-spread function through forward modeling. We present a point source catalog down to a flux limit of 4 Jy. We detect polarization from only one of the sources, PMN J0351-2744, at a level of 1.8% {+-} 0.4%, whereas the remaining sources have a polarization fraction below 2%. Compared to a reported average value of 7% at 1.4 GHz, the polarization fraction of compact sources significantly decreases at low frequencies. We find a wealth of diffuse polarized emission across a large area of the survey with a maximum peak of {approx}13 K, primarily with positive rotation measure values smaller than +10 rad m{sup -2}. The small values observed indicate that the emission is likely to have a local origin (closer than a few hundred parsecs). There is a large sky area at {alpha} {>=} 2{sup h}30{sup m} where the diffuse polarized emission rms is fainter than 1 K. Within this area of low Galactic polarization we characterize the foreground properties in a cold sky patch at ({alpha}, {delta}) = (4{sup h}, -27. Degree-Sign 6) in terms of three-dimensional power spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimating Residual Uncertainties for Lunar Irradiance Measurements Due to Imaging Acquisition Parameters

Lunar calibration is a commonly used method to track a climate satellite sensor's long-term radiometric stability. We present a modeling approach to examine the satellite sensor lunar observation uncertainties due to several important aspects related to the lunar image acquisition by the satellite sensor: lunar pixel shift, point spread function (PSF), lunar orientation, pitch, and oversampling rates. Our analyses can be summarized as follows: (1) The sensor observed lunar irradiance can vary due to small lunar pixel shift if the PSF is less than ideal; (2) During lunar calibration, an unstable oversampling rate due to spacecraft control will result in errors in observed lunar irradiance. A drift in oversampling rate would result in a bias in observed lunar irradiance and a random variation in oversampling rate would cause random error in lunar irradiance. Increasing the overall oversampling rates can reduce random error in observed lunar irradiance but would not change the biases in the observation; (3) Furthermore, the biases can vary when the Moon is observed at different orientations. Our results show impacts on observed lunar irradiance are on the order of 0.1 percent, which is a significant part of the overall uncertainty for a lunar irradiance measurement of a climate satellite sensor.

Oversampling↗

Mitigation of the Brighter-Fatter Effect in the LSST Camera

Thick, fully depleted charge-coupled devices (CCDs) are known to exhibit non-linear behavior at high signal levels due to the dynamic behavior of charges collecting in the potential wells of pixels, called the brighter-fatter effect (BFE). This particularly impacts bright calibration stars, which appear larger than their intrinsic shape, creating a flux-dependent point-spread function (PSF) that if left unmitigated, could make up a large fraction of the error budget in Stage IV weak-lensing (WL) surveys such as the Legacy Survey of Space and Time (LSST). In this paper, we analyze image measurements of flat fields and artificial stars taken at different illumination levels with the LSST Camera (LSSTCam) at SLAC National Accelerator Laboratory in order to quantify this effect in the LSST Camera before and after a previously introduced correction technique. We observe that the BFE evolves anisotropically as a function of flux due to higher-order BFEs, which violates the fundamental assumption of this correction method. We then introduce a new sampling method based on a physically motivated model to account these higher-order terms in the correction, and then we test the modified correction on both datasets. We find that the new method corrects the effect in flat fields better than it corrects the effect in artificial stars which we conclude is the result of a unmodeled curl component of the deflection field by the correction. We use these results to define a new metric for the full-well capacity of our sensors and advise image processing strategies to further limit the impact of the effect on LSST WL science pathways.

47 OTHER INSTRUMENTATION↗

High-resolution broadband spectroscopy using externally dispersed interferometry at the Hale telescope: Part 2, photon noise theory

High-resolution broadband spectroscopy at near-infrared (NIR) wavelengths (950 to 2450 nm) has been performed using externally dispersed interferometry (EDI) at the Hale telescope at Mt. Palomar, with the TEDI interferometer mounted within the central hole of the 200-in. primary mirror in series with the comounted TripleSpec NIR echelle spectrograph. These are the first multidelay EDI demonstrations on starlight. We demonstrated very high (10×) resolution boost and dramatic (20× or more) robustness to point spread function wavelength drifts in the native spectrograph. Data analysis, results, and instrument noise are described in a companion paper (part 1). This part 2 describes theoretical photon limited and readout noise limited behaviors, using simulated spectra and instrument model with noise added at the detector. We show that a single interferometer delay can be used to reduce the high frequency noise at the original resolution (1× boost case), and that except for delays much smaller than the native response peak half width, the fringing and nonfringing noises act uncorrelated and add in quadrature. This is due to the frequency shifting of the noise due to the heterodyning effect. We find a sum rule for the noise variance for multiple delays. The multiple delay EDI using a Gaussian distribution of exposure times has noise-to-signal ratio for photon-limited noise similar to a classical spectrograph with reduced slitwidth and reduced flux, proportional to the square root of resolution boost achieved, but without the focal spot limitation and pixel spacing Nyquist limitations. At low boost (~1×) EDI has ~1.4× smaller noise than conventional, and at >10× boost, EDI has ~1.4× larger noise than conventional. Readout noise is minimized by the use of three or four steps instead of 10 of TEDI. Net noise grows as step phases change from symmetrical arrangement with wavenumber across the band. As a result, for three (or four) steps, we calculate a multiplicative bandwidth of 1.8:1 (2.3:1), sufficient to handle the visible band (400 to 700 nm, 1.8:1) and most of TripleSpec (2.6:1).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗