Deformable 3D/3D CT-to-digital-tomosynthesis image registration in image-guided bronchoscopy interventions

Traditional navigational bronchoscopy procedures rely on preprocedural computed tomography (CT) and intraoperative chest radiography and cone-beam CT (CBCT) to biopsy peripheral lung lesions. This navigational approach is challenging due to the projective nature of radiography, and the high radiation dose, long imaging time, and large footprints of CBCT. Digital tomosynthesis (DTS) is considered an attractive alternative combining the advantages of radiography and CBCT. Only the depth resolution cannot match a full CBCT image due to the limited angle acquisition. To address this issue, preoperative CT is a good auxiliary in guiding bronchoscopy interventions. Nevertheless, CT-to-body divergence caused by anatomic changes and respiratory motion, hinders the effective use of CT imaging. To mitigate CT-to-body divergence, we propose a novel deformable 3D/3D CT-to-DTS registration algorithm employing a multistage, multiresolution approach and using affine and elastic B-spline transformation models with bone and lung mask images. A multiresolution strategy with a Gaussian image pyramid and a multigrid strategy within the B-spline model are applied. The normalized correlation coefficient is included in the cost function for the affine model and a multimetric weighted cost function is used for the B-spline model, with weights determined heuristically. Tested on simulated and real patient bronchoscopy data, the algorithm yields promising results. Assessed qualitatively by visual inspection and quantitatively by computing the Dice coefficient (DC) and the average symmetric surface distance (ASSD), the algorithm achieves mean DC of 0.82±0.05 and 0.74±0.05, and mean ASSD of 0.65±0.29mm and 0.93±0.43mm for simulated and real data, respectively. This algorithm lays the groundwork for CT-aided intraoperative DTS imaging in image-guided bronchoscopy interventions with future studies focusing on automated metric weight setting.

AAPM Task Group Report 238: 3D C-arms with volumetric imaging capability

This report reviews the image acquisition and reconstruction characteristics of C-arm Cone Beam Computed Tomography (C-arm CBCT) systems and provides guidance on quality control of C-arm systems with this volumetric imaging capability. The concepts of 3D image reconstruction, geometric calibration, image quality, and dosimetry covered in this report are also pertinent to CBCT for Image-Guided Radiation Therapy (IGRT). However, IGRT systems introduce a number of additional considerations, such as geometric alignment of the imaging at treatment isocenter, which are beyond the scope of the charge to the task group and the report. Section 1 provides an introduction to C-arm CBCT systems and reviews a variety of clinical applications. Section 2 briefly presents nomenclature specific or unique to these systems. A short review of C-arm fluoroscopy quality control (QC) in relation to 3D C-arm imaging is given in Section 3. Section 4 discusses system calibration, including geometric calibration and uniformity calibration. A review of the unique approaches and challenges to 3D reconstruction of data sets acquired by C-arm CBCT systems is give in Section 5. Sections 6 and 7 go in greater depth to address the performance assessment of C-arm CBCT units. First, Section 6 describes testing approaches and phantoms that may be used to evaluate image quality (spatial resolution and image noise and artifacts) and identifies several factors that affect image quality. Section 7 describes both free-in-air and in-phantom approaches to evaluating radiation dose indices. The methodologies described for assessing image quality and radiation dose may be used for annual constancy assessment and comparisons among different systems to help medical physicists determine when a system is not operating as expected. Baseline measurements taken either at installation or after a full preventative maintenance service call can also provide valuable data to help determine whether the performance of the system is acceptable. Collecting image quality and radiation dose data on existing phantoms used for CT image quality and radiation dose assessment, or on newly developed phantoms, will inform the development of performance criteria and standards. Phantom images are also useful for identifying and evaluating artifacts. In particular, comparing baseline data with those from current phantom images can reveal the need for system calibration before image artifacts are detected in clinical practice. Examples of artifacts are provided in Sections 4, 5, and 6.

Bayesian optimization to design a novel x-ray shaping device

PURPOSE

In radiation therapy, x-ray dose must be precisely sculpted to the tumor, while simultaneously avoiding surrounding organs at risk. This requires modulation of x-ray intensity in space and/or time. Typically, this is achieved using a multi leaf collimator (MLC)-a complex mechatronic device comprising over one hundred individually powered tungsten 'leaves' that move in or out of the radiation field as required. Here, an all-electronic x-ray collimation concept with no moving parts is presented, termed "SPHINX": Scanning Pencil-beam High-speed Intensity-modulated X-ray source. SPHINX utilizes a spatially distributed bremsstrahlung target and collimator array in conjunction with magnetic scanning of a high energy electron beam to generate a plurality of small x-ray "beamlets."

METHODS

A simulation framework was developed in Topas Monte Carlo incorporating a phase space electron source, transport through user defined magnetic fields, bremsstrahlung x-ray production, transport through a SPHINX collimator, and dose in water. This framework was completely parametric, meaning a simulation could be built and run for any supplied geometric parameters. This functionality was coupled with Bayesian optimization to find the best parameter set based on an objective function which included terms to maximize dose rate for a user defined beamlet width while constraining inter-channel cross talk and electron contamination. Designs for beamlet widths of 5, 7, and 10 mm² were generated. Each optimization was run for 300 iterations and took approximately 40 h on a 24-core computer. For the optimized 7-mm model, a simulation of all beamlets in water was carried out including a linear scanning magnet calibration simulation. Finally, a back-of-envelope dose rate formalism was developed and used to estimate dose rate under various conditions.

RESULTS

The optimized 5-, 7-, and 10-mm models had beamlet widths of 5.1 , 7.2 , and 10.1 mm² and dose rates of 3574, 6351, and 10 015 Gy/C, respectively. The reduction in dose rate for smaller beamlet widths is a result of both increased collimation and source occlusion. For the simulation of all beamlets in water, the scanning magnet calibration reduced the offset between the collimator channels and beam centroids from 2.9 ±1.9 mm to 0.01 ±0.03 mm. A slight reduction in dose rate of approximately 2% per degree of scanning angle was observed. Based on a back-of-envelope dose rate formalism, SPHINX in conjunction with next-generation linear accelerators has the potential to achieve substantially higher dose rates than conventional MLC-based delivery, with delivery of an intensity modulated 100 × 100 mm² field achievable in 0.9 to 10.6 s depending on the beamlet widths used.

CONCLUSIONS

Bayesian optimization was coupled with Monte Carlo modeling to generate SPHINX geometries for various beamlet widths. A complete Monte Carlo simulation for one of these designs was developed, including electron beam transport of all beamlets through scanning magnets, x-ray production and collimation, and dose in water. These results demonstrate that SPHINX is a promising candidate for sculpting radiation dose with no moving parts, and has the potential to vastly improve both the speed and robustness of radiotherapy delivery. A multi-beam SPHINX system may be a candidate for delivering magavoltage FLASH RT in humans.

Science and practice of imaging physics through 50 years of SPIE Medical Imaging conferences

Purpose: For 50 years now, SPIE Medical Imaging (MI) conferences have been the premier forum for disseminating and sharing new ideas, technologies, and concepts on the physics of MI. Approach: Our overarching objective is to demonstrate and highlight the major trajectories of imaging physics and how they are informed by the community and science present and presented at SPIE MI conferences from its inception to now. Results: These contributions range from the development of image science, image quality metrology, and image reconstruction to digital x-ray detectors that have revolutionized MI modalities including radiography, mammography, fluoroscopy, tomosynthesis, and computed tomography (CT). Recent advances in detector technology such as photon-counting detectors continue to enable new capabilities in MI. Conclusion: As we celebrate the past 50 years, we are also excited about what the next 50 years of SPIE MI will bring to the physics of MI.

Flat-panel conebeam CT in the clinic: history and current state

