Simulation of a new respiratory phase sorting method for 4D-imaging using optical surface information towards precision radiotherapy

Background: Respiratory signal detection is critical for 4-dimensional (4D) imaging. This study proposes and evaluates a novel phase sorting method using optical surface imaging (OSI), aiming to improve the precision of radiotherapy. Method: Based on 4D Extended Cardiac-Torso (XCAT) digital phantom, OSI in point cloud format was generated from the body segmentation, and image projections were simulated using the geometries of Varian 4D kV cone-beam-CT (CBCT). Respiratory signals were extracted respectively from the segmented diaphragm image (reference method) and OSI respectively, where Gaussian Mixture Model and Principal Component Analysis (PCA) were used for image registration and dimension reduction respectively. Breathing frequencies were compared using Fast-Fourier-Transform. Consistency of 4DCBCT images reconstructed using Maximum Likelihood Expectation Maximization algorithm was also evaluated quantitatively, where high consistency can be suggested by lower Root-Mean-Square-Error (RMSE), Structural-Similarity-Index (SSIM) value closer to 1, and larger Peak-Signal-To-Noise-Ratio (PSNR) respectively. Results: High consistency of breathing frequencies was observed between the diaphragm-based (0.232 Hz) and OSI-based (0.251 Hz) signals, with a slight discrepancy of 0.019Hz. Using end of expiration (EOE) and end of inspiration (EOI) phases as examples, the mean±1SD values of the 80 transverse, 100 coronal and 120 sagittal planes were 0.967, 0,972, 0.974 (SSIM); 1.657 ± 0.368, 1.464 ± 0.104, 1.479 ± 0.297 (RMSE); and 40.501 ± 1.737, 41.532 ± 1.464, 41.553 ± 1.910 (PSNR) for the EOE; and 0.969, 0.973, 0.973 (SSIM); 1.686 ± 0.278, 1.422 ± 0.089, 1.489 ± 0.238 (RMSE); and 40.535 ± 1.539, 41.605 ± 0.534, 41.401 ± 1.496 (PSNR) for EOI respectively. Conclusions: This work proposed and evaluated a novel respiratory phase sorting approach for 4D imaging using optical surface signals, which can potentially be applied to precision radiotherapy. Its potential advantages were non-ionizing, non-invasive, non-contact, and more compatible with various anatomic regions and treatment/imaging systems

Optimization strategies for respiratory motion management in stereotactic body radiation therapy

