Showing papers on "Imaging phantom published in 2010"

4D XCAT phantom for multimodality imaging research

Abstract: Purpose: The authors develop the 4D extended cardiac-torso (XCAT) phantom for multimodality imaging research. Methods: Highly detailed whole-body anatomies for the adult male and female were defined in the XCAT using nonuniform rational B-spline (NURBS) and subdivision surfaces based on segmentation of the Visible Male and Female anatomical datasets from the National Library of Medicine as well as patient datasets. Using the flexibility of these surfaces, the Visible Human anatomies were transformed to match body measurements and organ volumes for a 50th percentile (height and weight) male and female. The desired body measurements for the models were obtained using the PEOPLESIZE program that contains anthropometric dimensions categorized from 1st to the 99th percentile for US adults. The desired organ volumes were determined from ICRP Publication 89 [ICRP, ‘‘Basic anatomical and physiological data for use in radiological protection: reference values,” ICRP Publication 89 (International Commission on Radiological Protection, New York, NY, 2002)]. The male and female anatomies serve as standard templates upon which anatomical variations may be modeled in the XCAT through user-defined parameters. Parametrized models for the cardiac and respiratory motions were also incorporated into the XCAT based on high-resolution cardiac- and respiratory-gated multislice CT data. To demonstrate the usefulness of the phantom, the authors show example simulation studies in PET, SPECT, and CT using publicly available simulation packages. Results: As demonstrated in the pilot studies, the 4D XCAT (which includes thousands of anatomical structures) can produce realistic imaging data when combined with accurate models of the imaging process. With the flexibility of the NURBS surface primitives, any number of different anatomies, cardiac or respiratory motions or patterns, and spatial resolutions can be simulated to perform imaging research. Conclusions: With the ability to produce realistic, predictive 3D and 4D imaging data from populations of normal and abnormal patients under various imaging parameters, the authors conclude that the XCAT provides an important tool in imaging research to evaluate and improve imaging devices and techniques. In the field of x-ray CT, the phantom may also provide the necessary foundation with which to optimize clinical CT applications in terms of image quality versus radiation dose, an area of research that is becoming more significant with the growing use of CT.

SPIRiT: Iterative self-consistent parallel imaging reconstruction from arbitrary k-space.

Abstract: A new approach to autocalibrating, coil-by-coil parallel imaging reconstruction, is presented. It is a generalized reconstruction framework based on self-consistency. The reconstruction problem is formulated as an optimization that yields the most consistent solution with the calibration and acquisition data. The approach is general and can accurately reconstruct images from arbitrary k-space sampling patterns. The formulation can flexibly incorporate additional image priors such as off-resonance correction and regularization terms that appear in compressed sensing. Several iterative strategies to solve the posed reconstruction problem in both image and k-space domain are presented. These are based on a projection over convex sets and conjugate gradient algorithms. Phantom and in vivo studies demonstrate efficient reconstructions from undersampled Cartesian and spiral trajectories. Reconstructions that include off-resonance correction and nonlinear l(1)-wavelet regularization are also demonstrated.

Quantitative susceptibility map reconstruction from MR phase data using bayesian regularization: validation and application to brain imaging.

Abstract: The diagnosis of many neurologic diseases benefits from the ability to quantitatively assess iron in the brain. Paramagnetic iron modifies the magnetic susceptibility causing magnetic field inhomogeneity in MRI. The local field can be mapped using the MR signal phase, which is discarded in a typical image reconstruction. The calculation of the susceptibility from the measured magnetic field is an ill-posed inverse problem. In this work, a bayesian regularization approach that adds spatial priors from the MR magnitude image is formulated for susceptibility imaging. Priors include background regions of known zero susceptibility and edge information from the magnitude image. Simulation and phantom validation experiments demonstrated accurate susceptibility maps free of artifacts. The ability to characterize iron content in brain hemorrhage was demonstrated on patients with cavernous hemangioma. Additionally, multiple structures within the brain can be clearly visualized and characterized. The technique introduces a new quantitative contrast in MRI that is directly linked to iron in the brain.

A review of tissue substitutes for ultrasound imaging.

Abstract: The characterization and calibration of ultrasound imaging systems requires tissue-mimicking phantoms with known acoustic properties, dimensions and internal features. Tissue phantoms are available commercially for a range of medical applications. However, commercial phantoms may not be suitable in ultrasound system design or for evaluation of novel imaging techniques. It is often desirable to have the ability to tailor acoustic properties and phantom configurations for specific applications. A multitude of tissue-mimicking materials and phantoms are described in the literature that have been created using a variety of materials and preparation techniques and that have modeled a range of biological systems. This paper reviews ultrasound tissue-mimicking materials and phantom fabrication techniques that have been developed over the past four decades, and describes the benefits and disadvantages of the processes. Both soft tissue and hard tissue substitutes are explored.

MRI-Based Attenuation Correction for PET/MRI Using Ultrashort Echo Time Sequences