Research into conebeam CT concepts began as soon as the first clinical single-slice CT scanner was conceived. Early implementations of conebeam CT in the 1980s focused on high-contrast applications where concurrent high resolution ( < 200    μ m ), for visualization of small contrast-filled vessels, bones, or teeth, was an imaging requirement that could not be met by the contemporaneous CT scanners. However, the use of nonlinear imagers, e.g., x-ray image intensifiers, limited the clinical utility of the earliest diagnostic conebeam CT systems. The development of consumer-electronics large-area displays provided a technical foundation that was leveraged in the 1990s to first produce large-area digital x-ray detectors for use in radiography and then compact flat panels suitable for high-resolution and high-frame-rate conebeam CT. In this review, we show the concurrent evolution of digital flat panel (DFP) technology and clinical conebeam CT. We give a brief summary of conebeam CT reconstruction, followed by a brief review of the correction approaches for DFP-specific artifacts. The historical development and current status of flat-panel conebeam CT in four clinical areas-breast, fixed C-arm, image-guided radiation therapy, and extremity/head-is presented. Advances in DFP technology over the past two decades have led to improved visualization of high-contrast, high-resolution clinical tasks, and image quality now approaches the soft-tissue contrast resolution that is the standard in clinical CT. Future technical developments in DFPs will enable an even broader range of clinical applications; research in the arena of flat-panel CT shows no signs of slowing down.

Robust learning-based X-ray image denoising - potential pitfalls, their analysis and solutions

PURPOSE

Since guidance based on X-ray imaging is an integral part of interventional procedures, continuous efforts are taken towards reducing the exposure of patients and clinical staff to ionizing radiation. Even though a reduction in the X-ray dose may lower associated radiation risks, it is likely to impair the quality of the acquired images, potentially making it more difficult for physicians to carry out their procedures.

METHOD

We present a robust learning-based denoising strategy involving model- based simulations of low-dose X-ray images during the training phase. The method also utilizes a data-driven normalization step - based on an X-ray imaging model - to stabilize the mixed signal-dependent noise associated with X-ray images. We thoroughly analyze the method's sensitivity to a mismatch in dose levels used for training and application. We also study the impact of differing noise models used when training for low and very low-dose X-ray images on the denoising results.

RESULTS

A quantitative and qualitative analysis based on acquired phantom and clinical data has shown that the proposed learning-based strategy is stable across different dose levels and yields excellent denoising results, if an accurate noise model is applied. We also found that there can be severe artifacts when the noise characteristics of the training images are significantly different from those in the actual images to be processed. This problem can be especially acute at very low dose levels. During a thorough analysis of our experimental results, we further discovered that viewing the results from the perspective of denoising via thresholding of sub-band co efficients can be very beneficial to get a better understanding of the proposed learning-based denoising strategy.

CONCLUSION

The proposed learning-based denoising strategy provides scope for significant X-ray dose reduction without the loss of important image information if the characteristics of noise is accurately accounted for during the training ph.

Rigid and Non-rigid Motion Compensation in Weight-bearing CBCT of the Knee using Simulated Inertial Measurements

OBJECTIVE

Involuntary subject motion is the main source of artifacts in weight-bearing cone-beam CT of the knee. To achieve image quality for clinical diagnosis, the motion needs to be compensated. We propose to use inertial measurement units (IMUs) attached to the leg for motion estimation.

METHODS

We perform a simulation study using real motion recorded with an optical tracking system. Three IMU-based correction approaches are evaluated, namely rigid motion correction, non-rigid 2D projection deformation and non-rigid 3D dynamic reconstruction. We present an initialization process based on the system geometry. With an IMU noise simulation, we investigate the applicability of the proposed methods in real applications.

RESULTS

All proposed IMU-based approaches correct motion at least as good as a state-of-the-art marker-based approach. The structural similarity index and the root mean squared error between motion-free and motion corrected volumes are improved by 24-35% and 78-85%, respectively, compared with the uncorrected case. The noise analysis shows that the noise levels of commercially available IMUs need to be improved by a factor of 10⁵ which is currently only achieved by specialized hardware not robust enough for the application.

CONCLUSION

Our simulation study confirms the feasibility of this novel approach and defines improvements necessary for a real application.

SIGNIFICANCE

The presented work lays the foundation for IMU-based motion compensation in cone-beam CT of the knee and creates valuable insights for future developments.

X-Ray Scatter Estimation Using Deep Splines

X-ray scatter compensation is a very desirable technique in flat-panel X-ray imaging and cone-beam computed tomography. State-of-the-art U-net based scatter removal approaches yielded promising results. However, as there are no physics' constraints applied to the output of the U-Net, it cannot be ruled out that it yields spurious results. Unfortunately, in the context of medical imaging, those may be misleading and could lead to wrong conclusions. To overcome this problem, we propose to embed B-splines as a known operator into neural networks. This inherently constrains their predictions to well-behaved and smooth functions. In a study using synthetic head and thorax data as well as real thorax phantom data, we found that our approach performed on par with U-net when comparing both algorithms based on quantitative performance metrics. However, our approach not only reduces runtime and parameter complexity, but we also found it much more robust to unseen noise levels. While the U-net responded with visible artifacts, the proposed approach preserved the X-ray signal's frequency characteristics.

Special Section Guest Editorial: Computed tomography (CT) at 50 years

Guest editors Patrick La Riviere, Rebecca Fahrig, and Norbert Pelc introduce the JMI Special Section Celebrating X-Ray Computed Tomography at 50.

3D Non-Rigid Alignment of Low-Dose Scans Allows to Correct for Saturation in Lower Extremity Cone-Beam CT

Detector saturation in cone-beam computed tomography occurs when an object of highly varying shape and material composition is imaged using an automatic exposure control (AEC) system. When imaging a subject's knees, high beam energy ensures the visibility of internal structures but leads to overexposure in less dense border regions. In this work, we propose to use an additional low-dose scan to correct the saturation artifacts of AEC scans. Overexposed pixels are identified in the projection images of the AEC scan using histogram-based thresholding. The saturation-free pixels from the AEC scan are combined with the skin border pixels of the low-dose scan prior to volumetric reconstruction. To compensate for patient motion between the two scans, a 3D non-rigid alignment of the projection images in a backward-forward-projection process based on fiducial marker positions is proposed. On numerical simulations, the projection combination improved the structural similarity index measure from 0.883 to 0.999. Further evaluations were performed on two in vivo subject knee acquisitions, one without and one with motion between the AEC and low-dose scans. Saturation-free reference images were acquired using a beam attenuator. The proposed method could qualitatively restore the information of peripheral tissue structures. Applying the 3D non-rigid alignment made it possible to use the projection images with inter-scan subject motion for projection image combination. The increase in radiation exposure due to the additional low-dose scan was found to be negligibly low. The presented methods allow simple but effective correction of saturation artifacts.

Deep action learning enables robust 3D segmentation of body organs in various CT and MRI images

In this study, we propose a novel point cloud based 3D registration and segmentation framework using reinforcement learning. An artificial agent, implemented as a distinct actor based on value networks, is trained to predict the optimal piece-wise linear transformation of a point cloud for the joint tasks of registration and segmentation. The actor network estimates a set of plausible actions and the value network aims to select the optimal action for the current observation. Point-wise features that comprise spatial positions (and surface normal vectors in the case of structured meshes), and their corresponding image features, are used to encode the observation and represent the underlying 3D volume. The actor and value networks are applied iteratively to estimate a sequence of transformations that enable accurate delineation of object boundaries. The proposed approach was extensively evaluated in both segmentation and registration tasks using a variety of challenging clinical datasets. Our method has fewer trainable parameters and lower computational complexity compared to the 3D U-Net, and it is independent of the volume resolution. We show that the proposed method is applicable to mono- and multi-modal segmentation tasks, achieving significant improvements over the state-of-the-art for the latter. The flexibility of the proposed framework is further demonstrated for a multi-modal registration application. As we learn to predict actions rather than a target, the proposed method is more robust compared to the 3D U-Net when dealing with previously unseen datasets, acquired using different protocols or modalities. As a result, the proposed method provides a promising multi-purpose segmentation and registration framework, particular in the context of image-guided interventions.

Magnetic modeling of actively shielded rotating MRI magnets in the presence of environmental steel