Various challenges arise during the treatment of lung tumors with stereotactic body radiation therapy (SBRT), which is a form of hypofractionated high precision conformal radiation therapy delivered to small targets. The dose is applied in only a few fractions and respiratory organ and tumor motion is a source of uncertainty additional to interfractional set-up errors. Respiratory organ and tumor motion is highly patient-specific and it affects the whole radiotherapy treatment chain. In this thesis, motion management techniques for SBRT are evaluated and improved in a clinical setting. A clinical need for improvement has been present at the LMU university hospital for each issue addressed in this thesis: Initially, the usage of respiratory correlated computed tomography (4DCT), which is vital for SBRT treatment, was seen as impractical and prone to uncertainties in the data reconstruction in its current form. Therefore, the 4DCT reconstruction workflow has been improved to minimize these potential error sources. Secondly, treatment planning for tumors affected by respiratory motion was evaluated and subsequently improved. Finally, the treatment technique of respiratory gating was implemented at the clinic, which led to the need of evaluating the respiratory gating characteristics of the novel system configuration. At first, the 4DCT reconstruction workflow used in clinical practice was investigated, as in the presence of respiratory motion the knowledge of tumor position over time is essential in SBRT treatments. Using 4DCT, the full motion range of the individual tumor can be determined. However, certain 4DCT reconstruction methods can under- or overestimate tumor motion due to limitations in the data acquisition scheme and due to the incorrect sorting of certain X-ray computed tomography (CT) image slices into different respiratory phases. As the regular clinical workflow of cycle-based sorting (CBS) without maximum inspiration detection (and therefore no clear starting point for the individual breathing cycles) seemed to be affected by these potential errors, the usage of CBS with correct maximum detection and another sorting algorithm of the respiration states, so-called local amplitude-based sorting (LAS), both have been implemented for a reduction of image artifacts and improved 4DCT quality. The three phase binning algorithms have been investigated in a phantom study (using 10 different breathing waveforms) and in a patient study (with 10 different patients). The mis-representation of the tumor volume was reduced in both implemented sorting algorithms compared to the previously used CBS approach (without correct maximum detection) in the phantom and the patient study. The clinical recommendation was the use of CBS with improved maximum detection, as too many manual interventions would be needed for the LAS workflow. Secondly, a combination of the actual patient breathing trace during treatment, the log files generated by the linear accelerator (LINAC), and Monte Carlo (MC) four-dimensional (4D) dose calculations for each individual fraction was implemented as a 4D dose evaluation tool. This workflow was tested in a clinical environment for SBRT treatment planning on multiple CT datasets featuring: a native free-breathing 3DCT, an average intensity projection (AIP) as well as a maximum intensity projection (MIP), both obtained from the patient's 4DCT, and density overrides (DOs) in a 3DCT. This study has been carried out for 5 SBRT patients for three-dimensional conformal radiation therapy (3D-CRT) and volumetric modulated arc therapy (VMAT) treatment plans. The dose has been recalculated on each 4DCT breathing phase according the the patient's breathing waveform and accumulated to the gross tumor volume (GTV) at the end-of-exhale (EOE) breathing phase using deformable image registration. Even though the least differences in planned and recalculated dose were found for AIP and MIP treatment planning, the results indicate a strong dependency on individual tumor motion due to the variability of breathing motion in general, and on tumor size. The combination of the patient's individual breathing trace during each SBRT fraction with 4D MC dose calculation based on the LINAC log file information leads to a good approximation of actual dose delivery. Finally, in order to ensure precise and accurate treatment for respiratory gating techniques, the technical characteristics of the LINAC in combination with a breathing motion monitoring system as s surrogate for tumor motion have to be identified. The dose delivery accuracy and the latency of a surface imaging system in connection with a modern medical LINAC were investigated using a dynamic breathing motion phantom. The dosimetric evaluation has been carried out using a static 2D-diode array. The measurement of the dose difference between gated and ungated radiation delivery was found to be below 1% (for clinical relevant gating levels of about 30%). The beam-on latency, or time delay, determined using radiographic films was found to be up to 851 ms±100 ms. With these known parameters, an adjustment of the pre-selected gating level or the internal target volume (ITV) margins could be made. With the highly patient-specific character of respiratory motion, lung SBRT faces many additional challenges besides the specific issues addressed in this thesis. However, the findings of this thesis have improved clinical workflows at the Department of Radiation Oncology of the LMU University hospital. In a future perspective, a workflow using evaluation of the actual 4D dose in combination with accurate 4DCT image acquisition and specialized treatment delivery (such as respiratory gating) has the potential for a safe further reduction of treatment margins and increased sparing of organs-at-risk (OARs) in SBRT without compromising tumor dose targeting accuracy

Optimization strategies for respiratory motion management in stereotactic body radiation therapy

