This PDF file contains the front matter associated with SPIE Proceedings Volume 9783, including the Title Page, Copyright information, Table of Contents, and Conference Committee listing.

A conventional 3D/4D digital subtraction angiogram (DSA) requires two rotational acquisitions (mask and fill) to compute the log-subtracted projections that are used to reconstruct a 3D/4D volume. Since all of the vascular information is contained in the fill acquisition, it is hypothesized that it is possible to reduce the x-ray dose of the mask acquisition substantially and still obtain subtracted projections adequate to reconstruct a 3D/4D volume with noise level comparable to a full dose acquisition. A full dose mask and fill acquisition were acquired from a clinical study to provide a known full dose reference reconstruction. Gaussian noise was added to the mask acquisition to simulate a mask acquisition acquired at 10% relative dose. Noise in the low-dose mask projections was reduced with a weighted edge preserving (WEP) filter designed to preserve bony edges while suppressing noise. 2D log-subtracted projections were computed from the filtered low-dose mask and full-dose fill projections, and then 3D/4D-DSA reconstruction algorithms were applied. Additional bilateral filtering was applied to the 3D volumes. The signal-to-noise ratio measured in the filtered 3D/4D-DSA volumes was compared to the full dose case. The average ratio of filtered low-dose SNR to full-dose SNR was 1.07 for the 3D-DSA and 1.05 for the 4D-DSA, indicating the method is a feasible approach to restoring SNR in DSA scans acquired with a low-dose mask. The method was also tested in a phantom study with full dose fill and 22% dose mask.

Multiplanar reconstruction (MPR) in digital breast tomosynthesis (DBT) allows tomographic images to be portrayed in various orientations. We have conducted research to determine the resolution of tomosynthesis MPR. We built a phantom that houses a star test pattern to measure resolution. This phantom provides three rotational degrees of freedom. The design consists of two hemispheres with longitudinal and latitudinal grooves that reference angular increments. When joined together, the hemispheres form a dome that sits inside a cylindrical encasement. The cylindrical encasement contains reference notches to match the longitudinal and latitudinal grooves that guide the phantom’s rotations. With this design, any orientation of the star-pattern can be analyzed. Images of the star-pattern were acquired using a DBT mammography system at the Hospital of the University of Pennsylvania. Images taken were reconstructed and analyzed by two different methods. First, the maximum visible frequency (in line pairs per millimeter) of the star test pattern was measured. Then, the contrast was calculated at a fixed spatial frequency. These analyses confirm that resolution decreases with tilt relative to the breast support. They also confirm that resolution in tomosynthesis MPR is dependent on object orientation. Current results verify that the existence of super-resolution depends on the orientation of the frequency; the direction parallel to x-ray tube motion shows super-resolution. In conclusion, this study demonstrates that the direction of the spatial frequency relative to the motion of the x-ray tube is a determinant of resolution in MPR for DBT.

We have developed a clinically ready first generation stationary breast tomosynthesis system (s-DBT). In the s-DBT system, focal spot blur associated with x-ray source motion is completely eliminated, allowing for rapid acquisition of projection images over a larger angular span without changing the acquisition time.

In the phantom studies the 1st generation s‐DBT system has demonstrated 30% higher spatial resolution than the corresponding continuous motion DBT systems. The system is currently being evaluated for its diagnostic performance in 100 patient clinical evaluation against FFDM. Initial results indicate that the s­‐DBT system can produce increased lesion conspicuity and comparable MC visibility. However due to x­‐ray flux limitations, certain large size patients have to be excluded. Recent studies have shown that increasing the angular span beyond 30° can be beneficial for enhanced depth resolution. We report the preliminary characterization of the 2nd generation s-­DBT system with a new CNT x-­ray source array, increased tube flux and a larger angular span. Increasing x‐ray tube flux allows for a larger patient population and dual energy imaging. Results indicate that the system delivers more than twice the flux, allowing for imaging of all size patients with acquisition time of 2­‐4 seconds. A 7° increase in angular span over 1st generation decreased the ASF by 37%. Additionally, the 2nd generation s‐DBT system utilizing a specific AFVR reconstruction method resulted in a 92% increase in the in plane resolution over CM DBT system, and a 37% increase in spatial resolution over the 1st generation s-‐DBT system.

The coronary artery calcium score (CACS) measures the buildup of calcium on the coronary artery wall and has been shown to be an important predictor of the risk of coronary artery diseases (CAD). Currently CACS is measured using CT, though the relatively high cost and high radiation dose has limited its adoption as a routine screening procedure. Digital Chest Tomosynthesis (DCT), a low dose and low cost alternative to CT, and has been shown to achieve 90% of sensitivity of CT in lung disease screening. However commercial DCT requires long scanning time and cannot be adapted for high resolution gated cardiac imaging, necessary for CACS. The stationary DCT system (s- DCT), developed in our lab, has the potential to significantly shorten the scanning time and enables high resolution cardiac gated imaging. Here we report the preliminary results of using s-DCT to estimate the CACS. A phantom heart model was developed and scanned by the s-DCT system and a clinical CT in a phantom model with realistic coronary calcifications. The adapted fan-beam volume reconstruction (AFVR) method, developed specifically for stationary tomosynthesis systems, is used to obtain high resolution tomosynthesis images. A trained cardiologist segmented out the calcifications and the CACS was obtained. We observed a strong correlation between the tomosynthesis derived CACS and CT CACS (r² = 0.88). Our results shows s-DCT imaging has the potential to estimate CACS, thus providing a possible low cost and low dose imaging protocol for screening and monitoring CAD.

The aim of this study was to compare the detection of microcalcification clusters by human observers in breast images using 2D-mammography and narrow (15°/15 projections) and wide (50°/25 projections) angle digital breast tomosynthesis (DBT). Simulated microcalcification clusters with a range of microcalcification diameters (125 μm-275 μm) were inserted into 6 cm thick simulated compressed breasts. Breast images were produced with and without inserted microcalcification clusters using a set of image modelling tools, which were developed to represent clinical imaging by mammography and tomosynthesis. Commercially available software was used for image processing and image reconstruction. The images were then used in a series of 4-alternative forced choice (4AFC) human observer experiments conducted for signal detection with the microcalcification clusters as targets. The minimum detectable calcification diameter was found for each imaging modality: (i) 2D-mammography: 164±5 μm (ii) narrow angle DBT: 210±5 μm, (iii) wide angle DBT: 255±4 μm. A statistically significant difference was found between the minimum detectable calcification diameters that can be detected by the three imaging modalities. Furthermore, it was found that there was not a statistically significant difference between the results of the five observers that participated in this study. In conclusion, this study presents a method that quantifies the threshold diameter required for microcalcification detection, using high resolution, realistic images with observers, for the comparison of DBT geometries with 2D-mammography. 2Dmammography can visualise smaller detail diameter than both DBT imaging modalities and narrow-angle DBT can visualise a smaller detail diameter than wide-angle DBT.

Dual-energy contrast-enhanced digital mammography (DE CE-DM) uses an iodinated contrast agent to image the perfusion and vasculature of the breast. DE images are obtained by a weighted logarithmic subtraction of the high-energy (HE) and low-energy (LE) image pairs. We hypothesized that the optimal DE subtraction weighting factor is thickness-dependent, and developed a method for determining breast tissue composition and thickness in DE CE-DM. Phantoms were constructed using uniform blocks of 100% glandular-equivalent and 100% adipose-equivalent material. The thickness of the phantoms ranged from 3 to 8 cm, in 1 cm increments. For a given thickness, the glandular-adipose composition of the phantom was varied using different combinations of blocks. The logarithmic LE and logarithmic HE signal intensities were measured; they decrease linearly with increasing glandularity for a given thickness. The signals decrease with increasing phantom thickness and the x-ray signal decreases linearly with thickness for a given glandularity. As the thickness increases, the attenuation difference per additional glandular block decreases, indicating beam hardening. From the calibration mapping, we have demonstrated that we can predict percent glandular tissue and thickness when given two distinct signal intensities. Our results facilitate the subtraction of tissue at the boundaries of the breast, and aid in discriminating between contrast agent uptake in glandular tissue and subtraction artifacts.

