Search NASA⌕ Search

SEARCH · Search NASA

Results for “Image Reconstruction”

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 109 records · Page 6

Image reconstruction by parametric cubic convolution

Cubic convolution, which has been discussed by Rifman and McKinnon (1974), was originally developed for the reconstruction of Landsat digital images. In the present investigation, the reconstruction properties of the one-parameter family of cubic convolution interpolation functions are considered and thee image degradation associated with reasonable choices of this parameter is analyzed. With the aid of an analysis in the frequency domain it is demonstrated that in an image-independent sense there is an optimal value for this parameter. The optimal value is not the standard value commonly referenced in the literature. It is also demonstrated that in an image-dependent sense, cubic convolution can be adapted to any class of images characterized by a common energy spectrum.

Park, S. K.↗

Efficient Image Reconstruction Method for a Millimeter-wave Shoe Scanner

The Pacific Northwest National Laboratory (PNNL) has recently developed an active 3D microwave/millimeter-wave shoe scanner. This system is designed to detect threats concealed within the soles of common footwear. The system operates over the 10-40 GHz frequency range for excellent penetration of common materials. The system presents several challenges for efficient image reconstruction, including the thick dielectric window, multi-row linear arrays, and focusing close to the antenna elements. In this paper, an efficient back-projection reconstruction algorithm is presented that overcomes these challenges. Experimental imaging results are shown that demonstrate high-resolution imaging performance for this new scanner.

Sheen, David M.↗

Relevance Vector Machine Image Reconstruction Algorithm for Electrical Capacitance Tomography With Explicit Uncertainty Estimates

We present a Relevance Vector Machine (RVM) based algorithm for electrical capacitance tomography (ECT) applications that can concurrently provide image reconstruction results and uncertainty estimates about the reconstruction. To illustrate the RVM operation in ECT, we simulate typical ECT scenarios, making explicit the connection between the reconstructed pixel values and the corresponding uncertainty estimates in each case. We compare the RVM reconstruction performance with that of the Iterative Landweber Method (ILM) and the least absolute shrinkage and selection operator (LASSO) in all the considered scenarios. The results show that, in addition to the key advantage of providing uncertainty measures, RVM can achieve similar reconstruction results with either lower or similar computational complexity.

42 ENGINEERING↗

Indigo: A Domain-Specific Language for Fast, Portable Image Reconstruction

© 2018 IEEE. Linear operators used in iterative methods like conjugate gradient have typically been implemented either as ''matrix-driven'' subroutines backed by explicit sparse or dense matrices, or as ''matrix-free'' subroutines that implement specific linear operations directly (e.g. FFTs). The matrix-driven approach is generally more portable because it can target widely-Available BLAS libraries, but it can be inefficient in terms of time and space complexity. In contrast, the matrix-free approach is more performant because it leverages structure in operations, but it requires each operator be re-implemented on each new platform. To increase performance and portability, we propose a hybrid approach that represents linear operators as expression trees. Leaf nodes in the tree are either matrix-free or matrix-driven operators, and interior nodes represent mathematical compositions (sums, products, transposes) or structural compositions (stacks, block diagonals, etc.) of the leaf operators. This representation enables expert-guided reordering and fusion transformations that can improve performance or reduce memory pressure. We implement our approach in a domain-specific language called Indigo. We assess Indigo on image reconstruction problems arising in four application areas: magnetic resonance imaging, ptychography, magnetic particle imaging, and fluorescent microscopy. We give performance results from vendor BLAS libraries, and we introduce specializations to Sparse BLAS routines that achieve near-Roofline performance on multi-core, many-core, and GPU systems.

Driscoll, Michael↗

Extracellular vesicles of calcifying turkey leg tendon characterized by immunocytochemistry and high voltage electron microscopic tomography and 3-D graphic image reconstruction

