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 199 records · Page 11

Characterization of scatter in digital mammography from use of Monte Carlo simulations and comparison to physical measurements

Purpose: Monte Carlo simulations were performed with the goal of verifying previously published physical measurements characterizing scatter as a function of apparent thickness. A secondary goal was to provide a way of determining what effect tissue glandularity might have on the scatter characteristics of breast tissue. The overall reason for characterizing mammography scatter in this research is the application of these data to an image processing-based scatter-correction program. Methods: MCNPX was used to simulate scatter from an infinitesimal pencil beam using typical mammography geometries and techniques. The spreading of the pencil beam was characterized by two parameters: mean radial extent (MRE) and scatter fraction (SF). The SF and MRE were found as functions of target, filter, tube potential, phantom thickness, and the presence or absence of a grid. The SF was determined by separating scatter and primary by the angle of incidence on the detector, then finding the ratio of the measured scatter to the total number of detected events. The accuracy of the MRE was determined by placing ring-shaped tallies around the impulse and fitting those data to the point-spread function (PSF) equation using the value for MRE derived from the physical measurements. The goodness-of-fit was determined for each data set as a means of assessing the accuracy of the physical MRE data. The effect of breast glandularity on the SF, MRE, and apparent tissue thickness was also considered for a limited number of techniques. Results: The agreement between the physical measurements and the results of the Monte Carlo simulations was assessed. With a grid, the SFs ranged from 0.065 to 0.089, with absolute differences between the measured and simulated SFs averaging 0.02. Without a grid, the range was 0.28–0.51, with absolute differences averaging −0.01. The goodness-of-fit values comparing the Monte Carlo data to the PSF from the physical measurements ranged from 0.96 to 1.00 with a grid and 0.65 to 0.86 without a grid. Analysis of the data suggested that the nongrid data could be better described by a biexponential function than the single exponential used here. The simulations assessing the effect of breast composition on SF and MRE showed only a slight impact on these quantities. When compared to a mix of 50% glandular/50% adipose tissue, the impact of substituting adipose or glandular breast compositions on the apparent thickness of the tissue was about 5%. Conclusions: The findings show agreement between the physical measurements published previously and the Monte Carlo simulations presented here; the resulting data can therefore be used more confidently for an application such as image processing-based scatter correction. The findings also suggest that breast composition does not have a major impact on the scatter characteristics of breast tissue. Application of the scatter data to the development of a scatter-correction software program can be simplified by ignoring the variations in density among breast tissues.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Hubble Space Telescope Images of Nearby Luminous Quasars: Results for Eight Quasars and Tests of the Detection Sensitivity - 2

Observations with the Wide-Field Camera of the Hubble Space Telescope (HST) are presented for eight intrinsically luminous quasars with redshifts between 0.16 and 0.29. These observations, when combined with a similar HST study of the quasar PKS 2349-014, show that luminous nearby quasars exist in a variety of environments. Seven companion galaxies brighter than M(V) = 16.5 (H(sub 0) = 100 km s(sup -1) Mpc(sup -1), Omega(sub 0) = 1.0) lie within a projected distance of 25 kpc of the quasars; three of the companions are located closer than 3'' (6 kpc projected distance) from the quasars, well within the volume that would be enclosed by a typical L* host galaxy. The observed association of quasars and companion galaxies is statistically significant and may he an important element in the luminous-quasar phenomenon. Apparent host galaxies are detected for three of the quasars: PG 1116+215, 3C 273, and PG 1444+407; the hosts have an average absolute magnitude of about 0.6 mag brighter than L*. The agreement between the previously published major-axis directions in ground-based images and in the present HST images of 3C 273 and PG 1444+407 constitutes important evidence supporting the reality of these candidate host galaxies. Upper limits are placed on the visual-band brightnesses of representative galactic hosts for all the quasars. These limits are established by placing galaxy images obtained with HST underneath the quasars and measuring at what faintness level the known galaxies are detected. On average, the HST spirals would have been detected if they were as faint as 1 mag below L*, and the early-type galaxies could have been detected down to a brightness level of about L*, where L* is the Schechter characteristic luminosity of field galaxies. Smooth, featureless galaxy models (exponential disks or de Vaucouleurs profiles) are fitted to the residual light after a best-fitting point source is subtracted from the quasar images. The results show that smooth spiral galaxies brighter than, on average, about L*, would have been detected. These upper limits, or possible detections, are consistent with, for example, the eight luminous quasars studied in this paper, occurring in host galaxies that have a Shechter luminosity function with a lower cutoff in the range 0.01-0.1 L*. Tests are performed to determine if our failure to detect, in some cases, luminous host galaxies could be an artifact caused by our analysis procedures. These tests include comparing the measured point-spread function (PSF) for our HST observations with the PSFs used in previous ground-based studies of host galaxies, measuring the fluctuations in the sky signals that were subtracted from the quasar images, evaluating empirically the effects of using different stellar PSFs in the analysis, carrying out the subtraction of the stellar (nuclear) source in different ways, creating and analyzing artificial active galactic nuclei (AGNs) with known surface brightnesses, and fitting the observed quasar light to an analytic model that includes a host galaxy.