Rotating MRI systems could enable novel integrated medical devices such as MRI-Linacs, MRI-xray-angiography systems, and MRI-proton therapy systems. This work aimed to investigate the feasibility of rotating actively shielded superconducting MRI magnets in the presence of environmental steel-in particular, construction steel in the floor of the installation site. Two magnets were investigated: a 1.0 T split bore magnet, and a 1.5 T closed bore magnet. Each magnet was scaled to emulate field strengths of 0.5, 1.0, and 1.5 T. Finite Element Modeling was used to simulate these magnets in the presence of a 3 × 4 m steel plate located 1250 mm or 1400 mm below the isocenter. There are two possible rotation directions: around the longitudinal (z) axis or around the transverse (x) axis. Each model was solved for rotation angles between 0 and 360° in 30° intervals around each of these axes. For each simulation, a 300 mm DSV was extracted and decomposed into spherical harmonics. For the closed-bore magnet, total induced perturbation for the zero degree rotation angle was 223, 432, and 562 μT peak-to-peak (pk-pk) for the 0.5, 1.0, and 1.5 T models respectively (steel at 1250 mm). For the split-bore magnet, the same numbers were 1477, 16747, and 1766 μT. The substantially higher perturbation for the split-bore magnet can be traced to its larger fringe field. For rotation around the z-axis, total perturbation does not change as a function of angle but is exchanged between different harmonics. For rotation around the x-axis, total perturbation is different at each rotation angle. For the closed bore magnet, maximum perturbations occurred for a 90° rotation around the transverse axis. For the split-bore magnet, the opposite was observed, with the same 90° rotation yielding total perturbation lower than the conventional position. In all cases, at least 95% of the total perturbation was composed of 1st and 2nd order harmonics. The presence of environmental steel poses a major challenge to the realization of an actively shielded rotating superconducting MRI system, requiring some novel form of shimming. Possible shimming strategies are discussed at length.

XDose: toward online cross-validation of experimental and computational X-ray dose estimation

PURPOSE

As the spectrum of X-ray procedures has increased both for diagnostic and for interventional cases, more attention is paid to X-ray dose management. While the medical benefit to the patient outweighs the risk of radiation injuries in almost all cases, reproducible studies on organ dose values help to plan preventive measures helping both patient as well as staff. Dose studies are either carried out retrospectively, experimentally using anthropomorphic phantoms, or computationally. When performed experimentally, it is helpful to combine them with simulations validating the measurements. In this paper, we show how such a dose simulation method, carried out together with actual X-ray experiments, can be realized to obtain reliable organ dose values efficiently.

METHODS

A Monte Carlo simulation technique was developed combining down-sampling and super-resolution techniques for accelerated processing accompanying X-ray dose measurements. The target volume is down-sampled using the statistical mode first. The estimated dose distribution is then up-sampled using guided filtering and the high-resolution target volume as guidance image. Second, we present a comparison of dose estimates calculated with our Monte Carlo code experimentally obtained values for an anthropomorphic phantom using metal oxide semiconductor field effect transistor dosimeters.

RESULTS

We reconstructed high-resolution dose distributions from coarse ones (down-sampling factor 2 to 16) with error rates ranging from 1.62 % to 4.91 %. Using down-sampled target volumes further reduced the computation time by 30 % to 60 %. Comparison of measured results to simulated dose values demonstrated high agreement with an average percentage error of under [Formula: see text] for all measurement points.

CONCLUSIONS

Our results indicate that Monte Carlo methods can be accelerated hardware-independently and still yield reliable results. This facilitates empirical dose studies that make use of online Monte Carlo simulations to easily cross-validate dose estimates on-site.

Data-driven estimation of noise variance stabilization parameters for low-dose x-ray images

PURPOSE

Denoising x-ray images corrupted by signal-dependent mixed noise is usually approached either by considering noise statistics directly or by using noise variance stabilization (NVS) techniques. An advantage of the latter is that the noise variance can be stabilized to a known constant throughout the image, facilitating the application of denoising algorithms designed for the removal of additive Gaussian noise. A well-performing NVS is the generalized Anscombe transform (GAT). To calculate the GAT, the system gain as well as the variance of electronic noise are required. Unfortunately, these parameters are difficult to predict from the x-ray tube settings in clinical practice, because the system gain observed at the detector depends on the beam hardening caused by the patient.

MATERIALS AND METHODS

We propose a data-driven method for estimating the parameters required to carry out an NVS using the GAT. It utilizes the energy compaction property of the discrete cosine transform to obtain the NVS parameters using a robust regression approach relying on a linear Poisson-Gaussian model. The method has been experimentally validated with respect to beam hardening as well as denoising performance for different dose and scatter levels.

RESULTS

Across a range of low-dose x-ray settings, the proposed robust regression approach has estimated both system gain and electronic noise level with an average error of only 4.2%. When used to perform a GAT followed by the denoising of low-dose x-ray images, performance gains of 5% for peak-signal-to-noise ratio and 4% for structural similarity index can be obtained.

CONCLUSION

The parameters needed to calculate the GAT can be estimated efficiently and robustly using a data-driven approach. The improved parameter estimation method facilitates a more accurate GAT-based NVS and, hence, better denoising of low-dose x-ray images when algorithms designed for additive Gaussian noise are applied.

Automatic Orientation Estimation of Inertial Sensors in C-Arm CT Projections

To obtain CT images of the knee joint in a more lifelike position, data acquisition can be performed with patients in standing rather than in lying position. However, in that situation, people tend to show involuntary motion. One possibility to compensate for that motion is the use of Inertial Measurement Units, that capture the accelerations during the scan. For this purpose, their local coordinate system needs to be known. An estimation based on the SIFT algorithm was implemented and compared to an existing approach that uses the Fast Radial Symmetry transform and to expert labels for evaluation. The SIFT method showed to be superior to the existing approach as it could extract stable feature points from the projections that were used to estimate the three-dimensional coordinate system in a reliable manner. The final algorithm achieved a mean euclidean distance of 2.61 mm between the calculated position of the origin and the assumed ground truth by the expert labels.

Physics-based spectral compensation algorithm for x-ray CT with primary modulator

X-ray computed tomography (CT) scatter correction using primary modulator has been continuously developed over the past years, with progress in improving the performance of scatter correction. In this work, we further advance the primary modulator technique towards practical applications where the spectral nonuniformity caused by the modulator continues to be a challenging problem. A physics-based spectral compensation algorithm is proposed to adaptively correct for the spectral nonuniformity, and hence to reduce the resultant ring artifacts on reconstructed CT images. First, an initial spectrum of the CT system without the primary modulator is modeled using an understanding of x-ray CT physics, and optimized by an expectation maximization method; then, the optimized estimation of the initial spectrum is utilized to adaptively calculate the effective modulator thickness from measured transmissions of the primary modulator at each detector element, leading to a set of new spectra that is able to capture the nonuniform spectral distribution of the primary modulator; finally, using the modulator-modeled spectrum, a beam hardening mapping function is generated and beam hardening correction is applied to CT projections. A CatPhan600 phantom and an anthropomorphic thorax phantom were scanned with three different primary modulators to evaluate the approach. For the Catphan phantom, the spectral compensation algorithm efficiently removes the ring (and band) artifacts that otherwise dominate the reconstructed CT image. For the three modulators with nominal copper thickness of 52.5, 105 and 210 [Formula: see text]m, our method reduces the CT number nonuniformity from 147.9, 436.2 and 696.4 Hounsfield units (HU) to 14.6, 26.2 and 13.6 HU, respectively, close to that of the reference image (i.e. 7.5 HU). For the thorax phantom, the ring artifacts are also suppressed significantly on the transaxial image; on the sagittal image, the alternating black-and-white patterns are largely removed, with the CT number nonuniformity being reduced from 282.0 HU to 38.5 HU.

Registration of vascular structures using a hybrid mixture model

PURPOSE