To gain insight into the structure and possible function of extracellular vesicles in certain calcifying vertebrate tissues, normally mineralizing leg tendons from the domestic turkey, Meleagris gallopavo, have been studied in two separate investigations, one concerning the electron microscopic immunolocalization of the 66 kDa phosphoprotein, osteopontin, and the other detailing the organization and distribution of mineral crystals associated with the vesicles as determined by high voltage microscopic tomography and 3-D graphic image reconstruction. Immunolabeling shows that osteopontin is related to extracellular vesicles of the tendon in the sense that its initial presence appears coincident with the development of mineral associated with the vesicle loci. By high voltage electron microscopy and 3-D imaging techniques, mineral crystals are found to consist of small irregularly shaped particles somewhat randomly oriented throughout individual vesicles sites. Their appearance is different from that found for the mineral observed within calcifying tendon collagen, and their 3-D disposition is not regularly ordered. Possible spatial and temporal relationships of vesicles, osteopontin, mineral, and collagen are being examined further by these approaches.

NASA Discipline Musculoskeletal↗

SU-E-I-41: Non-Cartesian MR Image Reconstruction with Integrated Gradient Non-Linearity Correction

Purpose: Nonlinearities in the spatial encoding gradients of MRI systems cause geometric distortion in images. Typically, this is retrospectively corrected via image-domain interpolation (a.k.a., “gradwarp”) albeit with a loss of spatial resolution. For non-Cartesian MRI, the latter problem is exaggerated by noise and undersampling artifact. In this study, we describe a novel correction strategy that accounts for gradient nonlinearities during — rather than after — non-Cartesian MRI reconstruction, and demonstrate that this approach mitigates the resolution loss that can occur with standard methods. Methods: To test the proposed method, the American College of Radiology (ACR) quality control phantom was scanned on at 1.5 T (General Electric, v16.0, “zoom” gradient) using a 1.6x undersampled 3D non- Cartesian Shells trajectory (GRE, FOV=24 cm3, 120 shells, 16552 shots, 512 readout, matrix=2403). Image reconstruction was first performed via standard k-space density-compensated gridding and retrospectively corrected via cubic spline interpolation. Image reconstruction was then separately performed using a k-space and image-domain densitycompensated type-3 non-uniform fast Fourier transform (NUFFT), which provides a direct mapping between non-Cartesian k-space samples and warped image space voxel locations. Thus, no separate distortion correction procedure is needed for the proposed approach. The gradient distortion field was determined using vendor provided calibration data. Results: Phantom scan results show that both processing approaches successfully correct geometric distortion. However, visual inspection of the ACR phantom spatial resolution inserts shows that the proposed strategy preserves the resolution of the nominal (uncorrected) reconstruction while “gradwarp” imparts marked spatial blurring (especially for the 1.0 and 1.1 mm inserts) and thus resolution loss. Conclusion: We've presented a novel reconstruction strategy for non-Cartesian MRI that corrects for gradient nonlinearities during — rather than after — reconstruction, and thus better preserves image resolution than traditional interpolation-based methods. This approach is expected to be especially advantageous when imaging with non-standard magnet geometries. NIH RR018898; NIH EB10065.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Initial Refinement and Optimization of the Image Reconstruction (U-MBIR) Algorithm to Optimize its Performance in Detecting Damage and Flaws

Concrete is a critical component of nuclear power plants; thus, its safety and reliability must be thoroughly examined throughout the life cycle of the structural system. During the life cycle of this infrastructure, monitoring of the concrete for signs of degradation should be performed through nondestructive evaluation (NDE). Ultrasonic measurements have been an industry standard for both surface and subsurface inspections. As such, Oak Ridge National Laboratory (ORNL) is developing advanced image reconstruction algorithms to overcome the limitations of traditional ultrasonic NDE methodologies. The results and discussion presented herein summarize the current state of the ultrasonic model–based iterative reconstruction (U-MBIR) algorithm developed at ORNL.

36 MATERIALS SCIENCE↗

Improved proton computed tomography by dual modality image reconstruction