Abstract: One of the challenges in PET/MRI is the derivation of an attenuation map to correct the PET image for attenuation. Different methods have been suggested for deriving the attenuation map from an MR image. Because the low signal intensity of cortical bone on images acquired with conventional MRI sequences makes it difficult to detect this tissue type, these methods rely on some sort of anatomic precondition to predict the attenuation map, raising the question of whether these methods will be usable in the clinic when patients may exhibit anatomic abnormalities. Methods: We propose the use of the transverse relaxation rate, derived from images acquired with an ultrashort echo time sequence to classify the voxels into 1 of 3 tissue classes (bone, soft tissue, or air), without making any assumptions on patient anatomy. Each voxel is assigned a linear attenuation coefficient corresponding to its tissue class. A reference CT scan is used to determine the voxel-by-voxel accuracy of the proposed method. The overall accuracy of the MRI-based attenuation correction is evaluated using a method that takes into account the nonlocal effects of attenuation correction. Results: As a proof of concept, the head of a pig was used as a phantom for imaging. The new method yielded a correct tissue classification in 90% of the voxels. Five human brain PET/CT and MRI datasets were also processed, yielding slightly worse voxel-by-voxel performance, compared to a CT-derived attenuation map. The PET datasets were reconstructed using the segmented MRI attenuation map derived with the new method, and the resulting images were compared with segmented CT-based attenuation correction. An average error of around 5% was found in the brain. Conclusion: The feasibility of using the transverse relaxation rate map derived from ultrashort echo time MR images for the estimation of the attenuation map was shown on phantom and clinical brain data. The results indicate that the new method, compared with CT-based attenuation correction, yields clinically acceptable errors. The proposed method does not make any assumptions about patient anatomy and could therefore also be used in cases in which anatomic abnormalities are present.

Robust Super-Resolution Volume Reconstruction From Slice Acquisitions: Application to Fetal Brain MRI

Abstract: Fast magnetic resonance imaging slice acquisition techniques such as single shot fast spin echo are routinely used in the presence of uncontrollable motion These techniques are widely used for fetal magnetic resonance imaging (MRI) and MRI of moving subjects and organs Although high-quality slices are frequently acquired by these techniques, inter-slice motion leads to severe motion artifacts that are apparent in out-of-plane views Slice sequential acquisitions do not enable 3-D volume representation In this study, we have developed a novel technique based on a slice acquisition model, which enables the reconstruction of a volumetric image from multiple-scan slice acquisitions The super-resolution volume reconstruction is formulated as an inverse problem of finding the underlying structure generating the acquired slices We have developed a robust M-estimation solution which minimizes a robust error norm function between the model-generated slices and the acquired slices The accuracy and robustness of this novel technique has been quantitatively assessed through simulations with digital brain phantom images as well as high-resolution newborn images We also report here successful application of our new technique for the reconstruction of volumetric fetal brain MRI from clinically acquired data

Deriving Hounsfield units using grey levels in cone beam computed tomography

Abstract: Objectives: An in vitro study was performed to investigate the relationship between grey levels in dental cone beam CT (CBCT) and Hounsfield units (HU) in CBCT scanners. Methods: A phantom containing 8 different materials of known composition and density was imaged with 11 different dental CBCT scanners and 2 medical CT scanners. The phantom was scanned under three conditions: phantom alone and phantom in a small and large water container. The reconstructed data were exported as Digital Imaging and Communications in Medicine (DICOM) and analysed with On Demand 3D® by Cybermed, Seoul, Korea. The relationship between grey levels and linear attenuation coefficients was investigated. Results: It was demonstrated that a linear relationship between the grey levels and the attenuation coefficients of each of the materials exists at some “effective” energy. From the linear regression equation of the reference materials, attenuation coefficients were obtained for each of the materials and CT numbers in HU were der...

An Assessment of the Impact of Incorporating Time-of-Flight Information into Clinical PET/CT Imaging

Abstract: The introduction of fast scintillators with good stopping power for 511-keV photons has renewed interest in time-of-flight (TOF) PET. The ability to measure the difference between the arrival times of a pair of photons originating from positron annihilation improves the image signal-to-noise ratio (SNR). The level of improvement depends upon the extent and distribution of the positron activity and the time resolution of the PET scanner. While specific estimates can be made for phantom imaging, the impact of TOF PET is more difficult to quantify in clinical situations. The results presented here quantify the benefit of TOF in a challenging phantom experiment and then assess both qualitatively and quantitatively the impact of incorporating TOF information into the reconstruction of clinical studies. A clear correlation between patient body mass index and gain in SNR was observed in this study involving 100 oncology patient studies, with a gain due to TOF ranging from 1.1 to 1.8, which is consistent with the 590-ps time resolution of the TOF PET scanner. The visual comparison of TOF and non-TOF images performed by two nuclear medicine physicians confirmed the advantages of incorporating TOF into the reconstruction, advantages that include better definition of small lesions and image details, improved uniformity, and noise reduction.

Quantitative Accuracy of Clinical 99mTc SPECT/CT Using Ordered-Subset Expectation Maximization with 3-Dimensional Resolution Recovery, Attenuation, and Scatter Correction