Morphological changes to anatomy resulting from invasive surgical procedures or pathology, typically alter the surrounding vasculature. This makes it useful as a descriptor for feature-driven image registration in various clinical applications. However, registration of vasculature remains challenging, as vessels often differ in size and shape, and may even miss branches, due to surgical interventions or pathological changes. Furthermore, existing vessel registration methods are typically designed for a specific application. To address this limitation, we propose a generic vessel registration approach useful for a variety of clinical applications, involving different anatomical regions.

METHODS

A probabilistic registration framework based on a hybrid mixture model, with a refinement mechanism to identify missing branches (denoted as HdMM+) during vasculature matching, is introduced. Vascular structures are represented as 6-dimensional hybrid point sets comprising spatial positions and centerline orientations, using Student's t-distributions to model the former and Watson distributions for the latter.

RESULTS

The proposed framework is evaluated for intraoperative brain shift compensation, and monitoring changes in pulmonary vasculature resulting from chronic lung disease. Registration accuracy is validated using both synthetic and patient data. Our results demonstrate, HdMM+ is able to reduce more than [Formula: see text] of the initial error for both applications, and outperforms the state-of-the-art point-based registration methods such as coherent point drift and Student's t-distribution mixture model, in terms of mean surface distance, modified Hausdorff distance, Dice and Jaccard scores.

CONCLUSION

The proposed registration framework models complex vascular structures using a hybrid representation of vessel centerlines, and accommodates intricate variations in vascular morphology. Furthermore, it is generic and flexible in its design, enabling its use in a variety of clinical applications.

Pulse pileup analysis for a double-sided silicon strip detector using variable pulse shapes

Due to pulse pileup, photon counting detectors (PCDs) suffer from count loss and energy distortion when operating in high count rate environments. In this paper, we studied the pulse pileup of a double-sided silicon strip detector (DSSSD) to evaluate its potential application in a mammography system. We analyzed the pulse pileup using pulses of varied shapes, where the shape of the pulse depends on the location of photon interaction within the detector. To obtain the shaped pulses, first, transient currents for photons interacting at different locations were simulated using a Technology Computer-Aided Design (TCAD) software. Next, the currents were shaped by a CR-RC2 shaping circuit, calculated using Simulink. After obtaining these pulses, both the different orders of pileup and the energy spectrum were calculated by taking into account the following two factors: 1) spatial distribution of photon interactions within the detector, and 2) time interval distribution between successive photons under a given photon flux. We found that for a DSSSD with thickness of 300 μm, pitch of 25 μm and strip length of 1 cm, under a bias voltage of 50 V, the variable pulse shape model predicts the fraction free of pileup can be > 90 % under a photon flux of 3.75 Mcps/mm2. The double-sided silicon-strip detector is a promising candidate for digital mammography applications.

Preliminary results of DSA denoising based on a weighted low-rank approach using an advanced neurovascular replication system

PURPOSE

2D digital subtraction angiography (DSA) has become an important technique for interventional neuroradiology tasks, such as detection and subsequent treatment of aneurysms. In order to provide high-quality DSA images, usually undiluted contrast agent and a high X-ray dose are used. The iodinated contrast agent puts a burden on the patients' kidneys while the use of high-dose X-rays expose both patients and medical staff to a considerable amount of radiation. Unfortunately, reducing either the X-ray dose or the contrast agent concentration usually results in a sacrifice of image quality.

MATERIALS AND METHODS

To denoise a frame, the proposed spatiotemporal denoising method utilizes the low-rank nature of a spatially aligned temporal sequence where variation is introduced by the flow of contrast agent through a vessel tree of interest. That is, a constrained weighted rank-1 approximation of the stack comprising the frame to be denoised and its temporal neighbors is computed where the weights are used to prevent the contribution of non-similar pixels toward the low-rank approximation. The method has been evaluated using a vascular flow phantom emulating cranial arteries into which contrast agent can be manually injected (Vascular Simulations Replicator, Vascular Simulations, Stony Brook NY, USA). For the evaluation, image sequences acquired at different dose levels as well as different contrast agent concentrations have been used.

RESULTS

Qualitative and quantitative analyses have shown that with the proposed approach, the dose and the concentration of the contrast agent could both be reduced by about 75%, while maintaining the required image quality. Most importantly, it has been observed that the DSA images obtained using the proposed method have the closest resemblance to typical DSA images, i.e., they preserve the typical image characteristics best.

CONCLUSION

Using the proposed denoising approach, it is possible to improve the image quality of low-dose DSA images. This improvement could enable both a reduction in contrast agent and radiation dose when acquiring DSA images, thereby benefiting patients as well as clinicians. Since the resulting images are free from artifacts and as the inherent characteristics of the images are also preserved, the proposed method seems to be well suited for clinical images as well.

A photon recycling approach to the denoising of ultra-low dose X-ray sequences

PURPOSE

Clinical procedures that make use of fluoroscopy may expose patients as well as the clinical staff (throughout their career) to non-negligible doses of radiation. The potential consequences of such exposures fall under two categories, namely stochastic (mostly cancer) and deterministic risks (skin injury). According to the "as low as reasonably achievable" principle, the radiation dose can be lowered only if the necessary image quality can be maintained.

METHODS

Our work improves upon the existing patch-based denoising algorithms by utilizing a more sophisticated noise model to exploit non-local self-similarity better and this in turn improves the performance of low-rank approximation. The novelty of the proposed approach lies in its properly designed and parameterized noise model and the elimination of initial estimates. This reduces the computational cost significantly.

RESULTS

The algorithm has been evaluated on 500 clinical images (7 patients, 20 sequences, 3 clinical sites), taken at ultra-low dose levels, i.e. 50% of the standard low dose level, during electrophysiology procedures. An average improvement in the contrast-to-noise ratio (CNR) by a factor of around 3.5 has been found. This is associated with an image quality achieved at around 12 (square of 3.5) times the ultra-low dose level. Qualitative evaluation by X-ray image quality experts suggests that the method produces denoised images that comply with the required image quality criteria.

CONCLUSION

The results are consistent with the number of patches used, and they demonstrate that it is possible to use motion estimation techniques and "recycle" photons from previous frames to improve the image quality of the current frame. Our results are comparable in terms of CNR to Video Block Matching 3D-a state-of-the-art denoising method. But qualitative analysis by experts confirms that the denoised ultra-low dose X-ray images obtained using our method are more realistic with respect to appearance.

Thermal limits on MV x-ray production by bremsstrahlung targets in the context of novel linear accelerators

PURPOSE

To study the impact of target geometrical and linac operational parameters, such as target material and thickness, electron beam size, repetition rate, and mean current on the ability of the radiotherapy treatment head to deliver high-dose-rate x-ray irradiation in the context of novel linear accelerators capable of higher repetition rates/duty cycle than conventional clinical linacs.

METHODS

The depth dose in a water phantom without a flattening filter and heat deposition in an x-ray target by 10 MeV pulsed electron beams were calculated using the Monte-Carlo code MCNPX, and the transient temperature behavior of the target was simulated by ANSYS. Several parameters that affect both the dose distribution and temperature behavior were investigated. The target was tungsten with a thickness ranging from 0 to 3 mm and a copper heat remover layer. An electron beam with full width at half maximum (FWHM) between 0 and3 mm and mean current of 0.05-2 mA was used as the primary beam at repetition rates of 100, 200, 400, and 800 Hz.

RESULTS

For a 10 MeV electron beam with FWHM of 1 mm, pulse length of 5 μs, by using a thin tungsten target with thickness of 0.2 mm instead of 1 mm, and by employing a high repetition rate of 800 Hz instead of 100 Hz, the maximum dose rate delivered can increase two times from 0.57 to 1.16 Gy/s. In this simple model, the limiting factor on dose rate is the copper heat remover's softening temperature, which was considered to be 500°C in our study.

CONCLUSIONS

A high dose rate can be obtained by employing thin targets together with high repetition rate electron beams enabled by novel linac designs, whereas the benefit of thin targets is marginal at conventional repetition rates. Next generation linacs used to increase dose rate need different target designs compared to conventional linacs.

Motion compensation for cone-beam CT using Fourier consistency conditions