Virtual clinical trials are an emergent approach for the rapid evaluation and comparison of various breast imaging technologies and techniques using computer-based modeling tools. A fundamental requirement of this approach for mammography is the use of realistic looking breast anatomy in the studies to produce clinically relevant results. In this work, a biologically inspired approach has been used to simulate realistic synthetic breast phantom blocks for use in virtual clinical trials. A variety of high and low frequency features (including Cooper’s ligaments, blood vessels and glandular tissue) have been extracted from clinical digital breast tomosynthesis images and used to simulate synthetic breast blocks. The appearance of the phantom blocks was validated by presenting a selection of simulated 2D and DBT images interleaved with real images to a team of experienced readers for rating using an ROC paradigm. The average areas under the curve for 2D and DBT images were 0.53±.04 and 0.55±.07 respectively; errors are the standard errors of the mean. The values indicate that the observers had difficulty in differentiating the real images from simulated images. The statistical properties of simulated images of the phantom blocks were evaluated by means of power spectrum analysis. The power spectrum curves for real and simulated images closely match and overlap indicating good agreement.

An anthropomorphic digital breast phantom has been developed with the goal of generating random voxelized breast models that capture the anatomic variability observed in vivo. This is a new phantom and is not based on existing digital breast phantoms or segmentation of patient images. It has been designed at the outset to be modality agnostic (i.e., suitable for use in modeling x-ray based imaging systems, magnetic resonance imaging, and potentially other imaging systems) and open source so that users may freely modify the phantom to suit a particular study. In this work we describe the modeling techniques that have been developed, the capabilities and novel features of this phantom, and study simulated images produced from it. Starting from a base quadric, a series of deformations are performed to create a breast with a particular volume and shape. Initial glandular compartments are generated using a Voronoi technique and a ductal tree structure with terminal duct lobular units is grown from the nipple into each compartment. An additional step involving the creation of fat and glandular lobules using a Perlin noise function is performed to create more realistic glandular/fat tissue interfaces and generate a Cooper’s ligament network. A vascular tree is grown from the chest muscle into the breast tissue. Breast compression is performed using a neo-Hookean elasticity model. We show simulated mammographic and T₁-weighted MRI images and study properties of these images.

Spectral imaging is the acquisition of multiple images of an object at different energy spectra. In mammography, dual-energy imaging (spectral imaging with two energy levels) has been investigated for several applications, in particular material decomposition, which allows for quantitative analysis of breast composition and quantitative contrast-enhanced imaging. Material decomposition with dual-energy imaging is based on the assumption that there are two dominant photon interaction effects that determine linear attenuation: the photoelectric effect and Compton scattering. This assumption limits the number of basis materials, i.e. the number of materials that are possible to differentiate between, to two. However, Rayleigh scattering may account for more than 10% of the linear attenuation in the mammography energy range. In this work, we show that a modified version of a scanning multi-slit spectral photon-counting mammography system is able to acquire three images at different spectra and can be used for triple-energy imaging. We further show that triple-energy imaging in combination with the efficient scatter rejection of the system enables measurement of Rayleigh scattering, which adds an additional energy dependency to the linear attenuation and enables material decomposition with three basis materials. Three available basis materials have the potential to improve virtually all applications of spectral imaging.

Mammography is currently the standard imaging modality used to screen women for breast abnormalities and, as a result, it is a tool of great importance for the early detection of breast cancer. Physical phantoms are commonly used as surrogates of breast tissue to evaluate some aspects of the performance of mammography systems. However, most phantoms do not reproduce the anatomic heterogeneity of real breasts. New fabrication technologies, such as 3D printing, have created the opportunity to build more complex, anatomically realistic breast phantoms that could potentially assist in the evaluation of mammography systems. The primary objective of this work is to present a simple, easily reproducible methodology to design and print 3D objects that replicate the attenuation profile observed in real 2D mammograms. The secondary objective is to evaluate the capabilities and limitations of the competing 3D printing technologies, and characterize the x-ray properties of the different materials they use. Printable phantoms can be created using the open-source code introduced in this work, which processes a raw mammography image to estimate the amount of x-ray attenuation at each pixel, and outputs a triangle mesh object that encodes the observed attenuation map. The conversion from the observed pixel gray value to a column of printed material with equivalent attenuation requires certain assumptions and knowledge of multiple imaging system parameters, such as x-ray energy spectrum, source-to-object distance, compressed breast thickness, and average breast material attenuation. A detailed description of the new software, a characterization of the printed materials using x-ray spectroscopy, and an evaluation of the realism of the sample printed phantoms are presented.

Breast ultrasound tomography is an emerging imaging modality to reconstruct the sound speed, density, and ultrasound attenuation of the breast in addition to ultrasound reflection/beamforming images for breast cancer detection and characterization. We recently designed and manufactured a new synthetic-aperture breast ultrasound tomography prototype with two parallel transducer arrays consisting of a total of 768 transducer elements. The transducer arrays are translated vertically to scan the breast in a warm water tank from the chest wall/axillary region to the nipple region to acquire ultrasound transmission and reflection data for whole-breast ultrasound tomography imaging. The distance of these two ultrasound transducer arrays is adjustable for scanning breasts with different sizes. We use our breast ultrasound tomography prototype to acquire phantom and in vivo patient ultrasound data to study its feasibility for breast imaging. We apply our recently developed ultrasound imaging and tomography algorithms to ultrasound data acquired using our breast ultrasound tomography system. Our in vivo patient imaging results demonstrate that our breast ultrasound tomography can detect breast lesions shown on clinical ultrasound and mammographic images.

In this article, we present an easy-to-implement Multi-energy CT scanning strategy and a corresponding reconstruction method, which facilitate spectral CT imaging by improving the data efficiency the number-of-energy- channel fold without introducing visible limited-angle artifacts caused by reducing projection views. Leveraging the structure coherence at different energies, we first pre-reconstruct a prior structure information image using projection data from all energy channels. Then, we perform a k-means clustering on the prior image to generate a sparse dictionary representation for the image, which severs as a structure information constraint. We com- bine this constraint with conventional compressed sensing method and proposed a new model which we referred as Joint Clustering Prior and Sparsity Regularization (CPSR). CPSR is a convex problem and we solve it by Alternating Direction Method of Multipliers (ADMM). We verify our CPSR reconstruction method with a numerical simulation experiment. A dental phantom with complicate structures of teeth and soft tissues is used. X-ray beams from three spectra of different peak energies (120kVp, 90kVp, 60kVp) irradiate the phantom to form tri-energy projections. Projection data covering only 75◦ from each energy spectrum are collected for reconstruction. Independent reconstruction for each energy will cause severe limited-angle artifacts even with the help of compressed sensing approaches. Our CPSR provides us with images free of the limited-angle artifact. All edge details are well preserved in our experimental study.

