Performance studies of four-dimensional cone beam computed tomography

Four-dimensional cone beam computed tomography (4DCBCT) has been proposed to characterize the breathing motion of tumors before radiotherapy treatment. However, when the acquired cone beam projection data are retrospectively gated into several respiratory phases, the available data to reconstruct each phase is under-sampled and thus causes streaking artifacts in the reconstructed images. To solve the under-sampling problem and improve image quality in 4DCBCT, various methods have been developed. This paper presents performance studies of three different 4DCBCT methods based on different reconstruction algorithms. The aims of this paper are to study (1) the relationship between the accuracy of the extracted motion trajectories and the data acquisition time of a 4DCBCT scan and (2) the relationship between the accuracy of the extracted motion trajectories and the number of phase bins used to sort projection data. These aims will be applied to three different 4DCBCT methods: conventional filtered backprojection reconstruction (FBP), FBP with McKinnon-Bates correction (MB) and prior image constrained compressed sensing (PICCS) reconstruction. A hybrid phantom consisting of realistic chest anatomy and a moving elliptical object with known 3D motion trajectories was constructed by superimposing the analytical projection data of the moving object to the simulated projection data from a chest CT volume dataset. CBCT scans with gantry rotation times from 1 to 4 min were simulated, and the generated projection data were sorted into 5, 10 and 20 phase bins before different methods were used to reconstruct 4D images. The motion trajectories of the moving object were extracted using a fast free-form deformable registration algorithm. The root mean square errors (RMSE) of the extracted motion trajectories were evaluated for all simulated cases to quantitatively study the performance. The results demonstrate (1) longer acquisition times result in more accurate motion delineation for each method; (2) ten or more phase bins are necessary in 4DCBCT to ensure sufficient temporal resolution in tumor motion and (3) to achieve the same performance as FBP-4DCBCT with a 4 min data acquisition time, MB-4DCBCT and PICCS-4DCBCT need about 2- and 1 min data acquisition times, respectively.     

A method to correct for spectral artifacts in optical-CT dosimetry

The recent emergence of radiochromic dosimeters with low inherent light-scattering presents the possibility of fast 3D dosimetry using broad-beam optical computed tomography (optical-CT). Current broad beam scanners typically employ either a single or a planar array of light-emitting diodes (LED) for the light source. The spectrum of light from LED sources is polychromatic and this, in combination with the non-uniform spectral absorption of the dosimeter, can introduce spectral artifacts arising from preferential absorption of photons at the peak absorption wavelengths in the dosimeter. Spectral artifacts can lead to large errors in the reconstructed attenuation coefficients, and hence dose measurement. This work presents an analytic method for correcting for spectral artifacts which can be applied if the spectral characteristics of the light source, absorbing dosimeter, and imaging detector are known or can be measured. The method is implemented here for a PRESAGE? dosimeter scanned with the DLOS telecentric scanner (Duke Large field-of-view Optical-CT Scanner). Emission and absorption profiles were measured with a commercial spectrometer and spectrophotometer, respectively. Simulations are presented that show spectral changes can introduce errors of 8% for moderately attenuating samples where spectral artifacts are less pronounced. The correction is evaluated by application to a 16 cm diameter PRESAGE? cylindrical dosimeter irradiated along the axis with two partially overlapping 6 × 6 cm fields of different doses. The resulting stepped dose distribution facilitates evaluation of the correction as each step had different spectral contributions. The spectral artifact correction was found to accurately correct the reconstructed coefficients to within ～1.5%, improved from ～7.5%, for normalized dose distributions. In conclusion, for situations where spectral artifacts cannot be removed by physical filters, the method shown here is an effective correction. Physical filters may be less viable if they introduce strong sensitivity to Schlieren bands in the dosimeters.     

Patient dose from kilovoltage cone beam computed tomography imaging in radiation therapy

