An adaptive spot placement method on Cartesian grid for pencil beam scanning proton therapy

TEAM: Triangular-mEsh Adaptive and Multiscale proton spot generation method

BACKGROUND

Proton therapy is preferred for its dose conformality to spare normal tissues and organs-at-risk (OAR) via Bragg peaks with negligible exit dose. However, proton dose conformality can be further optimized: (1) the spot placement is based on the structured (e.g., Cartesian) grid, which may not offer conformal shaping to complex tumor targets; (2) the spot sampling pattern is uniform, which may be insufficient at the tumor boundary to provide the sharp dose falloff, and at the same time may be redundant at the tumor interior to provide the uniform dose coverage, for example, due to multiple Coulomb scattering (MCS); and (3) the lateral spot penumbra increases with respect to the depth due to MCS, which blurs the lateral dose falloff. On the other hand, while (1) the deliverable spots are subject to the minimum-monitor-unit (MMU) constraint, and (2) the dose rate is proportional to the MMU threshold, the current spot sampling method is sensitive to the MMU threshold and can fail to provide satisfactory plan quality for a large MMU threshold (i.e., high-dose-rate delivery).

PURPOSE

This work will develop a novel Triangular-mEsh-based Adaptive and Multiscale (TEAM) proton spot generation method to address these issues for optimizing proton dose conformality and plan delivery efficiency.

METHODS

Compared to the standard clinically-used spot placement method, three key elements of TEAM are as follows: (1) a triangular mesh instead of a structured grid: the triangular mesh is geometrically more conformal to complex target shapes and therefore more efficient and accurate for dose shaping inside and around the target; (2) adaptive sampling instead of uniform sampling: the adaptive sampling consists of relatively dense sampling at the tumor boundary to create the sharp dose falloff, which is more accurate, and coarse sampling at the tumor interior to uniformly cover the target, which is more efficient; and (3) depth-dependent sampling instead of depth-independent sampling: the depth-dependent sampling is used to compensate for MCS, that is, with increasingly dense sampling at the tumor boundary to improve dose shaping accuracy, and increasingly coarse sampling at the tumor interior to improve dose shaping efficiency, as the depth increases. In the TEAM method the spot locations are generated for each energy layer and layer-by-layer in the multiscale fashion; and then the spot weights are derived by solving the IMPT problem of dose-volume planning objectives, MMU constraints, and robustness optimization with respect to range and setup uncertainties.

RESULTS

Compared to the standard clinically-used spot placement method UNIFORM, TEAM achieved (1) better plan quality using <60% number of spots of UNIFORM; (2) better robustness to the number of spots; (3) better robustness to a large MMU threshold. Furthermore, TEAM provided better plan quality with fewer spots than other adaptive methods (Cartesian-grid or triangular-mesh).

CONCLUSIONS

A novel triangular-mesh-based proton spot placement method called TEAM is proposed, and it is demonstrated to improve plan quality, robustness to the number of spots, and robustness to the MMU threshold, compared to the clinically-used spot placement method and other adaptive methods.

A novel fast robust optimization algorithm for intensity-modulated proton therapy with minimum monitor unit constraint

BACKGROUND

Intensity-modulated proton therapy (IMPT) optimizes spot intensities and position, providing better conformability. However, the successful application of IMPT is dependent upon addressing the challenges posed by range and setup uncertainties. In order to address the uncertainties in IMPT, robust optimization is essential.

PURPOSE

This study aims to develop a novel fast algorithm for robust optimization of IMPT with minimum monitor unit (MU) constraint.

METHODS AND MATERIALS

The study formulates a robust optimization problem and proposes a novel, fast algorithm based on the alternating direction method of multipliers (ADMM) framework. This algorithm enables distributed computation and parallel processing. Ten clinical cases were used as test scenarios to evaluate the performance of the proposed approach. The robust optimization method (RBO-NEW) was compared with plans that only consider nominal optimization using CTV (NMO-CTV) without handling uncertainties and PTV (NMO-PTV) to handle the uncertainties, as well as with conventional robust-optimized plans (RBO-CONV). Dosimetric metrics, including D95, homogeneity index, and Dmean, were used to evaluate the dose distribution quality. The area under the root-mean-square dose (RMSD)-volume histogram curves (AUC) and dose-volume histogram (DVH) bands were used to evaluate the robustness of the treatment plan. Optimization time cost was also assessed to measure computational efficiency.