Abstract: We present a calibration method of a clinical SPECT/CT device for quantitative 99mTc SPECT. We use a commercially available reconstruction package including ordered-subset expectation maximization (OSEM) with depth-dependent 3-dimensional resolution recovery (OSEM-3D), CT-based attenuation correction, and scatter correction. We validated the method in phantom studies and applied it to images from patients injected with 99mTc-diphosponate. Methods: The following 3 steps were performed to derive absolute quantitative values from SPECT reconstructed images. In step 1, we used simulations to characterize the SPECT/CT system and derive emission recovery values for various imaging parameter settings. We simulated spheres of varying diameters and focused on the dependencies of activity estimation errors on structure size and position, pixel size, count density, and reconstruction parameters. In step 2, we cross-calibrated our clinical SPECT/CT system with the well counter using a large cylinder phantom. This step provided the mapping from image counts to kBq/mL. And in step 3, correction factors from steps 1 and 2 were applied to reconstructed images. We used a cylinder phantom with variable-sized spheres for verification of the method. For in vivo validation, SPECT/CT datasets from 16 patients undergoing 99mTc-diphosponate SPECT/CT examinations of the pelvis including the bladder were acquired. The radioactivity concentration in the patients9 urine served as the gold standard. Mean quantitative accuracy and SEs were calculated. Results: In the phantom experiments, the mean accuracy in quantifying radioactivity concentration in absolute terms was within 3.6% (SE, 8.0%), with a 95% confidence interval between −19.4% and +12.2%. In the patient studies, the mean accuracy was within 1.1% (SE, 8.4%), with a 95% confidence interval between −15.4% and +17.5%. Conclusion: Current commercially available SPECT/CT technology using OSEM-3D reconstruction, scatter correction, and CT-based attenuation correction allows quantification of 99mTc radioactivity concentration in absolute terms within 3.6% in phantoms and 1.1% in patients with a focus on the bladder. This opens up the opportunity of SPECT quantitation entering the routine clinical arena. Still, the imprecision caused by unavoidable measurement errors is a dominant factor for absolute quantitation in a clinical setup.

Nonlinear Regularization for Per Voxel Estimation of Magnetic Susceptibility Distributions From MRI Field Maps

Abstract: Magnetic susceptibility is an important physical property of tissues, and can be used as a contrast mechanism in magnetic resonance imaging (MRI). Recently, targeting contrast agents by conjugation with signaling molecules and labeling stem cells with contrast agents have become feasible. These contrast agents are strongly paramagnetic, and the ability to quantify magnetic susceptibility could allow accurate measurement of signaling and cell localization. Presented here is a technique to estimate arbitrary magnetic susceptibility distributions by solving an ill-posed inversion problem from field maps obtained in an MRI scanner. Two regularization strategies are considered: conventional Tikhonov regularization and a sparsity promoting nonlinear regularization using the l 1 norm. Proof of concept is demonstrated using numerical simulations, phantoms, and in a stroke model in a rat. Initial experience indicates that the nonlinear regularization better suppresses noise and streaking artifacts common in susceptibility estimation.

GPU-based Low Dose CT Reconstruction via Edge-preserving Total Variation Regularization

Abstract: High radiation dose in CT scans increases a lifetime risk of cancer and has become a major clinical concern. Recently, iterative reconstruction algorithms with Total Variation (TV) regularization have been developed to reconstruct CT images from highly undersampled data acquired at low mAs levels in order to reduce the imaging dose. Nonetheless, TV regularization may lead to over-smoothed images and lost edge information. To solve this problem, in this work we develop an iterative CT reconstruction algorithm with edge-preserving TV regularization to reconstruct CT images from highly undersampled data obtained at low mAs levels. The CT image is reconstructed by minimizing an energy consisting of an edge-preserving TV norm and a data fidelity term posed by the x-ray projections. The edge-preserving TV term is proposed to preferentially perform smoothing only on non-edge part of the image in order to avoid over-smoothing, which is realized by introducing a penalty weight to the original total variation norm. Our iterative algorithm is implemented on GPU to improve its speed. We test our reconstruction algorithm on a digital NCAT phantom, a physical chest phantom, and a Catphan phantom. Reconstruction results from a conventional FBP algorithm and a TV regularization method without edge preserving penalty are also presented for comparison purpose. The experimental results illustrate that both TV-based algorithm and our edge-preserving TV algorithm outperform the conventional FBP algorithm in suppressing the streaking artifacts and image noise under the low dose context. Our edge-preserving algorithm is superior to the TV-based algorithm in that it can preserve more information of fine structures and therefore maintain acceptable spatial resolution.

Microwave Radar-Based Differential Breast Cancer Imaging: Imaging in Homogeneous Breast Phantoms and Low Contrast Scenarios

Abstract: This paper presents an improved antenna array for radar-based breast cancer imaging. The improvement was achieved by increasing the number of antennas in the array to 31 elements, as well as by improving the antenna design itself. Using an experimental setup, with homogeneous curved breast phantoms, we have demonstrated substantial imaging improvement with the new antenna array. The new system is also able to detect 7 mm-diameter tumor phantoms in any location within the breast, even as close as 4 mm from the skin layer. Additionally, we have shown good imaging results in low-contrast scenarios, where the dielectric contrast between tumor and normal tissue was reduced to 2:1. Presented results clearly demonstrate the large impact of antenna's characteristics on imaging performance.

Reducing the Radiation Dose for CT Colonography Using Adaptive Statistical Iterative Reconstruction: A Pilot Study

Abstract: OBJECTIVE. The purpose of our study was to evaluate the feasibility of preserving image quality during CT colonography (CTC) using a reduced radiation dose with adaptive statistical iterative reconstruction (ASIR).MATERIALS AND METHODS. A proven colon phantom was imaged at standard dose settings (50 mAs) and at reduced doses (10–40 mAs) using six different ASIR levels (0–100%). We assessed 2D and 3D image quality and noise to determine the optimal dose and ASIR setting. Eighteen patients were then scanned with a standard CTC dose (50 mAs) in the supine position and at a reduced dose of 25 mAs with 40% ASIR in the prone position. Three radiologists blinded to the scanning techniques assessed 2D and 3D image quality and noise at three different colon locations. A score difference of ≥ 1 was considered clinically important. Actual noise measures were compared between the standard-dose and low-dose acquisitions.RESULTS. The phantom study showed image noise reduction that correlated with a higher percentage of...

Real-Time 2-D Temperature Imaging Using Ultrasound