Kilovoltage cone-beam computerized tomography (kV-CBCT) systems integrated into the gantry of linear accelerators can be used to acquire high-resolution volumetric images of the patient in the treatment position. Using on-line software and hardware, patient position can be determined accurately with a high degree of precision and, subsequently, set-up parameters can be adjusted to deliver the intended treatment. While the patient dose due to a single volumetric imaging acquisition is small compared to the therapy dose, repeated and daily image guidance procedures can lead to substantial dose to normal tissue. The dosimetric properties of a clinical CBCT system have been studied on an Elekta linear accelerator (Synergy RP, XVI system) and additional measurements performed on a laboratory system with identical geometry. Dose measurements were performed with an ion chamber and MOSFET detectors at the center, periphery, and surface of 30 and 16-cm-diam cylindrical shaped water phantoms, as a function of x-ray energy and longitudinal field-of-view (FOV) settings of 5,10,15, and 26 cm. The measurements were performed for full 360 degrees CBCT acquisition as well as for half-rotation scans for 120 kVp beams using the 30-cm-diam phantom. The dose at the center and surface of the body phantom were determined to be 1.6 and 2.3 cGy for a typical imaging protocol, using full rotation scan, with a technique setting of 120 kVp and 660 mAs. The results of our measurements have been presented in terms of a dose conversion factor fCBCT, expressed in cGy/R. These factors depend on beam quality and phantom size as well as on scan geometry and can be utilized to estimate dose for any arbitrary mAs setting and reference exposure rate of the x-ray tube at standard distance. The results demonstrate the opportunity to manipulate the scanning parameters to reduce the dose to the patient by employing lower energy (kVp) beams, smaller FOV, or by using half-rotation scan.     

牙种植中数字化曲面断层与锥体束CT应用的对比

背景：影像学检查可及时了解牙种植位置信息,但目前常用的数字化颌骨曲面断层对判定骨的高度和宽度存在误差。 目的：比较数字化曲面断层与锥体束CT在牙种植中的应用价值。 方法：选取要求行种植治疗修复缺失牙,牙周状况良好,符合种植术适应证的患者53例(65颗),种植前以一直径为3 mm的钢球为标记物置于患者待种植区,锥体束CT及普通数字化曲面断层分别测量三维重建图像及曲面断层片中标记物的直径,并与真实值进行比较,分析二者骨整合检出率差异。 结果与结论：相比于数字化曲面断层,口腔种植前锥体束CT测量的小钢球直径与真实值差异较小(P < 0.01)。种植后3个月,锥体束CT和普通数字化曲面断层得出的未发生骨结合的检出率接近(P > 0.05)。说明与普通数字化曲面断层相比,锥体束CT可以更准确的评估种植前牙槽骨的骨量,并可在种植后检查中更清晰的反映种植体周围骨质情况。     

Patient dose and image quality from mega-voltage cone beam computed tomography imaging

The evolution of ever more conformal radiation delivery techniques makes the subject of accurate localization of increasing importance in radiotherapy. Several systems can be utilized including kilo-voltage and mega-voltage cone-beam computed tomography (MV-CBCT), CT on rail or helical tomography. One of the attractive aspects of mega-voltage cone-beam CT is that it uses the therapy beam along with an electronic portal imaging device to image the patient prior to the delivery of treatment. However, the use of a photon beam energy in the mega-voltage range for volumetric imaging degrades the image quality and increases the patient radiation dose. To optimize image quality and patient dose in MV-CBCT imaging procedures, a series of dose measurements in cylindrical and anthropomorphic phantoms using an ionization chamber, radiographic films, and thermoluminescent dosimeters was performed. Furthermore, the dependence of the contrast to noise ratio and spatial resolution of the image upon the dose delivered for a 20-cm-diam cylindrical phantom was evaluated. Depending on the anatomical site and patient thickness, we found that the minimum dose deposited in the irradiated volume was 5-9 cGy and the maximum dose was between 9 and 17 cGy for our clinical MV-CBCT imaging protocols. Results also demonstrated that for high contrast areas such as bony anatomy, low doses are sufficient for image registration and visualization of the three-dimensional boundaries between soft tissue and bony structures. However, as the difference in tissue density decreased, the dose required to identify soft tissue boundaries increased. Finally, the dose delivered by MV-CBCT was simulated using a treatment planning system (TPS), thereby allowing the incorporation of MV-CBCT dose in the treatment planning process. The TPS-calculated doses agreed well with measurements for a wide range of imaging protocols.     

Exact filtered backprojection reconstruction for dynamic pitch helical cone beam computed tomography