Various challenges arise during the treatment of lung tumors with stereotactic body radiation therapy (SBRT), which is a form of hypofractionated high precision conformal radiation therapy delivered to small targets. The dose is applied in only a few fractions and respiratory organ and tumor motion is a source of uncertainty additional to interfractional set-up errors. Respiratory organ and tumor motion is highly patient-specific and it affects the whole radiotherapy treatment chain. In this thesis, motion management techniques for SBRT are evaluated and improved in a clinical setting. A clinical need for improvement has been present at the LMU university hospital for each issue addressed in this thesis: Initially, the usage of respiratory correlated computed tomography (4DCT), which is vital for SBRT treatment, was seen as impractical and prone to uncertainties in the data reconstruction in its current form. Therefore, the 4DCT reconstruction workflow has been improved to minimize these potential error sources. Secondly, treatment planning for tumors affected by respiratory motion was evaluated and subsequently improved. Finally, the treatment technique of respiratory gating was implemented at the clinic, which led to the need of evaluating the respiratory gating characteristics of the novel system configuration. At first, the 4DCT reconstruction workflow used in clinical practice was investigated, as in the presence of respiratory motion the knowledge of tumor position over time is essential in SBRT treatments. Using 4DCT, the full motion range of the individual tumor can be determined. However, certain 4DCT reconstruction methods can under- or overestimate tumor motion due to limitations in the data acquisition scheme and due to the incorrect sorting of certain X-ray computed tomography (CT) image slices into different respiratory phases. As the regular clinical workflow of cycle-based sorting (CBS) without maximum inspiration detection (and therefore no clear starting point for the individual breathing cycles) seemed to be affected by these potential errors, the usage of CBS with correct maximum detection and another sorting algorithm of the respiration states, so-called local amplitude-based sorting (LAS), both have been implemented for a reduction of image artifacts and improved 4DCT quality. The three phase binning algorithms have been investigated in a phantom study (using 10 different breathing waveforms) and in a patient study (with 10 different patients). The mis-representation of the tumor volume was reduced in both implemented sorting algorithms compared to the previously used CBS approach (without correct maximum detection) in the phantom and the patient study. The clinical recommendation was the use of CBS with improved maximum detection, as too many manual interventions would be needed for the LAS workflow. Secondly, a combination of the actual patient breathing trace during treatment, the log files generated by the linear accelerator (LINAC), and Monte Carlo (MC) four-dimensional (4D) dose calculations for each individual fraction was implemented as a 4D dose evaluation tool. This workflow was tested in a clinical environment for SBRT treatment planning on multiple CT datasets featuring: a native free-breathing 3DCT, an average intensity projection (AIP) as well as a maximum intensity projection (MIP), both obtained from the patient's 4DCT, and density overrides (DOs) in a 3DCT. This study has been carried out for 5 SBRT patients for three-dimensional conformal radiation therapy (3D-CRT) and volumetric modulated arc therapy (VMAT) treatment plans. The dose has been recalculated on each 4DCT breathing phase according the the patient's breathing waveform and accumulated to the gross tumor volume (GTV) at the end-of-exhale (EOE) breathing phase using deformable image registration. Even though the least differences in planned and recalculated dose were found for AIP and MIP treatment planning, the results indicate a strong dependency on individual tumor motion due to the variability of breathing motion in general, and on tumor size. The combination of the patient's individual breathing trace during each SBRT fraction with 4D MC dose calculation based on the LINAC log file information leads to a good approximation of actual dose delivery. Finally, in order to ensure precise and accurate treatment for respiratory gating techniques, the technical characteristics of the LINAC in combination with a breathing motion monitoring system as s surrogate for tumor motion have to be identified. The dose delivery accuracy and the latency of a surface imaging system in connection with a modern medical LINAC were investigated using a dynamic breathing motion phantom. The dosimetric evaluation has been carried out using a static 2D-diode array. The measurement of the dose difference between gated and ungated radiation delivery was found to be below 1% (for clinical relevant gating levels of about 30%). The beam-on latency, or time delay, determined using radiographic films was found to be up to 851 ms±100 ms. With these known parameters, an adjustment of the pre-selected gating level or the internal target volume (ITV) margins could be made. With the highly patient-specific character of respiratory motion, lung SBRT faces many additional challenges besides the specific issues addressed in this thesis. However, the findings of this thesis have improved clinical workflows at the Department of Radiation Oncology of the LMU University hospital. In a future perspective, a workflow using evaluation of the actual 4D dose in combination with accurate 4DCT image acquisition and specialized treatment delivery (such as respiratory gating) has the potential for a safe further reduction of treatment margins and increased sparing of organs-at-risk (OARs) in SBRT without compromising tumor dose targeting accuracy

Novel Method for More Efficient Optimizing the Knowledge-Based Planning: Specific Voxels of each Structure Influenced by Dominant Beamlets (SVSIDB)