Compared to conventional computed tomography (CT), dual energy CT allows for improved material decomposition by conducting measurements at two distinct energy spectra. Since radiation exposure is a major concern in clinical CT, there is a need for tools to reduce the noise level in images while preserving diagnostic information. One way to achieve this goal is the application of image-based denoising algorithms after an analytical reconstruction has been performed. We have developed a modified dictionary denoising algorithm for dual energy CT aimed at exploiting the high spatial correlation between between images obtained from different energy spectra. Both the low-and high energy image are partitioned into small patches which are subsequently normalized. Combined patches with improved signal-to-noise ratio are formed by a weighted addition of corresponding normalized patches from both images. Assuming that corresponding low-and high energy image patches are related by a linear transformation, the signal in both patches is added coherently while noise is neglected. Conventional dictionary denoising is then performed on the combined patches. Compared to conventional dictionary denoising and bilateral filtering, our algorithm achieved superior performance in terms of qualitative and quantitative image quality measures. We demonstrate, in simulation studies, that this approach can produce 2d-histograms of the high- and low-energy reconstruction which are characterized by significantly improved material features and separation. Moreover, in comparison to other approaches that attempt denoising without simultaneously using both energy signals, superior similarity to the ground truth can be found with our proposed algorithm.

In this work, we provide an initial characterization of a novel twin robotic X-ray system. This system is equipped with two motor-driven telescopic arms carrying X-ray tube and flat-panel detector, respectively. 2D radiographs and fluoroscopic image sequences can be obtained from different viewing angles. Projection data for 3D cone-beam CT reconstruction can be acquired during simultaneous movement of the arms along dedicated scanning trajectories. We provide an initial evaluation of the 3D image quality based on phantom scans and clinical images. Furthermore, initial evaluation of patient dose is conducted. The results show that the system delivers high image quality for a range of medical applications. In particular, high spatial resolution enables adequate visualization of bone structures. This system allows 3D X-ray scanning of patients in standing and weight-bearing position. It could enable new 2D/3D imaging workflows in musculoskeletal imaging and improve diagnosis of musculoskeletal disorders.

We propose an adapted method of our previously published five-dimensional (5D) motion compensation (MoCo) algorithm1, developed for micro-CT imaging of small animals, to provide for the first time motion artifact-free 5D cone-beam CT (CBCT) images from a conventional flat detector-based CBCT scan of clinical patients. Image quality of retrospectively respiratory- and cardiac-gated volumes from flat detector CBCT scans is deteriorated by severe sparse projection artifacts. These artifacts further complicate motion estimation, as it is required for MoCo image reconstruction. For high quality 5D CBCT images at the same x-ray dose and the same number of projections as todays 3D CBCT we developed a double MoCo approach based on motion vector fields (MVFs) for respiratory and cardiac motion. In a first step our already published four-dimensional (4D) artifact-specific cyclic motion-compensation (acMoCo) approach is applied to compensate for the respiratory patient motion. With this information a cyclic phase-gated deformable heart registration algorithm is applied to the respiratory motion-compensated 4D CBCT data, thus resulting in cardiac MVFs. We apply these MVFs on double-gated images and thereby respiratory and cardiac motion-compensated 5D CBCT images are obtained. Our 5D MoCo approach processing patient data acquired with the TrueBeam 4D CBCT system (Varian Medical Systems). Our double MoCo approach turned out to be very efficient and removed nearly all streak artifacts due to making use of 100% of the projection data for each reconstructed frame. The 5D MoCo patient data show fine details and no motion blurring, even in regions close to the heart where motion is fastest.

Mobile and compact C-arm systems are routinely used in interventional procedures for fluoroscopic CT imaging. The mechanical requirements guarantee for a maximum of flexibility and mobility but restrict the mechanical rotation range (e.g. 165°) and the lateral size of the field of measurement (FOM), typically about 160 mm. Recently, the rotate-plus-shift trajectory for the acquisition of complete datasets from 180° minus fan-angle has been published.^{1, 2} Here, we combine the rotate-plus-shift trajectory with a shifted detector approach for a fully motorized C-arm system. As the isocenter in non-centric C-arms can be freely chosen, the shifted detector can be equally well absorbed with an offset of the C parallel to the transaxial detector direction. The typical rotation range of 360° used in shifted detector trajectories is replaced by a double rotate-plus-shift scan requiring a rotation range of at least 180° minus fan-angle. The trajectory increasing the diameter of the FOM by up to a factor of two is presented and the practical application of variations with an asymmetric FOM is shown. For image reconstruction we use our modified FDK algorithm that is equipped with a generalized redundancy weight. The presented trajectory can increase the applicability and flexibility of C-arm systems and has the potential to perform intra-operative large volume control or overview scans and thus reduce the patient’s risk.

Striped ratio grids are a new concept for scatter management in cone-beam CT. These grids are a modification of conventional anti-scatter grids and consist of stripes which alternate between high grid ratio and low grid ratio. Such a grid is related to existing hardware concepts for scatter estimation such as blocker-based methods or primary modulation, but rather than modulating the primary, the striped ratio grid modulates the scatter. The transitions between adjacent stripes can be used to estimate and subtract the remaining scatter. However, these transitions could be contaminated by variation in the primary radiation. We describe a simple nonlinear image processing algorithm to estimate scatter, and proceed to validate the striped ratio grid on experimental data of a pelvic phantom. The striped ratio grid is emulated by combining data from two scans with different grids. Preliminary results are encouraging and show a significant reduction of scatter artifact.

Purpose: Cone-beam CT (CBCT) of the extremities provides high spatial resolution, but its quantitative accuracy may be challenged by involuntary sub-mm patient motion that cannot be eliminated with simple means of external immobilization. We investigate a two-step iterative motion compensation based on a multi-component metric of image sharpness. Methods: Motion is considered with respect to locally rigid motion within a particular region of interest, and the method supports application to multiple locally rigid regions. Motion is estimated by maximizing a cost function with three components: a gradient metric encouraging image sharpness, an entropy term that favors high contrast and penalizes streaks, and a penalty term encouraging smooth motion. Motion compensation involved initial coarse estimation of gross motion followed by estimation of fine-scale displacements using high resolution reconstructions. The method was evaluated in simulations with synthetic motion (1–4 mm) applied to a wrist volume obtained on a CMOS-based CBCT testbench. Structural similarity index (SSIM) quantified the agreement between motion-compensated and static data. The algorithm was also tested on a motion contaminated patient scan from dedicated extremities CBCT. Results: Excellent correction was achieved for the investigated range of displacements, indicated by good visual agreement with the static data. 10-15% improvement in SSIM was attained for 2-4 mm motions. The compensation was robust against increasing motion (4% decrease in SSIM across the investigated range, compared to 14% with no compensation). Consistent performance was achieved across a range of noise levels. Significant mitigation of artifacts was shown in patient data. Conclusion: The results indicate feasibility of image-based motion correction in extremities CBCT without the need for a priori motion models, external trackers, or fiducials.

C-arm-based cone-beam CT (CBCT) systems with flat-panel detectors are suitable for diagnostic knee imaging due to their potentially flexible selection of CT trajectories and wide volumetric beam coverage. In knee CT imaging, over-exposure artifacts can occur because of limitations in the dynamic range of the flat panel detectors present on most CBCT systems. We developed a straightforward but effective method for correction and detection of over-exposure for an Automatic Exposure Control (AEC)-enabled standard knee scan incorporating a prior low dose scan. The radiation dose associated with the low dose scan was negligible (0.0042mSv, 2.8% increase) which was enabled by partially sampling the projection images considering the geometry of the knees and lowering the dose further to be able to just see the skin-air interface. We combined the line integrals from the AEC and low dose scans after detecting over-exposed regions by comparing the line profiles of the two scans detector row-wise. The combined line integrals were reconstructed into a volumetric image using filtered back projection. We evaluated our method using in vivo human subject knee data. The proposed method effectively corrected and detected over-exposure, and thus recovered the visibility of exterior tissues (e.g., the shape and density of the patella, and the patellar tendon), incorporating a prior low dose scan with a negligible increase in radiation exposure.