Bahcall, John N.↗

Hypercalibration: A Pan-STARRS1-Based Recalibration of the Sloan Digital Sky Survey Photometry

In this paper, we present a recalibration of the Sloan Digital Sky Survey (SDSS) photometry with new flat fields and zero points derived from Pan-STARRS1. Using point-spread function (PSF) photometry of 60 million stars with 16 < r < 20, we derive a model of amplifier gain and flat-field corrections with per-run rms residuals of 3 millimagnitudes (mmag) in griz bands and 15 mmag in u band. The new photometric zero points are adjusted to leave the median in the Galactic north unchanged for compatibility with previous SDSS work. We also identify transient non-photometric periods in SDSS ("contrails") based on photometric deviations co-temporal in SDSS bands. Finally, the recalibrated stellar PSF photometry of SDSS and PS1 has an rms difference of {9, 7, 7, 8} mmag in griz, respectively, when averaged over 15' regions.

79 ASTRONOMY AND ASTROPHYSICS↗

A new approach to remove interorder background in high dispersion IUE images

A two dimensional Chebyshev interpolating scheme, followed by modeling of point spread functions as Voigt profiles, is proposed for final archiving IUESIPS to remove the interorder overlap from the smoothed background fluxes of high dispersion International Ultraviolet Explorer (IUE) images. Tests of the first and last stages of the algorithm suggest that it can give a considerably better zero point than the old IUESIPS algorithm which samples local interorder background.

Smith, Myron A.↗

Scale effects: HCMM data simulation. Usage of filtering techniques for scaling-up simulations

Image reduction used to simulate increase in altitude of an acquisition platform is equivalent to data smoothing, and can be achieved either by neighborhood averaging or by filtering techniques. The averaging approach is limited for accurate simulation. A filtering method is described which was based on the hypothesis that all changes due to altitude increase can be represented by a point spread function. Determination of the scale function and factor are discussed as well as implementation of the filtering. Filtering can be either in the spatial or frequency domain. In the spatial domain, filtering consists of the convolution of the image with the weights mask, and then of the declination of the points according to the appropriates scale factor. A simulation of an aircraft day image in the infrared channel is examined.

Digennaro, V.↗

SMASH: Survey of the MAgellanic Stellar History

The Large and Small Magellanic Clouds are unique local laboratories for studying the formation and evolution of small galaxies in exquisite detail. The Survey of the MAgellanic Stellar History (SMASH) is an NOAO community Dark Energy Camera (DECam) survey of the Clouds mapping 480 deg{sup 2} (distributed over ∼2400 square degrees at ∼20% filling factor) to ∼24th mag in ugriz . The primary goals of SMASH are to identify low surface brightness stellar populations associated with the stellar halos and tidal debris of the Clouds, and to derive spatially resolved star formation histories. Here, we present a summary of the survey, its data reduction, and a description of the first public Data Release (DR1). The SMASH DECam data have been reduced with a combination of the NOAO Community Pipeline, the PHOTRED automated point-spread-function photometry pipeline, and custom calibration software. The astrometric precision is ∼15 mas and the accuracy is ∼2 mas with respect to the Gaia reference frame. The photometric precision is ∼0.5%–0.7% in griz and ∼1% in u with a calibration accuracy of ∼1.3% in all bands. The median 5 σ point source depths in ugriz are 23.9, 24.8, 24.5, 24.2, and 23.5 mag. The SMASH data have already been used to discover the Hydra II Milky Way satellite, the SMASH 1 old globular cluster likely associated with the LMC, and extended stellar populations around the LMC out to R ∼ 18.4 kpc. SMASH DR1 contains measurements of ∼100 million objects distributed in 61 fields. A prototype version of the NOAO Data Lab provides data access and exploration tools.