RESULTS

The results demonstrated that the RBO plans exhibited better plan quality and robustness than the NMO plans, with RBO-NEW showing superior computational efficiency and plan quality compared to RBO-CONV. Specifically, statistical analysis results indicated that RBO-NEW was able to reduce the computational time from389.70 ± 207.40 $389.70\pm 207.40$ to228.60 ± 123.67 $228.60\pm 123.67$ s (p < 0.01 $p<0.01$ ) and reduce the mean organ-at-risk (OAR) dose from9.38 ± 12.80 $9.38\pm 12.80$ % of the prescription dose to9.07 ± 12.39 $9.07\pm 12.39$ % of the prescription dose (p < 0.05 $p<0.05$ ) compared to RBO-CONV.

CONCLUSION

This study introduces a novel fast robust optimization algorithm for IMPT treatment planning with minimum MU constraint. Such an algorithm is not only able to enhance the plan's robustness and computational efficiency without compromising OAR sparing but also able to improve treatment plan quality and reliability.

Cardinality-constrained plan-quality and delivery-time optimization method for proton therapy

BACKGROUND

While minimizing plan delivery time is beneficial for proton therapy in terms of motion management, patient comfort, and treatment throughput, it often poses a tradeoff with optimizing plan quality. A key component of plan delivery time is the energy switching time, which is approximately proportional to the number of energy layers, that is, the cardinality.

PURPOSE

This work aims to develop a novel optimization method that can efficiently compute the pareto surface between plan quality and energy layer cardinality, for the planner to navigate through this quality-and-efficiency tradeoff and select the appropriate plan of a balanced tradeoff.

METHODS

A new IMPT method CARD is proposed that (1) explicitly incorporates the minimization of energy layer cardinality as an optimization objective, and (2) automatically generates a set of plans sequentially with a descending order in number of energy layers. The energy layer cardinality is penalized through the l1,0-norm regularization with an upper bound, and the upper bound is monotonically decreased to compute a series of treatment plans with gradually decreased energy layer cardinality on the quality-and-efficiency pareto surface. For any given treatment plan, the plan optimality is enforced using dose-volume planning objectives and the plan deliverability is imposed through minimum-monitor-unit (MMU) constraints, with optimization solution algorithm based on iterative convex relaxation.

RESULTS

The new method CARD was validated in comparison with the benchmark plan of all energy layers (P0), and a state-of-the-art method called MMSEL, using prostate, head-and-neck (HN), lung, pancreas, liver and brain cases. While labor-intensive and time-consuming manual parameter tuning was needed for MMSEL to generate plans of predefined energy layer cardinality, CARD automatically and efficiently computed all plans with sequentially decreasing predefined energy layer cardinality all at once. With the acceptable plan quality (i.e., no more than 110% of total optimization objective value from P0), CARD achieved the reduction of number of energy layers to 52% (from 77 to 40), 48% (from 135 to 65), 59% (from 85 to 50), 67% (from 52 to 35), 80% (from 50 to 40), and 30% (from 66 to 20), for prostate, HN, lung, pancreas, liver, and brain cases, respectively, compared to P0, with overall better plan quality than MMSEL. Moreover, due to the nonconvexity of the MMU constraint, CARD provided the similar or even smaller optimization objective than P0, at the same time with fewer number of energy layers, that is, 55 versus 77, 85 versus 135, 45 versus 52, and 25 versus 66 for prostate, HN, pancreas, and brain cases, respectively.

CONCLUSIONS

We have developed a novel optimization algorithm CARD that can efficiently and automatically compute a series of treatment plans of any given energy layer sequentially, which allows the planner to navigate through the plan-quality and energy-layer-cardinality tradeoff and select the appropriate plan of a balanced tradeoff.

Biological optimization for hybrid proton-photon radiotherapy