In cone-beam CT, involuntary patient motion and inaccurate or irreproducible scanner motion substantially degrades image quality. To avoid artifacts this motion needs to be estimated and compensated during image reconstruction. In previous work we showed that Fourier consistency conditions (FCC) can be used in fan-beam CT to estimate motion in the sinogram domain. This work extends the FCC to [Formula: see text] cone-beam CT. We derive an efficient cost function to compensate for [Formula: see text] motion using [Formula: see text] detector translations. The extended FCC method have been tested with five translational motion patterns, using a challenging numerical phantom. We evaluated the root-mean-square-error and the structural-similarity-index between motion corrected and motion-free reconstructions. Additionally, we computed the mean-absolute-difference (MAD) between the estimated and the ground-truth motion. The practical applicability of the method is demonstrated by application to respiratory motion estimation in rotational angiography, but also to motion correction for weight-bearing imaging of knees. Where the latter makes use of a specifically modified FCC version which is robust to axial truncation. The results show a great reduction of motion artifacts. Accurate estimation results were achieved with a maximum MAD value of 708 μm and 1184 μm for motion along the vertical and horizontal detector direction, respectively. The image quality of reconstructions obtained with the proposed method is close to that of motion corrected reconstructions based on the ground-truth motion. Simulations using noise-free and noisy data demonstrate that FCC are robust to noise. Even high-frequency motion was accurately estimated leading to a considerable reduction of streaking artifacts. The method is purely image-based and therefore independent of any auxiliary data.

Comparison of Different Approaches for Measuring Tibial Cartilage Thickness

Osteoarthritis is a degenerative disease affecting bones and cartilage especially in the human knee. In this context, cartilage thickness is an indicator for knee cartilage health. Thickness measurements are performed on medical images acquired in-vivo. Currently, there is no standard method agreed upon that defines a distance measure in articular cartilage. In this work, we present a comparison of different methods commonly used in literature. These methods are based on nearest neighbors, surface normal vectors, local thickness and potential field lines. All approaches were applied to manual segmentations of tibia and lateral and medial tibial cartilage performed by experienced raters. The underlying data were contrast agent-enhanced cone-beam C-arm CT reconstructions of one healthy subject's knee. The subject was scanned three times, once in supine position and two times in a standing weight-bearing position. A comparison of the resulting thickness maps shows similar distributions and high correlation coefficients between the approaches above 0.90. The nearest neighbor method results on average in the lowest cartilage thickness values, while the local thickness approach assigns the highest values. We showed that the different methods agree in their thickness distribution. The results will be used for a future evaluation of cartilage change under weight-bearing conditions.

Intraoperative Imaging Modalities and Compensation for Brain Shift in Tumor Resection Surgery

Intraoperative brain shift during neurosurgical procedures is a well-known phenomenon caused by gravity, tissue manipulation, tumor size, loss of cerebrospinal fluid (CSF), and use of medication. For the use of image-guided systems, this phenomenon greatly affects the accuracy of the guidance. During the last several decades, researchers have investigated how to overcome this problem. The purpose of this paper is to present a review of publications concerning different aspects of intraoperative brain shift especially in a tumor resection surgery such as intraoperative imaging systems, quantification, measurement, modeling, and registration techniques. Clinical experience of using intraoperative imaging modalities, details about registration, and modeling methods in connection with brain shift in tumor resection surgery are the focuses of this review. In total, 126 papers regarding this topic are analyzed in a comprehensive summary and are categorized according to fourteen criteria. The result of the categorization is presented in an interactive web tool. The consequences from the categorization and trends in the future are discussed at the end of this work.

Development and testing of a database of NIH research funding of AAPM members: A report from the AAPM Working Group for the Development of a Research Database (WGDRD)

PURPOSE

To produce and maintain a database of National Institutes of Health (NIH) funding of the American Association of Physicists in Medicine (AAPM) members, to perform a top-level analysis of these data, and to make these data (hereafter referred to as the AAPM research database) available for the use of the AAPM and its members.

METHODS

NIH-funded research dating back to 1985 is available for public download through the NIH exporter website, and AAPM membership information dating back to 2002 was supplied by the AAPM. To link these two sources of data, a data mining algorithm was developed in Matlab. The false-positive rate was manually estimated based on a random sample of 100 records, and the false-negative rate was assessed by comparing against 99 member-supplied PI_ID numbers. The AAPM research database was queried to produce an analysis of trends and demographics in research funding dating from 2002 to 2015.

RESULTS

A total of 566 PI_ID numbers were matched to AAPM members. False-positive and -negative rates were respectively 4% (95% CI: 1-10%, N = 100) and 10% (95% CI: 5-18%, N = 99). Based on analysis of the AAPM research database, in 2015 the NIH awarded $USD 110M to members of the AAPM. The four NIH institutes which historically awarded the most funding to AAPM members were the National Cancer Institute, National Institute of Biomedical Imaging and Bioengineering, National Heart Lung and Blood Institute, and National Institute of Neurological Disorders and Stroke. In 2015, over 85% of the total NIH research funding awarded to AAPM members was via these institutes, representing 1.1% of their combined budget. In the same year, 2.0% of AAPM members received NIH funding for a total of $116M, which is lower than the historic mean of $120M (in 2015 USD).

CONCLUSIONS

A database of NIH-funded research awarded to AAPM members has been developed and tested using a data mining approach, and a top-level analysis of funding trends has been performed. Current funding of AAPM members is lower than the historic mean. The database will be maintained by members of the Working group for the development of a research database (WGDRD) on an annual basis, and is available to the AAPM, its committees, working groups, and members for download through the AAPM electronic content website. A wide range of questions regarding financial and demographic funding trends can be addressed by these data. This report has been approved for publication by the AAPM Science Council.

Improved targeting accuracy of lung tumor biopsies with scanning-beam digital x-ray tomosynthesis image guidance

PURPOSE

Electromagnetic navigation bronchoscopy (ENB) provides improved targeting accuracy during transbronchial biopsies of suspicious nodules. The greatest weakness of ENB-based guidance is the registration divergence that exists between the planning CT, acquired days or weeks before the intervention, and the patient on the table on the day of the intervention. Augmenting ENB guidance with real-time tomosynthesis imaging during the intervention could mitigate the divergence and further improve the yield of ENB-guided transbronchial biopsies. The real-time tomosynthesis prototype, the scanning-beam digital x-ray (SBDX) system, does not currently display images reconstructed by the iterative algorithm that was developed for this lung imaging application. A protocol using fiducial markers was therefore implemented to permit evaluation of potential improvements that would be provided by the SBDX system in a clinical setting.

METHODS

Ten 7 mm lesions (5 per side) were injected into the periphery of each of four preserved pig lungs. The lungs were then placed in a vacuum chamber that permitted simulation of realistic motion and deformation due to breathing. Standard clinical CT scans of the pig lung phantoms were acquired and reconstructed with isotropic resolution of 0.625 mm. Standard ENB-guided biopsy procedures including target identification, path planning, CT-to-lung registration and navigation to the lesion were carried out, and a fiducial marker was placed at the location at which a biopsy would have been acquired. The channel-to-target distance provided by the ENB system prior to fiducial placement was noted. The lung phantoms were then imaged using the SBDX system, and using high-resolution conebeam CT. The distance between the fiducial marker tip and the lesion was measured in SBDX images and in the gold-standard conebeam-CT images. The channel-to-target divergence predicted by the ENB system and measured in the SBDX images was compared to the gold standard to determine if improved targeting accuracy could be achieved using SBDX image guidance.

RESULTS

As expected, the ENB system showed poorer targeting accuracy for small peripheral nodules. Only 20 nodules of the 40 injected could be adequately reached using ENB guidance alone. The SBDX system was capable of visualizing these small lesions, and measured fiducial-to-target distances on SBDX agreed well with measurements in gold-standard conebeam-CT images (p = 0.0001). The correlation between gold-standard conebeam-CT distances and predicted fiducial-to-target distances provided by the ENB system was poor (p = 0.72), primarily due to inaccurate ENB CT-to-body registration and movement due to breathing.