79 ASTRONOMY AND ASTROPHYSICS↗

LANDSAT-D investigations in snow hydrology

The atmospheric radiative transfer calculation program (ATARD) and its supporting programs (setting up atmospheric profile, making Mie tables and an exponential-sum-fitting table) were completed. More sophisticated treatment of aerosol scattering (including angular phase function or asymmetric factor) and multichannel analysis of results from ATRAD are being developed. Some progress was made on a Monte Carlo program for examining two dimensional effects, specifically a surface boundary condition that varies across a scene. The MONTE program combines ATRAD and the Monte Carlo method together to produce an atmospheric point spread function. Currently the procedure passes monochromatic tests and the results are reasonable.

Dozier, J.↗

SYNTHESIZED SPECTRA OF OPTICALLY THIN EMISSION LINES PRODUCED BY THE BIFROST STELLAR ATMOSPHERE CODE, INCLUDING NONEQUILIBRIUM IONIZATION EFFECTS: A STUDY OF THE INTENSITY, NONTHERMAL LINE WIDTHS, AND DOPPLER SHIFTS

In recent years realistic 3D numerical models of the solar atmosphere have become available. The models attempt to recreate the solar atmosphere and mimic observations in the best way, in order to make it possible to couple complicated observations with physical properties such as the temperatures, densities, velocities, and magnetic fields. We here present a study of synthetic spectra created using the Bifrost code in order to assess how well they fit with previously taken solar data. A study of the synthetic intensity, nonthermal line widths, Doppler shifts, and correlations between any two of these three components of the spectra first assuming statistical equilibrium is made, followed by a report on some of the effects nonequilibrium ionization will have on the synthesized spectra. We find that the synthetic intensities compare well with the observations. The synthetic observations depend on the assumed resolution and point-spread function (PSF) of the instrument, and we find a large effect on the results, especially for intensity and nonthermal line width. The Doppler shifts produce the reported persistent redshifts for the transition region (TR) lines and blueshifts for the upper TR and corona lines. The nonthermal line widths reproduce the well-known turnoff point around (2–3) × 10{sup 5} K, but with much lower values than those observed. The nonthermal line widths tend to increase with decreasing assumed instrumental resolution, also when nonequilibrium ionization is included. Correlations between the nonthermal line width of any two TR line studies as reported by Chae et al. are reproduced, while the correlations of intensity to line width are reproduced only after applying a PSF to the data. Doppler shift correlations reported by Doschek for the TR lines and correlations of Doppler shift to nonthermal line width of the Fe xii{sub 19.5} line reported by Doschek et al. are reproduced.

79 ASTRONOMY AND ASTROPHYSICS↗

Limb brightening on Uranus - The visible spectrum

Results are presented for pinhole and slit narrowband photoelectric area-scanning photometry of the Uranus disk at eight central wavelengths ranging from 5600 to 7500 A. The observations were undertaken to detect limb brightening on the planet and to study the nature of its disk profile as a function of geometrical albedo. A coarse analysis of the results shows that the disk must exhibit limb brightening in the two CH4 bands at 6190 and 7300 A relative to all other wavebands investigated. Except for the two CH4 bands, the shapes and widths of the images are found to be essentially identical; point-spread-function data for all wavebands indicate no significant variation in the representative Gaussian l/e width with wavelength. A fine analysis of the data reveals absolute limb brightening in the 7300-A CH4 band and absolute limb darkening in the adjacent wavebands. These results are shown to be consistent with the hypothesis that clouds are absent from the visible atmosphere of Uranus.

Price, M. J.↗

Optimum constrained image restoration filters

The research described centered on development of an optimum image restoration filter (IRF) minimizing the radius of gyration of the corrected or composite system point-spread function (P-SF) subject to contraints, and reducing 2-dimensional spatial smearing or blurring of an image. The constraints are imposed on the radius of gyration of the IRF P-SF, the total restored image noise power, and the shape of the composite system frequency spectrum. The image degradation corresponds to mapping many points from the original image into a single resolution element. The P-SF is obtained as solution to a set of simultaneous differential equations obeying nonlinear integral constraints. Truncation errors due to edge effects are controlled by constraining the radius of gyration of the IRF P-SF. An iterative technique suppresses sidelobes of the composite system P-SF.

