Search

Purpose. Fast kV-switching (FKS) and dual-layer flat-panel detector (DL-FPD) technologies have been actively studied as promising dual-energy spectral imaging solutions for FPD-based cone-beam computed tomography (CT). However, cone-beam CT (CBCT) spectral imaging is known to face challenges in obtaining accurate and robust material discrimination performance. That is because the energy separation by either FKS or DL-FPD, alone, is still limited, along with apparently unpaired signal levels in the effective low- and high-energy projections in real applications, not to mention the x-ray scatter in cone-beam scan which will make the material decomposition almost impossible if no correction is applied. To further improve CBCT spectral imaging capability, this work aims to promote a source-detector joint multi-energy spectral imaging solution which takes advantages of both FKS and DL-FPD, and to conduct a feasibility study on the first tabletop CBCT system with the joint spectral imaging capability developed. Methods. For CBCT, development of multi-energy spectral imaging can be jointly realized by using an x-ray source with a generator whose kilo-voltages can alternate in tens of Hertz (i.e. FKS), and a DL-FPD whose top- and bottom-layer projections corresponds to different effective energy levels. Thanks to the complimentary characteristics inherent in FKS and DL-FPD, the overall energy separation will be significantly better when compared with FKS or DL-FPD alone, and the x-ray photon detection efficiency will be also improved when compared with FKS alone. In this work, a noise performance analysis using the Cramér–Rao lower bound (CRLB) method is conducted. The CRLB for basis material after a projection-domain material decomposition is derived, followed by a set of numerical calculations of CRLBs, for the FKS, the DL-FPD and the joint solution, respectively. To compensate for the slightly angular mismatch between low- and high- projections in FKS, a dual-domain projection completion scheme is implemented. Afterwards material decomposition from the complete projection data is carried out by using the maximum-likelihood method, followed by reconstruction of basis material and virtual monochromatic images (VMI). In this work, the first FKS and DL-FPD jointly enabled multi-energy tabletop CBCT system, to the best of our knowledge, has been developed in our laboratory. To evaluate its spectral imaging performance, a set of physics experiments are conducted, where the multi-energy and head phantoms are scanned using the 80/105/130 kVp switching pairs and projection data are collected using a prototype DL-FPD, whose both top and bottom layer of panels are composed of 550 μ m of cesium iodine (CsI) scintillators with no intermediate metal filter in-between. Results. The numerical simulations show that the joint spectral imaging solution can lead to a significant improvement in energy separation and lower noise levels in most of material decomposition cases. The physics experiments confirmed the feasibility and superiority of the joint spectral imaging, whose CNRs in the selected regions of interest of the multi-energy phantom were boosted by an average improvement of 21.9%, 20.4% for water basis images and 32.8%, 62.8% for iodine images when compared with that of the FKS and DL-FPD, respectively. For the head phantom case, the joint spectral imaging can effectively reduce the streaking artifacts as well, and the standard deviation in the selected regions of interest are reduced by an average decrement of 19.5% and 8.1% for VMI when compared with that of the FKS and DL-FPD, respectively. Conclusions. A feasibility study of the joint spectral imaging solution for CBCT by utilizing both the FKS and DL-FPD was conducted, with the first tabletop CBCT system having such a capability being developed, which exhibits improved CNR and is more effective in avoiding streaking artifacts as expected. [ABSTRACT FROM AUTHOR]

Objective. The aim of this study was to reconstruct volumetric computed tomography (CT) images in real-time from ultra-sparse two-dimensional x-ray projections, facilitating easier navigation and positioning during image-guided radiation therapy. Approach. Our approach leverages a voxel-sapce-searching Transformer model to overcome the limitations of conventional CT reconstruction techniques, which require extensive x-ray projections and lead to high radiation doses and equipment constraints. Main results. The proposed XTransCT algorithm demonstrated superior performance in terms of image quality, structural accuracy, and generalizability across different datasets, including a hospital set of 50 patients, the large-scale public LIDC-IDRI dataset, and the LNDb dataset for cross-validation. Notably, the algorithm achieved an approximately 300% improvement in reconstruction speed, with a rate of 44 ms per 3D image reconstruction compared to former 3D convolution-based methods. Significance. The XTransCT architecture has the potential to impact clinical practice by providing high-quality CT images faster and with substantially reduced radiation exposure for patients. The model's generalizability suggests it has the potential applicable in various healthcare settings. [ABSTRACT FROM AUTHOR]

Objective. To demonstrate that complete cone beam CT (CBCT) scans from both MV-energy and kV-energy LINAC sources can reduce metal artifacts in radiotherapy guidance, while maintaining standard-of-care x-ray doses levels. Approach. MV-CBCT and kV-CBCT scans are acquired at half normal dose. The impact of lowered dose on MV-CBCT data quality is mitigated by the use of a 4-layer MV-imager prototype and reduced LINAC energy settings (2.5 MV) to improve photon capture. Additionally, the MV-CBCT is used to determine the 3D position and pose of metal implants, which in turn is used to guide model-based poly-energetic correction and interleaving of the kV-CBCT and MV-CBCT data. Certain edge-preserving regularization steps incorporated into the model-based correction algorithm further reduce MV data noise. Main results. The method was tested in digital phantoms and a real pelvis phantom with large 2.5″ spherical inserts, emulating hip replacements of different materials. The proposed method demonstrated an appealing compromise between the high contrast of kV-CBCT and low artifact content of MV-CBCT. Contrast-to-noise improved 3-fold compared to MV-CBCT with a clinical 1-layer architecture at matched dose (37 mGy) and edge blur levels. Visual delineation of the bladder and prostate improved noteably over kV- or MV-CBCT alone. Significance. The proposed method demonstrates that a full MV-CBCT scan can be combined with kV-CBCT to reduce metal artifacts without resorting to complicated beam collimation strategies to limit the MV-CBCT dose contribution. Additionally, significant improvements in CNR can be achieved as compared to metal artifact reduction through current clinical MV-CBCT practices. [ABSTRACT FROM AUTHOR]

Cone-beam computed tomography (CBCT) is widely used in dental imaging, small animal imaging, radiotherapy, and non-destructive industrial inspection. The quality of CBCT images depends on the precise knowledge of the CBCT system's alignment. We introduce a distinct procedure, 'precision alignment loop (PAL)', to calibrate any CBCT system with a circular trajectory. We describe the calibration procedure by using a line-beads phantom, and how PAL determines the misalignments from a CBCT system. PAL also yields the uncertainties in the simulated calibration to give an estimate of the errors in the misalignments. From the analytical simulations, PAL can precisely obtain the source-to-rotation axis distance, and the geometric center G, 'the point in z -axis meets the detector', where the z -axis is coincident with the line from the x-ray source that intersects the axis of the rotation orthogonally. The uncertainties of three misalignment angles of the detector are within ±0.05°, which is close to ±0.04° for the results of Yang et al (2017 Med. Phys. 44 1692), but our method is easy and simple to implement. Our distinct procedure, on the other hand, yields the calibration of a micro-CT system and an example of reconstructed images, showing our calibration method for the CBCT system to be simple, precise, and accurate. [ABSTRACT FROM AUTHOR]

Objective. To create two non-coplanar, stereotactic ablative radiotherapy (SABR) lung patient treatment plans compliant with the radiation therapy oncology group (RTOG) 0813 dosimetric criteria using a simple, isocentric, therapy with kilovoltage arcs (SITKA) system designed to provide low cost external radiotherapy treatments for low- and middle-income countries (LMICs). Approach. A treatment machine design has been proposed featuring a 320 kVp x-ray tube mounted on a gantry. A deep learning cone-beam CT (CBCT) to synthetic CT (sCT) method was employed to remove the additional cost of planning CTs. A novel inverse treatment planning approach using GPU backprojection was used to create a highly non-coplanar treatment plan with circular beam shapes generated by an iris collimator. Treatments were planned and simulated using the TOPAS Monte Carlo (MC) code for two lung patients. Dose distributions were compared to 6 MV volumetric modulated arc therapy (VMAT) planned in Eclipse on the same cases for a Truebeam linac as well as obeying the RTOG 0813 protocols for lung SABR treatments with a prescribed dose of 50 Gy. Main results. The low-cost SITKA treatments were compliant with all RTOG 0813 dosimetric criteria. SITKA treatments showed, on average, a 6.7 and 4.9 Gy reduction of the maximum dose in soft tissue organs at risk (OARs) as compared to VMAT, for the two patients respectively. This was accompanied by a small increase in the mean dose of 0.17 and 0.30 Gy in soft tissue OARs. Significance. The proposed SITKA system offers a maximally low-cost, effective alternative to conventional radiotherapy systems for lung cancer patients, particularly in low-income countries. The system's non-coplanar, isocentric approach, coupled with the deep learning CBCT to sCT and GPU backprojection-based inverse treatment planning, offers lower maximum doses in OARs and comparable conformity to VMAT plans at a fraction of the cost of conventional radiotherapy. [ABSTRACT FROM AUTHOR]

Objective. This study aims at investigating a novel super resolution CBCT imaging approach with a dual-layer flat panel detector (DL-FPD). Approach. With DL-FPD, the low-energy and high-energy projections acquired from the top and bottom detector layers contain over-sampled spatial information, from which super-resolution CT images can be reconstructed. A simple mathematical model is proposed to explain the signal formation procedure in DL-FPD, and a dedicated recurrent neural network, named suRi-Net, is developed based upon the above imaging model to nonlinearly retrieve the high-resolution dual-energy information. Physical benchtop experiments are conducted to validate the performance of this newly developed super-resolution CBCT imaging method. Main Results. The results demonstrate that the proposed suRi-Net can accurately retrieve high spatial resolution information from the low-energy and high-energy projections of low spatial resolution. Quantitatively, the spatial resolution of the reconstructed CBCT images from the top and bottom detector layers is increased by about 45% and 54%, respectively. Significance. In the future, suRi-Net will provide a new approach to perform high spatial resolution dual-energy imaging in DL-FPD-based CBCT systems. [ABSTRACT FROM AUTHOR]