Abstract: We have previously introduced methods for noninvasive estimation of temperature change using diagnostic ultrasound. The basic principle was validated both in vitro and in vivo by several groups worldwide. Some limitations remain, however, that have prevented these methods from being adopted in monitoring and guidance of minimally invasive thermal therapies, e.g., RF ablation and high-intensity-focused ultrasound (HIFU). In this letter, we present first results from a real-time system for 2-D imaging of temperature change using pulse-echo ultrasound. The front end of the system is a commercially available scanner equipped with a research interface, which allows the control of imaging sequence and access to the RF data in real time. A high-frame-rate 2-D RF acquisition mode, M2D, is used to capture the transients of tissue motion/deformations in response to pulsed HIFU. The M2D RF data is streamlined to the back end of the system, where a 2-D temperature imaging algorithm based on speckle tracking is implemented on a graphics processing unit. The real-time images of temperature change are computed on the same spatial and temporal grid of the M2D RF data, i.e., no decimation. Verification of the algorithm was performed by monitoring localized HIFU-induced heating of a tissue-mimicking elastography phantom. These results clearly demonstrate the repeatability and sensitivity of the algorithm. Furthermore, we present in vitro results demonstrating the possible use of this algorithm for imaging changes in tissue parameters due to HIFU-induced lesions. These results clearly demonstrate the value of the real-time data streaming and processing in monitoring, and guidance of minimally invasive thermotherapy.

Cerenkov luminescence tomography for small-animal imaging

Abstract: Cerenkov radiation is a well-known phenomenon in which optical photons are emitted by charged particles moving faster than the speed of light in a medium. We have observed Cerenkov photons emitted from β-emitting radiotracers such as F18-fluorodeoxyglucose using a sensitive CCD camera. Phantom and in vivo mouse imaging experiments have demonstrated that surface measurements of the emitted Cerenkov optical photons could be used to reconstruct the radiotracer activity distribution inside an object by modeling the optical photon propagation with the diffusion equation and reconstructing the optical emission source distribution iteratively with a preconditioned conjugate gradient method.

Optimal radiofrequency and gradient spoiling for improved accuracy of T1 and B1 measurements using fast steady-state techniques.

Abstract: Variable flip angle T(1) mapping and actual flip-angle imaging B(1) mapping are widely used quantitative MRI methods employing radiofrequency spoiled gradient-echo pulse sequences. Incomplete elimination of the transverse magnetization in these sequences has been found to be a critical source of T(1) and B(1) measurement errors. In this study, comprehensive theoretical analysis of spoiling-related errors in variable flip angle and actual flip-angle imaging methods was performed using the combined isochromat summation and diffusion propagator model and validated by phantom experiments. The key theoretical conclusion is that correct interpretation of spoiling phenomena in fast gradient-echo sequences requires accurate consideration of the diffusion effect. A general strategy for improvement of T(1) and B(1) measurement accuracy was proposed based on the strong spoiling regimen, where diffusion-modulated spatial averaging of isochromats becomes a dominant factor determining magnetization evolution. Practical implementation of strongly spoiled variable flip angle and actual flip-angle imaging techniques requires sufficiently large spoiling gradient areas (A(G)) in combination with optimal radiofrequency phase increments (phi(0)). Optimal regimens providing <2% relative T(1) and B(1) measurement errors in a variety of tissues were theoretically derived for prospective in vivo variable flip angle (pulse repetition time = 15-20 ms, A(G) = 280-450 mT.ms/m, phi(0) = 169 degrees) and actual flip-angle imaging (pulse repetition time(1)/pulse repetition time(2) = 20/100 ms, A(G1)/A(G2) = 450/2250 mT.ms/m, phi(0) = 39 degrees) applications based on 25 mT/m maximal available gradient strength.

Empirical beam hardening correction (EBHC) for CT.

Abstract: Purpose: Due to x-ray beam polychromaticity and scattered radiation, attenuation measurements tend to be underestimated. Cupping and beam hardening artifacts become apparent in the reconstructed CT images. If only one material such as water, for example, is present, these artifacts can be reduced by precorrecting the rawdata. Higher order beam hardening artifacts, as they result when a mixture of materials such as water and bone, or water and bone and iodine is present, require an iterative beam hardening correction where the image is segmented into different materials and those are forward projected to obtain new rawdata. Typically, the forward projection must correctly model the beam polychromaticity and account for all physical effects, including the energy dependence of the assumed materials in the patient, the detector response, and others. We propose a new algorithm that does not require any knowledge about spectra or attenuation coefficients and that does not need to be calibrated. The proposed method corrects beam hardening in single energy CT data. Methods: The only a priori knowledge entering EBHC is the segmentation of the object into different materials. Materials other than water are segmented from the original image, e.g., by using simple thresholding. Then, a (monochromatic) forward projectionmore » of these other materials is performed. The measured rawdata and the forward projected material-specific rawdata are monomially combined (e.g., multiplied or squared) and reconstructed to yield a set of correction volumes. These are then linearly combined and added to the original volume. The combination weights are determined to maximize the flatness of the new and corrected volume. EBHC is evaluated using data acquired with a modern cone-beam dual-source spiral CT scanner (Somatom Definition Flash, Siemens Healthcare, Forchheim, Germany), with a modern dual-source micro-CT scanner (TomoScope Synergy Twin, CT Imaging GmbH, Erlangen, Germany), and with a modern C-arm CT scanner (Axiom Artis dTA, Siemens Healthcare, Forchheim, Germany). A large variety of phantom, small animal, and patient data were used to demonstrate the data and system independence of EBHC. Results: Although no physics apart from the initial segmentation procedure enter the correction process, beam hardening artifacts were significantly reduced by EBHC. The image quality for clinical CT, micro-CT, and C-arm CT was highly improved. Only in the case of C-arm CT, where high scatter levels and calibration errors occur, the relative improvement was smaller. Conclusions: The empirical beam hardening correction is an interesting alternative to conventional iterative higher order beam hardening correction algorithms. It does not tend to over- or undercorrect the data. Apart from the segmentation step, EBHC does not require assumptions on the spectra or on the type of material involved. Potentially, it can therefore be applied to any CT image.« less