Purpose: Proton computed tomography (CT) is a promising image modality for improving the stopping power estimates and dose calculations for particle therapy. However, the finite range of about 33 cm of water of most commercial proton therapy systems limits the sites that can be scanned from a full 360° rotation. In this paper the authors propose a method to overcome the problem using a dual modality reconstruction (DMR) combining the proton data with a cone-beam x-ray prior. Methods: A Catphan 600 phantom was scanned using a cone beam x-ray CT scanner. A digital replica of the phantom was created in the Monte Carlo code Geant4 and a 360° proton CT scan was simulated, storing the entrance and exit position and momentum vector of every proton. Proton CT images were reconstructed using a varying number of angles from the scan. The proton CT images were reconstructed using a constrained nonlinear conjugate gradient algorithm, minimizing total variation and the x-ray CT prior while remaining consistent with the proton projection data. The proton histories were reconstructed along curved cubic-spline paths. Results: The spatial resolution of the cone beam CT prior was retained for the fully sampled case and the 90° interval case, with the MTF = 0.5 (modulation transfer function) ranging from 5.22 to 5.65 linepairs/cm. In the 45° interval case, the MTF = 0.5 dropped to 3.91 linepairs/cm For the fully sampled DMR, the maximal root mean square (RMS) error was 0.006 in units of relative stopping power. For the limited angle cases the maximal RMS error was 0.18, an almost five-fold improvement over the cone beam CT estimate. Conclusions: Dual modality reconstruction yields the high spatial resolution of cone beam x-ray CT while maintaining the improved stopping power estimation of proton CT. In the case of limited angles, the use of prior image proton CT greatly improves the resolution and stopping power estimate, but does not fully achieve the quality of a 360° proton CT scan.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Tomographic Image Reconstruction for the Parallel-Slit Ring Collimator Fast Neutron Emission Tomography System

For the past three years, Oak Ridge National Laboratory (ORNL) has been developing a passive fast-neutron emission tomography capability. The goal of this development is the ability to quantify the neutron source strength of individual fuel pins (rods) in spent nuclear fuel assemblies. Such a system could be used to measure the burnup of each fuel pin in a spent fuel assembly in order to take burnup credit when loading dry storage casks or to count individual fuel pins in spent fuel assemblies for safeguards purposes. At present, a laboratory prototype imager is under construction. The purpose of this prototype is to demonstrate imaging capability sufficient to resolve individual fuel pins in spent fuel assemblies. This report documents the development of the iterative reconstruction code used to perform tomographic image reconstruction, the imager response calculation used by the reconstruction code, and the results of reconstructions of simulated tomographic imaging measurements for the prototype imager design.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Image reconstruction of shielded mixed-oxide fuel using a dual-particle imaging system

The dual-particle imaging system being developed at the University of Michigan was used at the Joint Research Centre in Ispra, Italy for measurements on samples of special nuclear material. A 1,150-g mixed-oxide (MOX) fuel sample was measured with various shielding configurations to determine how the presence of lead and/or polyethylene shielding degrades the system's ability to localize a source by simultaneous neutron and photon imaging. Three two-hour measurements were taken with the source shielded by an: a) 8-mm lead sheath and 5.1-cm lead bricks, b) 8-mm lead sheath and 6.5-cm polyethylene bricks, and c) 8-mm lead sheath, 5.1-cm lead bricks, and 6.5-cm polyethylene bricks. The 8-mm lead sheath was used in all cases to reduce the measured photon count rate. A “bare” measurement was also made by using only the 8-mm lead sheath, but the unexpected presence of additional sources has rendered the measurement unsuitable for comparison. Here, the resulting images show that the dual-particle imaging system is able to accurately localize the MOX canister in the presence of intervening material.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Sparsity-based photoacoustic image reconstruction with a linear array transducer and direct measurement of the forward model

Photoacoustic (PA) imaging is an emerging imaging technique for many clinical applications. One of the challenges posed by clinical translation is that imaging systems often rely on a finite-aperture transducer rather than a full tomography system. This results in imaging artifacts arising from an underdetermined reconstruction of the initial pressure distribution (IPD). Furthermore, clinical applications often require deep imaging, resulting in a low-signal-to-noise ratio for the acquired signal because of strong light attenuation in tissue. Conventional approaches to reconstruct the IPD, such as back projection and time-reversal, do not adequately suppress the artifacts and noise. We propose a sparsity-based optimization approach that improves the reconstruction of IPD in PA imaging with a linear array ultrasound transducer. In simulation studies, the forward model matrix was measured from k-Wave simulations, and the approach was applied to reconstruct simulated point objects and the Shepp–Logan phantom. The results were compared with the conventional back projection, time-reversal approach, frequency-domain reconstruction, and the iterative least-squares approaches. In experimental studies, the forward model of our imaging system is directly measured by scanning a graphite point source through the imaging field of view. Experimental images of graphite inclusions in tissue-mimicking phantoms are reconstructed and compared with the back projection and iterative least-squares approaches. Altogether these results show that our proposed optimization approach can leverage the sparsity of the PA images to improve the reconstruction of the IPD and outperform the existing popular reconstruction approaches.