Objective. Metal artifact reduction (MAR) has been a key issue in CT imaging. Recently, MAR methods based on deep learning have achieved promising results. However, when deploying deep learning-based MAR in real-world clinical scenarios, two prominent challenges arise. One limitation is the lack of paired training data in real applications, which limits the practicality of supervised methods. Another limitation is that image-domain methods suitable for more application scenarios are inadequate in performance while end-to-end approaches with better performance are only applicable to fan-beam CT due to large memory consumption. Approach. We propose a novel image-domain MAR method based on the generative adversarial network with variable constraints (MARGANVAC) to improve MAR performance. The proposed variable constraint is a kind of time-varying cost function that can relax the fidelity constraint at the beginning and gradually strengthen the fidelity constraint as the training progresses. To better deploy our image-domain supervised method into practical scenarios, we develop a transfer method to mimic the real metal artifacts by first extracting the real metal traces and then adding them to artifact-free images to generate paired training data. Main results. The effectiveness of the proposed method is validated in simulated fan-beam experiments and real cone-beam experiments. All quantitative and qualitative results demonstrate that the proposed method achieves superior performance compared with the competing methods. Significance. The MARGANVAC model proposed in this paper is an image-domain model that can be conveniently applied to various scenarios such as fan beam and cone beam CT. At the same time, its performance is on par with the cutting-edge dual-domain MAR approaches. In addition, the metal artifact transfer method proposed in this paper can easily generate paired data with real artifact features, which can be better used for model training in real scenarios. [ABSTRACT FROM AUTHOR]

In vivo dosimetry (IVD) is an important tool in external beam radiotherapy (EBRT) to detect major errors by assessing differences between expected and delivered dose and to record the received dose by individual patients. Also, in intraoperative radiation therapy (IORT), IVD is highly relevant to register the delivered dose. This is especially relevant in low-risk breast cancer patients since a high dose of IORT is delivered in a single fraction. In contrast to EBRT, online treatment planning based on intraoperative imaging is only under development for IORT. Up to date, two commercial treatment planning systems proposed intraoperative ultrasound or in-room cone-beam CT for real-time IORT planning. This makes IVD even more important because of the possibility for real-time treatment adaptation. Here, we summarize recent developments and applications of IVD methods for IORT in clinical practice, highlighting important contributions and identifying specific challenges such as a treatment planning system for IORT. HDR brachytherapy as a delivery technique was not considered. We add IVD for ultrahigh dose rate (FLASH) radiotherapy that promises to improve the treatment efficacy, when compared to conventional radiotherapy by limiting the rate of toxicity while maintaining similar tumour control probabilities. To date, FLASH IORT is not yet in clinical use. [ABSTRACT FROM AUTHOR]

Objective. To compare a not adapted (NA) robust planning strategy with three fully automated online adaptive proton therapy (OAPT) workflows based on the same optimization method: dose mimicking (DM). The added clinical value and limitations of the OAPT methods are investigated for head and neck cancer (HNC) patients. Approach. The three OAPT strategies aimed at compensating for inter-fractional anatomical changes by mimiking different dose distributions on corrected cone beam CT images (corrCBCTs). Order by complexity, the OAPTs were: (1) online adaptive dose restoration (OADR) where the approved clinical dose on the planning-CT (pCT) was mimicked, (2) online adaptation using DM of the deformed clinical dose from the pCT to corrCBCTs (OADEF), and (3) online adaptation applying DM to a predicted dose on corrCBCTs (OAML). Adaptation was only applied in fractions where the target coverage criteria were not met (D98% < 95% of the prescribed dose). For 10 HNC patients, the accumulated dose distributions over the 35 fractions were calculated for NA, OADR, OADEF, and OAML. Main results. Higher target coverage was observed for all OAPT strategies compared to no adaptation. OADEF and OAML outperformed both NA and OADR and were comparable in terms of target coverage to initial clinical plans. However, only OAML provided comparable NTCP values to those from the clinical dose without statistically significant differences. When the NA initial plan was evaluated on corrCBCTs, 51% of fractions needed adaptation. The adaptation rate decreased significantly to 25% when the last adapted plan with OADR was selected for delivery, to 16% with OADEF, and to 21% with OAML. The reduction was even greater when the best plan among previously generated adapted plans (instead of the last one) was selected. Significance. The implemented OAPT strategies provided superior target coverage compared to no adaptation, higher OAR sparing, and fewer required adaptations. [ABSTRACT FROM AUTHOR]

Objective. CBCTs in image-guided radiotherapy provide crucial anatomy information for patient setup and plan evaluation. Longitudinal CBCT image registration could quantify the inter-fractional anatomic changes, e.g. tumor shrinkage, and daily OAR variation throughout the course of treatment. The purpose of this study is to propose an unsupervised deep learning-based CBCT-CBCT deformable image registration which enables quantitative anatomic variation analysis. Approach. The proposed deformable registration workflow consists of training and inference stages that share the same feed-forward path through a spatial transformation-based network (STN). The STN consists of a global generative adversarial network (GlobalGAN) and a local GAN (LocalGAN) to predict the coarse- and fine-scale motions, respectively. The network was trained by minimizing the image similarity loss and the deformable vector field (DVF) regularization loss without the supervision of ground truth DVFs. During the inference stage, patches of local DVF were predicted by the trained LocalGAN and fused to form a whole-image DVF. The local whole-image DVF was subsequently combined with the GlobalGAN generated DVF to obtain the final DVF. The proposed method was evaluated using 100 fractional CBCTs from 20 abdominal cancer patients in the experiments and 105 fractional CBCTs from a cohort of 21 different abdominal cancer patients in a holdout test. Main Results. Qualitatively, the registration results show good alignment between the deformed CBCT images and the target CBCT image. Quantitatively, the average target registration error calculated on the fiducial markers and manually identified landmarks was 1.91 ± 1.18 mm. The average mean absolute error, normalized cross correlation between the deformed CBCT and target CBCT were 33.42 ± 7.48 HU, 0.94 ± 0.04, respectively. Significance. In summary, an unsupervised deep learning-based CBCT-CBCT registration method is proposed and its feasibility and performance in fractionated image-guided radiotherapy is investigated. This promising registration method could provide fast and accurate longitudinal CBCT alignment to facilitate inter-fractional anatomic changes analysis and prediction. [ABSTRACT FROM AUTHOR]

Objective. Existing dual-layer flat panel detectors (DL-FPDs) use a thin scintillator layer to preferentially detect low-energy x-rays, followed by a permanent Cu filter to absorb residual low-energy x-rays, and finally, a thicker scintillator layer to preferentially detect high-energy x-rays. The image outputs of the two scintillator layers can be jointly processed for dual-energy (DE) planar and cone-beam CT imaging. In clinical practice, a given FPD is often used for not only DE imaging but also routine single-energy (SE) imaging. With the permanent Cu layer, the total x-ray absorption is unsatisfactory for SE imaging since more than 30% of x-rays can be lost in the Cu layer. The purpose of this work was to demonstrate the benefits of using a removable filter material in DL-FPDs for SE and DE imaging applications. Approach. The proposed detector contains a removable filter between the two scintillator layers. The filter can be either a chamber filled with a liquid high- Z eff material or a removable solid filter. When DE imaging is not clinically indicated, the DL-FPD can switch to a high-efficiency SE imaging mode by retracting the filter from the inter-scintillator space. For commonly available filter materials (iodine, gadolinium, and Cu), their optimal area densities were theoretically calculated for both water-bone decomposition and water-iodine decomposition DE imaging tasks. Preliminary experimental studies were also performed to compare the SE performance of the proposed DL-FPD with the existing DL-FPD with the permanent Cu filter and study the stability of the liquid filter on a rotating gantry. Main results. The optimal filter material was found to be an iodine solution (approximately 375 mg cm−2). With this liquid filter in place, the proposed DL-FPD has equivalent or better DE imaging performance compared with the existing DL-FPD with the Cu filter. When the filter is removed from the inter-scintillator space for SE imaging, the total x-ray absorption efficiency of the proposed DL-FPD ranges from 73% (100 kVp) to 54% (140 kVp), compared with 51% (100 kVp) to 41% (140 kVp) for the existing DL-FPD with a permanent 1 mm Cu filter. Significance. The removable filter provides a boost to the total x-ray absorption efficiency of DL-FPDs for SE imaging without compromising DE imaging. This can facilitate the adoption of DL-FPDs in clinical x-ray imaging systems that usually perform more SE imaging procedures than DE imaging series. [ABSTRACT FROM AUTHOR]

Purpose. Target localization in pulmonary interventions (e.g. transbronchial biopsy of a lung nodule) is challenged by deformable motion and may benefit from fluoroscopic overlay of the target to provide accurate guidance. We present and evaluate a 3D–2D image registration method for fluoroscopic overlay in the presence of tissue deformation using a multi-resolution/multi-scale (MRMS) framework with an objective function that drives registration primarily by soft-tissue image gradients. Methods. The MRMS method registers 3D cone-beam CT to 2D fluoroscopy without gating of respiratory phase by coarse-to-fine resampling and global-to-local rescaling about target regions-of-interest. A variation of the gradient orientation ( GO) similarity metric (denoted G O ′ ) was developed to downweight bone gradients and drive registration via soft-tissue gradients. Performance was evaluated in terms of projection distance error at isocenter (PDEiso). Phantom studies determined nominal algorithm parameters and capture range. Preclinical studies used a freshly deceased, ventilated porcine specimen to evaluate performance in the presence of real tissue deformation and a broad range of 3D–2D image mismatch. Results. Nominal algorithm parameters were identified that provided robust performance over a broad range of motion (0–20 mm), including an adaptive parameter selection technique to accommodate unknown mismatch in respiratory phase. The G O ′ metric yielded median PDEiso = 1.2 mm, compared to 6.2 mm for conventional GO. Preclinical studies with real lung deformation demonstrated median PDEiso = 1.3 mm with MRMS + G O ′ registration, compared to 2.2 mm with a conventional transform. Runtime was 26 s and can be reduced to 2.5 s given a prior registration within ∼5 mm as initialization. Conclusions. MRMS registration via soft-tissue gradients achieved accurate fluoroscopic overlay in the presence of deformable lung motion. By driving registration via soft-tissue image gradients, the method avoided false local minima presented by bones and was robust to a wide range of motion magnitude. [ABSTRACT FROM AUTHOR]