Riemer, T. E.↗

ZFIRE: The Evolution of the Stellar Mass Tully–Fisher Relation to Redshift ∼2.2

Using observations made with MOSFIRE on Keck I as part of the ZFIRE survey, we present the stellar mass Tully–Fisher relation at 2.0<z<2.5. The sample was drawn from a stellar-mass-limited, K{sub s}-band-selected catalog from ZFOURGE over the CANDELS area in the COSMOS field. We model the shear of the Hα emission line to derive rotational velocities at 2.2× the scale radius of an exponential disk (V{sub 2.2}). We correct for the blurring effect of a 2D point-spread function (PSF) and the fact that the MOSFIRE PSF is better approximated by a Moffat than a Gaussian, which is more typically assumed for natural seeing. We find for the Tully–Fisher relation at 2.0<z<2.5 that logV{sub 2.2}=(2.18±0.051)+(0.193 ± 0.108)(logM/M{sub ⊙}−10) and infer an evolution of the zero-point of ΔM/M{sub ⊙}=−0.25±0.16 dex or ΔM/M{sub ⊙}=−0.39±0.21 dex compared to z = 0 when adopting a fixed slope of 0.29 or 1/4.5, respectively. We also derive the alternative kinematic estimator S{sub 0.5}, with a best-fit relation logS{sub 0.5}=(2.06±0.032) + (0.211±0.086) (logM/M{sub ⊙}−10), and infer an evolution of ΔM/M{sub ⊙} =−0.45±0.13 dex compared to z<1.2 if we adopt a fixed slope. We investigate and review various systematics, such as PSF effects, projection effects, systematics related to stellar mass derivation, selection biases, and slope. We find that discrepancies between the various literature values are reduced when taking these into account. Our observations correspond well with the gradual evolution predicted by semianalytic models.

79 ASTRONOMY AND ASTROPHYSICS↗

MISR Level 1B1 Local Mode Radiance Data (MIB1LM_V1)

The results of two types of processing are included in this product. First, the Radiance Scaling operation converts the camera's digital number output to a measure of energy incident on the front optical surface. The measurement is expressed in units called radiance (energy per unit area, wavelength, and solid angle) as defined by the International Standard (SI). Second, Radiance Conditioning modifies the radiances to remove instrument-dependent effects. Specifically, image sharpening is applied, and focal-plane scattering is removed. Additionally, all radiances are adjusted to remove slight spectral sensitivity differences among the 1504 detector elements of each spectral band and each camera. In addition to the Level 1B1 radiometric product for MISR's Global Mode imagery, there is a separate Level 1B1 product for each high-resolution Local Mode scene. The Radiometric Product contains spectral radiances for all MISR channels (four spectral bands and nine cameras). Each radiance value represents the incident radiance averaged over the sensor's total band response. Processing includes both radiance scaling and conditioning steps. Radiance scaling converts the Level 1A data from digital counts to radiances using coefficients derived in combination with the On-Board Calibrator (OBC) and vicarious calibrations. The OBC contains Spectralon calibration panels which are deployed monthly and reflect sunlight into the cameras. The OBC detector standards then measure this reflected light to provide the calibration. Vicarious field campaigns are conducted less frequently but provide an independent methodology useful for reducing systematic errors. Radiance conditioning removes undesirable instrument effects. Image enhancement is provided by deconvolving the scene with the sensor's point-spread-function. Additionally, in-band scaling adjusts the reported radiances to correspond to a nominal band response profile. This frees the Level 2 software from the need to correct for detector element non-uniformities. No out-of-band correction is done for this product, nor are the data geometrically corrected or resampled at this point. In summary, the Level 1B1 Product contains the Data Numbers (DNs) radiometrically-scaled to radiances with no geometric resampling. [Location=GLOBAL] [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=275 m for red band only; Longitude_Resolution=275 m for red band only; Horizontal_Resolution_Range=250 km - < 500 km or approximately 2.5 degrees - < 5.0 degrees; Temporal_Resolution=about 15 orbits/day].

NEAR IR IMAGERY↗

MISR Level 1B1 Radiance Data (MI1B1_V1)