CONCLUSIONS

The SBDX system permits visualization of small lung nodules, as well as accurate measurement of channel-to-target distances. Combined use of ENB with SBDX real-time image guidance could improve accuracy and yield of biopsies, particularly of those lesions located in the periphery of the lung.

Interventional dual-energy imaging-Feasibility of rapid kV-switching on a C-arm CT system

PURPOSE

In the last years, dual-energy CT imaging has shown clinical value, thanks to its ability to differentiate materials based on their atomic number and to exploit different properties of images acquired at two different energies. C-arm CT systems are used to guide procedures in the interventional suite. Until now, there are no commercially available systems that employ dual-energy material decomposition. This paper explores the feasibility of implementing a fast kV-switching technique on a clinically available angiographic system for acquiring dual-energy C-arm CT images.

METHODS

As an initial proof of concept, a fast kV-switching approach was implemented on an angiographic C-arm system and the peak tube voltage during 3D rotational scans was measured. The tube voltage measurements during fast kV-switching scans were compared to corresponding measurements on kV-constant scans. Additionally, to prove stability of the requested exposure parameters, the accuracy of the delivered tube current and pulse width were also recorded and compared. In a first phantom experiment, the voxel intensity values of the individual tube voltage components of the fast kV-switching scans were compared to their corresponding kV-constant scans. The same phantom was used for a simple material decomposition between different iodine concentrations and pure water using a fast kV-switching protocol of 81 and 125 kV. In the last experiment, the same kV-switching protocol as in the phantom scan was used in an in vivo pig study to demonstrate the clinical feasibility.

RESULTS

During rapid kV-switching acquisitions, the measured tube voltage of the x-ray tube during fast switching scans has an absolute deviation of 0.23 ± 0.13 kV compared to the measured tube voltage produced during kV-constant acquisitions. The stability of the peak tube voltage over different scan requests was about 0.10 kV for the low and 0.46 for the high energy kV-switching scans and less than 0.1 kV for kV-constant scans, indicating slightly lower stability for kV-switching scans. The tube current resulted in a relative deviation of -1.6% for the low and 6.6% overestimation for the high tube voltage of the kV-switching scans compared to the kV-constant scans. The pulse width showed no deviation for the longer pulse width and only minor deviations (0.02 ± 0.02 ms) for the shorter pulse widths compared to the kV-constant scans. The phantom experiment using different iodine concentrations showed an accurate correlation (R² > 0.99) between the extracted intensity values in the kV-switching and kV-constant reconstructed volumes, and allows for an automatic differentiation between contrast concentration down to 10% (350 mg/ml iodine) and pure water under low-noise conditions. Preliminary results of iodine and soft tissue separation showed also promising results in the first in vivo pig study.

CONCLUSIONS

The feasibility of dual-energy imaging using a fast kV-switching method on an angiographic C-arm CT system was investigated. Direct measurements of beam quality in the x-ray field demonstrate the stability of the kV-switching method. Phantom and in vivo experiments showed that images did not deviate from those of corresponding kV-constant scans. All performed experiments confirmed the capability of performing fast kV-switching scans on a clinically available C-arm CT system. More complex material decomposition tasks and postprocessing steps will be part of future investigations.

Performance of a clinical gridded electron gun in magnetic fields: Implications for MRI-linac therapy

PURPOSE

MRI-linac therapy is a rapidly growing field, and requires that conventional linear accelerators are operated with the fringe field of MRI magnets. One of the most sensitive accelerator components is the electron gun, which serves as the source of the beam. The purpose of this work was to develop a validated finite element model (FEM) model of a clinical triode (or gridded) electron gun, based on accurate geometric and electrical measurements, and to characterize the performance of this gun in magnetic fields.

METHODS

The geometry of a Varian electron gun was measured using 3D laser scanning and digital calipers. The electric potentials and emission current of these guns were measured directly from six dose matched true beam linacs for the 6X, 10X, and 15X modes of operation. Based on these measurements, a finite element model (FEM) of the gun was developed using the commercial software opera/scala. The performance of the FEM model in magnetic fields was characterized using parallel fields ranging from 0 to 200 G in the in-line direction, and 0-35 G in the perpendicular direction.

RESULTS

The FEM model matched the average measured emission current to within 5% across all three modes of operation. Different high voltage settings are used for the different modes; the 6X, 10X, and 15X modes have an average high voltage setting of 15, 10, and 11 kV. Due to these differences, different operating modes show different sensitivities in magnetic fields. For in line fields, the first current loss occurs at 40, 20, and 30 G for each mode. This is a much greater sensitivity than has previously been observed. For perpendicular fields, first beam loss occurred at 8, 5, and 5 G and total beam loss at 27, 22, and 20 G.

CONCLUSIONS

A validated FEM model of a clinical triode electron gun has been developed based on accurate geometric and electrical measurements. Three different operating modes were simulated, with a maximum mean error of 5%. This gun shows greater sensitivity to in-line magnetic fields than previously presented models, and different operating modes show different sensitivity.

Virtual scatter modulation for X-ray CT scatter correction using primary modulator

A new scatter estimation algorithm with a concept of virtual scatter modulation for X-ray scatter correction using primary modulator is proposed to reduce the aliasing errors in the estimated scatter. Virtual scatter modulation can be realized through dividing the measured primary-modulated image by the measured modulation function. After the division, the aggravation of the aliasing of primary due to the non-uniformity of the modulation function is largely transferred to that of scatter. Since scatter in general has less high frequencies than primary does, the aggravation of its aliasing is expected to be weaker, and therefore the overall aliasing can be reduced. A CatPhan©600 phantom and an anthropomorphic thorax phantom are scanned on a tabletop X-ray cone-beam computed tomography system to validate our proposed algorithm. On the Catphan phantom, the oscillations that are clearly observed in the central region of the Catphan scatter profile estimated using the original primary-modulation algorithm, are mostly eliminated with the proposed scatter modulation algorithm, leading to less residual artifacts and better CT number uniformity in the reconstructed image. Compared with 38.9 HU of CT nonuniformity in a selected uniform region when the primary-modulation algorithm is used, the new algorithm significantly reduces it to 4.5 HU, reaching the same level of uniformity as the ground truth reference. On the thorax phantom, overall better CT number uniformity is also achieved.

Marker-free motion correction in weight-bearing cone-beam CT of the knee joint

PURPOSE

To allow for a purely image-based motion estimation and compensation in weight-bearing cone-beam computed tomography of the knee joint.

METHODS

Weight-bearing imaging of the knee joint in a standing position poses additional requirements for the image reconstruction algorithm. In contrast to supine scans, patient motion needs to be estimated and compensated. The authors propose a method that is based on 2D/3D registration of left and right femur and tibia segmented from a prior, motion-free reconstruction acquired in supine position. Each segmented bone is first roughly aligned to the motion-corrupted reconstruction of a scan in standing or squatting position. Subsequently, a rigid 2D/3D registration is performed for each bone to each of K projection images, estimating 6 × 4 × K motion parameters. The motion of individual bones is combined into global motion fields using thin-plate-spline extrapolation. These can be incorporated into a motion-compensated reconstruction in the backprojection step. The authors performed visual and quantitative comparisons between a state-of-the-art marker-based (MB) method and two variants of the proposed method using gradient correlation (GC) and normalized gradient information (NGI) as similarity measure for the 2D/3D registration.

RESULTS

The authors evaluated their method on four acquisitions under different squatting positions of the same patient. All methods showed substantial improvement in image quality compared to the uncorrected reconstructions. Compared to NGI and MB, the GC method showed increased streaking artifacts due to misregistrations in lateral projection images. NGI and MB showed comparable image quality at the bone regions. Because the markers are attached to the skin, the MB method performed better at the surface of the legs where the authors observed slight streaking of the NGI and GC methods. For a quantitative evaluation, the authors computed the universal quality index (UQI) for all bone regions with respect to the motion-free reconstruction. The authors quantitative evaluation over regions around the bones yielded a mean UQI of 18.4 for no correction, 53.3 and 56.1 for the proposed method using GC and NGI, respectively, and 53.7 for the MB reference approach. In contrast to the authors registration-based corrections, the MB reference method caused slight nonrigid deformations at bone outlines when compared to a motion-free reference scan.