Objective. Periodic respiratory motion and inter-fraction variations are sources of geometric uncertainty in stereotactic body radiation therapy (SBRT) of pulmonary lesions. This study extensively evaluates and validates the separate and combined dosimetric effect of both factors using 4D-CT and daily 4D-cone beam CT (CBCT) dose accumulation scenarios. Approach. A first cohort of twenty early stage or metastatic disease lung cancer patients were retrospectively selected to evaluate each scenario. The planned-dose (3DRef) was optimized on a 3D mid-position CT. To estimate the dosimetric impact of respiratory motion (4DRef), inter-fractional variations (3DAcc) and the combined effect of both factors (4DAcc), three dose accumulation scenarios based on 4D-CT, daily mid-cone beam CT (CBCT) position and 4D-CBCT were implemented via CT-CT/CT-CBCT deformable image registration (DIR) techniques. Each scenario was compared to 3DRef. A separate cohort of ten lung SBRT patients was selected to validate DIR techniques. The distance discordance metric (DDM) was implemented per voxel and per patient for tumor and organs at risk (OARs), and the dosimetric impact for CT-CBCT DIR geometric errors was calculated. Main results. Median and interquartile range (IQR) of the dose difference per voxel were 0.05/2.69 Gy and −0.12/2.68 Gy for 3 D Acc − 3 D Ref and 4 D Acc − 3 D Ref. For 4 D Ref − 3 D Ref the IQR was considerably smaller −0.15/0.78 Gy. These findings were confirmed by dose volume histogram parameters calculated in tumor and OARs. For CT-CT/CT-CBCT DIR validation, DDM (95th percentile) was highest for heart (6.26 mm)/spinal cord (8.00 mm), and below 3 mm for tumor and the rest of OARs. The dosimetric impact of CT-CBCT DIR errors was below 2 Gy for tumor and OARs. Significance. The dosimetric impact of inter-fraction variations were shown to dominate those of periodic respiration in SBRT for pulmonary lesions. Therefore, treatment evaluation and dose-effect studies would benefit more from dose accumulation focusing on day-to-day changes then those that focus on respiratory motion. [ABSTRACT FROM AUTHOR]

Objective. Existing clinical C-arm interventional systems use scintillator-based energy-integrating flat panel detectors (FPDs) to generate cone-beam CT (CBCT) images. Despite its volumetric coverage, FPD-CBCT does not provide sufficient low-contrast detectability desired for certain interventional procedures. The purpose of this work was to develop a C-arm photon counting detector (PCD) CT system with a step-and-shoot data acquisition method to further improve the tomographic imaging performance of interventional systems. Approach. As a proof-of-concept, a cadmium telluride-based 51 cm Ă— 0.6 cm PCD was mounted in front of a FPD in an Artis Zee biplane system. A total of 10 C-arm sweeps (5 forward and 5 backward) were prescribed. A motorized patient table prototype was synchronized with the C-arm system such that it translates the object by a designated distance during the sub-second rest time in between gantry sweeps. To evaluate whether this multi-sweep step-and-shoot acquisition strategy can generate high-quality and volumetric PCD-CT images without geometric distortion artifacts, experiments were performed using physical phantoms, a human cadaver head, and an in vivo swine subject. Comparison with FPD-CT was made under matched narrow beam collimation and radiation dose conditions. Main results. Compared with FPD-CT images, PCD-CT images had lower noise and improved visualization of low-contrast lesion models, as well as improved visibility of small iodinated blood vessels. Fine structures were visualized more clearly by the PCD-CT than the highest-available resolution provided by FPD-CBCT and MDCT. No perceivable geometric distortion artifacts were observed in the multi-planar PCD-CT images. Significance. This work is the first demonstration of the feasibility of high-quality and multi-planar (volumetric) PCD-CT imaging with a rotating C-arm gantry. [ABSTRACT FROM AUTHOR]

Objective. Image guidance and precise irradiation are fundamental to ensure the reliability of small animal oncology studies. Accurate positioning of the animal and the in-beam monitoring of the delivered radio-therapeutic treatment necessitate several imaging modalities. In the particular context of proton therapy with a pulsed beam, information on the delivered dose can be retrieved by monitoring the thermoacoustic waves resulting from the brief and local energy deposition induced by a proton beam (ionoacoustics). The objective of this work was to fabricate a multimodal phantom (x-ray, proton, ultrasound, and ionoacoustics) allowing for sufficient imaging contrast for all the modalities. Approach. The phantom anatomical parts were extracted from mouse computed tomography scans and printed using polylactic acid (organs) and a granite/polylactic acid composite (skeleton). The anatomical pieces were encapsulated in silicone rubber to ensure long term stability. The phantom was imaged using x-ray cone-beam computed tomography, proton radiography, ultrasound imaging, and monitoring of a 20 MeV pulsed proton beam using ionoacoustics. Main results. The anatomical parts could be visualized in all the imaging modalities validating the phantom capability to be used for multimodal imaging. Ultrasound images were simulated from the x-ray cone-beam computed tomography and co-registered with ultrasound images obtained before the phantom irradiation and low-resolution ultrasound images of the mouse phantom in the irradiation position, co-registered with ionoacoustic measurements. The latter confirmed the irradiation of a tumor surrogate for which the reconstructed range was found to be in reasonable agreement with the expectation. Significance. This study reports on a realistic small animal phantom which can be used to investigate ionoacoustic range (or dose) verification together with ultrasound, x-ray, and proton imaging. The co-registration between ionoacoustic reconstructions of the impinging proton beam and x-ray imaging is assessed for the first time in a pre-clinical scenario. [ABSTRACT FROM AUTHOR]