For grating interferometer-based x-ray differential phase contrast (DPC) imaging systems, their noise performance is strongly dependent on both the visibility of the interference fringe pattern and the total number of photons used to acquire and extract the DPC signal. For a given interferometer, it is usually designed to work at a specific x-ray energy, therefore any deviation from the designed energy may result in certain visibility loss. In this work, a single photon counting detector (PCD) was incorporated into a DPC imaging system, which enabled photons with energies close to the designed operation energy of the interferometer to be selectively used for DPC signal extraction. This approach led to significant boost in the fringe visibility, but it also discarded x-ray photons with other energies incident on the detector and might result in degradations of the overall radiation dose efficiency of the DPC imaging systems. This work presents a novel singular value decomposition (SVD)-based method to leverage the entire spectrum of x-ray photons detected by the PCD, enabling both fringe visibility improvement and reduction in image noise. As evidenced by the results of experimental phantom studies, the contrast-to-noise ratio of the final DPC images could be effectively improved by the proposed method.

One of the most exciting aspects of the grating based x-ray differential phase contrast (DPC) acquisition method is the concurrent generation of the so-called dark field (DF) signal, along with the classical absorption signal and the novel DPC signal. The DF signal is associated with local distribution of small angle scatterers in an image object, while the absorption signal and DPC signal are often used to characterize the relatively uniform structure of the image object. Besides the endogenous image contrast, exogenous contrast media are often used in x-ray imaging to locally enhance the image signal. This paper proposes a potential contrast medium for DF signal enhancement: microbubbles (MBs). MBs have already been developed for clinical use in ultrasound imaging, and recent experimental studies have shown that MBs may also enhance the DF signal, although it remained unclear how the physical characteristics of the MBs quantitatively impact the DF signal. In this paper, a systematic study was performed to investigate the quantitative relationships between the DF signal and the following properties of MBs: size, concentration, shell thickness, size uniformity, and whether gold nanoparticles were attached. The experimental results demonstrated that, an increased MB size (about 4 microns) may generate a stronger DF signal for our DPC imaging system; additionally, a moderately increased shell thickness and the use of gold nanoparticles on the shell surface also resulted in further enhancement of the DF signal. These findings may provide critical information needed for using MBs as the contrast agent of x-ray DF imaging.

X-ray grating-based phase-contrast Talbot-Lau interferometry is a promising imaging technology that has the potential to raise soft tissue contrast in comparison to conventional attenuation-based imaging. Additionally, it is sensitive to attenuation, refraction and scattering of the radiation and thus provides complementary and otherwise inaccessible information due to the dark-field image, which shows the sub-pixel size granularity of the measured object. Until recent progress the method has been mainly limited to photon energies below 40 keV. Scaling the method to photon energies that are sufficient to pass large and spacious objects represents a challenging task. This is caused by increasing demands regarding the fabrication process of the gratings and the broad spectra that come along with the use of polychromatic X-ray sources operated at high acceleration voltages. We designed a setup that is capable to reach high visibilities in the range from 50 to 120 kV. Therefore, spacious and dense parts of the human body with high attenuation can be measured, such as a human knee. The authors will show investigations on the resulting attenuation, differential phase-contrast and dark-field images. The images experimentally show that X-ray grating-based phase-contrast radiography is feasible with highly absorbing parts of the human body containing massive bones.

X-ray phase contrast imaging has been investigated during the last two decades for potential benefits in soft tissue imaging. Long imaging time, high radiation dose and general measurement complexity involving motion of x-ray optical components have prevented the clinical translation of these methods. In all existing popular phase contrast imaging methods, multiple measurements per projection angle involving motion of optical components are required to achieve quantitatively accurate estimation of absorption, phase and differential phase. Recently we proposed an algorithmic approach to use spectral detection data in a phase contrast imaging setup to obtain absorption, phase and differential phase in a single-step. Our generic approach has been shown via simulations in all three types of phase contrast imaging: propagation, coded aperture and grating interferometry. While other groups have used spectral detector in phase contrast imaging setups, our proposed method is unique in outlining an approach to use this spectral data to simplify phase contrast imaging. In this abstract we show the first experimental proof of our single-shot phase retrieval using a Medipix3 photon counting detector in an edge illumination aperture (also referred to as coded aperture) phase contrast set up as well as for a free space propagation setup. Our preliminary results validate our new transport equation for edge illumination PCI and our spectral phase retrieval algorithm for both PCI methods being investigated. Comparison with simulations also point to excellent performance of Medipix3 built-in charge sharing correction mechanism.

Flat-panel cone-beam CT (FP-CBCT) is a promising imaging modality, partly due to its potential for high spatial resolution reconstructions in relatively compact scanners. Despite this potential, FP-CBCT can face difficulty resolving important fine scale structures (e.g, trabecular details in dedicated extremities scanners and microcalcifications in dedicated CBCT mammography). Model-based methods offer one opportunity to improve high-resolution performance without any hardware changes. Previous work, based on a linearized forward model, demonstrated improved performance when both system blur and spatial correlations characteristics of FP-CBCT systems are modeled. Unfortunately, the linearized model relies on a staged processing approach that complicates tuning parameter selection and can limit the finest achievable spatial resolution. In this work, we present an alternative scheme that leverages a full nonlinear forward model with both system blur and spatially correlated noise. A likelihood-based objective function is derived from this forward model and we derive an iterative optimization algorithm for its solution. The proposed approach is evaluated in simulation studies using a digital extremities phantom and resolution-noise trade-offs are quantitatively evaluated. The correlated nonlinear model outperformed both the uncorrelated nonlinear model and the staged linearized technique with up to a 86% reduction in variance at matched spatial resolution. Additionally, the nonlinear models could achieve finer spatial resolution (correlated: 0.10 mm, uncorrelated: 0.11 mm) than the linear correlated model (0.15 mm), and traditional FDK (0.40 mm). This suggests the proposed nonlinear approach may be an important tool in improving performance for high-resolution clinical applications.

We present a new algorithm that allows for raw data-based automated cardiac and respiratory intrinsic gating in cone-beam CT scans. It can be summarized in three steps: First, a median filter is applied to an initially reconstructed volume. The forward projection of this volume contains less motion information and is subtracted from the original projections. This results in new raw data that contain only moving and not static anatomy like bones, that would otherwise impede the cardiac or respiratory signal acquisition. All further steps are applied to these modified raw data. Second, the raw data are cropped to a region of interest (ROI). The ROI in the raw data is determined by the forward projection of a binary volume of interest (VOI) that includes the diaphragm for respiratory gating and most of the edge of the heart for cardiac gating. Third, the mean gray value in this ROI is calculated for every projection and the respiratory/cardiac signal is acquired using a bandpass filter. Steps two and three are carried out simultaneously for 64 or 1440 overlapping VOI inside the body for the respiratory or cardiac signal respectively. The signals acquired from each ROI are compared and the most consistent one is chosen as the desired cardiac or respiratory motion signal. Consistency is assessed by the standard deviation of the time between two maxima. The robustness and efficiency of the method is evaluated using simulated and measured patient data by computing the standard deviation of the mean signal difference between the ground truth and the intrinsic signal.