The results of two types of processing are included in this product. First, the Radiance Scaling operation converts the camera's digital number output to a measure of energy incident on the front optical surface. The measurement is expressed in units called radiance (energy per unit area, wavelength, and solid angle) as defined by the International Standard (SI). Second, Radiance Conditioning modifies the radiances to remove instrument-dependent effects. Specifically, image sharpening is applied, and focal-plane scattering is removed. Additionally, all radiances are adjusted to remove slight spectral sensitivity differences among the 1504 detector elements of each spectral band and each camera. The Radiometric Product contains spectral radiances for all MISR channels (four spectral bands and nine cameras). Each radiance value represents the incident radiance averaged over the sensor's total band response. Processing includes both radiance scaling and conditioning steps. Radiance scaling converts the Level 1A data from digital counts to radiances using coefficients derived in combination with the On-Board Calibrator (OBC) and vicarious calibrations. The OBC contains Spectralon calibration panels which are deployed monthly and reflect sunlight into the cameras. The OBC detector standards then measure this reflected light to provide the calibration. Vicarious field campaigns are conducted less frequently but provide an independent methodology useful for reducing systematic errors. Radiance conditioning removes undesirable instrument effects. Image enhancement is provided by deconvolving the scene with the sensor's point-spread-function. Additionally, in-band scaling adjusts the reported radiances to correspond to a nominal band response profile. This frees the Level 2 software from the need to correct for detector element non-uniformities. No out-of-band correction is done for this product, nor are the data geometrically corrected or resampled at this point. In summary, the Level 1B1 Product contains the Data Numbers (DNs) radiometrically-scaled to radiances with no geometric resampling. [Location=GLOBAL] [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=275 m for red band only; Longitude_Resolution=275 m for red band only; Horizontal_Resolution_Range=250 meters - < 500 meters; Temporal_Resolution=about 15 orbits/day; Temporal_Resolution_Range=Daily - < Weekly, about 15 orbits/day, Daily - < Weekly].

EOSDIS↗

MISR Level 1B1 Radiance Data (MI1B1_V2)

The results of two types of processing are included in this product. First, the Radiance Scaling operation converts the camera's digital number output to a measure of energy incident on the front optical surface. The measurement is expressed in units called radiance (energy per unit area, wavelength, and solid angle) as defined by the International Standard (SI). Second, Radiance Conditioning modifies the radiances to remove instrument-dependent effects. Specifically, image sharpening is applied, and focal-plane scattering is removed. Additionally, all radiances are adjusted to remove slight spectral sensitivity differences among the 1504 detector elements of each spectral band and each camera. The Radiometric Product contains spectral radiances for all MISR channels (four spectral bands and nine cameras). Each radiance value represents the incident radiance averaged over the sensor's total band response. Processing includes both radiance scaling and conditioning steps. Radiance scaling converts the Level 1A data from digital counts to radiances using coefficients derived in combination with the On-Board Calibrator (OBC) and vicarious calibrations. The OBC contains Spectralon calibration panels which are deployed monthly and reflect sunlight into the cameras. The OBC detector standards then measure this reflected light to provide the calibration. Vicarious field campaigns are conducted less frequently but provide an independent methodology useful for reducing systematic errors. Radiance conditioning removes undesirable instrument effects. Image enhancement is provided by deconvolving the scene with the sensor's point-spread-function. Additionally, in-band scaling adjusts the reported radiances to correspond to a nominal band response profile. This frees the Level 2 software from the need to correct for detector element non-uniformities. No out-of-band correction is done for this product, nor are the data geometrically corrected or resampled at this point. In summary, the Level 1B1 Product contains the Data Numbers (DNs) radiometrically-scaled to radiances with no geometric resampling. [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=1.1 km; Longitude_Resolution=1.1 km; Temporal_Resolution=about 15 orbits/day; Temporal_Resolution_Range=about 15 orbits/day].

AM-1↗

MISR Level 1B1 Local Mode Radiance Data (MIB1LM_V2)