CONCLUSIONS

The authors showed that their method based on the NGI similarity measure yields reconstruction quality close to the MB reference method. In contrast to the MB method, the proposed method does not require any preparation prior to the examination which will improve the clinical workflow and patient comfort. Further, the authors found that the MB method causes small, nonrigid deformations at the bone outline which indicates that markers may not accurately reflect the internal motion close to the knee joint. Therefore, the authors believe that the proposed method is a promising alternative to MB motion management.

A novel electron accelerator for MRI-Linac radiotherapy

PURPOSE

MRI guided radiotherapy is a rapidly growing field; however, current electron accelerators are not designed to operate in the magnetic fringe fields of MRI scanners. As such, current MRI-Linac systems require magnetic shielding, which can degrade MR image quality and limit system flexibility. The purpose of this work was to develop and test a novel medical electron accelerator concept which is inherently robust to operation within magnetic fields for in-line MRI-Linac systems.

METHODS

Computational simulations were utilized to model the accelerator, including the thermionic emission process, the electromagnetic fields within the accelerating structure, and resulting particle trajectories through these fields. The spatial and energy characteristics of the electron beam were quantified at the accelerator target and compared to published data for conventional accelerators. The model was then coupled to the fields from a simulated 1 T superconducting magnet and solved for cathode to isocenter distances between 1.0 and 2.4 m; the impact on the electron beam was quantified.

RESULTS

For the zero field solution, the average current at the target was 146.3 mA, with a median energy of 5.8 MeV (interquartile spread of 0.1 MeV), and a spot size diameter of 1.5 mm full-width-tenth-maximum. Such an electron beam is suitable for therapy, comparing favorably to published data for conventional systems. The simulated accelerator showed increased robustness to operation in in-line magnetic fields, with a maximum current loss of 3% compared to 85% for a conventional system in the same magnetic fields.

CONCLUSIONS

Computational simulations suggest that replacing conventional DC electron sources with a RF based source could be used to develop medical electron accelerators which are robust to operation in in-line magnetic fields. This would enable the development of MRI-Linac systems with no magnetic shielding around the Linac and reduce the requirements for optimization of magnetic fringe field, simplify design of the high-field magnet, and increase system flexibility.

Reproducibility of Parenchymal Blood Volume Measurements Using an Angiographic C-arm CT System

RATIONALE AND OBJECTIVES

Intra-procedural measurement of hepatic perfusion following liver embolization continues to be a challenge. Blood volume imaging before and after interventional procedures would allow identifying the treatment end point or even allow predicting treatment outcome. Recent liver oncology studies showed the feasibility of parenchymal blood volume (PBV) imaging using an angiographic C-arm system. This study was done to evaluate the reproducibility of PBV measurements using cone beam computed tomography (CBCT) before and after embolization of the liver in a swine model.

MATERIALS AND METHODS

CBCT imaging was performed before and after partial bland embolization of the left lobe of the liver in five adult pigs. Intra-arterial injection of iodinated contrast with a 6-second x-ray delay was used with a two-sweep 8-second rotation imaging protocol. Three acquisitions, each separated by 10 minutes to allow for contrast clearance, were obtained before and after embolization in each animal. Post-processing was carried out using dedicated software to generate three-dimensional (3D) PBV maps. Two region-of-interest measurements were placed on two views within the right and left lobe on each CBCT 3D PBV map. Variation in PBV for scans acquired within each animal was determined by the coefficient of variation and intraclass correlation. A Wilcoxon signed-rank test was used to test post-procedure reduction in PBV.

RESULTS

The CBCT PBV maps showed mean coefficients of variation of 7% (range: 2%-16%) and 25% (range:  13%-34%) for baseline and embolized PBV maps, respectively. The intraclass correlation for PBV measurements was 0.89, demonstrating high reproducibility, with measurable reduction in PBV displayed after embolization (P = 0.007).

CONCLUSIONS

Intra-procedural acquisition of 3D PBV maps before and after liver embolization using CBCT is highly reproducible and shows promising application for obtaining intra-procedural PBV maps during locoregional therapy.

Coping with real world data: Artifact reduction and denoising for motion-compensated cardiac C-arm CT

PURPOSE

Detailed analysis of cardiac motion would be helpful for supporting clinical workflow in the interventional suite. With an angiographic C-arm system, multiple heart phases can be reconstructed using electrocardiogram gating. However, the resulting angular undersampling is highly detrimental to the quality of the reconstructed images, especially in nonideal intraprocedural imaging conditions. Motion-compensated reconstruction has previously been shown to alleviate this problem, but it heavily relies on a preliminary reconstruction suitable for motion estimation. In this work, the authors propose a processing pipeline tailored to augment these initial images for the purpose of motion estimation and assess how it affects the final images after motion compensation.

METHODS

The following combination of simple, direct methods inspired by the core ideas of existing approaches proved beneficial: (a) Streak reduction by masking high-intensity components in projection domain after filtering. (b) Streak reduction by subtraction of estimated artifact volumes in reconstruction domain. (c) Denoising in spatial domain using a joint bilateral filter guided by an uncompensated reconstruction. (d) Denoising in temporal domain using an adaptive Gaussian smoothing based on a novel motion detection scheme.

RESULTS

Experiments on a numerical heart phantom yield a reduction of the relative root-mean-square error from 89.9% to 3.6% and an increase of correlation with the ground truth from 95.763% to 99.995% for the motion-compensated reconstruction when the authors' processing is applied to the initial images. In three clinical patient data sets, the signal-to-noise ratio measured in an ideally homogeneous region is increased by 37.7% on average. Overall visual appearance is improved notably and some anatomical features are more readily discernible.

CONCLUSIONS

The authors' findings suggest that the proposed sequence of steps provides a clear advantage over an arbitrary sequence of individual image enhancement methods and is fit to overcome the issue of lacking image quality in motion-compensated C-arm imaging of the heart. As for future work, the obtained results pave the way for investigating how accurately cardiac functional motion parameters can be determined with this modality.

Optimization of an on-board imaging system for extremely rapid radiation therapy

PURPOSE

Next-generation extremely rapid radiation therapy systems could mitigate the need for motion management, improve patient comfort during the treatment, and increase patient throughput for cost effectiveness. Such systems require an on-board imaging system that is competitively priced, fast, and of sufficiently high quality to allow good registration between the image taken on the day of treatment and the image taken the day of treatment planning. In this study, three different detectors for a custom on-board CT system were investigated to select the best design for integration with an extremely rapid radiation therapy system.

METHODS

Three different CT detectors are proposed: low-resolution (all 4×4 mm pixels), medium-resolution (a combination of 4×4 mm pixels and 2×2 mm pixels), and high-resolution (all 1×1 mm pixels). An in-house program was used to generate projection images of a numerical anthropomorphic phantom and to reconstruct the projections into CT datasets, henceforth called "realistic" images. Scatter was calculated using a separate Monte Carlo simulation, and the model included an antiscatter grid and bowtie filter. Diagnostic-quality images of the phantom were generated to represent the patient scan at the time of treatment planning. Commercial deformable registration software was used to register the diagnostic-quality scan to images produced by the various on-board detector configurations. The deformation fields were compared against a "gold standard" deformation field generated by registering initial and deformed images of the numerical phantoms that were used to make the diagnostic and treatment-day images. Registrations of on-board imaging system data were judged by the amount their deformation fields differed from the corresponding gold standard deformation fields--the smaller the difference, the better the system. To evaluate the registrations, the pointwise distance between gold standard and realistic registration deformation fields was computed.

RESULTS

By most global metrics (e.g., mean, median, and maximum pointwise distance), the high-resolution detector had the best performance but the medium-resolution detector was comparable. For all medium- and high-resolution detector registrations, mean error between the realistic and gold standard deformation fields was less than 4 mm. By pointwise metrics (e.g., tracking a small lesion), the high- and medium-resolution detectors performed similarly. For these detectors, the smallest error between the realistic and gold standard registrations was 0.6 mm and the largest error was 3.6 mm.