X-ray fluorescence computed tomography (XFCT) imaging of gold nanoparticle-loaded objects using 110 kVp x-rays

Abstract: A conventional x-ray fluorescence computed tomography (XFCT) technique requires monochromatic synchrotron x-rays to simultaneously determine the spatial distribution and concentration of various elements such as metals in a sample. However, the synchrotron-based XFCT technique appears to be unsuitable for in vivo imaging under a typical laboratory setting. In this study we demonstrated, for the first time to our knowledge, the possibility of performing XFCT imaging of a small animal-sized object containing gold nanoparticles (GNPs) at relatively low concentrations using polychromatic diagnostic energy range x-rays. Specifically, we created a phantom made of polymethyl methacrylate plastic containing two cylindrical columns filled with saline solution at 1 and 2 wt% GNPs, respectively, mimicking tumors/organs within a small animal. XFCT scanning of the phantom was then performed using microfocus 110 kVp x-ray beam and cadmium telluride (CdTe) x-ray detector under a pencil beam geometry after proper filtering of the x-ray beam and collimation of the detector. The reconstructed images clearly identified the locations of the two GNP-filled columns with different contrast levels directly proportional to gold concentration levels. On the other hand, the current pencil-beam implementation of XFCT is not yet practical for routine in vivo imaging tasks with GNPs, especially in terms of scanning time. Nevertheless, with the use of multiple detectors and a limited number of projections, it may still be used to image some objects smaller than the current phantom size. The current investigation suggests several modification strategies of the current XFCT setup, such as the adoption of the quasi-monochromatic cone/fan x-ray beam and XFCT-specific spatial filters or pinhole detector collimators, in order to establish the ultimate feasibility of a bench-top XFCT system for GNP-based preclinical molecular imaging applications.

Ultrasound-Guided Robot for Flexible Needle Steering

Abstract: The success rate of medical procedures involving needle insertion is often directly related to needle placement accuracy. Due to inherent limitations of commonly used freehand needle placement techniques, there is a need for a system providing for controlled needle steering for procedures that demand high positional accuracy. This paper describes a robotic system developed for flexible needle steering inside soft tissues under real-time ultrasound imaging. An inverse kinematics algorithm based on a virtual spring model is applied to calculate needle base manipulations required for the tip to follow a curved trajectory while avoiding physiological obstacles. The needle tip position is derived from ultrasound images and is used in calculations to minimize the tracking error, enabling a closed-loop needle insertion. In addition, as tissue stiffness is a necessary input to the control algorithm, a novel method to classify tissue stiffness from localized tissue displacements is proposed and shown to successfully distinguish between soft and stiff tissue. The system performance was experimentally verified by robotic manipulation of the needle base inside a phantom with layers of varying stiffnesses. The closed-loop experiment with updated tissue stiffness parameters demonstrated a needle-tip tracking error of ~ 1 mm and proved to be significantly more accurate than the freehand method.

Image-based point spread function implementation in a fully 3D OSEM reconstruction algorithm for PET

Abstract: The interest in positron emission tomography (PET) and particularly in hybrid integrated PET/CT systems has significantly increased in the last few years due to the improved quality of the obtained images. Nevertheless, one of the most important limits of the PET imaging technique is still its poor spatial resolution due to several physical factors originating both at the emission (e.g. positron range, photon non-collinearity) and at detection levels (e.g. scatter inside the scintillating crystals, finite dimensions of the crystals and depth of interaction). To improve the spatial resolution of the images, a possible way consists of measuring the point spread function (PSF) of the system and then accounting for it inside the reconstruction algorithm. In this work, the system response of the GE Discovery STE operating in 3D mode has been characterized by acquiring (22)Na point sources in different positions of the scanner field of view. An image-based model of the PSF was then obtained by fitting asymmetric two-dimensional Gaussians on the (22)Na images reconstructed with small pixel sizes. The PSF was then incorporated, at the image level, in a three-dimensional ordered subset maximum likelihood expectation maximization (OS-MLEM) reconstruction algorithm. A qualitative and quantitative validation of the algorithm accounting for the PSF has been performed on phantom and clinical data, showing improved spatial resolution, higher contrast and lower noise compared with the corresponding images obtained using the standard OS-MLEM algorithm.

Fast 3D chemical exchange saturation transfer (CEST) imaging of the human brain.