The results of two types of processing are included in this product. First, the Radiance Scaling operation converts the camera's digital number output to a measure of energy incident on the front optical surface. The measurement is expressed in units called radiance (energy per unit area, wavelength, and solid angle) as defined by the International Standard (SI). Second, Radiance Conditioning modifies the radiances to remove instrument-dependent effects. Specifically, image sharpening is applied, and focal-plane scattering is removed. Additionally, all radiances are adjusted to remove slight spectral sensitivity differences among the 1504 detector elements of each spectral band and each camera. In addition to the Level 1B1 radiometric product for MISR's Global Mode imagery, there is a separate Level 1B1 product for each high-resolution Local Mode scene. The Radiometric Product contains spectral radiances for all MISR channels (four spectral bands and nine cameras). Each radiance value represents the incident radiance averaged over the sensor's total band response. Processing includes both radiance scaling and conditioning steps. Radiance scaling converts the Level 1A data from digital counts to radiances using coefficients derived in combination with the On-Board Calibrator (OBC) and vicarious calibrations. The OBC contains Spectralon calibration panels which are deployed monthly and reflect sunlight into the cameras. The OBC detector standards then measure this reflected light to provide the calibration. Vicarious field campaigns are conducted less frequently but provide an independent methodology useful for reducing systematic errors. Radiance conditioning removes undesirable instrument effects. Image enhancement is provided by deconvolving the scene with the sensor's point-spread-function. Additionally, in-band scaling adjusts the reported radiances to correspond to a nominal band response profile. This frees the Level 2 software from the need to correct for detector element non-uniformities. No out-of-band correction is done for this product, nor are the data geometrically corrected or resampled at this point. In summary, the Level 1B1 Product contains the Data Numbers (DNs) radiometrically-scaled to radiances with no geometric resampling. [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=275 m for red band only; Longitude_Resolution=275 m for red band only; Temporal_Resolution=about 15 orbits/day].

AM-1↗

AN ELECTRON-TRACKING COMPTON TELESCOPE FOR A SURVEY OF THE DEEP UNIVERSE BY MeV GAMMA-RAYS

Photon imaging for MeV gammas has serious difficulties due to huge backgrounds and unclearness in images, which originate from incompleteness in determining the physical parameters of Compton scattering in detection, e.g., lack of the directional information of the recoil electrons. The recent major mission/instrument in the MeV band, Compton Gamma Ray Observatory/COMPTEL, which was Compton Camera (CC), detected a mere ∼30 persistent sources. It is in stark contrast with the ∼2000 sources in the GeV band. Here we report the performance of an Electron-Tracking Compton Camera (ETCC), and prove that it has a good potential to break through this stagnation in MeV gamma-ray astronomy. The ETCC provides all the parameters of Compton-scattering by measuring 3D recoil electron tracks; then the Scatter Plane Deviation (SPD) lost in CCs is recovered. The energy loss rate (dE/dx), which CCs cannot measure, is also obtained, and is found to be helpful to reduce the background under conditions similar to those in space. Accordingly, the significance in gamma detection is improved severalfold. On the other hand, SPD is essential to determine the point-spread function (PSF) quantitatively. The SPD resolution is improved close to the theoretical limit for multiple scattering of recoil electrons. With such a well-determined PSF, we demonstrate for the first time that it is possible to provide reliable sensitivity in Compton imaging without utilizing an optimization algorithm. As such, this study highlights the fundamental weak-points of CCs. In contrast we demonstrate the possibility of ETCC reaching the sensitivity below 1 × 10{sup −12} erg cm{sup −2} s{sup −1} at 1 MeV.

79 ASTRONOMY AND ASTROPHYSICS↗

Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation

Purpose: Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction. Methods: Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according toFessler [“Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography,” IEEE Trans. Image Process. 5(3), 493–506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map. Results: Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization “map” designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter. Conclusions: Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

HST image restoration: A comparison of pre- and post-servicing mission results

A variety of image restoration techniques (e.g., Wiener filter, Lucy-Richardson, MEM) have been applied quite successfully to the aberrated HST images. The HST servicing mission (scheduled for late 1993 or early 1994) will install a corrective optics system (COSTAR) for the Faint Object Camera and spectrographs and replace the Wide Field/Planetary Camera with a second generation instrument (WF/PC-II) having its own corrective elements. The image quality is expected to be improved substantially with these new instruments. What then is the role of image restoration for the HST in the long term? Through a series of numerical experiments using model point-spread functions for both aberrated and unaberrated optics, we find that substantial improvements in image resolution can be obtained for post-servicing mission data using the same or similar algorithms as being employed now to correct aberrated images. Included in our investigations are studies of the photometric integrity of the restoration algorithms and explicit models for HST pointing errors (spacecraft jitter).

Hanisch, R. J.↗