Purpose: Prompt and reliable detection of intracranial hemorrhage (ICH) has substantial clinical impact in diagnosis and treatment of stroke and traumatic brain injury. This paper describes the design, development, and preliminary performance characterization of a dedicated cone-beam CT (CBCT) head scanner prototype for imaging of acute ICH. Methods: A task-based image quality model was used to analyze the detectability index as a function of system configuration, and hardware design was guided by the results of this model-based optimization. A robust artifact correction pipeline was developed using GPU-accelerated Monte Carlo (MC) scatter simulation, beam hardening corrections, detector veiling glare, and lag deconvolution. An iterative penalized weighted least-squares (PWLS) reconstruction framework with weights adjusted for artifact-corrected projections was developed. Various bowtie filters were investigated for potential dose and image quality benefits, with a MC-based tool providing estimates of spatial dose distribution. Results: The initial prototype will feature a source-detector distance of 1000 mm and source-axis distance of 550 mm, a 43x43 cm² flat panel detector, and a 15° rotating anode x-ray source with 15 kW power and 0.6 focal spot size. Artifact correction reduced image nonuniformity by ~250 HU, and PWLS reconstruction with modified weights improved the contrast to noise ratio by 20%. Inclusion of a bowtie filter can potentially reduce dose by 50% and improve CNR by 25%. Conclusions: A dedicated CBCT system capable of imaging millimeter-scale acute ICH was designed. Preliminary findings support feasibility of point-of-care applications in TBI and stroke imaging, with clinical studies beginning on a prototype.

In this work, a newly developed reconstruction algorithm, Synchronized MultiArtifact Reduction with Tomographic RECONstruction (SMART-RECON), was applied to C-arm cone beam CT perfusion (CBCTP) imaging. This algorithm contains a special rank regularizer, designed to reduce limited-view artifacts associated with super- short scan reconstructions. As a result, high temporal sampling and temporal resolution image reconstructions were achieved using an interventional C-arm x-ray system. The algorithm was evaluated in terms of the fidelity of the dynamic contrast update curves and the accuracy of perfusion parameters through numerical simulation studies. Results shows that, not only were the dynamic curves accurately recovered (relative root mean square error ∈ [3%, 5%] compared with [13%, 22%] for FBP), but also the noise in the final perfusion maps was dramatically reduced. Compared with filtered backprojection, SMART-RECON generated CBCTP maps with much improved capability in differentiating lesions with perfusion deficits from the surrounding healthy brain tissues.

Currently, clinical acquisition of IV 3D-DSA requires two separate scans: one mask scan without contrast medium and a filled scan with contrast injection. Having two separate scans adds radiation dose to the patient and increases the likelihood of suffering inadvertent patient motion induced mis-registration and the associated mis-registraion artifacts in IV 3D-DSA images. In this paper, a new technique, SMART-RECON is introduced to generate IV 3D-DSA images from a single Cone Beam CT (CBCT) acquisition to eliminate the mask scan. Potential benefits of eliminating mask scan would be: (1) both radiation dose and scan time can be reduced by a factor of 2; (2) intra-sweep motion can be eliminated; (3) inter-sweep motion can be mitigated. Numerical simulations were used to validate the algorithm in terms of contrast recoverability and the ability to mitigate limited view artifacts.

When an automatic exposure control is introduced in C-arm cone beam CT data acquisition, the spectral inconsistencies between acquired projection data are exacerbated. As a result, conventional water/bone correction schemes are not as effective as in conventional diagnostic x-ray CT acquisitions with a fixed tube potential. In this paper, a new method was proposed to reconstruct several images with different degrees of spectral consistency and thus different levels of beam hardening artifacts. The new method relies neither on prior knowledge of the x-ray beam spectrum nor on prior compositional information of the imaging object. Numerical simulations were used to validate the algorithm.

We introduce a novel strategy for fluence field modulation (FFM) in x-ray CT using multiple aperture devices (MADs). MAD filters permit FFM by blocking or transmitting the x-ray beam on a fine (0.1-1 mm) scale. The filters have a number of potential advantages over other beam modulation strategies including the potential for a highly compact design, modest actuation speed and acceleration requirements, and spectrally neutral filtration due to their essentially binary action. In this work, we present the underlying MAD filtration concept including a design process to achieve a specific class of FFM patterns. A set of MAD filters is fabricated using a tungsten laser sintering process and integrated into an x-ray CT test bench. A characterization of the MAD filters is conducted and compared to traditional attenuating bowtie filters and the ability to flatten the fluence profile for a 32 cm acrylic phantom is demonstrated. MAD-filtered tomographic data was acquired on the CT test bench and reconstructed without artifacts associated with the MAD filter. These initial studies suggest that MAD-based FFM is appropriate for integration in clinical CT system to create patient-specific fluence field profile and reduce radiation exposures.

Physicians rely on CT Perfusion (CTP) images and quantitative image data, including cerebral blood flow, cerebral blood volume, and bolus arrival delay, to diagnose and treat stroke patients. However, the quantification of these metrics may vary depending on the computational method used. Therefore, we have developed a dynamic and realistic digital brain phantom upon which CTP scans can be simulated based on a set of ground truth scenarios. Building upon the previously developed 4D extended cardiac-torso (XCAT) phantom containing a highly detailed brain model, this work consisted of expanding the intricate vasculature by semi-automatically segmenting existing MRA data and fitting nonuniform rational B-spline surfaces to the new vessels. Using time attenuation curves input by the user as reference, the contrast enhancement in the vessels changes dynamically. At each time point, the iodine concentration in the arteries and veins is calculated from the curves and the material composition of the blood changes to reflect the expected values. CatSim, a CT system simulator, generates simulated data sets of this dynamic digital phantom which can be further analyzed to validate CTP studies and post-processing methods. The development of this dynamic and realistic digital phantom provides a valuable resource with which current uncertainties and controversies surrounding the quantitative computations generated from CTP data can be examined and resolved.

Quantitative computed tomography (CT) applications based on statistical iterative reconstruction algorithms require accurate models of the CT acquisition process, with a key component being the x-ray fan beam intensity. We present a method to experimentally determine the extra-focal radiation profile incident on individual CT detectors. Using a tungsten cylinder as a knife edge, a super-sampled signal was created from sinogram data, which traced the “occlusion” of the x-ray source as seen by a detector. By differentiating this signal and correcting for finite detector size and motion blur, the effective source profile can be recovered. Extra-focal scatter was found to be on the order of 1-3 percent of the focal beam intensity, with lower relative magnitude at the isocenter and increasing towards the edge of the fan beam, with its profile becoming asymmetric at large angles. The implications for reconstruction algorithms and QCT applications will be discussed.

Cerebral CT perfusion (CTP) imaging is playing an important role in the diagnosis and treatment of acute ischemic strokes. Meanwhile, the reliability of CTP-based ischemic lesion detection has been challenged due to the noisy appearance and low signal-to-noise ratio of CTP maps. To reduce noise and improve image quality, a rigorous study on the noise transfer properties of CTP systems is highly desirable to provide the needed scientific guidance. This paper concerns how noise in the CTP source images propagates to the final CTP maps. Both theoretical deviations and subsequent validation experiments demonstrated that, the noise level of background frames plays a dominant role in the noise of the cerebral blood volume (CBV) maps. This is in direct contradiction with the general belief that noise of non-background image frames is of greater importance in CTP imaging. The study found that when radiation doses delivered to the background frames and to all non-background frames are equal, lowest noise variance is achieved in the final CBV maps. This novel equality condition provides a practical means to optimize radiation dose delivery in CTP data acquisition: radiation exposures should be modulated between background frames and non-background frames so that the above equality condition is satisïnAed. For several typical CTP acquisition protocols, numerical simulations and in vivo canine experiment demonstrated that noise of CBV can be effectively reduced using the proposed exposure modulation method.