There is a huge problem and time-consuming computation to optimize the IMRT treatment plan. Extracting the optimized plan from the predicted 3D3 so-called optimizing the KBP is also involved in this challenge. Some algorithms and methods have been presented for clustering and down-sampling the voxels to make the problem smaller, in recent years. In the current research, a novel down-sampling method is presented for optimizing the knowledge-based planning more efficiently. The concept of SVSIDB and corresponding down-sampling algorithm are proposed under title of SMP-2. The algorithm has been run on the data of 30 patients from the Open-KBP dataset. For each patient, there are 19 sets of dose prediction data in this dataset. Therefore, a total of 570 KBP-optimizing problems have been solved by applying the QuadLin model in the CVX framework. Resulted plans are evaluated and compared regarding two main fields which are the quality of the treatment plan as well as the computation efficiency. Solve time is the evaluation criteria for the latter field i.e. computation efficiency. The results of the current study indicated a remarkable improvement in the computation efficiency. Accordingly, the proposed method, SMP-2, reduced the average solving time by 46% in comparison to the full-data QuadLin model. The results also show an up to 53% reduction in solve time along with up to 22% improvement in clinical criteria compared to the previous research. Evaluation of the research results indicated that the SVSIDB has not only reduced the solve time but also improved the quality of the treatment plans. This is a remarkable achievement of the proposed model compared to the previous research and confirmed the significant effectiveness of the SVSIDB method which has the potential of even more improvement of the computation efficiency.Comment: Keywords: Computing Efficiency; Voxels Down-Sampling; Radiotherapy Treatment Planning; Clustering. arXiv admin note: text overlap with arXiv:2206.0583

Evaluation of interplay and organ motion effects by means of 4D dose reconstruction and accumulation

PURPOSE: Pencil beam scanned proton therapy (PBS-PT) treatment quality might be compromised by interplay and motion effects. Via fraction-wise reconstruction of 4D dose distributions and dose accumulation, we assess the clinical relevance of motion related target dose degradation in thoracic cancer patients. METHODS AND MATERIALS: For the ten thoracic patients (Hodgkin lymphoma and non-small cell lung cancer) treated at our proton therapy facility, daily breathing pattern records, treatment delivery log-files and weekly repeated 4DCTs were collected. Patients exhibited point-max target motion of up to 20 mm. They received robustly optimized treatment plans, delivered with five-times rescanning in fractionated regimen. Treatment delivery records were used to reconstruct 4D dose distributions and the accumulated treatment course dose per patient. Fraction-wise target dose degradations were analyzed and the accumulated treatment course dose, representing an estimation of the delivered dose, was compared with the prescribed dose. RESULTS: No clinically relevant loss of target dose homogeneity was found in the fraction-wise reconstructed 4D dose distributions. Overall, in 97% of all reconstructed fraction doses, D98 remained within 5% from the prescription dose. The V95 of accumulated treatment course doses was higher than 99.7% for all ten patients. CONCLUSIONS: 4D dose reconstruction and accumulation enables the clinical estimation of actual exhibited interplay and motion effects. In the patients considered here, the loss of homogeneity caused by interplay and organ motion did not show systematic pattern and smeared out throughout the course of fractionated PBS-PT treatment. Dose degradation due to anatomical changes showed to be more severe and triggered treatment adaptations for five patients

Artificial Intelligence and Medicine

The introduction of artificial intelligence (AI) has resulted in numerous technological advancements in the medical profession and a radical transformation of the old medical model. Artificial intelligence in medicine consists mostly of machine learning, deep learning, expert systems, intelligent robotics, the internet of medical things, and other prevalent and new AI technology. The primary applications of AI in the medical industry are intelligent screening, intelligent diagnosis, risk prediction, and supplemental treatment. Presently, medical AI has achieved significant advances, and big data quality management, new technology empowerment innovation, multi-domain knowledge integration, and personalized medical decision-making will exhibit greater growth potential in the clinical arena