47 OTHER INSTRUMENTATION↗

Poster – 02: Positron Emission Tomography (PET) Imaging Reconstruction using higher order Scattered Photon Coincidences

PET images are affected by the presence of scattered photons. Incorrect scatter-correction may cause artifacts, particularly in 3D PET systems. Current scatter reconstruction methods do not distinguish between single and higher order scattered photons. A dual-scattered reconstruction method (GDS-MLEM) that is independent of the number of Compton scattering interactions and less sensitive to the need for high energy resolution detectors, is proposed. To avoid overcorrecting for scattered coincidences, the attenuation coefficient was calculated by integrating the differential Klein-Nishina cross-section over a restricted energy range, accounting only for scattered photons that were not detected. The optimum image can be selected by choosing an energy threshold which is the upper energy limit for the calculation of the cross-section and the lower limit for scattered photons in the reconstruction. Data was simulated using the GATE platform. 500,000 multiple scattered photon coincidences with perfect energy resolution were reconstructed using various methods. The GDS-MLEM algorithm had the highest confidence (98%) in locating the annihilation position and was capable of reconstructing the two largest hot regions. 100,000 photon coincidences, with a scatter fraction of 40%, were used to test the energy resolution dependence of different algorithms. With a 350–650 keV energy window and the restricted attenuation correction model, the GDS-MLEM algorithm was able to improve contrast recovery and reduce the noise by 7.56%–13.24% and 12.4%–24.03%, respectively. This approach is less sensitive to the energy resolution and shows promise if detector energy resolutions of 12% can be achieved.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Hybrid-dual-fourier tomographic algorithm for a fast three-dimensionial optical image reconstruction in turbid media

A reconstruction technique for reducing computation burden in the 3D image processes, wherein the reconstruction procedure comprises an inverse and a forward model. The inverse model uses a hybrid dual Fourier algorithm that combines a 2D Fourier inversion with a 1D matrix inversion to thereby provide high-speed inverse computations. The inverse algorithm uses a hybrid transfer to provide fast Fourier inversion for data of multiple sources and multiple detectors. The forward model is based on an analytical cumulant solution of a radiative transfer equation. The accurate analytical form of the solution to the radiative transfer equation provides an efficient formalism for fast computation of the forward model.

Alfano, Robert R.↗

Investigation of dose reduction in cardiac perfusion SPECT via optimization and choice of the image reconstruction strategy

Background: We investigated the extent to which the administered dose (activity) level can be reduced without sacrificing diagnostic accuracy for three reconstruction strategies for SPECT-myocardial perfusion imaging (MPI). Methods: We optimized the parameters of the three reconstruction strategies for perfusion-defect detection over a range of simulated administered dose levels using a set of hybrid studies (derived from 190 subjects) consisting of clinical SPECT-MPI data modified to contain realistic simulated lesions. The optimized strategies we considered are filtered backprojection (FBP) with no correction for degradations, ordered-subsets expectation-maximization (OS-EM) with attenuation correction (AC), scatter correction (SC), and resolution correction (RC), and OS-EM with scatter and resolution correction only. Each study was evaluated using a total perfusion deficit (TPD) score computed by the Quantitative Perfusion SPECT (QPS) software package. We conducted a receiver operating characteristics (ROC) study based on the TPD scores for each dose level and reconstruction strategy. Results: For FBP, the achieved optimum values of the area under the ROC curve (AUC) at 100%, 50%, 25%, and 12.5% of standard dose were 0.75, 0.74, 0.72, and 0.70, respectively, compared to 0.81, 0.79, 0.76, and 0.74 for OS-EM with AC–SC–RC and 0.78, 0.77, 0.74, 0.72 for OS-EM with SC–RC. Conclusions: Our results suggest that studies reconstructed by OS-EM with AC–SC–RC could possibly be reduced, on average, to 25% of the originally administered dose without causing diagnostic accuracy (AUC) to decrease below that of FBP.

62 RADIOLOGY AND NUCLEAR MEDICINE↗