The overall goal of this work is to develop a rapid, accurate and fully automated software tool to estimate patient-specific organ doses from computed tomography (CT) scans using a deterministic Boltzmann Transport Equation solver and automated CT segmentation algorithms. This work quantified the accuracy of organ dose estimates obtained by an automated segmentation algorithm. The investigated algorithm uses a combination of feature-based and atlas-based methods. A multiatlas approach was also investigated. We hypothesize that the auto-segmentation algorithm is sufficiently accurate to provide organ dose estimates since random errors at the organ boundaries will average out when computing the total organ dose. To test this hypothesis, twenty head-neck CT scans were expertly segmented into nine regions. A leave-one-out validation study was performed, where every case was automatically segmented with each of the remaining cases used as the expert atlas, resulting in nineteen automated segmentations for each of the twenty datasets. The segmented regions were applied to gold-standard Monte Carlo dose maps to estimate mean and peak organ doses. The results demonstrated that the fully automated segmentation algorithm estimated the mean organ dose to within 10% of the expert segmentation for regions other than the spinal canal, with median error for each organ region below 2%. In the spinal canal region, the median error was 7% across all data sets and atlases, with a maximum error of 20%. The error in peak organ dose was below 10% for all regions, with a median error below 4% for all organ regions. The multiple-case atlas reduced the variation in the dose estimates and additional improvements may be possible with more robust multi-atlas approaches. Overall, the results support potential feasibility of an automated segmentation algorithm to provide accurate organ dose estimates.

Energy-resolved photon-counting CT (PCCT) is promising for material decomposition with multi-contrast agents. However, corrections for non-idealities of PCCT detectors are required, which are still active research areas. In addition, PCCT is associated with very high cost due to lack of mass production. In this work, we proposed an alternative approach to performing multi-energy CT, which was achieved by acquiring triple or quadruple x-ray beam measurements on a dual-source CT scanner. This strategy was based on a “Twin Beam” design on a single-source scanner for dual-energy CT. Examples of beam filters and spectra for triple and quadruple x-ray beam were provided. Computer simulation studies were performed to evaluate the accuracy of material decomposition for multi-contrast mixtures using both tri-beam and quadruple-beam configurations. The proposed strategy can be readily implemented on a dual-source scanner, which may allow material decomposition of multi-contrast agents to be performed on clinical CT scanners with energy-integrating detector.

Respiration-correlated CT (4DCT) forms the basis of clinical 4D radiotherapy workflows for patients with thoracic and abdominal lesions. 4DCT image data, however, often suffers from motion artifacts due to unfulfilled assumptions during reconstruction and image/projection data sorting. In this work and focusing on low-pitch helical scanning protocols, two questionable assumptions are addressed: (1) the need for regular breathing patterns and (2) a constant correlation between the external breathing signal acquired for image/projection sorting and internal motion patterns. To counteract (1), a patient-specific upper breathing signal amplitude threshold is introduced to avoid artifacts due to unusual deep inspiration (helpful for both amplitude- and phase-based reconstruction). In addition, a projection data binning algorithm based on a statistical analysis of the patient's breathing signal is proposed to stabilize phase-based sorting. To further alleviate the need for (2), an image artifact metric is incorporated into and minimized during the reconstruction process. The optimized reconstruction is evaluated using 30 clinical 4DCT data sets and demonstrated to significantly reduce motion artifacts.

Detective Quantum Efficiency (DQE) is one of the most important image quality metrics for evaluating the spatial resolution performance of flat-panel x-ray detectors. In this work, we simulate the DQE of amorphous selenium (a-Se) xray detectors with a detailed Monte Carlo transport code (ARTEMIS) for modeling semiconductor-based direct x-ray detectors. The transport of electron-hole pairs is achieved with a spatiotemporal model that accounts for recombination and trapping of carriers and Coulombic effects of space charge and external applied electric field. A range of x-ray energies has been simulated from 10 to 100 keV. The DQE results can be used to study the spatial resolution characteristics of detectors at different energies.

Indirect conversion flat panel detectors (FPDs) based on amorphous silicon (a-Si) technology are widely used in digital X-ray imaging. In such FPDs a scintillator layer is used for converting X-rays into visible light photons. However, the lateral spread of these photons inside the scintillator layer reduces spatial resolution of the FPD. In this study, FPDs incorporating pixelated scintillators with a barrier rib structure were developed to limit lateral spread of light photons thereby improving spatial resolution. For the pixelated scintillator, a two-dimensional barrier rib structure was first manufactured on a substrate layer, coated with reflective materials, and filled to the rim with the scintillating material of gadolinium oxysulfide (GOS). Several scintillator samples were fabricated, with pitch size varying from 160 to 280 μm and rib height from 200 to 280 μm. The samples were directly coupled to an a-Si flat panel photodiode array with a pitch of 200 μm to convert optical photons to electronic signals. With the pixelated scintillator, the detector modulation transfer function was shown to improve significantly (by 94% at 2 cycle/mm) compared to a detector using an unstructured GOS layer. However, the prototype does show lower sensitivity due to the decrease in scintillator fill factor. The preliminary results demonstrated the feasibility of using the barrier-rib structure to improve the spatial resolution of FPDs. Such an improvement would greatly benefit nondestructive testing applications where the spatial resolution is the most important parameter. Further investigation will focus on improving the detector sensitivity and exploring its medical applications.

Digital Breast Tomosynthesis (DBT) requires excellent image quality in a dynamic mode at very low dose levels while Full Field Digital Mammography (FFDM) is a static imaging modality that requires high saturation dose levels. These opposing requirements can only be met by a dynamic detector with a high dynamic range. This paper will discuss a wafer-scale CMOS-based mammography detector with 49.5 μm pixels and a CsI scintillator. Excellent image quality is obtained for FFDM as well as DBT applications, comparing favorably with a-Se detectors that dominate the X-ray mammography market today. The typical dynamic range of a mammography detector is not high enough to accommodate both the low noise and the high saturation dose requirements for DBT and FFDM applications, respectively. An approach based on gain switching does not provide the signal-to-noise benefits in the low-dose DBT conditions. The solution to this is to add frame summing functionality to the detector. In one X-ray pulse several image frames will be acquired and summed. The requirements to implement this into a detector are low noise levels, high frame rates and low lag performance, all of which are unique characteristics of CMOS detectors. Results are presented to prove that excellent image quality is achieved, using a single detector for both DBT as well as FFDM dose conditions. This method of frame summing gave the opportunity to optimize the detector noise and saturation level for DBT applications, to achieve high DQE level at low dose, without compromising the FFDM performance.

Active matrix flat panel imagers (AMFPI) have become the dominant detector technology for digital radiography and fluoroscopy. For low dose imaging, electronic noise from the amorphous silicon thin film transistor (TFT) array degrades imaging performance. We have fabricated the first prototype solid-state AMFPI using a uniform layer of avalanche amorphous selenium (a-Se) photoconductor to amplify the signal to eliminate the effect of electronic noise. We have previously developed a large area solid-state avalanche a-Se sensor structure referred to as High Gain Avalanche Rushing Photoconductor (HARP) capable of achieving gains of 75. In this work we successfully deposited this HARP structure onto a 24 x 30 cm² TFT array with a pixel pitch of 85 μm. An electric field (E_Se) up to 105 Vμm^-1 was applied across the a-Se layer without breakdown. Using the HARP layer as a direct detector, an X-ray avalanche gain of 15 ± 3 was achieved at E_Se = 105 Vμm^-1. In indirect mode with a 150 μm thick structured CsI scintillator, an optical gain of 76 ± 5 was measured at E_Se = 105 Vμm^-1. Image quality at low dose increases with the avalanche gain until the electronic noise is overcome at a constant exposure level of 0.76 mR. We demonstrate the success of a solid-state HARP X-ray imager as well as the largest active area HARP sensor to date.