CONCLUSIONS

The medium-resolution CT detector was selected as the best for an extremely rapid radiation therapy system. In essentially all test cases, data from this detector produced a significantly better registration than data from the low-resolution detector and a comparable registration to data from the high-resolution detector. The medium-resolution detector provides an appropriate compromise between registration accuracy and system cost.

Practical dose point-based methods to characterize dose distribution in a stationary elliptical body phantom for a cone-beam C-arm CT system

PURPOSE

To propose new dose point measurement-based metrics to characterize the dose distributions and the mean dose from a single partial rotation of an automatic exposure control-enabled, C-arm-based, wide cone angle computed tomography system over a stationary, large, body-shaped phantom.

METHODS

A small 0.6 cm(3) ion chamber (IC) was used to measure the radiation dose in an elliptical body-shaped phantom made of tissue-equivalent material. The IC was placed at 23 well-distributed holes in the central and peripheral regions of the phantom and dose was recorded for six acquisition protocols with different combinations of minimum kVp (109 and 125 kVp) and z-collimator aperture (full: 22.2 cm; medium: 14.0 cm; small: 8.4 cm). Monte Carlo (MC) simulations were carried out to generate complete 2D dose distributions in the central plane (z = 0). The MC model was validated at the 23 dose points against IC experimental data. The planar dose distributions were then estimated using subsets of the point dose measurements using two proposed methods: (1) the proximity-based weighting method (method 1) and (2) the dose point surface fitting method (method 2). Twenty-eight different dose point distributions with six different point number cases (4, 5, 6, 7, 14, and 23 dose points) were evaluated to determine the optimal number of dose points and their placement in the phantom. The performances of the methods were determined by comparing their results with those of the validated MC simulations. The performances of the methods in the presence of measurement uncertainties were evaluated.

RESULTS

The 5-, 6-, and 7-point cases had differences below 2%, ranging from 1.0% to 1.7% for both methods, which is a performance comparable to that of the methods with a relatively large number of points, i.e., the 14- and 23-point cases. However, with the 4-point case, the performances of the two methods decreased sharply. Among the 4-, 5-, 6-, and 7-point cases, the 7-point case (1.0% [±0.6%] difference) and the 6-point case (0.7% [±0.6%] difference) performed best for method 1 and method 2, respectively. Moreover, method 2 demonstrated high-fidelity surface reconstruction with as few as 5 points, showing pixelwise absolute differences of 3.80 mGy (±0.32 mGy). Although the performance was shown to be sensitive to the phantom displacement from the isocenter, the performance changed by less than 2% for shifts up to 2 cm in the x- and y-axes in the central phantom plane.

CONCLUSIONS

With as few as five points, method 1 and method 2 were able to compute the mean dose with reasonable accuracy, demonstrating differences of 1.7% (±1.2%) and 1.3% (±1.0%), respectively. A larger number of points do not necessarily guarantee better performance of the methods; optimal choice of point placement is necessary. The performance of the methods is sensitive to the alignment of the center of the body phantom relative to the isocenter. In body applications where dose distributions are important, method 2 is a better choice than method 1, as it reconstructs the dose surface with high fidelity, using as few as five points.

Experimental validation of L-shell x-ray fluorescence computed tomography imaging: phantom study

Thanks to the current advances in nanoscience, molecular biochemistry, and x-ray detector technology, x-ray fluorescence computed tomography (XFCT) has been considered for molecular imaging of probes containing high atomic number elements, such as gold nanoparticles. The commonly used XFCT imaging performed with K-shell x rays appears to have insufficient imaging sensitivity to detect the low gold concentrations observed in small animal studies. Low energy fluorescence L-shell x rays have exhibited higher signal-to-background ratio and appeared as a promising XFCT mode with greatly enhanced sensitivity. The aim of this work was to experimentally demonstrate the feasibility of L-shell XFCT imaging and to assess its achievable sensitivity. We built an experimental L-shell XFCT imaging system consisting of a miniature x-ray tube and two spectrometers, a silicon drift detector (SDD), and a CdTe detector placed at [Formula: see text] with respect to the excitation beam. We imaged a 28-mm-diameter water phantom with 4-mm-diameter Eppendorf tubes containing gold solutions with concentrations of 0.06 to 0.1% Au. While all Au vials were detectable in the SDD L-shell XFCT image, none of the vials were visible in the CdTe L-shell XFCT image. The detectability limit of the presented L-shell XFCT SDD imaging setup was 0.007% Au, a concentration observed in small animal studies.

Proton-induced x-ray fluorescence CT imaging

PURPOSE

To demonstrate the feasibility of proton-induced x-ray fluorescence CT (pXFCT) imaging of gold in a small animal sized object by means of experiments and Monte Carlo (MC) simulations.

METHODS

First, proton-induced gold x-ray fluorescence (pXRF) was measured as a function of gold concentration. Vials of 2.2 cm in diameter filled with 0%-5% Au solutions were irradiated with a 220 MeV proton beam and x-ray fluorescence induced by the interaction of protons, and Au was detected with a 3 × 3 mm(2) CdTe detector placed at 90° with respect to the incident proton beam at a distance of 45 cm from the vials. Second, a 7-cm diameter water phantom containing three 2.2-diameter vials with 3%-5% Au solutions was imaged with a 7-mm FWHM 220 MeV proton beam in a first generation CT scanning geometry. X-rays scattered perpendicular to the incident proton beam were acquired with the CdTe detector placed at 45 cm from the phantom positioned on a translation/rotation stage. Twenty one translational steps spaced by 3 mm at each of 36 projection angles spaced by 10° were acquired, and pXFCT images of the phantom were reconstructed with filtered back projection. A simplified geometry of the experimental data acquisition setup was modeled with the MC TOPAS code, and simulation results were compared to the experimental data.

RESULTS

A linear relationship between gold pXRF and gold concentration was observed in both experimental and MC simulation data (R(2) > 0.99). All Au vials were apparent in the experimental and simulated pXFCT images. Specifically, the 3% Au vial was detectable in the experimental [contrast-to-noise ratio (CNR) = 5.8] and simulated (CNR = 11.5) pXFCT image. Due to fluorescence x-ray attenuation in the higher concentration vials, the 4% and 5% Au contrast were underestimated by 10% and 15%, respectively, in both the experimental and simulated pXFCT images.

CONCLUSIONS

Proton-induced x-ray fluorescence CT imaging of 3%-5% gold solutions in a small animal sized water phantom has been demonstrated for the first time by means of experiments and MC simulations.

Robust primary modulation-based scatter estimation for cone-beam CT

PURPOSE

Scattered radiation is one of the major problems facing image quality in flat detector cone-beam computed tomography (CBCT). Previously, a new scatter estimation and correction method using primary beam modulation has been proposed. The original image processing technique used a frequency-domain-based analysis, which proved to be sensitive to the accuracy of the modulator pattern both spatially and in amplitude as well as to the frequency of the modulation pattern. In addition, it cannot account for penumbra effects that occur, for example, due to the finite focal spot size and the scatter estimate can be degraded by high-frequency components of the primary image.

METHODS

In this paper, the authors present a new way to estimate the scatter using primary modulation. It is less sensitive to modulator nonidealities and most importantly can handle arbitrary modulator shapes and changes in modulator attenuation. The main idea is that the scatter estimation can be expressed as an optimization problem, which yields a separation of the scatter and the primary image. The method is evaluated using simulated and experimental CBCT data. The scattering properties of the modulator itself are analyzed using a Monte Carlo simulation.

RESULTS

All reconstructions show strong improvements of image quality. To quantify the results, all images are compared to reference images (ideal simulations and collimated scans).

CONCLUSIONS

The proposed modulator-based scatter reduction algorithm may open the field of flat detector-based imaging to become a quantitative modality. This may have significant impact on C-arm imaging and on image-guided radiation therapy.