Effect of detruncation on the accuracy of Monte Carlo-based scatter estimation in truncated CBCT.

PURPOSE: The purpose of this study is to investigate the necessity of detruncation for scatter estimation of truncated cone-beam CT (CBCT) data and to evaluate different detruncation algorithms. Scattered radiation results in some of the most severe artifacts in CT and depends strongly on the size and the shape of the scanned object. Especially in CBCT systems the large cone-angle and the small detector-to-isocenter distance lead to a large amount of scatter detected, resulting in cupping artifacts, streak artifacts, and inaccurate CT-values. If a small field of measurement (FOM) is used, as it is often the case in CBCT systems, data are truncated in longitudinal and lateral direction. Since only truncated data are available as input for the scatter estimation, the already challenging correction of scatter artifacts becomes even more difficult. METHODS: The following detruncation methods are compared and evaluated with respect to scatter estimation: constant detruncation, cosine detruncation, adaptive detruncation, and prior-based detruncation using anatomical data from a similar phantom or patient, also compared to the case where no detruncation was performed. Each of the resulting, detruncated reconstructions serve as input volume for a Monte Carlo (MC) scatter estimation and subsequent scatter correction. An evaluation is performed on a head simulation, measurements of a head phantom and a patient using a dental CBCT geometry with a FOM diameter of 11 cm. Additionally, a thorax phantom is measured to assess performance in a C-Arm geometry with a FOM of up to 20 cm. RESULTS: If scatter estimation is based on simple detruncation algorithms like a constant or a cosine detruncation scatter is estimated inaccurately, resulting in incorrect CT-values as well as streak artifacts in the corrected volume. For the dental CBCT phantom measurement CT-values for soft tissue were corrected from -204 HU (no scatter correction) to -87 HU (no detruncation), -218 HU (constant detruncation), -141 HU (cosine detruncation), -91 HU (adaptive detruncation), -34 HU (prior-based detruncation using a different prior) and -24 HU (prior-based detruncation using the identical prior) for a reference value of -26 HU measured in slit scan mode. In all cases the prior-based detruncation results in the best scatter correction, followed by the adaptive detruncation, as these algorithms provide a rather accurate model of high-density structures outside the FOM, compared to a simple constant or a cosine detruncation. CONCLUSIONS: Our contribution is twofold: first we give a comprehensive comparison of various detruncation methods for the purpose of scatter estimation. We find that the choice of the detruncation method has a significant influence on the quality of MC-based scatter correction. Simple or no detruncation is often insufficient for artifact removal and results in inaccurate CT-values. On the contrary, prior-based detruncation can achieve a high CT-value accuracy and nearly artifact-free volumes from truncated CBCT data when combined with other state-of-the-art artifact corrections. Secondly, we show that prior-based detruncation is effective even with data from a different patient or phantom. The fact that data completion does not require data from the same patient dramatically increases the applicability and usability of this scatter estimation.

Segmentation-free empirical beam hardening correction for CT.

PURPOSE: The polychromatic nature of the x-ray beams and their effects on the reconstructed image are often disregarded during standard image reconstruction. This leads to cupping and beam hardening artifacts inside the reconstructed volume. To correct for a general cupping, methods like water precorrection exist. They correct the hardening of the spectrum during the penetration of the measured object only for the major tissue class. In contrast, more complex artifacts like streaks between dense objects need other techniques of correction. If using only the information of one single energy scan, there are two types of corrections. The first one is a physical approach. Thereby, artifacts can be reproduced and corrected within the original reconstruction by using assumptions in a polychromatic forward projector. These assumptions could be the used spectrum, the detector response, the physical attenuation and scatter properties of the intersected materials. A second method is an empirical approach, which does not rely on much prior knowledge. This so-called empirical beam hardening correction (EBHC) and the previously mentioned physical-based technique are both relying on a segmentation of the present tissues inside the patient. The difficulty thereby is that beam hardening by itself, scatter, and other effects, which diminish the image quality also disturb the correct tissue classification and thereby reduce the accuracy of the two known classes of correction techniques. The herein proposed method works similar to the empirical beam hardening correction but does not require a tissue segmentation and therefore shows improvements on image data, which are highly degraded by noise and artifacts. Furthermore, the new algorithm is designed in a way that no additional calibration or parameter fitting is needed. METHODS: To overcome the segmentation of tissues, the authors propose a histogram deformation of their primary reconstructed CT image. This step is essential for the proposed algorithm to be segmentation-free (sf). This deformation leads to a nonlinear accentuation of higher CT-values. The original volume and the gray value deformed volume are monochromatically forward projected. The two projection sets are then monomially combined and reconstructed to generate sets of basis volumes which are used for correction. This is done by maximization of the image flatness due to adding additionally a weighted sum of these basis images. sfEBHC is evaluated on polychromatic simulations, phantom measurements, and patient data. The raw data sets were acquired by a dual source spiral CT scanner, a digital volume tomograph, and a dual source micro CT. Different phantom and patient data were used to illustrate the performance and wide range of usability of sfEBHC across different scanning scenarios. The artifact correction capabilities are compared to EBHC. RESULTS: All investigated cases show equal or improved image quality compared to the standard EBHC approach. The artifact correction is capable of correcting beam hardening artifacts for different scan parameters and scan scenarios. CONCLUSIONS: sfEBHC generates beam hardening-reduced images and is furthermore capable of dealing with images which are affected by high noise and strong artifacts. The algorithm can be used to recover structures which are hardly visible inside the beam hardening-affected regions.

Nonlinear quantum optomechanics via individual intrinsic two-level defects.

We propose to use the intrinsic two-level system (TLS) defect states found naturally in integrated optomechanical devices for exploring cavity QED-like phenomena with localized phonons. The Jaynes-Cummings-type interaction between TLS and mechanics can reach the strong coupling regime for existing nano-optomechanical systems, observable via clear signatures in the optomechanical output spectrum. These signatures persist even at finite temperature, and we derive an explicit expression for the temperature at which they vanish. Further, the ability to drive the defect with a microwave field allows for realization of phonon blockade, and the available controls are sufficient to deterministically prepare non-classical states of the mechanical resonator.

Electromagnetic interaction of split-ring resonators: The role of separation and relative orientation.

Extinction cross-section spectra of split-ring-resonator dimers have been measured at near-infrared frequencies with a sensitive spatial modulation technique. The resonance frequency of the dimer's coupled mode as well as its extinction cross-section and its quality factor depend on the relative orientation and separation of the two split-ring resonators. The findings can be interpreted in terms of electric and magnetic dipole-dipole interaction. Numerical calculations based on a Discontinuous Galerkin Time-Domain approach are in good agreement with the experiments and support our physical interpretation.

Discontinuous Galerkin time-domain computations of metallic nanostructures.

We apply the three-dimensional Discontinuous-Galerkin Time-Domain method to the investigation of the optical properties of bar- and V-shaped metallic nanostructures on dielectric substrates. A flexible finite element-like mesh together with an expansion into high-order basis functions allows for an accurate resolution of complex geometries and strong field gradients. In turn, this provides accurate results on the optical response of realistic structures. We study in detail the influence of particle size and shape on resonance frequencies as well as on scattering and absorption efficiencies. Beyond a critical size which determines the onset of the quasi-static limit we find significant deviations from the quasi-static theory. Furthermore, we investigate the influence of the excitation by comparing normal illumination and attenuated total internal reflection setups. Finally, we examine the possibility of coherently controlling the local field enhancement of V-structures via chirped pulses.