This work focuses on the determination of the radiation doses for a total sample of 916 children, categorised into four age groups (<1, 1â€"5, <5â€"10, <10â€"15 years) undergoing the most frequent paediatric CT scans performed in different scan facilities in Tunisia in order to establish the national diagnostic reference levels (DRLs). Dose evaluation concerned the dosimetric indicators: volume computed tomography dose index (CTDIvol) and doseâ€"length product (DLP). The different paediatric CT protocols and practices were also evaluated. The results show a large variation in doses between different radiology departments. For head scans, the respective DRLs for children aged <1, 1â€"5, 5â€"10 and 10â€"15 years were 26, 38, 51 and 51 mGy, respectively, for CTDIvol and 384, 664, 873 and 978 mGy cm, respectively, for DLP. For the chest, the equivalent respective DRLs were 8, 10, 12 and 15 mGy for CTDIvol and 118, 330, 442 and 526 mGy cm for DLP. For the abdomen, the respective DRLs were 9, 13, 19 and 18 mGy for CTDIvol and 353, 485, 592 and 1073 mGy cm for DLP. This study shows that the optimisation of paediatric CT procedures should be a priority, especially within regional hospitals. The implementation of corrective actions will take place after the initial DRLs. These actions, including recommendations and guidelines to good practice, should be a joint effort of all stakeholders, including health authorities, the radiation protection regulator, professional societies and universities. [ABSTRACT FROM AUTHOR]

Objective. Cone-Beam CT (CBCT) often results in severe image artifacts and inaccurate HU values, meaning poor quality CBCT images cannot be directly applied to dose calculation in radiotherapy. To overcome this, we propose a cycle-residual connection with a dilated convolution-consistent generative adversarial network (Cycle-RCDC-GAN). Approach. The cycle-consistent generative adversarial network (Cycle-GAN) was modified using a dilated convolution with different expansion rates to extract richer semantic features from input images. Thirty pelvic patients were used to investigate the effect of synthetic CT (sCT) from CBCT, and 55 head and neck patients were used to explore the generalizability of the model. Three generalizability experiments were performed and compared: the pelvis trained model was applied to the head and neck; the head and neck trained model was applied to the pelvis, and the two datasets were trained together. Main results. The mean absolute error (MAE), the root mean square error (RMSE), peak signal to noise ratio (PSNR), the structural similarity index (SSIM), and spatial nonuniformity (SNU) assessed the quality of the sCT generated from CBCT. Compared with CBCT images, the MAE improved from 28.81 to 18.48, RMSE from 85.66 to 69.50, SNU from 0.34 to 0.30, and PSNR from 31.61 to 33.07, while SSIM improved from 0.981 to 0.989. The sCT objective indicators of Cycle-RCDC-GAN were better than Cycle-GAN’s. The objective metrics for generalizability were also better than Cycle-GAN’s. Significance. Cycle-RCDC-GAN enhances CBCT image quality and has better generalizability than Cycle-GAN, which further promotes the application of CBCT in radiotherapy. [ABSTRACT FROM AUTHOR]

Objective. We develop a model-based optimization algorithm for ‘one-step’ dual-energy (DE) CT decomposition of three materials directly from projection measurements. Approach. Since the three-material problem is inherently undetermined, we incorporate the volume conservation principle (VCP) as a pair of equality and nonnegativity constraints into the objective function of the recently reported model-based material decomposition (MBMD). An optimization algorithm (constrained MBMD, CMBMD) is derived that utilizes voxel-wise separability to partition the volume into a VCP-constrained region solved using interior-point iterations, and an unconstrained region (air surrounding the object, where VCP is violated) solved with conventional two-material MBMD. Constrained MBMD (CMBMD) is validated in simulations and experiments in application to bone composition measurements in the presence of metal hardware using DE cone-beam CT (CBCT). A kV-switching protocol with non-coinciding low- and high-energy (LE and HE) projections was assumed. CMBMD with decomposed base materials of cortical bone, fat, and metal (titanium, Ti) is compared to MBMD with (i) fat-bone and (ii) fat-Ti bases. Main results. Three-material CMBMD exhibits a substantial reduction in metal artifacts relative to the two-material MBMD implementations. The accuracies of cortical bone volume fraction estimates are markedly improved using CMBMD, with ∼5â€"10Ă— lower normalized root mean squared error in simulations with anthropomorphic knee phantoms (depending on the complexity of the metal component) and ∼2â€"2.5Ă— lower in an experimental test-bench study. Significance. In conclusion, we demonstrated one-step three-material decomposition of DE CT using volume conservation as an optimization constraint. The proposed method might be applicable to DE applications such as bone marrow edema imaging (fat-bone-water decomposition) or multi-contrast imaging, especially on CT/CBCT systems that do not provide coinciding LE and HE ray paths required for conventional projection-domain DE decomposition. [ABSTRACT FROM AUTHOR]

Objective. Modern preclinical small animal radiation platforms utilize cone beam computerized tomography (CBCT) for image guidance and experiment planning purposes. The resolution of CBCT images is of particular importance for visualizing fine animal anatomical structures. One major cause of spatial resolution reduction is the finite size of the x-ray focal spot. In this work, we proposed a simple method to measure x-ray focal spot intensity map and a CBCT image domain deblurring model to mitigate the effect of focal spot-induced image blurring. Approach. We measured a projection image of a tungsten ball bearing using the flat panel detector of the CBCT platform. We built a forward blurring model of the projection image and derived the spot intensity map by deconvolving the measured projection image. Based on the measured spot intensity map, we derived a CBCT image domain blurring model for images reconstructed by the filtered backprojection algorithm. Based on this model, we computed image domain blurring kernel and improved the CBCT image resolution by deconvolving the CBCT image. Main results. We successfully measured the x-ray focal spot intensity map. The spot size characterized by full width at half maximum was ∼0.75 Ă— 0.55 mm2 at 40 kVp. We computed image domain convolution kernels caused by the x-ray focal spot. A simulation study on noiseless projections was performed to evaluate the spatial resolution improvement exclusively by the focal spot kernel, and the modulation transfer function (MTF) at 50% was increased from 1.40 to 1.65 mmâˆ'1 for in-plane images and 1.05â€"1.32 mmâˆ'1 for cross-plane images. Experimental studies on a CT insert phantom and a plastinated mouse phantom demonstrated improved spatial resolution after image domain deconvolution, as indicated by visually improved resolution of fine structures. MTF at 50% was improved from 1.00 to 1.12 mmâˆ'1 for in-plane direction and from 0.72 to 0.84 mmâˆ'1 for cross-plane direction. Significance. The proposed method to mitigate blurring caused by finite x-ray spot size and improve CBCT image resolution is simple and effective. [ABSTRACT FROM AUTHOR]

Objective. 4D-CBCT provides phase-resolved images valuable for radiomics analysis for outcome prediction throughout treatment courses. However, 4D-CBCT suffers from streak artifacts caused by under-sampling, which severely degrades the accuracy of radiomic features. Previously we developed group-patient-trained deep learning methods to enhance the 4D-CBCT quality for radiomics analysis, which was not optimized for individual patients. In this study, a patient-specific model was developed to further improve the accuracy of 4D-CBCT based radiomics analysis for individual patients. Approach. This patient-specific model was trained with intra-patient data. Specifically, patient planning 4D-CT was augmented through image translation, rotation, and deformation to generate 305 CT volumes from 10 volumes to simulate possible patient positions during the onboard image acquisition. 72 projections were simulated from 4D-CT for each phase and were used to reconstruct 4D-CBCT using FDK back-projection algorithm. The patient-specific model was trained using these 305 paired sets of patient-specific 4D-CT and 4D-CBCT data to enhance the 4D-CBCT image to match with 4D-CT images as ground truth. For model testing, 4D-CBCT were simulated from a separate set of 4D-CT scan images acquired from the same patient and were then enhanced by this patient-specific model. Radiomics features were then extracted from the testing 4D-CT, 4D-CBCT, and enhanced 4D-CBCT image sets for comparison. The patient-specific model was tested using 4 lung-SBRT patients’ data and compared with the performance of the group-based model. The impact of model dimensionality, region of interest (ROI) selection, and loss function on the model accuracy was also investigated. Main results. Compared with a group-based model, the patient-specific training model further improved the accuracy of radiomic features, especially for features with large errors in the group-based model. For example, the 3D whole-body and ROI loss-based patient-specific model reduces the errors of the first-order median feature by 83.67%, the wavelet LLL feature maximum by 91.98%, and the wavelet HLL skewness feature by 15.0% on average for the four patients tested. In addition, the patient-specific models with different dimensionality (2D versus 3D) or loss functions (L1 versus L1 + VGG + GAN) achieved comparable results for improving the radiomics accuracy. Using whole-body or whole-body+ROI L1 loss for the model achieved better results than using the ROI L1 loss alone as the loss function. Significance. This study demonstrated that the patient-specific model is more effective than the group-based model on improving the accuracy of the 4D-CBCT radiomic features analysis, which could potentially improve the precision for outcome prediction in radiotherapy. [ABSTRACT FROM AUTHOR]

Objective. Cone-beam CT (CBCT) in modern pre-clinical small-animal radiation research platforms provides volumetric images for image guidance and experiment planning purposes. In this work, we implemented multi-energy element-resolved (MEER) CBCT using three scans with different kVps on a SmART platform (Precision x-ray Inc.) to determine images of relative electron density (rED) and elemental composition (EC) that are needed for Monte Carlo-based radiation dose calculation. Approach. We performed comprehensive calibration tasks to achieve sufficient accuracy for this quantitative imaging purpose. For geometry calibration, we scanned a ball bearing phantom and used an analytical method together with an optimization approach to derive gantry angle specific geometry parameters. Intensity calibration and correction included the corrections for detector lag, glare, and beam hardening. The corrected CBCT projection images acquired at 30, 40, and 60 kVp in multiple scans were used to reconstruct CBCT images using the Feldkampâ€"Davisâ€"Kress reconstruction algorithm. After that, an optimization problem was solved to determine images of rED and EC. We demonstrated the effectiveness of our CBCT calibration steps by showing improvements in image quality and successful material decomposition in cases with a small animal CT calibration phantom and a plastinated mouse phantom. Main results. It was found that artifacts induced by geometry inaccuracy, detector lag, glare, and beam hardening were visually reduced. CT number mean errors were reduced from 19% to 5%. In the CT calibration phantom case, median errors in H, O, and Ca fractions for all the inserts were below 1%, 2%, and 4% respectively, and median error in rED was less than 5%. Compared to the standard approach deriving material type and rED via CT number conversion, our approach improved Monte Carlo simulation-based dose calculation accuracy in bone regions. Mean dose error was reduced from 47.5% to 10.9%. Significance. The MEER-CBCT implemented on an existing CBCT system of a small animal irradiation platform achieved accurate material decomposition and significantly improved Monte Carlo dose calculation accuracy. [ABSTRACT FROM AUTHOR]

Objective. In this study we developed an automatic method to predict tumour volume and shape in weeks 3 and 4 of radiotherapy (RT), using cone-beam computed tomography (CBCT) scans acquired up to week 2, allowing identification of large tumour changes. Approach. 240 non-small cell lung cancer (NSCLC) patients, treated with 55 Gy in 20 fractions, were collected. CBCTs were rigidly registered to the planning CT. Intensity values were extracted in each voxel of the planning target volume across all CBCT images from days 1, 2, 3, 7 and 14. For each patient and in each voxel, four regression models were fitted to voxel intensity; applying linear, Gaussian, quadratic and cubic methods. These models predicted the intensity value for each voxel in weeks 3 and 4, and the tumour volume found by thresholding. Each model was evaluated by computing the root mean square error in pixel value and structural similarity index metric (SSIM) for all patients. Finally, the sensitivity and specificity to predict a 30% change in volume were calculated for each model. Main results. The linear, Gaussian, quadratic and cubic models achieved a comparable similarity score, the average SSIM for all patients was 0.94, 0.94, 0.90, 0.83 in week 3, respectively. At week 3, a sensitivity of 84%, 53%, 90% and 88%, and specificity of 99%, 100%, 91% and 42% were observed for the linear, Gaussian, quadratic and cubic models respectively. Overall, the linear model performed best at predicting those patients that will benefit from RT adaptation. The linear model identified 21% and 23% of patients in our cohort with more than 30% tumour volume reduction to benefit from treatment adaptation in weeks 3 and 4 respectively. Significance. We have shown that it is feasible to predict the shape and volume of NSCLC tumours from routine CBCTs and effectively identify patients who will respond to treatment early. [ABSTRACT FROM AUTHOR]

Proton therapy treatment for lungs remains challenging as images enabling the detection of inter- and intra-fractional motion, which could be used for proton dose adaptation, are not readily available. 4D computed tomography (4DCT) provides high image quality but is rarely available in-room, while in-room 4D cone beam computed tomography (4DCBCT) suffers from image quality limitations stemming mostly from scatter detection. This study investigated the feasibility of using virtual 4D computed tomography (4DvCT) as a prior for a phase-per-phase scatter correction algorithm yielding a 4D scatter corrected cone beam computed tomography image (4DCBCTcor), which can be used for proton dose calculation. 4DCT and 4DCBCT scans of a porcine lung phantom, which generated reproducible ventilation, were acquired with matching breathing patterns. Diffeomorphic Morphons, a deformable image registration algorithm, was used to register the mid-position 4DCT to the mid-position 4DCBCT and yield a 4DvCT. The 4DCBCT was reconstructed using motion-aware reconstruction based on spatial and temporal regularization (MA-ROOSTER). Successively for each phase, digitally reconstructed radiographs of the 4DvCT, simulated without scatter, were exploited to correct scatter in the corresponding CBCT projections. The 4DCBCTcor was then reconstructed with MA-ROOSTER using the corrected CBCT projections and the same settings and deformation vector fields as those already used for reconstructing the 4DCBCT. The 4DCBCTcor and the 4DvCT were evaluated phase-by-phase, performing proton dose calculations and comparison to those of a ground truth 4DCT by means of dose-volume-histograms (DVH) and gamma pass-rates (PR). For accumulated doses, DVH parameters deviated by at most 1.7% in the 4DvCT and 2.0% in the 4DCBCTcor case. The gamma PR for a (2%, 2 mm) criterion with 10% threshold were at least 93.2% (4DvCT) and 94.2% (4DCBCTcor), respectively. The 4DCBCTcor technique enabled accurate proton dose calculation, which indicates the potential for applicability to clinical 4DCBCT scans. [ABSTRACT FROM AUTHOR]

The aim of the present study was to describe patient dose indicator levels during intra-arterial catheter (IAC) implantation for liver chemotherapy, and to determine factors affecting the dose indicators. Between January 2017 and January 2019, 61 patients with hepatic metastases from colorectal cancer were retrospectively included. Interventions were carried out in a standardised manner by three experienced radiologists on the same angiographic table without changes in protocol parameters. For each patient, clinical, radiological and dosimetry data were collected, including the air kerma area product (KAP), part of KAP due to the fluoroscopy and fluoroscopy time (FT), total kerma at the reference interventional point and peak skin dose (PSD). Local dose reference levels (RLs) were determined as the third quartile of the patient dose distributions. Univariate and multivariate analysis of factors affecting dose indicators was performed. The mean KAP was 111 Gy cm2, the mean reference point air kerma (Ka,r) was 648 mGy, the mean PSD was 613 mGy, and the mean FT was 3190 s (62% of the KAP). The mean cone beam computed tomography dose was 37.3 ± 11.8 Gy cm2, which accounted for 37% of the KAP. The RL could be proposed taking into account the third quartiles (KAP = 164.6 Gy cm2, Ka,r = 904.5 mGy, FT = 4011 s and standard deviation = 772.7 mGy). The factors affecting dose indicators were related to the patients (sex, cardiovascular risk factors, weight, body mass index), to the vascular anatomies (coeliac trunk angulation) and to the procedures (number of embolised arteries). This study allowed a better understanding of dose indicators and factors affecting these indicators during the implantation of IACs for hepatic chemotherapy, which is a long and difficult procedure. Local dose RLs were determined. Multicentre, multi-equipment studies are necessary. [ABSTRACT FROM AUTHOR]

Purpose. Although deep learning (DL) technique has been successfully used for computed tomography (CT) reconstruction, its implementation on cone-beam CT (CBCT) reconstruction is extremely challenging due to memory limitations. In this study, a novel DL technique is developed to resolve the memory issue, and its feasibility is demonstrated for CBCT reconstruction from sparsely sampled projection data. Methods. The novel geometry-guided deep learning (GDL) technique is composed of a GDL reconstruction module and a post-processing module. The GDL reconstruction module learns and performs projection-to-image domain transformation by replacing the traditional single fully connected layer with an array of small fully connected layers in the network architecture based on the projection geometry. The DL post-processing module further improves image quality after reconstruction. We demonstrated the feasibility and advantage of the model by comparing ground truth CBCT with CBCT images reconstructed using (1) GDL reconstruction module only, (2) GDL reconstruction module with DL post-processing module, (3) Feldkamp, Davis, and Kress (FDK) only, (4) FDK with DL post-processing module, (5) ray-tracing only, and (6) ray-tracing with DL post-processing module. The differences are quantified by peak-signal-to-noise ratio (PSNR), structural similarity index measure (SSIM), and root-mean-square error (RMSE). Results. CBCT images reconstructed with GDL show improvements in quantitative scores of PSNR, SSIM, and RMSE. Reconstruction time per image for all reconstruction methods are comparable. Compared to current DL methods using large fully connected layers, the estimated memory requirement using GDL is four orders of magnitude less, making DL CBCT reconstruction feasible. Conclusion. With much lower memory requirement compared to other existing networks, the GDL technique is demonstrated to be the first DL technique that can rapidly and accurately reconstruct CBCT images from sparsely sampled data. [ABSTRACT FROM AUTHOR]

Active matrix, flat-panel imagers (AMFPIs) suffer from decreased detective quantum efficiency under conditions of low dose per image frame (such as for digital breast tomosynthesis, fluoroscopy and cone-beam CT) due to low signal compared to the additive electronic noise. One way to address this challenge is to introduce a high-gain x-ray converter called particle-in-binder mercuric iodide (PIB HgI2) which exhibits 3–10 times higher x-ray sensitivity compared to that of a-Se and CsI:Tl converters employed in commercial AMFPI systems. However, a remaining challenge for practical implementation of PIB HgI2 is the high level of image lag, which is believed to largely originate from the trapping of holes. Towards addressing this challenge, this paper reports a theoretical investigation of the use of a Frisch grid structure embedded in the converter to suppress hole signal—which would be expected to reduce image lag. The grid acts as a third electrode sandwiched between a continuous top electrode and pixelated bottom electrodes having a 100 μm pitch. Signal properties of such a detector are investigated as a function of VDR (the ratio of the voltage difference between the electrodes in the region below the grid to that above the grid), grid pitch (the center-to-center distance between two neighboring grid wires) and RGRID (the ratio of grid wire width to grid pitch) for mammographic x-ray energies. The results show that smaller grid pitch suppresses hole signal to a higher degree (up to ∼96%) while a larger gap between grid wires and higher VDR provide minimally impeded electron transport. Examination of the tradeoff between maximizing electron signal and minimizing hole signal indicates that a grid design having a grid pitch of 20 μm with RGRID of 50% and 65% provides hole signal suppression of ∼93% and ∼95% for VDR of 1 and 3, respectively. [ABSTRACT FROM AUTHOR]

High cone-angle artifacts (HCAAs) appear frequently in circular cone-beam computed tomography (CBCT) images and can heavily affect diagnosis and treatment planning. To reduce HCAAs in CBCT scans, we propose a novel deep learning approach that reduces the three-dimensional (3D) nature of HCAAs to two-dimensional (2D) problems in an efficient way. Specifically, we exploit the relationship between HCAAs and the rotational scanning geometry by training a convolutional neural network (CNN) using image slices that were radially sampled from CBCT scans. We evaluated this novel approach using a dataset of input CBCT scans affected by HCAAs and high-quality artifact-free target CBCT scans. Two different CNN architectures were employed, namely U-Net and a mixed-scale dense CNN (MS-D Net). The artifact reduction performance of the proposed approach was compared to that of a Cartesian slice-based artifact reduction deep learning approach in which a CNN was trained to remove the HCAAs from Cartesian slices. In addition, all processed CBCT scans were segmented to investigate the impact of HCAAs reduction on the quality of CBCT image segmentation. We demonstrate that the proposed deep learning approach with geometry-aware dimension reduction greatly reduces HCAAs in CBCT scans and outperforms the Cartesian slice-based deep learning approach. Moreover, the proposed artifact reduction approach markedly improves the accuracy of the subsequent segmentation task compared to the Cartesian slice-based workflow. [ABSTRACT FROM AUTHOR]

This study aimed to investigate if a commercial, knowledge-based tool for radiotherapy planning could be used to estimate the amount of sparing in organs at risk (OARs) in the re-planning strategy for adaptive radiotherapy (ART). Eighty head and neck (HN) VMAT Pareto plans from our institute's database were used to train a knowledge-based planning (KBP) model. An evaluation set of another 20 HN patients was randomly selected. For each patient in the evaluation set, the planning computed tomography (CT) and 2 sets of on-board cone-beam CT, corresponding to the middle and second half of the radiotherapy treatment course, were extracted. The original plan was re-calculated on a daily deformed CT (delivered dose‐volume histogram (DVH)) and compared with the KBP DVH predictions and with the final KBP DVH after optimisation of the plan, which was performed on the same image sets. To evaluate the feasibility of this method, the range of KBP DVH uncertainties was compared with the gains obtained from re-planning. DVH differences and receiver operating characteristic (ROC) curve analysis were used for this purpose. On average, final KBP uncertainties were smaller than the gain in re-planning. Statistical tests confirmed significant differences between the two groups. ROC analysis showed KBP performance in terms of area under the curve values higher than 0.7, which confirmed a good accuracy in predicted values. Overall, for 48% of cases, KBP predicted a desirable outcome from re-planning, and the final dose confirmed an effective gain in 47% of cases. We have established a systematic workflow to identify effective OAR sparing in re-planning based on KBP predictions that can be implemented in an on-line, ART process. [ABSTRACT FROM AUTHOR]

Simultaneous acquisition of cone beam CT (CBCT) projections using both the kV and MV imagers of an image guided radiotherapy system reduces set-up scan times—a benefit to lung cancer radiation oncology patients—but increases noise in the 3D reconstruction. In this article, we present a kV‐MV scan time reduction technique that uses two noise-reducing measures to achieve superior performance. The first is a high-DQE multi-layer MV imager prototype. The second is a beam hardening correction algorithm which combines poly-energetic modeling with edge-preserving, regularized smoothing of the projections. Performance was tested in real acquisitions of the Catphan 604 and a thorax phantom. Percent noise was quantified from voxel values in a soft tissue volume of interest (VOI) while edge blur was quantified from a VOI straddling a boundary between air and soft material. Comparisons in noise/resolution performance trade-off were made between our proposed approach, a dose-equivalent kV-only scan, and a kV‐MV reconstruction technique previously published by Yin et al (2005 Med. Phys. 32 9). The proposed technique demonstrated lower noise as a function of spatial resolution than the baseline kV‐MV method, notably a 50% noise reduction at typical edge blur levels. Our proposed method also exhibited fainter non-uniformity artifacts and in some cases superior contrast. Overall, we find that the combination of a multi-layer MV imager, acquiring at a LINAC source energy of 2.5 MV, and a denoised beam hardening correction algorithm enables noise, resolution, and dose performance comparable to standard kV-imager only set-up CBCT, but with nearly half the gantry rotation time. [ABSTRACT FROM AUTHOR]

Purpose. Accurate localization and labeling of vertebrae in computed tomography (CT) is an important step toward more quantitative, automated diagnostic analysis and surgical planning. In this paper, we present a framework (called Ortho2D) for vertebral labeling in CT in a manner that is accurate and memory-efficient. Methods. Ortho2D uses two independent faster R-convolutional neural network networks to detect and classify vertebrae in orthogonal (sagittal and coronal) CT slices. The 2D detections are clustered in 3D to localize vertebrae centroids in the volumetric CT and classify the region (cervical, thoracic, lumbar, or sacral) and vertebral level. A post-process sorting method incorporates the confidence in network output to refine classifications and reduce outliers. Ortho2D was evaluated on a publicly available dataset containing 302 normal and pathological spine CT images with and without surgical instrumentation. Labeling accuracy and memory requirements were assessed in comparison to other recently reported methods. The memory efficiency of Ortho2D permitted extension to high-resolution CT to investigate the potential for further boosts to labeling performance. Results. Ortho2D achieved overall vertebrae detection accuracy of 97.1%, region identification accuracy of 94.3%, and individual vertebral level identification accuracy of 91.0%. The framework achieved 95.8% and 83.6% level identification accuracy in images without and with surgical instrumentation, respectively. Ortho2D met or exceeded the performance of previously reported 2D and 3D labeling methods and reduced memory consumption by a factor of ∼50 (at 1 mm voxel size) compared to a 3D U-Net, allowing extension to higher resolution datasets than normally afforded. The accuracy of level identification increased from 80.1% (for standard/low resolution CT) to 95.1% (for high-resolution CT). Conclusions. The Ortho2D method achieved vertebrae labeling performance that is comparable to other recently reported methods with significant reduction in memory consumption, permitting further performance boosts via application to high-resolution CT. [ABSTRACT FROM AUTHOR]

Patient movement affects image quality in oral and maxillofacial cone-beam computed tomography imaging. While many efforts are made to minimize the possibility of motion during a scan, relatively little attention has been given to motion correction after acquisition. We propose a novel method which can improve the image quality after an oral and maxillofacial scan. The proposed method is based on our previous work and is a retrospective motion estimation and motion compensation (ME/MC) approach that iteratively estimates and compensates for rigid pose change over time. During motion estimation, image update and motion update are performed alternately in a multi-resolution scheme to obtain the motion. We propose use of a feature-based motion update and patch-based image update in the iterative estimation process, to alleviate the effect of limited scan field of view on estimation. During motion compensation, a fine-resolution image reconstruction was performed with compensation for the estimated motion. The proposed ME/MC method was evaluated with simulations, phantom and patient studies. Two experts in dentomaxillofacial radiology assessed the diagnostic importance of the resulting motion artifact suppression. The quality of the reconstructed images was improved after motion compensation, and most of the image artifacts were eliminated. Quantitative analysis by comparison to a reference image and by calculation of a sharpness metric agreed with the qualitative observation. The results are promising, and further evaluation is required to assess the clinical value of the proposed method. [ABSTRACT FROM AUTHOR]

Adaptive-radiation-therapy (ART) is applied to account for anatomical variations observed over the treatment course. Daily or weekly cone-beam computed tomography (CBCT) is commonly used in clinic for patient positioning, but CBCT's inaccuracy in Hounsfield units (HU) prevents its application to dose calculation and treatment planning. Adaptive re-planning can be performed by deformably registering planning CT (pCT) to CBCT. However, scattering artifacts and noise in CBCT decrease the accuracy of deformable registration and induce uncertainty in treatment plan. Hence, generating from CBCT a synthetic CT (sCT) that has the same anatomical structure as CBCT but accurate HU values is desirable for ART. We proposed an unsupervised style-transfer-based approach to generate sCT based on CBCT and pCT. Unsupervised learning was desired because exactly matched CBCT and CT are rarely available, even when they are taken a few minutes apart. In the proposed model, CBCT and pCT are two inputs that provide anatomical structure and accurate HU information, respectively. The training objective function is designed to simultaneously minimize (1) contextual loss between sCT and CBCT to maintain the content and structure of CBCT in sCT and (2) style loss between sCT and pCT to achieve pCT-like image quality in sCT. We used CBCT and pCT images of 114 patients to train and validate the designed model, and another 29 independent patient cases to test the model's effectiveness. We quantitatively compared the resulting sCT with the original CBCT using the deformed same-day pCT as reference. Structure-similarity-index, peak-signal-to-noise-ratio, and mean-absolute-error in HU of sCT were 0.9723, 33.68, and 28.52, respectively, while those of CBCT were 0.9182, 29.67, and 49.90, respectively. We have demonstrated the effectiveness of the proposed model in using CBCT and pCT to synthesize CT-quality images. This model may permit using CBCT for advanced applications such as adaptive treatment planning. [ABSTRACT FROM AUTHOR]

Fixed-gantry radiation therapy has been proposed as a low-cost alternative to the conventional rotating-gantry radiation therapy, that may help meet the rising global treatment demand. Fixed-gantry systems require gravitational motion compensated reconstruction algorithms to produce cone-beam CT (CBCT) images of sufficient quality for image guidance. The aim of this work was to adapt and investigate five CBCT reconstruction algorithms for fixed-gantry CBCT images. The five algorithms investigated were Feldkamp–Davis–Kress (FDK), prior image constrained compressed sensing (PICCS), gravitational motion compensated FDK (GMCFDK), motion compensated PICCS (MCPICCS) (a novel CBCT reconstruction algorithm) and simultaneous motion estimation and iterative reconstruction (SMEIR). Fixed-gantry and rotating-gantry CBCT scans were acquired of 3 rabbits, with the rotating-gantry scans used as a reference. Projections were sorted into rotation bins, based on the angle of rotation of the rabbit during image acquisition. The algorithms were compared using the structural similarity index measure root mean square error, and reconstruction time. Evaluation of the reconstructed volumes showed that, when compared with the reference rotating-gantry volume, the conventional FDK algorithm did not accurately reconstruct fixed-gantry CBCT scans. Whilst the PICCS reconstruction algorithm reduced some motion artefacts, the motion estimation reconstruction methods (GMCFDK, MCPICCS and SMEIR) were able to greatly reduce the effect of motion artefacts on the reconstructed volumes. This finding was verified quantitatively, with GMCFDK, MCPICCS and SMEIR reconstructions having RMSE 17%–19% lower and SSIM 1% higher than a conventional FDK. However, all motion compensated fixed-gantry CBCT reconstructions had a 56%–61% higher RMSE and 1.5% lower SSIM than FDK reconstructions of conventional rotating-gantry CBCT scans. The results show that motion compensation is required to reduce motion artefacts for fixed-gantry CBCT reconstructions. This paper further demonstrates the feasibility of fixed-gantry CBCT scans, and the ability of CBCT reconstruction algorithms to compensate for motion due to horizontal rotation. [ABSTRACT FROM AUTHOR]

X-ray scatter remains a major physics challenge in volumetric computed tomography (CT), whose physical and statistical behaviors have been commonly leveraged in order to eliminate its impact on CT image quality. In this work, we conduct an in-depth derivation of how the scatter distribution and scatter to primary ratio (SPR) will change during the spectral correction, leading to an interesting finding on the property of scatter. Such a characterization of scatter's behavior provides an analytic approach of compensating for the SPR as well as approximating the change of scatter distribution after spectral correction, even though both of them might be significantly distorted as the linearization mapping function in spectral correction could vary a lot from one detector pixel to another. We conduct an evaluation of SPR compensations (SPRCs) on a Catphan phantom and an anthropomorphic chest phantom to validate the characteristics of scatter. In addition, this scatter property is also directly adopted into CT imaging using a spectral modulator with flying focal spot technology (SMFFS) as an example to demonstrate its potential in practical applications. For cone-beam CT (CBCT) scans at both 80 and 120 kVp, CT images with accurate CT numbers can be achieved after spectral correction followed by the appropriate SPRC based on our presented scatter property. In the case of the SMFFS based CBCT scan of the Catphan phantom at 120 kVp, after a scatter correction using an analytic algorithm derived from the scatter property, CT image quality was significantly improved, with the averaged root mean square error reduced from 297.9 to 6.5 Hounsfield units. [ABSTRACT FROM AUTHOR]

Acquiring CBCTs from a limited scan angle can help to reduce the imaging time, save the imaging dose, and allow continuous target localizations through arc-based treatments with high temporal resolution. However, insufficient scan angle sampling leads to severe distortions and artifacts in the reconstructed CBCT images, limiting their clinical applicability. 2D–3D deformable registration can map a prior fully-sampled CT/CBCT volume to estimate a new CBCT, based on limited-angle on-board cone-beam projections. The resulting CBCT images estimated by 2D–3D deformable registration can successfully suppress the distortions and artifacts, and reflect up-to-date patient anatomy. However, traditional iterative 2D–3D deformable registration algorithm is very computationally expensive and time-consuming, which takes hours to generate a high quality deformation vector field (DVF) and the CBCT. In this work, we developed an unsupervised, end-to-end, 2D–3D deformable registration framework using convolutional neural networks (2D3D-RegNet) to address the speed bottleneck of the conventional iterative 2D–3D deformable registration algorithm. The 2D3D-RegNet was able to solve the DVFs within 5 seconds for 90 orthogonally-arranged projections covering a combined 90° scan angle, with DVF accuracy superior to 3D–3D deformable registration, and on par with the conventional 2D–3D deformable registration algorithm. We also performed a preliminary robustness analysis of 2D3D-RegNet towards projection angular sampling frequency variations, as well as scan angle offsets. The synergy of 2D3D-RegNet with biomechanical modeling was also evaluated, and demonstrated that 2D3D-RegNet can function as a fast DVF solution core for further DVF refinement. [ABSTRACT FROM AUTHOR]

This study introduces and evaluates respiratory-correlated four-dimensional (4D) inverse geometry computed tomography (IGCT). The projection data of the IGCT were acquired in a single gantry rotation over 120 s. Three virtual phantoms—static Defrise, 4D Shepp–Logan, and 4D extended cardiac-torso (XCAT)—were used to obtain projection data for the IGCT and cone-beam computed tomography (CBCT). The projection acquisition parameters were determined to eliminate vacancies in the Radon space for an accurate rebinning process. Phase-based sorting was conducted within 10 phase bins, and the sorted projection data were binned into a cone beam geometry. Finally, Feldkamp–Davis–Kress reconstruction was conducted independently at each phase. The reconstructed images were compared using the structural similarity index measure (SSIM) and root mean square error (RMSE). The vertical profile of the Defrise reconstruction image was uniform, and the cone beam artefact was reduced in the IGCT image. Under an ideal projection acquisition condition, the mean coronal plane SSIMs of the Shepp–Logan and 4D XCAT phantoms were 0.899 and 0.706, respectively, which were higher than those of the CBCT (0.784 and 0.623, respectively). Similarly, the mean RMSEs of the coronal plane IGCT (0.036 and 0.158) exhibited an improvement over those of the CBCT (0.165 and 0.261, respectively). The mean standard deviations of the SSIM and RMSE were lower for IGCT than for CBCT. In particular, the SSIM and RMSE of the sagittal and coronal planes of the Shepp–Logan IGCT images were stable in all phase bins; however, those of the CBCT changed depending on the phase bins. Poor image quality was observed for IGCT under inappropriate conditions. This was caused by a vacancy in the Radon space, owing to an inappropriate scan setting. Overall, the proposed 4D IGCT exhibited better image quality than conventional CBCT. [ABSTRACT FROM AUTHOR]

Beam hardening in x-ray computed tomography (CT) is inevitable because of the polychromatic x-ray spectrum and energy-dependent attenuation coefficients of materials, leading to the underestimation of artifacts arising from projection data, especially on metal regions. State-of-the-art research on beam-hardening artifacts is based on a numerical method that recursively performs CT reconstruction, which leads to a heavy computational burden. To address this computational issue, we propose a constrained beam-hardening estimator that provides an efficient numerical solution via a linear combination of two images reconstructed only once during the entire process. The proposed estimator reflects the geometry of metal objects and physical characteristics of beam hardening during the transmission of polychromatic x-rays through a material. Most of the associated parameters are numerically obtained from an initial uncorrected CT image and forward projection transformation without additional optimization procedures. Only the unknown parameter related to beam-hardening artifacts is fine-tuned by linear optimization, which is performed only in the reconstruction image domain. The proposed approach was systematically assessed using numerical simulations and phantom data for qualitative and quantitative comparisons. Compared with existing sinogram inpainting-based and model-based approaches, the proposed scheme in conjunction with the constrained beam-hardening estimator not only provided improved image quality in areas surrounding the metal but also achieved fast beam-hardening correction owing to the analytical reconstruction structure. This work may have significant implications in improving dose calculation accuracy or target volume delineation for treatment planning in radiotherapy. [ABSTRACT FROM AUTHOR]

Automated segmentation of the esophagus is critical in image-guided/adaptive radiotherapy of lung cancer to minimize radiation-induced toxicities such as acute esophagitis. We have developed a semantic physics-based data augmentation method for segmenting the esophagus in both planning CT (pCT) and cone beam CT (CBCT) using 3D convolutional neural networks. One hundred and ninety-one cases with their pCTs and CBCTs from four independent datasets were used to train a modified 3D U-Net architecture and a multi-objective loss function specifically designed for soft-tissue organs such as the esophagus. Scatter artifacts and noises were extracted from week-1 CBCTs using a power-law adaptive histogram equalization method and induced to the corresponding pCT were reconstructed using CBCT reconstruction parameters. Moreover, we leveraged physics-based artifact induction in pCTs to drive the esophagus segmentation in real weekly CBCTs. Segmentations were evaluated using the geometric Dice coefficient and Hausdorff distance as well as dosimetrically using mean esophagus dose and D5cc. Due to the physics-based data augmentation, our model trained just on the synthetic CBCTs was robust and generalizable enough to also produce state-of-the-art results on the pCTs and CBCTs, achieving Dice overlaps of 0.81 and 0.74, respectively. It is concluded that our physics-based data augmentation spans the realistic noise/artifact spectrum across patient CBCT/pCT data and can generalize well across modalities, eventually improving the accuracy of treatment setup and response analysis. [ABSTRACT FROM AUTHOR]

Purpose. A radiotherapy system with a fixed treatment beam and a rotating patient positioning system could be smaller, more robust and more cost effective compared to conventional rotating gantry systems. However, patient rotation could cause anatomical deformation and compromise treatment delivery. In this work, we demonstrate an image-guided treatment workflow with a fixed beam prototype system that accounts for deformation during rotation to maintain dosimetric accuracy. Methods. The prototype system consists of an Elekta Synergy linac with the therapy beam orientated downward and a custom-built patient rotation system (PRS). A phantom that deforms with rotation was constructed and rotated within the PRS to quantify the performance of two image guidance techniques: motion compensated cone-beam CT (CBCT) for pre-treatment volumetric imaging and kilovoltage infraction monitoring (KIM) for real-time image guidance. The phantom was irradiated with a 3D conformal beam to evaluate the dosimetric accuracy of the workflow. Results. The motion compensated CBCT was used to verify pre-treatment position and the average calculated position was within −0.3 ± 1.1 mm of the phantom's ground truth position at 0°. KIM tracked the position of the target in real-time as the phantom was rotated and the average calculated position was within −0.2 ± 0.8 mm of the phantom's ground truth position. A 3D conformal treatment delivered on the prototype system with image guidance had a 3%/2 mm gamma pass rate of 96.3% compared to 98.6% delivered using a conventional rotating gantry linac. Conclusions. In this work, we have shown that image guidance can be used with fixed-beam treatment systems to measure and account for changes in target position in order to maintain dosimetric coverage during horizontal rotation. This treatment modality could provide a viable treatment option when there insufficient space for a conventional linear accelerator or where the cost is prohibitive. [ABSTRACT FROM AUTHOR]

The high conformality of intensity-modulated proton therapy (IMPT) dose distributions causes treatment plans to be sensitive to geometrical changes during the course of a fractionated treatment. This can be addressed using adaptive proton therapy (APT). One important question in APT is the frequency of adaptations performed during a fractionated treatment, which is related to the question whether plan adaptation has to be done online or offline. The purpose of this work is to investigate the impact of weekly and daily online IMPT plan adaptation on the treatment quality for head and neck patients. A cohort of ten head and neck patients with daily acquired cone-beam CT (CBCT) images was evaluated retrospectively. Dose tracking of the IMPT treatment was performed for three scenarios: base plan with no adaptation (BP), weekly online adaptation (OAW), and daily online adaptation (OAD). Both adaptation schemes used an in-house developed online APT workflow, performing Monte Carlo dose calculations on scatter-corrected CBCTs. IMPT plan adaptation was achieved by only tuning the weights of a subset of beamlets, based on deformable image registration from the planning CT to each CBCT. Although OAD mitigated random delivery errors more effectively than OAW on a fraction per fraction basis, both OAW and OAD achieved the clinical goals for all ten patients, while BP failed for six cases. In the high-risk CTV, accumulated values of D98% ranged between 97.15% and 99.73% of the prescription dose for OAD, with a median of 98.07%. For OAW, values between 95.02% and 99.26% were obtained, with a median of 97.61% of the prescription dose. Otherwise, the dose to most organs at risk was similar for all three scenarios. Globally, our results suggest that OAW could be used as an alternative approach to OAD for most patients in order to reduce the clinical workload. [ABSTRACT FROM AUTHOR]

Model-based iterative reconstruction (MBIR) for cone-beam CT (CBCT) offers better noise-resolution tradeoff and image quality than analytical methods for acquisition protocols with low x-ray dose or limited data, but with increased computational burden that poses a drawback to routine application in clinical scenarios. This work develops a comprehensive framework for acceleration of MBIR in the form of penalized weighted least squares optimized with ordered subsets separable quadratic surrogates. The optimization was scheduled on a set of stages forming a morphological pyramid varying in voxel size. Transition between stages was controlled with a convergence criterion based on the deviation between the mid-band noise power spectrum (NPS) measured on a homogeneous region of the evolving reconstruction and that expected for the converged image, computed with an analytical model that used projection data and the reconstruction parameters. A stochastic backprojector was developed by introducing a random perturbation to the sampling position of each voxel for each ray in the reconstruction within a voxel-based backprojector, breaking the deterministic pattern of sampling artifacts when combined with an unmatched Siddon forward projector. This fast, forward and backprojector pair were included into a multi-resolution reconstruction strategy to provide support for objects partially outside of the field of view. Acceleration from ordered subsets was combined with momentum accumulation stabilized with an adaptive technique that automatically resets the accumulated momentum when it diverges noticeably from the current iteration update. The framework was evaluated with CBCT data of a realistic abdomen phantom acquired on an imaging x-ray bench and with clinical CBCT data from an angiography robotic C-arm (Artis Zeego, Siemens Healthineers, Forchheim, Germany) acquired during a liver embolization procedure. Image fidelity was assessed with the structural similarity index (SSIM) computed with a converged reconstruction. The accelerated framework provided accurate reconstructions in 60 s (SSIM = 0.97) and as little as 27 s (SSIM = 0.94) for soft-tissue evaluation. The use of simple forward and backprojectors resulted in 9.3× acceleration. Accumulation of momentum provided extra ∼3.5× acceleration with stable convergence for 6–30 subsets. The NPS-driven morphological pyramid resulted in initial faster convergence, achieving similar SSIM with 1.5× lower runtime than the single-stage optimization. Acceleration of MBIR to provide reconstruction time compatible with clinical applications is feasible via architectures that integrate algorithmic acceleration with approaches to provide stable convergence, and optimization schedules that maximize convergence speed. [ABSTRACT FROM AUTHOR]

Image-guided therapies in the abdomen and pelvis are often hindered by motion artifacts in cone-beam CT (CBCT) arising from complex, non-periodic, deformable organ motion during long scan times (5–30 s). We propose a deformable image-based motion compensation method to address these challenges and improve CBCT guidance. Motion compensation is achieved by selecting a set of small regions of interest in the uncompensated image to minimize a cost function consisting of an autofocus objective and spatiotemporal regularization penalties. Motion trajectories are estimated using an iterative optimization algorithm (CMA-ES) and used to interpolate a 4D spatiotemporal motion vector field. The motion-compensated image is reconstructed using a modified filtered backprojection approach. Being image-based, the method does not require additional input besides the raw CBCT projection data and system geometry that are used for image reconstruction. Experimental studies investigated: (1) various autofocus objective functions, analyzed using a digital phantom with a range of sinusoidal motion magnitude (4, 8, 12, 16, 20 mm); (2) spatiotemporal regularization, studied using a CT dataset from The Cancer Imaging Archive with deformable sinusoidal motion of variable magnitude (10, 15, 20, 25 mm); and (3) performance in complex anatomy, evaluated in cadavers undergoing simple and complex motion imaged on a CBCT-capable mobile C-arm system (Cios Spin 3D, Siemens Healthineers, Forchheim, Germany). Gradient entropy was found to be the best autofocus objective for soft-tissue CBCT, increasing structural similarity (SSIM) by 42%–92% over the range of motion magnitudes investigated. The optimal temporal regularization strength was found to vary widely (0.5–5 mm−2) over the range of motion magnitudes investigated, whereas optimal spatial regularization strength was relatively constant (0.1). In cadaver studies, deformable motion compensation was shown to improve local SSIM by ∼17% for simple motion and ∼21% for complex motion and provided strong visual improvement of motion artifacts (reduction of blurring and streaks and improved visibility of soft-tissue edges). The studies demonstrate the robustness of deformable motion compensation to a range of motion magnitudes, frequencies, and other factors (e.g. truncation and scatter). [ABSTRACT FROM AUTHOR]

Organ-at-risk (OAR) delineation is a key step for cone-beam CT (CBCT) based adaptive radiotherapy planning that can be a time-consuming, labor-intensive, and subject-to-variability process. We aim to develop a fully automated approach aided by synthetic MRI for rapid and accurate CBCT multi-organ contouring in head-and-neck (HN) cancer patients. MRI has superb soft-tissue contrasts, while CBCT offers bony-structure contrasts. Using the complementary information provided by MRI and CBCT is expected to enable accurate multi-organ segmentation in HN cancer patients. In our proposed method, MR images are firstly synthesized using a pre-trained cycle-consistent generative adversarial network given CBCT. The features of CBCT and synthetic MRI (sMRI) are then extracted using dual pyramid networks for final delineation of organs. CBCT images and their corresponding manual contours were used as pairs to train and test the proposed model. Quantitative metrics including Dice similarity coefficient (DSC), Hausdorff distance 95% (HD95), mean surface distance, and residual mean square distance (RMS) were used to evaluate the proposed method. The proposed method was evaluated on a cohort of 65 HN cancer patients. CBCT images were collected from those patients who received proton therapy. Overall, DSC values of 0.87 ± 0.03, 0.79 ± 0.10/0.79 ± 0.11, 0.89 ± 0.08/0.89 ± 0.07, 0.90 ± 0.08, 0.75 ± 0.06/0.77 ± 0.06, 0.86 ± 0.13, 0.66 ± 0.14, 0.78 ± 0.05/0.77 ± 0.04, 0.96 ± 0.04, 0.89 ± 0.04/0.89 ± 0.04, 0.83 ± 0.02, and 0.84 ± 0.07 for commonly used OARs for treatment planning including brain stem, left/right cochlea, left/right eye, larynx, left/right lens, mandible, optic chiasm, left/right optic nerve, oral cavity, left/right parotid, pharynx, and spinal cord, respectively, were achieved. This study provides a rapid and accurate OAR auto-delineation approach, which can be used for adaptive radiation therapy. [ABSTRACT FROM AUTHOR]

Cone-beam computed tomography (CBCT) images acquired during volumetric modulated arc therapy (VMAT; ii-CBCT) can be used to calculate actual delivered doses (ADDs). However, such ii-CBCT images are degraded by scattered megavoltage x-rays (MV-scatters). We aimed to evaluate the dose calculation accuracy of the MV-scatter uncorrected or corrected ii-CBCT images acquired during VMAT deliveries. For MV-scatter correction on concurrent kilovoltage projections (PMVkV), projections consisting only of MV-scatters (PMVS) were acquired under the same MV beam parameters and gantry angles and subtracted from PMVkV (PMVScorr). In addition, the projections by kilovoltage beams were acquired for reference (PkV). The corresponding CBCT images were reconstructed using the Feldkamp–Davis–Kress algorithm (CBCTMVkV, CBCTMVScorr, and CBCTkV as reference). A multi-energy phantom with rods of various relative electron densities (REDs) was used to generate a CBCT-number-to-RED conversion table. First, CBCTkV was reconstructed. Then, the mean CBCT-numbers within each rod were extracted, and a reference table was generated. Concurrent kilovoltage imaging with various field sizes was also demonstrated, and CBCTMVkV and CBCTMVScorr were reconstructed. The extracted CBCT-numbers of each ii-CBCT image were converted into REDs using the reference table. Next, the absolute differences of RED between the ii-CBCT image and CBCTkV were calculated. Ten VMAT plans using a 10 MV flattening-filter-free beam were used for concurrent imaging of an anthropomorphic torso phantom. Moreover, an iterative reconstruction algorithm (IRA) was used for CBCTMVScorr. The plans were recalculated for the corresponding CBCTMVkV, CBCTMVScorr, CBCTMVScorr+IRA, and CBCTkV with the reference table. Finally, the doses were evaluated using 3D gamma analysis (1%/1 mm). The median difference ranges between CBCTMVkV/CBCTMVScorr and the reference values were −0.58 to −0.10/−0.03 to 0.00. The median gamma pass rates of the doses on CBCTMVkV, CBCTMVScorr, and CBCTMVScorr+IRA to the rate on CBCTkV were 70.4, 99.5, and 98.2%, respectively. CBCTMVScorr were comparable with CBCTkV for calculating the ADD from VMAT. [ABSTRACT FROM AUTHOR]

Monte Carlo simulation (MCS) is one of the most accurate computation methods for dose calculation and image formation in radiation therapy. However, the high computational complexity and long execution time of MCS limits its broad use. In this paper, we present a novel strategy to accelerate MCS using a graphic processing unit (GPU), and we demonstrate the application in mega-voltage (MV) cone-beam computed tomography (CBCT) simulation. A new framework that generates a series of MV projections from a single simulation run is designed specifically for MV-CBCT acquisition. A Geant4-based GPU code for photon simulation is incorporated into the framework for the simulation of photon transport through a phantom volume. The FastEPID method, which accelerates the simulation of MV images, is modified and integrated into the framework. The proposed GPU-based simulation strategy was tested for its accuracy and efficiency in a Catphan 604 phantom and an anthropomorphic pelvis phantom with beam energies at 2.5 MV, 6 MV, and 6 MV FFF. In all cases, the proposed GPU-based simulation demonstrated great simulation accuracy and excellent agreement with measurement and CPU-based simulation in terms of reconstructed image qualities. The MV-CBCT simulation was accelerated by factors of roughly 900–2300 using an NVIDIA Tesla V100 GPU card against a 2.5 GHz AMD Opteron™ Processor 6380. [ABSTRACT FROM AUTHOR]

Cone-beam computed tomography (CBCT)- and magnetic resonance (MR)-images allow a daily observation of patient anatomy but are not directly suited for accurate proton dose calculations. This can be overcome by creating synthetic CTs (sCT) using deep convolutional neural networks. In this study, we compared sCTs based on CBCTs and MRs for head and neck (H&N) cancer patients in terms of image quality and proton dose calculation accuracy. A dataset of 27 H&N-patients, treated with proton therapy (PT), containing planning CTs (pCTs), repeat CTs, CBCTs and MRs were used to train two neural networks to convert either CBCTs or MRs into sCTs. Image quality was quantified by calculating mean absolute error (MAE), mean error (ME) and Dice similarity coefficient (DSC) for bones. The dose evaluation consisted of a systematic non-clinical analysis and a clinical recalculation of actually used proton treatment plans. Gamma analysis was performed for non-clinical and clinical treatment plans. For clinical treatment plans also dose to targets and organs at risk (OARs) and normal tissue complication probabilities (NTCP) were compared. CBCT-based sCTs resulted in higher image quality with an average MAE of 40 ± 4 HU and a DSC of 0.95, while for MR-based sCTs a MAE of 65 ± 4 HU and a DSC of 0.89 was observed. Also in clinical proton dose calculations, sCTCBCT achieved higher average gamma pass ratios (2%/2 mm criteria) than sCTMR (96.1% vs. 93.3%). Dose-volume histograms for selected OARs and NTCP-values showed a very small difference between sCTCBCT and sCTMR and a high agreement with the reference pCT. CBCT- and MR-based sCTs have the potential to enable accurate proton dose calculations valuable for daily adaptive PT. Significant image quality differences were observed but did not affect proton dose calculation accuracy in a similar manner. Especially the recalculation of clinical treatment plans showed high agreement with the pCT for both sCTCBCT and sCTMR. [ABSTRACT FROM AUTHOR]

The purpose of this study was to compare various multisource configurations applied to cone beam CT (CBCT) using phantom imaging and Monte Carlo simulations. Image quality, scatter, and dose were evaluated in both overlapping (large cone angle) and collimated (small cone angle) configurations for CBCT. Four x-ray tube configurations were considered: traditional one source, three source overlapping, six source overlapping, and six source collimated. Image quality was evaluated on a prototype breast CT system using the following five phantoms: a Defrise phantom, a previously reported CBCT QA phantom (Corgi), a polyethylene cylinder, and two anthropomorphic phantoms (hand and knee). Scatter contamination and radiation dose were evaluated using Monte Carlo simulations of a voxelized polyethylene cylinder. The modulation of the Defrise phantom disks on average was 2.7X greater for the six source collimated configuration than the six source overlapping configuration. The data lost from cone beam artifact (spatial domain) and the null cone (frequency domain) in the overlapping configuration were completely recovered using the collimated configuration. The maximum scatter-to-primary ratio (SPR) for the overlapping configuration was 0.81 and the maximum SPR for the collimated configuration was 0.26. The average dose and maximum dose was 4X less in the collimated six source configuration when compared with the overlapping configurations. The maximum dose for the overlapping configurations (one, three & six) remained constant, but the average dose for the multisource (three & six source) overlapping configurations increased 25% when compared to the one source configuration. Use of a collimated multisource x-ray tube configuration was shown to provide significant improvements in image quality throughout the cone-beam geometry field-of-view, reduction in scatter contamination, and more efficient use of dose in comparison to both the traditional CBCT geometry with a single source and the overlapping multisource configurations. [ABSTRACT FROM AUTHOR]

Osteoporosis, characterized by reduced bone mineral density and micro-architectural degeneration, significantly enhances fracture-risk. There are several viable methods for trabecular bone micro-imaging, which widely vary in terms of technology, reconstruction principle, spatial resolution, and acquisition time. We have performed an excised cadaveric bone specimen study to evaluate different computed tomography (CT)-imaging modalities for trabecular bone micro-structural analysis. Excised cadaveric bone specimens from the distal radius were scanned using micro-CT and four in vivo CT imaging modalities: high-resolution peripheral quantitative computed tomography (HR-pQCT), dental cone beam CT (CBCT), whole-body multi-row detector CT (MDCT), and extremity CBCT. A new algorithm was developed to optimize soft thresholding parameters for individual in vivo CT modalities for computing quantitative bone volume fraction maps. Finally, agreement of trabecular bone micro-structural measures, derived from different in vivo CT imaging, with reference measures from micro-CT imaging was examined. Observed values of most trabecular measures, including trabecular bone volume, network area, transverse and plate-rod micro-structure, thickness, and spacing, for in vivo CT modalities were higher than their micro-CT-based reference values. In general, HR-pQCT-based trabecular bone measures were closer to their reference values as compared to other in vivo CT modalities. Despite large differences in observed values of measures among modalities, high linear correlation (r ε [0.94 0.99]) was found between micro-CT and in vivo CT-derived measures of trabecular bone volume, transverse and plate micro-structural volume, and network area. All HR-pQCT-derived trabecular measures, except the erosion index, showed high correlation (r ε [0.91 0.99]). The plate-width measure showed a higher correlation (r ε [0.72 0.91]) among in vivo and micro-CT modalities than its counterpart binary plate-rod characterization-based measure erosion index (r ε [0.65 0.81]). Although a strong correlation was observed between micro-structural measures from in vivo and micro-CT imaging, large shifts in their values for in vivo modalities warrant proper scanner calibration prior to adopting in multi-site and longitudinal studies. [ABSTRACT FROM AUTHOR]