Abstract: Chemical exchange saturation transfer magnetic resonance imaging can detect low-concentration compounds with exchangeable protons through saturation transfer to water. This technique is generally slow, as it requires acquisition of saturation images at multiple frequencies. In addition, multislice imaging is complicated by saturation effects differing from slice to slice because of relaxation losses. In this study, a fast three-dimensional chemical exchange saturation transfer imaging sequence is presented that allows whole-brain coverage for a frequency-dependent saturation spectrum (z-spectrum, 26 frequencies) in less than 10 min. The approach employs a three-dimensional gradient- and spin-echo readout using a prototype 32-channel phased-array coil, combined with two-dimensional sensitivity encoding accelerations. Results from a homogenous protein-containing phantom at 3T show that the sequence produced a uniform contrast across all slices. To show translational feasibility, scans were also performed on five healthy human subjects. Results for chemical exchange saturation transfer images at 3.5 ppm downfield of the water resonance, so-called amide proton transfer images, show that lipid signals are sufficiently suppressed and artifacts caused by B(0) inhomogeneity can be removed in postprocessing. The scan time and image quality of these in vivo results show that three-dimensional chemical exchange saturation transfer MRI using gradient- and spin-echo acquisition is feasible for whole-brain chemical exchange saturation transfer studies at 3T in a clinical time frame.

Optical phantoms of varying geometry based on thin building blocks with controlled optical properties

Abstract: Current innovations in optical imaging, measurement tech- niques, and data analysis algorithms express the need for reliable test- ing and comparison methods. We present the design and character- ization of silicone elastomer-based optical phantoms. Absorption is included by adding a green dye and scattering by adding TiO2 or SiO2 particles. Optical coherence tomography measurements demon- strate a linear dependence of the attenuation coefficient with scatterer concentration in the absence of absorbers. Optical transmission spec- troscopy of the nonscattering absorbing phantoms shows a linear con- centration dependent absorption coefficient. Both types of samples are stable over a period of 6 months. Confocal microscopy of the samples demonstrates a homogeneous distribution of the scatterers, albeit with some clustering. Based on layers with thicknesses as small as 50 m, we make multifaceted structures resembling flow channels, wavy skin-like structures, and a layered and curved phantom resem- bling the human retina. Finally, we demonstrate the ability to incor- porate gold nanoparticles within the phantoms. In conclusion, our phantoms are easy to make, are based on affordable materials, exhibit well-defined and controllable thickness, refractive index, absorption, and scattering coefficients, are homogeneous, and allow the incorpo- ration of novel types of nanoparticle contrast agents. We believe our phantoms fulfill many of the requirements for an "ideal" tissue phan- tom, and will be particularly suited for novel optical coherence to- mography applications. © 2010 Society of Photo-Optical Instrumentation Engineers. DOI: 10.1117/1.3369003

Shading correction for on-board cone-beam CT in radiation therapy using planning MDCT images.

Abstract: Purpose: Applications of cone-beam CT(CBCT) to image-guided radiation therapy (IGRT) are hampered by shading artifacts in the reconstructed images. These artifacts are mainly due to scatter contamination in the projections but also can result from uncorrected beam hardening effects as well as nonlinearities in responses of the amorphous silicon flat panel detectors. While currently, CBCT is mainly used to provide patient geometry information for treatment setup, more demanding applications requiring high-quality CBCTimages are under investigation. To tackle these challenges, many CBCT correction algorithms have been proposed; yet, a standard approach still remains unclear. In this work, we propose a shading correction method for CBCT that addresses artifacts from low-frequency projection errors. The method is consistent with the current workflow of radiation therapy. Methods: With much smaller inherent scatter signals and more accurate detectors, diagnostic multidetector CT (MDCT) provides high quality CTimages that are routinely used for radiation treatment planning. Using the MDCT image as “free” prior information, we first estimate the primary projections in the CBCT scan via forward projection of the spatially registered MDCT data. Since most of the CBCT shading artifacts stem from low-frequency errors in the projections such as scatter, these errors can be accurately estimated by low-pass filtering the difference between the estimated and raw CBCT projections. The error estimates are then subtracted from the raw CBCT projections. Our method is distinct from other published correction methods that use the MDCT image as a prior because it is projection-based and uses limited patient anatomical information from the MDCT image. The merit of CBCT-based treatment monitoring is therefore retained. Results: The proposed method is evaluated using two phantom studies on tabletop systems. On the Catphan©600 phantom, our approach reduces the reconstruction error from 348 Hounsfield unit (HU) without correction to 4 HU around the object center after correction, and from 375 HU to 17 HU in the high-contrast regions. In the selected regions of interest (ROIs), the average image contrast is increased by a factor of 3.3. When noise suppression is implemented, the proposed correction substantially improves the contrast-to-noise ratio(CNR) and therefore the visibility of low-contrast objects, as seen in a more challenging pelvis phantom study. Besides a significant improvement in image uniformity, a low-contrast object of ∼ 25 HU , which is otherwise buried in the shading artifacts, can be clearly identified after the proposed correction due to a CNR increase of 3.1. Compared to a kernel-based scatter correction method coupled with an analytical beam hardening correction, our approach also shows an overall improved performance with some residual artifacts. Conclusions: By providing effective shading correction, our approach has the potential to improve the accuracy of more advanced CBCT-based clinical applications for IGRT, such as tumor delineation and dose calculation.

Development and validation of segmentation and interpolation techniques in sinograms for metal artifact suppression in CT

