Masoud Zarepisheh
Associate Professor
Department of Medical Physics
Memorial Sloan Kettering Cancer Center
321 East 61st Street, New York, NY, 10065
| Email | Google Scholar | LinkedIn | Github |
About
I have finished my Ph.D. in applied mathematics, mainly focusing on the theoretical aspects and algorithm design of Mathematical Optimization techniques. I have later shifted my focus, during my postdocs, on the applications of mathematical modeling and optimization in Medical Physics, and particularly in radiotherapy cancer treatment planning. I am currently a faculty member at MSKCC where I am working closely with clinicians to leverage mathematical optimization techniques in real clinical environments.
What is Radiotherapy Cancer Treatment Planning?
Radiotherapy (aka radiation therapy) is one of the major modalities of cancer treatment that can be prescribed as the main treatment or can be given in combination with surgery and/or chemotherapy. Radiotherapy is the use of radiation beams, precisely directed to diseased cells, to fight cancer. The main objective of radiotherapy is to deliver radiation to the tumor to kill cancerous cells without harming the surrounding healthy tissue. Each cancer patient’s tumor is unique in terms of its shape, location, and the surrounding healthy tissue. The machine settings, including the beam angles and beam shapes, must be customized based on the patient’s image data and the physician’s dose prescription for both the tumor and nearby healthy tissues. This gives rise to large-scale, multi-criteria, and often non-convex optimization problems.
Want to learn more? Read review papers (Zarepisheh et al. 2021, Breedveld et al. 2019, Ehrgott et al. 2010) or watch YouTube videos (Edelman competition, Varian IMRT, Elekta VMAT).
Clinical Deployment of Automated Planning at MSKCC
In radiotherapy, the machine settings must be customized for each cancer patient and this requires extensive iterative manual tuning of the parameters. This process is: 1) time-consuming, 2) labor-intensive, and 3) the plan quality heavily dependens on the skill and experience of the planner. We have recently developed and clinically deployed an in-house automated treatment planning system called ECHO (Expedited Constrained Hierarchical Optimization). ECHO is seamlessly integrated into an FDA-approved commertial treatment planning system Eclipse and being used in our daily clinical routine with more than 10,000 treated to date (May, 2024). ECHO was nominated as a finalist for the prestigious Franz Edelman 2021 award (YouTube | Slides | Paper | Podcast | News).
The following table summarizes some of the technical chanllenges in this project and the optimization tools we have employed to tackle them.
Implementation Status | Challenge | Optimization Tools | Publication |
---|---|---|---|
Clinically Implemented | Multi-Criteria Challenge | Hierarchical Constrained Nonlinear Optimization | Zarepisheh et al. (2019), Zarepisheh et al. (2021) |
Clinically Implemented | Computational Complexity (Dense Influence Matrix) | Lagrangian Methods | Zarepisheh et al. (2019) |
Clinically Implemented | Non-convexity in Clinical Criteria (DVH) | Mixed Integer Programming, Convex Relaxtion | Mukherjee et al. (2020) |
Clinically Implemented | Non-convexity in Machine Parameters (VMAT) | Sequential Convex Programming | Dursun et al. (2021), Dursun et al. (2023) |
Research/Development | Beam Angle Selection (Combinatorial Optimization Problem) | Bayesian Optimization | Taasti et al. (2020) |
Research/Development | Uncertainty Management | Robust Optimization (P-Norm Function) | Taasti et al. (2020) |
Research/Development | Multi-Criteria Challenge | Deep Learning Predictions (Moment-Based Loss Functions) | Jhanwar et al. (2022) |
Research/Development | Increasing Computational Complexity With Uncertainty Managements | Distributed Optimization (ADMM) | Fu et al. (2022), Zarepisheh et al. (2018) |
Research/Development | Reducing IMRT Plan Complexity | Wavelet-induced smoothness | Tefagh, Zarepisheh (2023) |
Open-Source Python Package (PortPy)
PortPy is a codebase open-source package providing research-ready data and code to facilitate the development and clinical translation of radiotherapy treatment planning algorithms. PortPy includes benchmark data (pre-calculated data from an FDA-approved commercial treatment planning system Eclipse) and benchmark code (algorithms for optimizing beamlets, beam angles, apertures, and robust optimization).
PortPy is integrated with an open-source optimization package CVXPy which is a high-level programming language for optimization and provides out-of-the-box access to many commercial/open-source optimization engines (e.g., CPLEX, Mosek, Gorubi). PortPy is also integrated with a widely used open-source medical imaging package 3D-Slicer which offers advanced visualization tools.