We present an exact filtered backprojection reconstruction formula for helical cone beam computed tomography in which the pitch of the helix varies with time. We prove that the resulting algorithm, which is functionally identical to the constant pitch case, provides exact reconstruction provided that the projection of the helix onto the detector forms convex boundaries and that PI lines are unique. Furthermore, we demonstrate that both of these conditions are satisfied provided the sum of the translational velocity and the derivative of the translational acceleration does not change sign. As a special case, we show that gantry tilt can also be handled by our dynamic pitch formula. Simulation results demonstrate the resulting algorithm.     

Streaking artifacts reduction in four-dimensional cone-beam computed tomography

Cone-beam computed tomography (CBCT) using an "on-board" x-ray imaging device integrated into a radiation therapy system has recently been made available for patient positioning, target localization, and adaptive treatment planning. One of the challenges for gantry mounted image-guided radiation therapy (IGRT) systems is the slow acquisition of projections for cone-beam CT (CBCT), which makes them sensitive to any patient motion during the scans. Aiming at motion artifact reduction, four-dimensional CBCT (4D CBCT) techniques have been introduced, where a surrogate for the target's motion profile is utilized to sort the cone-beam data by respiratory phase. However, due to the limited gantry rotation speed and limited readout speed of the on-board imager, fewer than 100 projections are available for the image reconstruction at each respiratory phase. Thus, severe undersampling streaking artifacts plague 4D CBCT images. In this paper, the authors propose a simple scheme to significantly reduce the streaking artifacts. In this method, a prior image is first reconstructed using all available projections without gating, in which static structures are well reconstructed while moving objects are blurred. The undersampling streaking artifacts from static structures are estimated from this prior image volume and then can be removed from the phase images using gated reconstruction. The proposed method was validated using numerical simulations, experimental phantom data, and patient data. The fidelity of stationary and moving objects is maintained, while large gains in streak artifact reduction are observed. Using this technique one can reconstruct 4D CBCT datasets using no more projections than are acquired in a 60 s scan. At the same time, a temporal gating window as narrow as 100 ms was utilized. Compared to the conventional 4D CBCT reconstruction, streaking artifacts were reduced by 60% to 70%.     

A reconstruction method for equidistant fan beam differential phase contrast computed tomography

We report a reconstruction method, called a back-projection filtered (BPF) algorithm, for fan beam differential phase contrast computed tomography (DPC-CT) with equidistant geometrical configuration. This work comprises a numerical study of the algorithm and its experimental verification with a three-grating interferometer and an x-ray tube source. The numerical simulation and experimental results demonstrate that the proposed method can deal with several classes of truncated datasets. It could be of interest in future medical phase contrast imaging applications.     

Bilateral cleft lip and palate: A morphometric analysis of facial skeletal form using cone beam computed tomography

Bilateral cleft lip and palate (BCLP) is caused by a lack of merging of maxillary and nasal facial prominences during development and morphogenesis. BCLP is associated with congenital defects of the oronasal facial region that can impair ingestion, mastication, speech, and dentofacial development. Using cone beam computed tomography (CBCT) images, 7‐ to 18‐year old individuals born with BCLP (n = 15) and age‐ and sex‐matched controls (n = 15) were retrospectively assessed. Coordinate values of three‐dimensional facial skeletal anatomical landmarks (n = 32) were measured from each CBCT image. Data were evaluated using principal coordinates analysis (PCOORD) and Euclidean Distance Matrix Analysis (EDMA). PCOORD axes 1–3 explain approximately 45% of the morphological variation between samples, and specific patterns of morphological differences were associated with each axis. Approximately, 30% of facial skeletal measures significantly differ by confidence interval testing (α = 0.10) between samples. While significant form differences occur across the facial skeleton, strong patterns of differences are localized to the lateral and superioinferior aspects of the nasal aperture. In conclusion, the BCLP deformity significantly alters facial skeletal morphology of the midface and oronasal regions of the face, but morphological differences were also found in the upper facial skeleton and to a lesser extent, the lower facial skeleton. This pattern of strong differences in the oronasal region of the facial skeleton combined with differences across the rest of the facial complex underscores the idea that bones of the craniofacial skeleton are integrated. Clin. Anat. 28:584–592, 2015. © 2015 Wiley Periodicals, Inc.     