Abstract: Purpose: Metal prostheses cause artifacts in computed tomography(CT)images. The purpose of this work was to design an efficient and accurate metal segmentation in raw data to achieve artifact suppression and to improve CTimage quality for patients with metal hip or shoulder prostheses. Methods: The artifact suppression technique incorporates two steps: metal object segmentation in raw data and replacement of the segmented region by new values using an interpolation scheme, followed by addition of the scaled metal signal intensity. Segmentation of metal is performed directly in sinograms, making it efficient and different from current methods that perform segmentation in reconstructed images in combination with Radon transformations. Metal signal segmentation is achieved by using a Markov random field model (MRF). Three interpolation methods are applied and investigated. To provide a proof of concept, CT data of five patients with metal implants were included in the study, as well as CT data of a PMMA phantom with Teflon, PVC, and titanium inserts. Accuracy was determined quantitatively by comparing mean Hounsfield (HU) values and standard deviation (SD) as a measure of distortion in phantom images with titanium (original and suppressed) and without titanium insert. Qualitative improvement was assessed by comparing uncorrected clinical images with artifact suppressed images. Results: Artifacts in CT data of a phantom and five patients were automatically suppressed. The general visibility of structures clearly improved. In phantom images, the technique showed reduced SD close to the SD for the case where titanium was not inserted, indicating improved image quality. HU values in corrected images were different from expected values for all interpolation methods. Subtle differences between interpolation methods were found. Conclusions: The new artifact suppression design is efficient, for instance, in terms of preserving spatial resolution, as it is applied directly to original raw data. It successfully reduced artifacts in CTimages of five patients and in phantom images. Sophisticated interpolation methods are needed to obtain reliable HU values close to the prosthesis.

Real-time 3D target tracking in MRI guided focused ultrasound ablations in moving tissues.

Abstract: Magnetic resonance imaging-guided high intensity focused ultrasound is a promising method for the noninvasive ablation of pathological tissue in abdominal organs such as liver and kidney. Due to the high perfusion rates of these organs, sustained sonications are required to achieve a sufficiently high temperature elevation to induce necrosis. However, the constant displacement of the target due to the respiratory cycle render continuous ablations challenging, since dynamic repositioning of the focal point is required. This study demonstrates subsecond 3D high intensity focused ultrasound-beam steering under magnetic resonance-guidance for the real-time compensation of respiratory motion. The target is observed in 3D space by coupling rapid 2D magnetic resonance-imaging with prospective slice tracking based on pencil-beam navigator echoes. The magnetic resonance-data is processed in real-time by a computationally efficient reconstruction pipeline, which provides the position, the temperature and the thermal dose on-the-fly, and which feeds corrections into the high intensity focused ultrasound-ablator. The effect of the residual update latency is reduced by using a 3D Kalman-predictor for trajectory anticipation. The suggested method is characterized with phantom experiments and verified in vivo on porcine kidney. The results show that for update frequencies of more than 10 Hz and latencies of less then 114 msec, temperature elevations can be achieved, which are comparable to static experiments.

Block matching 3D random noise filtering for absorption optical projection tomography

Abstract: Absorption and emission optical projection tomography (OPT), alternatively referred to as optical computed tomography (optical-CT) and optical-emission computed tomography (optical-ECT), are recently developed three-dimensional imaging techniques with value for developmental biology and ex vivo gene expression studies. The techniques' principles are similar to the ones used for x-ray computed tomography and are based on the approximation of negligible light scattering in optically cleared samples. The optical clearing is achieved by a chemical procedure which aims at substituting the cellular fluids within the sample with a cell membranes' index matching solution. Once cleared the sample presents very low scattering and is then illuminated with a light collimated beam whose intensity is captured in transillumination mode by a CCD camera. Different projection images of the sample are subsequently obtained over a 360° full rotation, and a standard backprojection algorithm can be used in a similar fashion as for x-ray tomography in order to obtain absorption maps. Because not all biological samples present significant absorption contrast, it is not always possible to obtain projections with a good signal-to-noise ratio, a condition necessary to achieve high-quality tomographic reconstructions. Such is the case for example, for early stage's embryos. In this work we demonstrate how, through the use of a random noise removal algorithm, the image quality of the reconstructions can be considerably improved even when the noise is strongly present in the acquired projections. Specifically, we implemented a block matching 3D (BM3D) filter applying it separately on each acquired transillumination projection before performing a complete three-dimensional tomographical reconstruction. To test the efficiency of the adopted filtering scheme, a phantom and a real biological sample were processed. In both cases, the BM3D filter led to a signal-to-noise ratio increment of over 30 dB on severe noise-affected reconstructions revealing original—noise-hidden—image details. These results show the utility of the BM3D approach for OPT under typical conditions of very low light absorption, suggesting its implementation as an efficient alternative to other filtering schemes such as for example the median filter.

Anatomical background and generalized detectability in tomosynthesis and cone-beam CT

Abstract: Purpose: Anatomical background presents a major impediment to detectability in 2D radiography as well as 3D tomosynthesis and cone-beam CT (CBCT). This article incorporates theoretical and experimental analysis of anatomical background ''noise'' in cascaded systems analysis of 2D and 3D imaging performance to yield ''generalized'' metrics of noise-equivalent quanta (NEQ) and detectability index as a function of the orbital extent of the (circular arc) source-detector orbit. Methods: A physical phantom was designed based on principles of fractal self-similarity to exhibit power-law spectral density ({kappa}/f{sup {beta}}) comparable to various anatomical sites (e.g., breast and lung). Background power spectra [S{sub B}(f)] were computed as a function of source-detector orbital extent, including tomosynthesis ({approx}10 deg. - 180 deg.) and CBCT (180 deg. +fan to 360 deg.) under two acquisition schemes: (1) Constant angular separation between projections (variable dose) and (2) constant total number of projections (constant dose). The resulting S{sub B} was incorporated in the generalized NEQ, and detectability index was computed from 3D cascaded systems analysis for a variety of imaging tasks. Results: The phantom yielded power-law spectra within the expected spatial frequency range, quantifying the dependence of clutter magnitude ({kappa}) and correlation ({beta}) with increasing tomosynthesis angle. Incorporation of S{sub B}more » in the 3D NEQ provided a useful framework for analyzing the tradeoffs among anatomical, quantum, and electronic noise with dose and orbital extent. Distinct implications are posed for breast and chest tomosynthesis imaging system design--applications varying significantly in {kappa} and {beta}, and imaging task and, therefore, in optimal selection of orbital extent, number of projections, and dose. For example, low-frequency tasks (e.g., soft-tissue masses or nodules) tend to benefit from larger orbital extent and more fully 3D tomographic imaging, whereas high-frequency tasks (e.g., microcalcifications) require careful, application-specific selection of orbital extent and number of projections to minimize negative effects of quantum and electronic noise. Conclusions: The complex tradeoffs among anatomical background, quantum noise, and electronic noise in projection imaging, tomosynthesis, and CBCT can be described by generalized cascaded systems analysis, providing a useful framework for system design and optimization.« less

