- Short report
Reconstruction of the treatment area by use of sinogram in helical tomotherapy
Radiation Oncologyvolume 9, Article number: 252 (2014)
TomoTherapy (Accuray, USA) has an image-guided radiotherapy system with a megavoltage (MV) X-ray source and an on-board imaging device. This system allows one to acquire the delivery sinogram during the actual treatment, which partly includes information from the irradiated object. In this study, we try to develop image reconstruction during treatment with helical tomotherapy.
Sinogram data were acquired during helical tomotherapy delivery using an arc-shaped detector array that consists of 576 xenon-gas filled detector cells. In preprocessing, these were normalized with full air-scan data. A software program was developed that reconstructs 3D images during treatment with corrections as; (1) the regions outside the field were masked not to be added in the backprojection (a masking correction), and (2) each voxel of the reconstructed image was divided by the number of the beamlets passing through its voxel (a ray-passing correction).
The masking correction produced a reconstructed image, however, it contained streak artifacts. The ray-passing correction reduced this artifact. Although the SNR (the ratio of mean to standard deviation in a homogeneous region) and the contrast of the reconstructed image were slightly improved with the ray-passing correction, use of only the masking correction was sufficient for the visualization purpose.
The visualization of the treatment area was feasible by using the sinogram in helical tomotherapy. This proposed method would be useful in the treatment verification.
As radiotherapy is complex, treatment verification becomes significant. The evaluation of the absorbed dose in phantoms is strongly recommended for all patients having intensity-modulated radiotherapy (IMRT) . In addition, accuracy of the patient setup is more important in the IMRT than that in the conventional radiotherapy. Image-guided radiotherapy (IGRT) can entail correcting the patient position just prior to treatment by gathering information about anatomical locations during setup. IGRT can utilize various imaging technologies such as the portal images of the treatment beam -, magnetic resonance imaging , ultrasound ,, and computed tomography (CT) ,.
TomoTherapy® has an IGRT system with megavoltage (MV) X-ray source and an on-board imaging device . With the MV CT, it became feasible to perform efficient daily-3D registration of the patient position before each treatment delivery. This system also allows one to acquire the delivery sinogram during the actual treatment. The sinogram has often been used in the treatment verification ,, and one can come up with the visualization of treatment area from the sinogram. For conventional linear accelerators, in fact, CT reconstruction with portal images during rotational treatment such as a volumetric modulated arc therapy (VMAT) has been successfully performed ,. Also the MV CT reconstruction during treatment in helical TomoTherapy® delivery has been first tried in Ref. , where the insertion of full field-of-view (FOV) beamlets was cooperated with the treatment sinogram. In this note, we focus on the image reconstruction of treatment area without full FOV insertion.
For the preliminary arrangement, we developed a helical CT reconstruction algorithm that includes corrections for the heterogeneous beam profile and the geometrical disagreement between the X-ray source position and the detector curvature. Then, the feasibility of the reconstruction of the treatment area was examined with the treatment sinogram to a phantom, which includes the information of the irradiated part of the phantom.
Materials and methods
All data were obtained on a TomoTherapy® unit at Accuray in Madison, Wisconsin. The factory IMRT test plan was irradiated to TomoPhantom® (Accuray, USA) with density plugs. This test plan prescribes 10 Gy for two cylindrical targets (3 cm radius and 6 cm length) which are seen as blue and sky-blue circles in Figure 1(a). Sinogram data were acquired during helical TomoTherapy® deliveries using an arc-shaped detector array consisting of 576 xenon-gas filled detector cells, of which the data from middle 527 cells were used in the reconstruction. The data has the following acquisition properties: gantry rotation period 14 s, data sampling rate 300 Hz, starting view angle 229.371°, and couch travel distance over entire data sample 90.6 mm. There has 52,788 samples, and thus, the couch speed was 0.515 mm/s.
The dose distribution in the treatment planning system (TPS) and a part of the corresponding delivery sinogram are shown in Figures 1(a) and (b), respectively. The sinogram was normalized with full air-scan data to the correct heterogeneous beam profile and modulated beam intensity of TomoTherapy®.
The reconstruction was performed with an in-house program employing the filtered back projection (FBP) algorithm using Shepp-Logan filter. Because the source-to-isocenter distance (85.0 cm) is smaller than the detector radius of curvature (99.8 cm), our program converts the original data in each detector cell into the virtual one with the curvature corresponding source-to-isocenter geometry. Then, the virtual data with a constant cell-to-cell interval was created by linear interpolation.
In general, it is impossible to make a correct reconstruction with a limited FOV using FBP . As shown in Figure 1(b), the area blocked by binary multileaf collimator (MLC) in the sinogram has a lower X-ray intensity than that inside the FOV, so that the conventional reconstruction scheme does not successfully visualize the object. Instead this yields an unrealistically high-attenuation area outside the irradiated site in the object. For the visualization of the irradiated area, therefore, we employed two corrections. One was a masking correction, which masks the area outside the FOV so as not to include this area in the backprojection process. For this, sinogram normalized by full air scan, where β and γ denote the gantry and fan angles, respectively, can be expressed as,
where p’ means the threshold for masking region, and here, we employed p’ = 0.2. With this correction, the outside field is regarded as air and the boundary of masking region is discontinuous. Of course, this is not true, but it enhances the information from the irradiated area in the FBP reconstruction scheme.
The other correction was a ray-passing adjustment, which normalized each voxel of the reconstructed image to the number of the X-rays passing through the corresponding voxel. Namely, using the masking function,
the correction factor for ray-passing can be expressed as,
The backprojection generates stronger signals from the angles passing more X-rays. The ray-passing correction corrects this effect. The reconstructed region was controlled by β’ in Eq. (3). In this study, the area irradiated with more than 35% (R(γ, β’) = 0.35) and 55% (R(γ, β’) = 0.55) of the maximum was reconstructed and the other area was masked.
The contrast in the images was evaluated by the ratio of the signal in high-density regions to that in low-density regions in the object. The homogeneity was evaluated from the three regions that are composed of the same material. The signal-to-noise ratio (SNR) inside the region-of-interest (ROI) was also evaluated.
Results and discussion
Figures 2(a)-(f) show the reconstructed images using a full FOV, with no masking and no ray-passing corrections, with masking correction only, and with both corrections, respectively. Without masking and ray-passing corrections (Figure 2(b)), the visualization of the treatment area was poor. Use of the masking correction and the ray-passing correction improved the reconstructed image for the treatment site (Figures 2(c)-(e)). The masking correction made the image clearly visible, however, a streak artifact was introduced (Figure 2(c)). The ray-passing correction reduced this artifact (Figures 2(d) and (e)). In Figure 2(f), the image difference between Figures 2(a) and (e) is shown. It is found that the location and the size of inner plugs were well reproduced in Figure 2(f). On the other hand, the difference can depend on the density plug, because of the low contrast of Figure 2(e) in comparison with Figure 2(a).
Table 1 shows the quantitative evaluation of contrast, homogeneity, and SNR. The contrast index was calculated by the ratio of the mean value of regions A, C, and D to that of region B (see Figure 3). As expected, the contrast for the in-treatment images is considerably poorer than for the full FOV image. The homogeneity, as evaluated by the min-max values of the mean value in the regions A, C, and D has the same tendency. On the other hand, the SNR can be enhanced in the in-treatment images, presumably, due to the blurring effect.
In the analysis of image contrast and homogeneity, the ray-passing correction improved the image quality, but no dramatic change in visibility was yielded. Of course, this is not a general conclusion. One set of questions might be to further examine how well this type of technique works for different cases, such as different target sizes, different anatomical regions, and with different levels of leaf modulation.
Although a further study will be required, the present result encouraged us to develop the record-and-verify system with the reconstructed delivery area from the actual treatment. Although a further study will be required, the present result encouraged us to develop the record-and-verify system with the reconstructed delivery area from the actual treatment. Also one may be interested in the dose reconstruction using present method. The present method cannot be applied for the dose reconstruction directly. However, the dose reconstruction in each treatment session requires the information of the patient location during treatment, which can be provided by the present method. Thus, the development of the image reconstruction using the delivery sinogram would be a promising tool for in-vivo dosimetry as well as for verification of irradiated areas.
A reconstruction technique using the treatment sinogram has been developed for helical Tomotherapy. The improved visibility of structures in the reconstructed image makes this a promising tool for verifying relative anatomical positions during the course of a treatment.
ICRU Report 83: Prescribing, recording, and reporting photon-beam intensity-modulated radiation therapy (IMRT) J ICRU 2010, 10: 1-106. 10.1093/jicru/ndq003
Balter JM, Lam KL, Sandler HM, Littles JF, Bree RL, Ten Haken RK: Automated localization of the prostate at the time of treatment using implanted radiopaque markers: technical feasibility. Int J Radiat Oncol Biol Phys 1995, 33: 1281-1286. 10.1016/0360-3016(95)02083-7
Alasti H, Petric MP, Catton CN, Warde PR: Portal imaging for evaluation of daily on-line setup errors and off-line organ motion during conformal irradiation of carcinoma of the prostate. Int J Radiat Oncol Biol Phys 2001, 49: 869-884. 10.1016/S0360-3016(00)01446-2
Bergstrom P, Lofrot PO, Widmark A: High-precision conformal radiotherapy (HPCRT) of prostate cancer––a new technique for exact positioning of the prostate at the time of treatment. Int J Radiat Oncol Biol Phys 1998, 42: 305-311. 10.1016/S0360-3016(98)00229-6
Vigneault E, Pouliot J, Laverdiere J, Roy J, Dorion M: Electronic portal imaging device detection of radioopaque markers for the evaluation of prostate position during megavoltage irradiation: a clinical study. Int J Radiat Oncol Biol Phys 1997, 37: 205-212. 10.1016/S0360-3016(96)00341-0
Mah D, Freedman G, Milestone B, Hanlon A, Palacio E, Richardson T, Movsas B, Mitra R, Horwitz E, Hanks GE: Measurement of intrafractional prostate motion using magnetic resonance imaging. Int J Radiat Oncol Biol Phys 2002, 54: 568-575. 10.1016/S0360-3016(02)03008-0
Chandra A, Dong L, Huang E, Kuban DA, O’Neill L, Rosen I, Pollack A, Kuban DA, O’Neill L, Rosen I, Pollack A: Experience of ultrasound-based daily prostate localization. Int J Radiat Oncol Biol Phys 2003, 56: 436-447. 10.1016/S0360-3016(02)04612-6
Morr J, DiPetrillo T, Tsai JS, Engler M, Wazer DE: Implementation and utility of a daily ultrasound-based localization system with intensity-modulated radiotherapy for prostate cancer. Int J Radiat Oncol Biol Phys 2002, 53: 1124-1129. 10.1016/S0360-3016(02)02820-1
van Herk M, Bruce A, Kroes AP, Shouman T, Touw A, Lebesque JV: Quantification of organ motion during conformal radiotherapy of the prostate by three dimensional image registration. Int J Radiat Oncol Biol Phys 1995, 33: 1311-1320. 10.1016/0360-3016(95)00116-6
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-1349. 10.1016/S0360-3016(02)02884-5
Ruchala KJ, Olivera GH, Schloesser EA, Mackie TR: Megavoltage CT on a tomotherapy system. Phys Med Biol 1999, 44: 2597-2621. 10.1088/0031-9155/44/10/316
Sevillano D, Mínguez C, Sánchez A, Sánchez-Reyes A: Measurement and correction of leaf open times in helical tomotherapy. Med Phys 2012, 39: 6972-6980. 10.1118/1.4762565
Hashimoto M, Uematsu M, Ito M, Hama Y, Inomata T, Fujii M, Nishio T, Nakamura N, Nakagawa K: Investigation of the feasibility of a simple method for verifying the motion of a binary multileaf collimator synchronized with rotation of the gantry for helical tomotherapy. J Appl Clin Med Phys 2012, 13: 27-43.
Poludniowski G, Thomas MDR, Evans PM, Webb S: CT reconstruction from portal images acquired during volumetric-modulated arc therapy. Phys Med Biol 2011, 55: 5635-5651. 10.1088/0031-9155/55/19/002
Kida S, Saotome N, Masutani Y, Yamashita H, Ohtomo K, Nakagawa K, Sakumi A, Haga A: 4D-CBCT reconstruction using MV portal imaging during volumetric modulated arc therapy. Radiother Oncol 2011, 100: 380-385. 10.1016/j.radonc.2011.08.047
Ruchala KJ, Olivera GH, Kapatoes JM, Schloesser EA, Reckwerdt PJ, Mackie TR: Megavoltage CT image reconstruction during tomotherapy treatments. Phys Med Biol 2000, 45: 3545-3562. 10.1088/0031-9155/45/12/303
The authors would like to thank Tadashi Nakabayashi and Juki Hozumi (Accuray Japan) for a lot of efforts about preparing this present study. This work was supported by JSPS KAKENHI 24234567.
K.N. received the research grant from Accuray Japan.
KN conceived the idea. CM, KR, EC and DC (Accuray group) generated the factory plan and all data used in this manuscript. The reconstruction method was discussed among AH, SK, DS, MN, TM, and YM (Tokyo group). Based on the discussion, AH developed the program and analyzed the data obtained by Accuray group. AH drafted the manuscript. All authors read and approved the final manuscript.