The ability to improve visualization of structural information in digital radiography without increasing radiation exposures requires improved image quality across all spatial frequencies, especially at high frequencies. The detective quantum efficiency (DQE) as a function of spatial frequency quantifies image quality given by an x-ray detector. We present a method of increasing DQE at high spatial frequencies by improving the modulation transfer function (MTF) and reducing noise aliasing. The Apodized Aperature Pixel (AAP) design uses a detector with micro-elements to synthesize desired pixels and provide higher DQE than conventional detector designs. A cascaded system analysis (CSA) that incorporates x-ray interactions is used for comparison of the theoretical MTF, noise power spectrum (NPS), and DQE. Signal and noise transfer through the converter material is shown to consist of correlated an uncorrelated terms. The AAP design was shown to improve the DQE of both material types that have predominantly correlated transfer (such as CsI) and predominantly uncorrelated transfer (such as Se). Improvement in the MTF by 50% and the DQE by 100% at the sampling cut-off frequency is obtained when uncorrelated transfer is prevalent through the converter material. Optimizing high-frequency DQE results in improved image contrast and visualization of small structures and fine-detail.

Sequential CT studies present an excellent opportunity to apply prior-image-based reconstruction (PIBR) methods that leverage high-fidelity prior imaging studies to improve image quality and/or reduce x-ray exposure in subsequent studies. One major obstacle in using PIBR is that the initial and subsequent studies are often performed on different scanners (e.g. diagnostic CT followed by CBCT for interventional guidance); this results in mismatch in attenuation values due to hardware and software differences. While improved artifact correction techniques can potentially mitigate such differences, the correction is often incomplete. Here, we present an alternate strategy where the PIBR itself is used to mitigate these differences. We define a new penalty for the previously introduced PIBR called Reconstruction of Difference (RoD). RoD differs from many other PIBRs in that it reconstructs only changes in the anatomy (vs. reconstructing the current anatomy). Direct regularization of the difference image in RoD provides an opportunity to selectively penalize spatial frequencies of the difference image (e.g. low frequency differences associated with attenuation offsets and shading artifacts) without interfering with the variations in unchanged background image. We leverage this flexibility and introduce a novel regularization strategy using a generalized Fourier penalty within the RoD framework and develop the modified reconstruction algorithm. We evaluate the performance of the new approach in both simulation studies and in physical CBCT test-bench data. We find that generalized Fourier penalty can be highly effective in reducing low-frequency x-ray artifacts through selective suppression of spatial frequencies in the reconstructed difference image.

Until today, several software-based approaches to increase the temporal resolution in cardiac computed tomography by estimating motion vector fields (MVFs) have been developed. Thereunder, the majority are motion compensation algorithms, which estimate the MVFs employing a three-dimensional registration routine working on reconstructions of multiple cardiac phases.^{2, 6, 7, 12} We present an algorithm that requires nothing more than the data needed for a short scan reconstruction for motion estimation and motion-compensated reconstruction, which both are based on the reconstruction of volumes from a limited angular range.^{2, 3, 7, 8} Those partial angle reconstructions are centered at different time points during the short scan and have a temporal resolution of about 10ms each. The MVFs are estimated by a constrained cost function optimization routine employing a motion artifact measuring cost function. During optimization, the MVFs are applied directly by warping the partial angle reconstructions, and the motion compensation is established by simply adding the shifted images. In order to enforce smooth vector fields and keep the number of parameters low, the motion is modeled by a low degree polynomial. Furthermore, to find a good estimation of the MVFs even in phases with rapid cardiac motion, the constrained optimization is re-initialized multiple times. The algorithm is validated with the help of a simulation study and applied to patient data, where motion- compensated reconstructions are performed in various cardiac phases. We show that the image quality can be improved, also in more rapid cardiac phases due to re-initialization of the optimization routine.

Lack of access to projection data from patient CT scans is a major limitation for development and validation of new reconstruction algorithms. To meet this critical need, we are building a library of CT patient projection data in an open and vendor-neutral format, DICOM-CT-PD, which is an extended DICOM format that contains sinogram data, acquisition geometry, patient information, and pathology identification. The library consists of scans of various types, including head scans, chest scans, abdomen scans, electrocardiogram (ECG)-gated scans, and dual-energy scans. For each scan, three types of data are provided, including DICOM-CT-PD projection data at various dose levels, reconstructed CT images, and a free-form text file. Several instructional documents are provided to help the users extract information from DICOM-CT-PD files, including a dictionary file for the DICOM-CT-PD format, a DICOM-CT-PD reader, and a user manual. Radiologist detection performance based on the reconstructed CT images is also provided. So far 328 head cases, 228 chest cases, and 228 abdomen cases have been collected for potential inclusion. The final library will include a selection of 50 head, chest, and abdomen scans each from at least two different manufacturers, and a few ECG-gated scans and dual-source, dual-energy scans. It will be freely available to academic researchers, and is expected to greatly facilitate the development and validation of CT reconstruction algorithms.

Iterative reconstruction methods have become very popular and show the potential to reduce dose. We present a limit to the maximum dose reduction possible with new reconstruction algorithms obtained by analyzing the information content of the raw data, assuming the reconstruction algorithm does not have a priori knowledge about the object or correlations between pixels. This limit applies to the task of estimating the density of a lesion embedded in a known background object, where the shape of the lesion is known but its density is not. Under these conditions, the density of the lesion can be estimated directly from the raw data in an optimal manner. This optimal estimate will meet or outperform the performance of any reconstruction method operating on the raw data, under the condition that the reconstruction method does not introduce a priori information. The raw data bound can be compared to the lesion density estimate from FBP in order to produce a limit on the dose reduction possible from new reconstruction algorithms. The possible dose reduction from iterative reconstruction varies with the object, but for a lesion embedded in the center of a water cylinder, it is less than 40%. Additionally, comparisons between iterative reconstruction and filtered backprojection are sometimes confounded by the effect of through-slice blurring in the iterative reconstruction. We analyzed the magnitude of the variance reduction brought about by through-slice blurring on scanners from two different vendors and found it to range between 11% and 48%.

When the scan field of view (SFOV) of a CT system is not large enough to enclose the entire cross-section of a patient, or the patient needs to be intentionally positioned partially outside the SFOV for certain clinical CT scans, truncation artifacts are often observed in the reconstructed CT images. Conventional wisdom to reduce truncation artifacts is to complete the truncated projection data via data extrapolation with different a priori assumptions. This paper presents a novel truncation artifact reduction method that directly works in the CT image domain. Specifically, a discriminative dictionary that includes a sub-dictionary of truncation artifacts and a sub-dictionary of non-artifact image information was used to separate a truncation artifact-contaminated image into two sub-images, one with reduced truncation artifacts, and the other one containing only the truncation artifacts. Both experimental phantom and retrospective human subject studies have been performed to characterize the performance of the proposed truncation artifact reduction method.

We present a novel reconstruction for motion correction of non-cardiac organs. With non-cooperative patients or in emergency case, breathing motion or motion of the skull may compromise image quality. Our algorithm is based on the optimization of either motion artefact metrics or data-driven metrics. This approach was successfully applied in cardiac CTA [1]. While motion correction of the coronary vessels requires a local motion model, global motion models are sufficient for organs like the lung or the skull. The parameter vector for the global affine motion is estimated iteratively, using the open source optimization library NLOPT. The image is updated using motion compensated reconstruction in each of the iterations. Evaluation of the metric value, e.g. the image entropy, provides information for the next iteration loop. After reaching the fixed point of the iteration, the final motion parameters are used for a motion-compensated full quality reconstruction. In head imaging the motion model is based on translation and rotation, in thoracic imaging the rotation is replaced by non-isotropic scaling in all three dimensions. We demonstrate the efficiency of the method in thoracic imaging by evaluating PET-CT data from free-breathing patients. In neuro imaging, data from stroke patients showing skull tremor were analyzed. It was shown that motion artefacts can be largely reduced and spatial resolution was restored. In head imaging, similar results can be obtained using motion artefact metrics or data-driven metrics. In case of image-based metrics, the entropy of the image proved to be superior. Breathing motion could also be significantly reduced using entropy metric. However, in this case data driven metrics cannot be applied because the line integrals associated to the ROI of the lung have to be computed using the local ROI mechanism [2] It was shown that the lung signal is corrupted by signals originating from the complement of the lung. Thus a meaningful optimization of a data-driven cost function is not possible.