Objective.Hybrid proton-photon radiotherapy (RT) is a cancer treatment option to broaden access to proton RT. Additionally, with a refined treatment planning method, hybrid RT has the potential to offer superior plan quality compared to proton-only or photon-only RT, particularly in terms of target coverage and sparing organs-at-risk (OARs), when considering robustness to setup and range uncertainties. However, there is a concern regarding the underestimation of the biological effect of protons on OARs, especially those in close proximity to targets. This study seeks to develop a hybrid treatment planning method with biological dose optimization, suitable for clinical implementation on existing proton and photon machines, with each photon or proton treatment fraction delivering a uniform target dose.Approach.The proposed hybrid biological dose optimization method optimized proton and photon plan variables, along with the number of fractions for each modality, minimizing biological dose to the OARs and surrounding normal tissues. To mitigate underestimation of hot biological dose spots, proton biological dose was minimized within a ring structure surrounding the target. Hybrid plans were designed to be deliverable separately and robustly on existing proton and photon machines, with enforced uniform target dose constraints for the proton and photon fraction doses. A probabilistic formulation was utilized for robust optimization of setup and range uncertainties for protons and photons. The nonconvex optimization problem, arising from minimum monitor unit constraint and dose-volume histogram constraints, was solved using an iterative convex relaxation method.Main results.Hybrid planning with biological dose optimization effectively eliminated hot spots of biological dose, particularly in normal tissues surrounding the target, outperforming proton-only planning. It also provided superior overall plan quality and OAR sparing compared to proton-only or photon-only planning strategies.Significance.This study presents a novel hybrid biological treatment planning method capable of generating plans with reduced biological hot spots, superior plan quality to proton-only or photon-only plans, and clinical deliverability on existing proton and photon machines, separately and robustly.

Efficient contour-based annotation by iterative deep learning for organ segmentation from volumetric medical images

PURPOSE

Training deep neural networks usually require a large number of human-annotated data. For organ segmentation from volumetric medical images, human annotation is tedious and inefficient. To save human labour and to accelerate the training process, the strategy of annotation by iterative deep learning recently becomes popular in the research community. However, due to the lack of domain knowledge or efficient human-interaction tools, the current AID methods still suffer from long training time and high annotation burden.

METHODS

We develop a contour-based annotation by iterative deep learning (AID) algorithm which uses boundary representation instead of voxel labels to incorporate high-level organ shape knowledge. We propose a contour segmentation network with a multi-scale feature extraction backbone to improve the boundary detection accuracy. We also developed a contour-based human-intervention method to facilitate easy adjustments of organ boundaries. By combining the contour-based segmentation network and the contour-adjustment intervention method, our algorithm achieves fast few-shot learning and efficient human proofreading.

RESULTS

For validation, two human operators independently annotated four abdominal organs in computed tomography (CT) images using our method and two compared methods, i.e. a traditional contour-interpolation method and a state-of-the-art (SOTA) convolutional network (CNN) method based on voxel label representation. Compared to these methods, our approach considerably saved annotation time and reduced inter-rater variabilities. Our contour detection network also outperforms the SOTA nnU-Net in producing anatomically plausible organ shape with only a small training set.

CONCLUSION

Taking advantage of the boundary shape prior and the contour representation, our method is more efficient, more accurate and less prone to inter-operator variability than the SOTA AID methods for organ segmentation from volumetric medical images. The good shape learning ability and flexible boundary adjustment function make it suitable for fast annotation of organ structures with regular shape.

The emerging role of photoacoustic imaging in clinical oncology

Clinical oncology can benefit substantially from imaging technologies that reveal physiological characteristics with multiscale observations. Complementing conventional imaging modalities, photoacoustic imaging (PAI) offers rapid imaging (for example, cross-sectional imaging in real time or whole-breast scanning in 10-15 s), scalably high levels of spatial resolution, safe operation and adaptable configurations. Most importantly, this novel imaging modality provides informative optical contrast that reveals details on anatomical, functional, molecular and histological features. In this Review, we describe the current state of development of PAI and the emerging roles of this technology in cancer screening, diagnosis and therapy. We comment on the performance of cutting-edge photoacoustic platforms, and discuss their clinical applications and utility in various clinical studies. Notably, the clinical translation of PAI is accelerating in the areas of macroscopic and mesoscopic imaging for patients with breast or skin cancers, as well as in microscopic imaging for histopathology. We also highlight the potential of future developments in technological capabilities and their clinical implications, which we anticipate will lead to PAI becoming a desirable and widely used imaging modality in oncological research and practice.