Cone-beam CT reconstruction for non-periodic organ motion using time-ordered chain graph model
© The Author(s) 2017
Received: 22 February 2017
Accepted: 23 August 2017
Published: 4 September 2017
The purpose of this study is to introduce the new concept of a four-dimensional (4D) cone-beam computed tomography (CBCT) reconstruction approach for non-periodic organ motion in cooperation with the time-ordered chain graph model (TCGM) and to compare it with previously developed methods such as total variation-based compressed sensing (TVCS) and prior-image constrained compressed sensing (PICCS).
Materials and Methods
Our proposed reconstruction is based on a model including the constraint originating from the images of neighboring time phases. Namely, the reconstructed time-series images depend on each other in this TCGM scheme, and the time-ordered images are concurrently reconstructed in the iterative reconstruction approach. In this study, iterative reconstruction with the TCGM was carried out with 90° projection ranges. The images reconstructed by the TCGM were compared with the images reconstructed by TVCS (200° projection ranges) and PICCS (90° projection ranges).
Two kinds of projection data sets–an elliptic-cylindrical digital phantom and two clinical patients’ data–were used. For the digital phantom, an air sphere was contained and virtually moved along the longitudinal axis by 3 cm/30 s and 3 cm/60 s; the temporal resolution was evaluated by measuring the penumbral width of the air sphere. The clinical feasibility of the non-periodic time-ordered 4D CBCT image reconstruction was examined with the patient data in the pelvic region.
In the evaluation of the digital-phantom reconstruction, the penumbral widths of the TCGM yielded the narrowest result; the results obtained by PICCS and TCGM using 90° projection ranges were 2.8% and 18.2% for 3 cm/30 s, and 5.0% and 23.1% for 3 cm/60 s narrower than that of TVCS using 200° projection ranges. This suggests that the TCGM has a better temporal resolution, whereas PICCS seems similar to TVCS. These reconstruction methods were also compared using patients’ projection data sets. Although all three reconstruction results showed motion related to rectal gas or stool, the result obtained by the TCGM was visibly clearer with less blurring.
The TCGM is a feasible approach to visualize non-periodic organ motion. The digital-phantom results demonstrated that the proposed method provides 4D image series with a better temporal resolution compared to TVCS and PICCS. The clinical patients’ results also showed that the present method enables us to visualize motion related to rectal gas and flatus in the rectum.
The volumetric imaging technique using a cone-beam computed tomography (CBCT) device mounted on the gantry of a radiotherapy linear accelerator is crucial in both image-guided radiotherapy (IGRT) and adaptive radiotherapy (ART). Since first introduced by Jaffray et al. with a volumetric reconstruction algorithm widely known as the Feldkamp-Davis-Kress (FDK) method proposed by Feldkamp et al., [1, 2] it has enabled the highly accurate positioning of radiotherapy patients [3–5]. Recently, CBCT images have been expected to provide image sets as the material of a retrospective analysis for margin evaluation, anatomical deformation, and dose distribution [6–8]. Online four-dimensional (4D) CBCT reconstruction has already been used in the treatment of tumor sites with periodic motion, such as a lung tumor, and it has also provided an accurate setup that visually considers the intrafractional motion of tumors and risk organs [9–12]. Especially, in the scenario of adaptive radiotherapy, 4D CBCT and in-treatment CBCT imaging technique has been expected to enable understanding 4D dose distribution with better temporal resolution .
For the purpose of 4D CBCT image reconstruction accompanied by periodic organ motion such as respiration and cardiac activity, information related to the respiratory or cardiac motion signal is necessary to classify projections into several phase groups. In the case of respiratory motion, the measurement of a signal synchronized with projection images is categorized by two methods–the use of an external respiratory monitoring system and image-based respiratory phase recognition.
As an example of the latter case, Zijp et al. proposed the method known as the Amsterdam Shroud, which produces a 1D projection in the craniocaudal (CC) direction after applying a CC derivative filter,[14, 15] and this has been used in a commercial linac-integrated CBCT device, the X-ray Volumetric Imaging (XVI) System mounted on Elekta Synergy (Elekta, Crawley, UK). Regarding cardiac activity, Lauzier et al. demonstrated 4D image reconstruction using electrocardiogram (ECG) signals in order to sort projections into several cardiac phases .
In contrast to the various studies on volumetric 4D CBCT imaging with periodic motion, few studies have dealt with volumetric imaging with non-periodic intrafractional motion or deformation using CBCT in radiotherapy. In the modalities outside CBCT, volumetric imaging methods available for non-periodic motion monitoring have been introduced using MRI and transperineal ultrasound, but the former has required quite a huge system combined with MRI and linac and the latter has been limited to monitor organs around the anorectal area [17, 18]. Regarding CT and CBCT as an 4D volumetric imaging modality, the temporal resolution and continuous time-ordered image reconstruction have been intensely studied in the literature in the concept of short-scan CT image reconstruction [19–21]. In this approach, the temporal resolution is directly related to the angular range, and a shorter range surely provides better temporal resolution. However, it might cause a degradation in the image quality and artifacts if the range is shorter than 180° plus the fan angle . On the other hand, Pang and Rowlands introduced “just-in-time tomography,” which reconstructs a digital tomosynthesis image from the projections of the cone-beam acquisition geometry . The method was actually a time-ordered imaging approach, but the created tomosynthesis images were still two-dimensional.
To overcome the above limitations, iterative reconstruction approaches have been proposed and are still actively studied [16, 23–30]. In particular, compressed sensing (CS)-based iterative reconstruction approaches such as total-variation-based compressed sensing (TVCS) and prior-image constrained compressed sensing (PICCS) have enabled image reconstruction with a limited number of projections and have strengthened the reliability of 4D CBCT images [16, 23, 24, 31]. The key idea of CS-based image reconstruction method is that, the sparsified image can be reconstructed from an undersampled projection data set instead of directly reconstructing a target image . The sparsified image, which is obtained by the application of sparsifying transform, contains significantly fewer image pixels which have significant pixel values. As the sparsifying transform, the ℓ 1-norm of the local spatial gradient is widely used for image reconstruction, which is called the total variation (TV). TVCS employs the TV transform of the target image itself, and PICCS employs the one of the difference image between the target and the prior image .
The non-periodic-motion correlated image reconstruction approach is based on a method that bunches projections into several time-phase groups. In contrast to the periodic motion case, there is no need to measure the motion signals for phase binning. After classifying the projections with the constant projection-angle interval, the images can be independently reconstructed in each time-phase group. In principle, however, the reconstructed time-ordered image sets are correlated with the prior image, which is the image prepared in advance of the iterative process, and the image sets of previous or subsequent time phases in a series of continuous 4D images. In the present study, we explicitly regard the continuous 4D reconstructed image series as sequential data and express them in the manner of a chain graph model. This is formulated as the regularization term of the statistical image reconstruction . In this study, we introduce a time-ordered 4D CBCT reconstruction method with the constraint of the time-ordered chain graph model (TCGM) and compare images reconstructed by previously introduced methods–TVCS and PICCS.
Materials and Methods
The projection data sets used for image reconstruction will be also introduced in this section, which include those of a digital phantom with the virtual motion of an air sphere and those of two clinical patients’ data in the pelvic region.
Brief review of statistical image reconstruction
gives reconstructed images. Hereafter the term lnP(μ ∗ ) in Eq. (8) is regarded as the regularization or constraint term, λ R(μ ∗ ), discussed in the next section. λ is the regularization parameter that describes the trade off between the regularization function term and the data fidelity term.
Brief review of the constraint term used in total-variation-based compressed sensing (TVCS) and prior-image constrained compressed sensing (PICCS)
where ε is a small constant that ensures that the total variation is differentiable at the origin. The TV operator Ψ is known as a sparsifying transform, and its ℓ 1-norm is to be minimized in compressed sensing.
where α and β are the weights of the TV and PIC terms, respectively. In the case of image reconstruction using PICCS, both constraint terms are combined and used as a weight; the values of 0.09 and 0.91 were respectively proposed for the TV and PIC terms by Chen et al.  TVCS reconstruction is implemented with the values of 1.0 and 0.0 for α and β, respectively.
Time-ordered chain graph model (TCGM)
As shown in Fig. 2, μ t and y t are the latent variables and observed data, respectively. Assuming the continuous measurement of projections, y t can be regarded as one projection angle datum for each t-th time phase and regarded as a group of several continuous projection angle data when assuming volumetric image reconstruction.
Figure 2 demonstrates the concept of the TCGM, in which the state of the t-th time phase located between the (t−1)-th and (t+1)-th time phases is constrained by the image objects of the previous time phase μ t−1 and the next subsequent time phase μ t+1. In all time phases, the images are renewed during the iteration step so that the constraints can be also dynamically changed. That is, the converged images of all time phases are concurrently obtained.
The first term on the right hand side of Eq. (13) is the distribution of the spatial difference between the previous and current phases, whereas the second term is that between the current and subsequent phases. In the case of the first and the last phases, the adjacent phase exists only in one side and Δ μ T is formed as doubled TV form of one-sided subtracted image.
The parameter set (α,β,γ) in Eq. (12) controls the relative weight of the three sparsity-promoting terms in the objective function. In this study, (α,β,γ) was chosen to be (1.0,0.0,0.0),(0.1,0.9,0.0), and (0.1,0.0,0.9), for TVCS, PICCS, and the TCGM, respectively. TVCS is used to create the initial images in PICCS and the TCGM. In addition, TVCS-reconstructed images are used as the prior-image constraint in PICCS.
Projection data sets for image reconstruction
Two patients’ projection data sets acquired for IGRT using the XVI system (version 4.2) mounted on the Synergy linac gantry were prepared. The data were acquired as pre-treatment CBCT imaging for the purpose of patient setup for prostate IMRT, and the chosen data sets contained the motion of rectal gas during projection acquisition. The acquisition time was approximately 120 s (696 projections) for patient A and 60 s (345 projections) for patient B.
The XVI system has three types of collimators, S, M, and L, which allow for the reconstruction of small, medium, and large fields of view (FOVs), respectively. With S collimator cassettes, the center of the kV radiation field is in line with the central axis of the kV X-ray source. With M and L collimator cassettes the kV radiation field is offset in the cross-line direction by 11.5 and 19 cm, respectively, at the flat-panel detector (FPD). The patients’ data sets in the present study were collected in the M FOV mode with an offset location of the FPD unit of 11.5 cm, whereas the geometry of the digital phantom projection assumed the XVI system with S FOV mode. The M FOV mode is frequently used for the projection acquisition of the abdominal or pelvic region, and the projections contain truncation [36–38]. Therefore, before the initial reconstruction using TVCS, the projection data were extended using the following procedure, which consists of three steps. First, the normal three-dimensional (3D) reconstruction containing the whole pelvis using all of the 360° projections was implemented. Second, virtual projections with 800 ×512 pixels were produced by the reprojection process of the reconstructed 3D image. Third, “extended" mosaic projections were created by combining real projections with “reproduced" peripheral area from virtual projections. Thus, the size of the extended 2D projection image was 800 ×512 for the clinical patient cases. This expansion allows for reconstruction of the whole body with 180° plus the fan angle, even in TVCS.
Image reconstruction was performed for nine continuous time phases with equal intervals in all TVCS, PICCS, and TCGM cases. An image matrix size of 256 ×256×60 with voxel size of 1 mm for the digital phantom reconstruction, and 400 ×400×60 with the same voxel size for patient data reconstruction was used. For TVCS, the projection data range was 180° + the fan angle = 200°; therefore, the projection data used in a certain phase were partly used in the reconstruction of the adjacent phases. For PICCS and the TCGM, the projection data range used was reduced by 90°. The projection data were still partly overlapped; however, the range was drastically decreased.
For the evaluation of our time-ordered reconstruction algorithms, the penumbral width of the moving air sphere along the CC direction was evaluated using the above digital phantom; reconstructed images of the digital phantom contained blurring of the air sphere due to longitudinal motion. The penumbral width of the longitudinal profile of the pixel values were compared among four reconstructions–normal 3D, TVCS (200° projection range), PICCS (90° projection range), and the TCGM (90° projection range). In this study, the penumbral width was defined as 10-90% penumbra of the pixel-inverse profile of the moving air cavity to compare the width of the time window.
Results of digital-phantom image reconstruction
Penumbral widths W penumbra of the moving sphere along the longitudinal axis
3 cm/60 s
3 cm/30 s
Results of patients’ pelvic image reconstruction
Consistency evaluation of pixel values in 4D reconstructed images
All 4D CBCT techniques previously proposed for IGRT devices were based on the assumption that the motion of the object of interest was periodic, such as lung motion due to respiration. With this assumption, the motion signal (e.g., respiratory signal) was used to classify the projection images into several motion-phase bins so that the volumetric-image series were reconstructed using projection images within each bin. These reconstructed images admittedly included the dimension of time. However, they were still “averaged” 4D images over the gantry rotation and could not represent non-periodic organ motion, including the baseline shift in the lung-tumor location, which is obtained as time progresses.
To obtain the time-ordered images, the projection data should be classified by time-ordered phase bins. This idea would work even in the FDK reconstruction algorithm if a rapid scanning system was available.
In the IGRT device available in our institution, however, the gantry rotation speed was limited by 6.0°/s, though C-arm angiographic system has much faster options of the rotation speed [39, 40] and linac gantry speeds faster than 6.0°/s are currently possible under specific conditions according to IEC 60601-2-1. Therefore, the projection angle range in each time phase has to be as narrow as possible –much less than 180° + the fan angle– to ensure a high time resolution in the reconstructed images. The FDK algorithm could not work now, which requires a range of at least 180° + the fan angle, owing to the Shannon-Nyquist sampling theorem. A possible approach is to relax this requirement, and compressed sensing with prior information is a prime candidate. So far, an FDK reconstruction image with full projection data was used as the prior information. On the other hand, the time-ordered images should be mutually correlated. Therefore, it is natural to introduce the time-chain model as the prior information, as in the proposed method.
In this study, the advantage of the TCGM over TVCS and PICCS was confirmed by the results of the 10-90% penumbral width measurement in the digital phantom; the temporal window width becomes narrower in PICCS and TCGM, and TCGM was narrower than PICCS. The reason why the TVCS result is the widest is obvious; simply because resulted by the widest projection range of 200° was used in this study in comparison with those of PICCS and the TCGM, which were both 90°. It is noted that a range less than 200° (= 180° + the fan angle) cannot be applied in TVCS for successful image reconstruction of 40-cm FOV without any prior information, since projections of at least 180° range for all voxels in the FOV of 40 cm are needed . Although the time-ordered images in the present study were reconstructed by extended projection image creation even in TVCS, TVCS has a disadvantage regarding the time resolution. On the other hand, the difference between PICCS and the TCGM is caused by the difference in the prior images used in the constraint terms. In the case of PICCS, an iteratively reconstructed image is always compared and constrained with the initial image, which is reconstructed using TVCS with a 200° projection range. In contrast, the TCGM uses the previous and subsequent images, which are being reconstructed in the same iterative process using 90° range of projections. Each 90° range overlaps with each other, which means that the related projection range is inherently narrower than that of PICCS.
The time resolution could be improved in both PICCS and the TCGM by increasing the number of time phases and narrowing the corresponding projection range, although there would be a trade-off regarding the image quality. We presented the nine-phase results with 90° projection ranges, and the results of clinical patients’ pelvic images were reasonably good; they showed rectal gas and flatus moving in the rectum. In the case of the TCGM, the images were reconstructed with fewer or weaker artifacts due to rectal gas, as shown in Fig. 9.
It is inherently possible to provide time-ordered 4D image series with a better temporal resolution with the presented TCGM method. The demonstrated results were reconstructed simply using just a 90° projection range and the weights of 0.1 and 0.9 for the TV and TCGM terms, and finding the optimized parameters is planned for future work. The behavior of the present method might be sensitive to the choice of weights. For the weight, one can introduce prior distributions and maximize probability function marginalized with respect to these hyper parameters to reconstruct images. More or less, this method definitely enables a reduction in the temporal window and provides image series of non-periodic time-ordered motion. The method should be used to understand the 4D dose distribution with in-treatment CBCT acquisition in adaptive radiotherapy scenarios [41–43].
The remarkable aspect of the present method is that the constraint term is based on the series of time-ordered statuses, which means that the t-th status is located between the (t+1)-th and (t−1)-th one and the constraint terms purely originate from the time-ordered change in the statuses. The present method can be improved by incorporating a deformation vector field to demonstrate organ motion and deformation using a more essential expression of motion.
It might be a limitation to reconstruct image volumes accompanied by sudden change and relatively fast motions, such as respiration, cough and swallowing, as the result of digital phantom study of double speed motion indicates. An initial motivation of the present method was to capture peristaltic activity of lower intestinal tract. So application of this method to the periodic motion like respiration might be resulting larger blurring compared to the 4D reconstruction from bunched projections using phase recognition introduced by the previous studies, [14, 15, 41] though base-line shift of diaphragm will probably be detected by the present method. Regarding sudden changes like cough and swallowing, it may be difficult to capture the time series of the motion using the present method. However the situation before and after the sudden change, e.g. status before and after cough or swallowing, can be reconstructed without strong motion artifacts caused by those sudden changes. Those images will probably have enough image quality in order for irradiated 4D dose reconstruction and evaluation of tumor shrinkage. Another limitation of the present method might be inconsistency of pixel values through all time phases. Pixel values, especially in peripheral areas of reconstructed FOV, might be fluctuated regarding projection ranges used in the reconstruction of particular phase. In the previous studies the methods to avoid those inconsistencies in dose reconstruction on CBCT images were already proposed, such as ROI mapping approach and WAB method . Those method can be applied to the dose reconstruction with image sets reconstructed using the present method.
A new concept of 4D CBCT reconstruction called the TCGM has been proposed, and reconstructed images obtained with this method have been compared to those obtained by TVCS and PICCS. The digital-phantom results demonstrated that the proposed method can provide 4D image series with a better temporal resolution compared to the two other methods. The clinical patients’ results show that the present method is applicable to the visualization of the motion of rectal gas and flatus in the rectum.
This work was partially supported by Grant-in-Aids from the Japan Society for the Promotion of Science (JSPS) KAKENHI JP Scientific Research (C), Grant Number 15K08691 and Research Activity Start-up, Grant Number 16H07429.
Availability of data and material
The data will not be shared, because the ethics committee does not allow sharing of the data.
MN, AH, JK and TM conceived the method. MN, AH, TM, SH and KN discussed the detailed algorithm and designed the project. MN and AH developed image reconstruction program. MN, AH, YM, SH and SK generated projection data of the digital phantom and analyzed the data. MN and AH drafted the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
The present study is ethically approved by institutional review board in the University of Tokyo Hospital. The reference number is 3372. Written informed consent was obtained from all patients whose data were used in this study.
Consent for publication
This is included in above consent.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Jaffray DA, Siewerdsen JH, Wong JW, Martinez AA. Flat-panel cone-beam computed tomography for image-guided radiation therapy. Int J Radiat Oncol Biol Phys. 2002; 53:1337–49.View ArticlePubMedGoogle Scholar
- Feldkamp LA, Davis LC, Kress JW. Practical cone-beam algorithm. J Opt Soc Am A. 1984; 1(6):612–9.View ArticleGoogle Scholar
- Boydev C, Taleb-Ahmed A, Derraz F, Peyrodie L, Thiran JP, Pasquier D. Development of CBCT-based prostate setup correction strategies and impact of rectal distension. Radiat Oncol. 2015; 10(1):83.View ArticlePubMedPubMed CentralGoogle Scholar
- Yao L, Zhu L, Wang J, Liu L, Zhou S, Jiang S, Cao Q, Qu A, Tian S.Positioning accuracy during VMAT of gynecologic malignancies and the resulting dosimetric impact by a 6-degree-of-freedom couch in combination with daily kilovoltage cone beam computed tomography. Radiat Oncol. 2015; 10(1):104.View ArticlePubMedPubMed CentralGoogle Scholar
- Oechsner M, Chizzali B, Devecka M, Combs SE, Wilkens JJ, Duma MN. Registration uncertainties between 3D cone beam computed tomography and different reference CT datasets in lung stereotactic body radiation therapy. Radiat Oncol. 2016; 11(1):142.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakagawa K, Shiraishi K, Nakagawa K, Shiraishi K, Kida S, Haga A, Yamamoto K, Saegusa S, Terahara A, Itoh S, et al. First report on prostate displacements immediately before and after treatment relative to the position during vmat delivery. Acta Oncologica. 2009; 48(8):1206–8.View ArticlePubMedGoogle Scholar
- Yang Y, Schreibmann E, Li T, Wang C, Xing L. Evaluation of on-board kv cone beam ct (cbct)-based dose calculation. Phys Med Biol. 2007; 52(3):685.View ArticlePubMedGoogle Scholar
- Létourneau D, Wong JW, Oldham M, Gulam M, Watt L, Jaffray DA, Siewerdsen JH, Martinez AA. Cone-beam-ct guided radiation therapy: technical implementation. Radiother Oncol. 2005; 75(3):279–86.View ArticlePubMedGoogle Scholar
- Nakagawa K, Haga A, Kida S, Masutani Y, Yamashita H, Takahashi W, Sakumi A, Saotome N, Shiraki T, Ohtomo K, Iwai Y. Yoda K. 4D registration and 4D verification of lung tumor position for stereotactic volumetric modulated arc therapy using respiratory-correlated cone-beam CT. J Radiat Res. 2013; 54(1):152–6.View ArticlePubMedGoogle Scholar
- Li J, Harrison A, Yu Y, Xiao Y, Werner-Wasik M, Lu B. Evaluation of Elekta 4D cone beam CT-based automatic image registration for radiation treatment of lung cancer. Br J Radiol. 2015; 88(1053):1–6.View ArticleGoogle Scholar
- Jin P, van Wieringen N, Maarten C, Hulshof CM, Bel A, Alderliesten T. 4D cone-beam CT imaging for guidance in radiation therapy: setup verification by use of implanted fiducial markers. In: Proc. SPIE Med. Imaging, vol. 9786.2016. p. 1–10.Google Scholar
- Santoso AP, Song KH, Qin Y, Gardner SJ, Liu C, Chetty IJ, Movsas B, Ajlouni M, Wen N. Evaluation of gantry speed on image quality and imaging dose for 4D cone-beam CT acquisition. Radiat Oncol. 2016; 11:98.View ArticlePubMedPubMed CentralGoogle Scholar
- Cai W, Dhou S, Cifter F, Myronakis M, Hurwitz MH, Christopher L. Williams, Ross I Berbeco, Joao Seco, and John H Lewis. 4D cone beam CT-based dose assessment for SBRT lung cancer treatment. Phys Med Biol. 2016; 61(2):554–68.View ArticlePubMedGoogle Scholar
- Zijp L, Sonke J-J, van Herk M. Extraction of the respiratory signal from sequential thorax cone-beam x-ray images. Int Conf Use Comput Radiat Ther. 2014; 1:507–9.Google Scholar
- Sonke JJ, Zijp L, Remeijer P, van Herk M. Respiratory correlated cone beam ct. Med Phys. 2005; 32(4):1176–86.View ArticlePubMedGoogle Scholar
- Lauzier PT, Tang J, Chen GH. Time-resolved cardiac interventional cone-beam CT reconstruction from fully truncated projections using the prior image constrained compressed sensing (PICCS) algorithm. Phys Med Biol. 2012; 57(9):2461–76.View ArticlePubMedPubMed CentralGoogle Scholar
- Stam MK, Crijns SPM, Zonnenberg BA, Barendrecht MM, van Vulpen M, Lagendijk JJW, Raaymakers BW. Navigators for motion detection during real-time MRI-guided radiotherapy. Phys Med Biol. 2012; 57(21):6797–805.View ArticlePubMedGoogle Scholar
- Trivedi A, Ashikaga T, Hard D, Archambault J, Lachaine M, Cooper DT, Wallace III JH. Development of 3-dimensional transperineal ultrasound for image guided radiation therapy of the prostate: Early evaluations of feasibility and use for inter- and intrafractional prostate localization. Practical Radiat Oncol. 2017; 7(1):e27–33.View ArticleGoogle Scholar
- Parker DL. Optimal short scan convolution reconstruction for fan beam ct. Med Phys. 1982; 9(2):254–7.View ArticlePubMedGoogle Scholar
- Noo F, Defrise M, Clackdoyle R, Kudo H. Image reconstruction from fan-beam projections on less than a short scan. Phys Med Biol. 2002; 47(14):2525–46.View ArticlePubMedGoogle Scholar
- Liu Y, Liu H, Wang Y, Wang G. Half-scan cone-beam CT fluoroscopy with multiple x-ray sources. Med Phys. 2001; 28(7):1466.View ArticlePubMedGoogle Scholar
- Pang G, Rowlands JA. Just-in-time tomography (jitt): a new concept for image-guided radiation therapy. Phys Med Biol. 2005; 50(21):N323.View ArticlePubMedGoogle Scholar
- Chen GH, Tang J, Leng S. Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. Med Phys. 2008; 35(2):660–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Sidky EY, Pan X. Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization. Phys Med Biol. 2008; 53(17):4777–807.View ArticlePubMedPubMed CentralGoogle Scholar
- Tang J, Nett BE, Chen GH. Performance comparison between total variation (TV)-based compressed sensing and statistical iterative reconstruction algorithms. Phys Med Biol. 2009; 54(19):5781–804.View ArticlePubMedPubMed CentralGoogle Scholar
- Bian J, Siewerdsen JH, Han X, Sidky EY, Prince JL, Pelizzari CA, Pan X. Evaluation of sparse-view reconstruction from flat-panel-detector cone-beam CT. Phys Med Biol. 2010; 55:6575–99.View ArticlePubMedPubMed CentralGoogle Scholar
- Jia X, Dong B, Lou Y, Jiang SB. GPU-based iterative cone-beam CT reconstruction using tight frame regularization. Phys Med Biol. 2011; 56(13):3787–807.View ArticlePubMedGoogle Scholar
- Yan H, Wang X, Shi F, Bai T, Folkerts M, Cervino L, Jiang SB, Jia X. Towards the clinical implementation of iterative low-dose cone-beam CT reconstruction in image-guided radiation therapy: Cone/ring artifact correction and multiple GPU implementation. Med Phys. 2014; 41(11):111912.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen GH, Theriault-Lauzier P, Student Member, IEEE, Tang J, Nett B, Leng S, Zambelli J, Qi Z, Bevins N, Raval A, Reeder S, Rowley H. Time-Resolved Interventional Cardiac C-arm Cone-Beam CT: An Application of the PICCS Algorithm. IEEE Trans. Med. Imaging. 2012; 31(4):907–23.View ArticlePubMedGoogle Scholar
- Chen GH, Yinsheng L. Synchronized multiartifact reduction with tomographic reconstruction (SMART-RECON): A statistical model based iterative image reconstruction method to eliminate limited-view artifacts and to mitigate the temporal-average artifacts in time-resolved CT. Med Phys. 2015; 42(8):4698–707.View ArticlePubMedPubMed CentralGoogle Scholar
- Candės EJ, Romberg J, Tao T. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information. IEEE Trans Inf Theory. 2006; 52(2):489–509.View ArticleGoogle Scholar
- Koller D, Friedman N. Probabilistic graphical models: principles and techniques. Cambridge: MIT Press; 2009.Google Scholar
- Hanson KM, Wecksung GW. Bayesian approach to limited-angle reconstruction in computed tomography. J Opt Soc Am. 1983; 73(11):1501.View ArticleGoogle Scholar
- Green PJ. Bayesian reconstructions from emission tomography data using a modified EM algorithm. IEEE Trans. Med. Imaging. 1990; 9(1):84–93.View ArticlePubMedGoogle Scholar
- Lange K, Fessler JA. Globally convergent algorithms for maximum a posteriori transmission tomography. IEEE Trans. Image Process. 1995; 4(10):1430–8.View ArticlePubMedGoogle Scholar
- Lehmann J, Perks J, Semon S, Harse R, Purdy A. Commissioning experience with cone-beam computed tomography for image-guided radiation therapy. 2007; 8(3):21–36.Google Scholar
- Oelfke U, Tu̇cking T, Nill S, Seeber A, Hesse B, Huber P, Thilmann C. Linac-integrated kV-cone beam CT: technical features and first applications. Med Dosim. 2006; 31(1):62–70.View ArticlePubMedGoogle Scholar
- Hansis E, Bredno J, Sowards-Emmerd D, Shao L. Iterative reconstruction for circular cone-beam CT with an offset flat-panel Detector. IEEE Nucl Sci Symp Med Imaging Conf. 2010; 1:2228–31.View ArticleGoogle Scholar
- Sykes JR, Amer A, Czajka J, Moore CJ. A feasibility study for image guided radiotherapy using low dose, high speed, cone beam X-ray volumetric imaging. Radiother Oncol. 2005; 77(1):45–52.View ArticlePubMedGoogle Scholar
- Bai M, Liu B, Mu H, Liu X, Jiang Y. The comparison of radiation dose between C-arm flat-detector CT (DynaCT) and multi-slice CT (MSCT). A phantom study. Eur J Radiol. 2012; 81(11):3577–80.View ArticlePubMedGoogle Scholar
- Kida S, Masutani Y, Yamashita H, Imae T, Matsuura T, Saotome N, Ohtomo K, Nakagawa K, Haga A. In-treatment 4D cone-beam CT with image-based respiratory phase recognition. Radiol Phys Technol. 2012; 5(2):138–47.View ArticlePubMedGoogle Scholar
- Nakagawa K, Haga A, Shiraishi K, Yamashita H, Igaki H, Terahara A, Ohtomo K, Saegusa S, Shiraki T, Oritate T, Yoda K. First clinical cone-beam CT imaging during volumetric modulated arc therapy. Radiother Oncol. 2009; 90(3):422–3.View ArticlePubMedGoogle Scholar
- Nakagawa K, Yamashita H, Shiraishi K, Igaki H, Terahara A, Nakamura N, Ohtomo K, Saegusa S, Shiraki T, Oritate T, Yoda K. Verification of in-treatment tumor position using kilovoltage cone-beam computed tomography: a preliminary study. Int J Radiat Oncol Biol Phys. 2007; 69(4):970–3.View ArticlePubMedGoogle Scholar
- Fotina I, Hopfgartner J, Stock M, Steininger T, Lutgendorf-Caucig C, Georg D, Yan H, Tian Z, Shao Y, Jiang SB, Jia X. Feasibility of CBCT-based dose calculation: Comparative analysis of HU adjustment techniques. Radiother Oncol. 2012; 104(2):249–56.View ArticlePubMedGoogle Scholar