Photon-counting CT (PCCT) may yield potential value for many clinical applications due to its relative immunity to electronic noise, increased geometric efficiency relative to current scintillating detectors, and the ability to resolve energy information about the detected photons. However, there are a large number of parameters that require optimization, particularly the energy thresholds configuration. Fast and accurate estimation of signal and noise in PCCT can benefit the optimization of acquisition parameters for specific diagnostic tasks. Based on the acquisition parameters and detector response of our research PCCT system, we derived mathematical models for both signal and noise. The signal model took the tube spectrum, beam filtration, object attenuation, water beam hardening, and detector response into account. The noise model considered the relationship between noise and radiation dose, as well as the propagation of noise as threshold data are subtracted to yield energy bin data. To determine the absolute noise value, a noise look-up table (LUT) was acquired using a limited number of calibration scans. The noise estimation algorithm then used the noise LUT to estimate noise for scans with a variety of combination of energy thresholds, dose levels, and object attenuation. Validation of the estimation algorithms was performed on our whole-body research PCCT system using semianthropomorphic water phantoms and solutions of calcium and iodine. The algorithms achieved accurate estimation of signal and noise for a variety of scanning parameter combinations. The proposed method can be used to optimize energy thresholds configuration for many clinical applications of PCCT.

The energy resolving capabilities of Photon Counting Detectors (PCD) in Computed Tomography (CT) facilitate energy-sensitive measurements. The provided image-information can be processed with Dual Energy and Multi Energy algorithms. A research PCD-CT firstly allows acquiring images with a close to clinical configuration of both the X-ray tube and the CT-detector. In this study, two algorithms (Material Decomposition and Virtual Non-Contrast-imaging (VNC)) are applied on a data set acquired from an anesthetized rabbit scanned using the PCD-CT system. Two contrast agents (CA) are applied: A gadolinium (Gd) based CA used to enhance contrasts for vascular imaging, and xenon (Xe) and air as a CA used to evaluate local ventilation of the animal's lung. Four different images are generated: a) A VNC image, suppressing any traces of the injected Gd imitating a native scan, b) a VNC image with a Gd-image as an overlay, where contrast enhancements in the vascular system are highlighted using colored labels, c) another VNC image with a Xe-image as an overlay, and d) a 3D rendered image of the animal's lung, filled with Xe, indicating local ventilation characteristics. All images are generated from two images based on energy bin information. It is shown that a modified version of a commercially available dual energy software framework is capable of providing images with diagnostic value obtained from the research PCD-CT system.

Shot noise processes are omnipresent in physics and many of their properties have been extensively studied in the past, including the particular problem of level crossing of shot noise. Energy-sensitive, photon-counting detectors using comparators to discriminate pulse-heights are currently heavily investigated for medical applications, e.g. for x-ray computed tomography and x-ray mammography. Surprisingly, no mention of the close relation between the two topics can be found in the literature on photon-counting detectors. In this paper, we point out the close analogy between level crossing of shot noise and the problem of determining count rates of photon- counting detectors subject to pulse pile-up. The latter is very relevant for obtaining precise forward models for photon-counting detectors operated under conditions of very high x-ray flux employed in clinical x-ray computed tomography. Although several attempts have been made to provide reasonably accurate, approximative models for the registered number of counts in x-ray detectors under conditions of high flux and arbitrary x-ray spectra, see, e.g., no exact, analytic solution is given in the literature for general continuous pulse shapes. In this paper we present such a solution for arbitrary response functions, x-ray spectra and continuous pulse shapes based on a result from the theory of level crossing. We briefly outline the theory of level crossing including the famous Rice theorem and translate from the language of level crossing to the language of photon-counting detection.

A high-resolution (HR) data collection mode has been introduced to a whole-body, research photon-counting-detector CT system installed in our laboratory. In this mode, 64 rows of 0.45 mm x 0.45 mm detector pixels were used, which corresponded to a pixel size of 0.25 mm x 0.25 mm at the iso-center. Spatial resolution of this HR mode was quantified by measuring the MTF from a scan of a 50 micron wire phantom. An anthropomorphic lung phantom, cadaveric swine lung, temporal bone and heart specimens were scanned using the HR mode, and image quality was subjectively assessed by two experienced radiologists. High spatial resolution of the HR mode was evidenced by the MTF measurement, with 15 lp/cm and 20 lp/cm at 10% and 2% modulation. Images from anthropomorphic phantom and cadaveric specimens showed clear delineation of small structures, such as lung vessels, lung nodules, temporal bone structures, and coronary arteries. Temporal bone images showed critical anatomy (i.e. stapes superstructure) that was clearly visible in the PCD system. These results demonstrated the potential application of this imaging mode in lung, temporal bone, and vascular imaging. Other clinical applications that require high spatial resolution, such as musculoskeletal imaging, may also benefit from this high resolution mode.

With many attractive attributes, photon counting detectors with many energy bins are being considered for clinical CT systems. In practice, a large amount of projection data acquired for multiple energy bins must be transferred in real time through slip rings and data storage subsystems, causing a bandwidth bottleneck problem. The higher resolution of these detectors and the need for faster acquisition additionally contribute to this issue. In this work, we introduce a new approach to lossless compression, specifically for projection data from photon counting detectors, by utilizing the dependencies in the multi-energy data. The proposed predictor estimates the value of a projection data sample as a weighted average of its neighboring samples and an approximation from other energy bins, and the prediction residuals are then encoded. Context modeling using three or four quantized local gradients is also employed to detect edge characteristics of the data. Using three simulated phantoms including a head phantom, compression of 2.3:1-2.4:1 was achieved. The proposed predictor using zero, three, and four gradient contexts was compared to JPEG-LS and the ideal predictor (noiseless projection data). Among our proposed predictors, three-gradient context is preferred with a compression ratio from Golomb coding 7% higher than JPEG-LS and only 3% lower than the ideal predictor. In encoder efficiency, the Golomb code with the proposed three-gradient contexts has higher compression than block floating point. We also propose a lossy compression scheme, which quantizes the prediction residuals with scalar uniform quantization using quantization boundaries that limit the ratio of quantization error variance to quantum noise variance. Applying our proposed predictor with three-gradient context, the lossy compression achieved a compression ratio of 3.3:1 but inserted a 2.1% standard deviation of error compared to that of quantum noise in reconstructed images. From the initial simulation results, the proposed algorithm shows good control over the bits needed to represent multienergy projection data.

There exists interest in designing a PET system with reduced detectors due to cost concerns, while not significantly compromising the PET utility. Recently developed optimization-based algorithms, which have demonstrated the potential clinical utility in image reconstruction from sparse CT data, may be used for enabling such design of innovative PET systems. In this work, we investigate a PET configuration with reduced number of detectors, and carry out preliminary studies from patient data collected by use of such sparse-PET configuration. We consider an optimization problem combining Kullback-Leibler (KL) data fidelity with an image TV constraint, and solve it by using a primal-dual optimization algorithm developed by Chambolle and Pock. Results show that advanced algorithms may enable the design of innovative PET configurations with reduced number of detectors, while yielding potential practical PET utilities.