Correction for beam hardening in computed tomography

Corrections for beam-hardening artifacts in computed tomography can be made by using a model which assumes that water and bone mineral are the only constituents of tissue. With this model, a correction factor for the measured transmission values can be calculated such that the reconstructed attenuation coefficients have values corresponding to a monoenergetic source of known energy. Systematic errors in the uncorrected attenuation coefficients, which may be 5%, can be reduced to less than 1% if corrected transmission values are used.     

Techniques to alleviate the effects of view aliasing artifacts in computed tomography.

Due to practical limitations in data acquisition, 3-D computed tomography systems must attempt to provide rapid reconstructions of acceptable quality from a limited number of views. The use of convolution backprojection (CBP) for image reconstruction from an inadequate number of projections, results in view aliasing artifacts. In this paper we investigate different post-processing methods of alleviating the effects of view aliasing artifacts. Two distinct methods and their variants are considered. The first, termed the intermediate view reprojection (IVR) method, involves estimating a set of intermediate views by reprojection, followed by a reconstruction using the augmented set of views. The second, termed the error-correction (EC) method, incorporates a correction on the initial reconstruction based on the projection-domain error. Suitable modifications and variants of the above methods are indicated. Of the methods discussed, the IVR method is simple, tends to reduce the effects of artifacts with less susceptibility to secondary effects, and is applicable to region-of-interest reconstructions.     

Position of the greater palatine foramen: an anatomical study through cone beam computed tomography images

Purpose

The block anesthesia of the greater palatine foramen (GPF) is largely used in minor oral surgeries, periodontics and general dentistry. Furthermore, the area of the GPF serves as a donor of soft tissue graft. So, the aim of this study was to evaluate the position and characteristics of the GPF in Brazilian patients using cone beam computed tomography (CBCT) providing anatomical information for the greater palatine nerve block anesthesia and indicate site to collect palatal donor tissue.

Methods

Fifty CBCT exams of Brazilian patients with a mean age of 35.8 years (27 male/23 female) were evaluated. All patients had erupted first, second and third upper molars. A total of 100 GPF were evaluated bilaterally. The GPFs were assessed regarding position, diameter and distances to the midline maxillary suture (MMS) and to alveolar ridge (AR). Guidelines were drawn in the CBCT axial image depicting all molar interproximal surfaces, bilaterally. The guidelines were located between first, second and third molar and in the center of the second and third, performing five guidelines in each side. These guidelines and the molars were landmarks to assess the GPF anatomic position.

Results

From the 100 GPF analyzed, 92 were located in the third molar region (24 male/22 female). The 92 GPF were distributed as 47 in the left side and 45 in the right side. The average GPF diameter and the distance to both the AR and the MMS were 3.1 mm; 7.9 and 15.3 mm, respectively.

Conclusions

Within the limits of this study, we concluded that the in Brazilian patients studied, the GPF location was more closely related to third molar. Therefore, whenever the third molar is erupted, it could be used as landmark for successful GPN block anesthesia. Moreover, harvesting palatal mucosa graft around the third molar should be done cautiously to prevent damage to the GPF vascular-nerve complex.     

Anatomical study of the human discomallear ligament using cone beam computed tomography imaging and morphological observations

In the present study, the human discomallear ligament (DML) was observed in structures at both macroscopic and cone beam computed tomography levels. Assessments were made regarding the distribution of calcitonin-gene-related peptide (CGRP), protein gene-product (PGP) 9.5, and substance P (SP) of the DML based on immunohistochemical analyses of the anatomical properties of jaw movements using 27 Japanese human cadavers (mean, 79.3 +/- 8.6 years; male, 74.9 +/- 8.0; female, 82.8 +/- 7.5). The DML of the anterior region was connected to the TMJ disc. The DML of the posterior region was attached to both the head and the anterior process of the malleus through the petrotympanic fissure, which formed a narrow channel. The structure of the petrotympanic fissure through the DML was attached to the malleus, and this structure was associated with the mobility of the malleus. In the anterior and posterior parts of the disc-associated connective tissue of the DML, CGRP-, PGP9.5- and SP-positive nerve fibers were located around numerous blood vessels, a condition which may be correlated with chronic pain syndrines disorders and the auditory system.