The imaging performance of a LaBr3-based PET scanner

Abstract: A prototype time-of-flight (TOF) PET scanner based on cerium-doped lanthanum bromide [LaBr3 (5% Ce)] has been developed. LaBr3 has a high light output, excellent energy resolution and fast timing properties that have been predicted to lead to good image quality. Intrinsic performance measurements of spatial resolution, sensitivity and scatter fraction demonstrate good conventional PET performance; the results agree with previous simulation studies. Phantom measurements show the excellent image quality achievable with the prototype system. Phantom measurements and corresponding simulations show a faster and more uniform convergence rate, as well as more uniform quantification, for TOF reconstruction of the data, which have 375 ps intrinsic timing resolution, compared to non-TOF images. Measurements and simulations of a hot and cold sphere phantom show that the 7% energy resolution helps to mitigate residual errors in the scatter estimate because a high energy threshold (>480 keV) can be used to restrict the amount of scatter accepted without a loss of true events. Preliminary results with incorporation of a model of detector blurring in the iterative reconstruction algorithm not only show improved contrast recovery but also point out the importance of an accurate resolution model of the tails of LaBr3's point spread function. The LaBr3 TOF-PET scanner demonstrated the impact of superior timing and energy resolutions on image quality.

Evaluation of the attenuation properties of MR equipment for its use in a whole-body PET/MR scanner

Abstract: The combination of magnetic resonance imaging (MR) and positron emission tomography (PET) scanners can provide a powerful tool for clinical diagnosis and investigation. Among the challenges of developing a combined scanner, obtaining attenuation maps for PET reconstruction is of critical importance. This requires accounting for the presence of MR hardware in the field of view. The attenuation introduced by this hardware cannot be obtained from MR data. We propose the creation of attenuation models of MR hardware, to be registered into the MR-based attenuation map prior to PET reconstruction. Two steps were followed to assess the viability of this method. First, transmission and emission measurements were performed on MR components (RF coils and medical probes). The severity of the artifacts in the reconstructed PET images was evaluated. Secondly, a high-exposure computed tomography (CT) scan was used to obtain a model of a head coil. This model was registered into the attenuation map of PET/CT scans of a uniform phantom fitted with the coil. The resulting PET images were compared to the PET/CT reconstruction in the absence of coils. The artifacts introduced by misregistration of the model were studied. The transmission scans revealed 17% count loss due to the presence of head and neck coils in the field of view. Important sources of attenuation were found in the lock, signal cables and connectors. However, the worst source of attenuation was the casing between both coils. None of the measured medical probes introduced a significant amount of attenuation. Concerning the attenuation model of the head coil, reconstructed PET images with model-based correction were comparable to the reference PET/CT reconstruction. However, inaccuracies greater than 1-2 mm in the axial positioning of the model led to important artifacts. In conclusion, the results show that model-based attenuation correction is possible. Using a high-exposure scan to create an attenuation model of the coils has been proved feasible. However, adequate registration of the model is mandatory.

Patient-specific radiation dose and cancer risk estimation in CT: part I. development and validation of a Monte Carlo program.

Abstract: Purpose: Radiation-dose awareness and optimization in CT can greatly benefit from a dose-reporting system that provides dose and risk estimates specific to each patient and each CT examination. As the first step toward patient-specific dose and risk estimation, this article aimed to develop a method for accurately assessing radiation dose from CT examinations. Methods: A Monte Carlo program was developed to model a CT system (LightSpeed VCT, GE Healthcare). The geometry of the system, the energy spectra of the x-ray source, the three-dimensional geometry of the bowtie filters, and the trajectories of source motions during axial and helical scans were explicitly modeled. To validate the accuracy of the program, a cylindrical phantom was built to enable dose measurements at seven different radial distances from its central axis. Simulated radial dose distributions in the cylindrical phantom were validated against ion chamber measurements for single axial scans at all combinations of tube potential and bowtie filter settings. The accuracy of the program was further validated using two anthropomorphic phantoms (a pediatric one-year-old phantom and an adult female phantom). Computer models of the two phantoms were created based on their CT data and were voxelized for input into the Monte Carlo program. Simulated dose at various organ locations was compared against measurements made with thermoluminescent dosimetry chips for both single axial and helical scans. Results: For the cylindrical phantom, simulations differed from measurements by −4.8% to 2.2%. For the two anthropomorphic phantoms, the discrepancies between simulations and measurements ranged between (−8.1%, 8.1%) and (−17.2%, 13.0%) for the single axial scans and the helical scans, respectively. Conclusions: The authors developed an accurate Monte Carlo program for assessing radiation dose from CT examinations. When combined with computer models of actual patients, the program can provide accurate dose estimates for specific patients.