IP Library Granted Patent US 12,239,853
Granted Patent B2
US 12,239,853 · App. 18/157,792 · Granted Mar 4, 2025

Systems and methods for modeling radiation source

Inventor: Zhi Shi (Shanghai, CN)
Assignee: SHANGHAI UNITED IMAGING HEALTHCARE CO., LTD.
A61N5/1075A61N5/1031A61N5/1071A61N2005/1034A61N2005/1035A61N2005/1076A61N2005/1089
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,239,853
App. No.
18/157,792
Granted
Mar 4, 2025
Kind
B2
Abstract

Systems and methods for determining a target multi-source model of a radiation source corresponding to an energy spectrum is provided. The systems may obtain an initial multi-source model of the radiation source, which includes an initial phase space file that includes information of a plurality of simulated particles of a plurality of energy levels. The systems may estimate, based on the initial phase space file, a plurality of component PDD curves corresponding to the plurality of energy levels. The systems may obtain a measured PDD curve corresponding to radiation of the energy spectrum. For each energy level, the systems may determine, based on the plurality of component PDD curves and the measured PDD curve, a weight for the each energy level. The systems may further determine the target multi-source model of the radiation source based at least in part on the initial multi-source model and the weights.

Claims (72)

1. A system, comprising:

a storage device storing a set of instructions for modeling a radiation source configured to emit radiation of an energy spectrum that includes a plurality of energy levels; and

at least one processor in communication with the storage device, wherein when executing the set of instructions, the at least one processor is configured to cause the system to perform operations including:

obtaining an initial multi-source model of the radiation source, wherein the initial multi-source model includes an initial phase space file that includes information of a plurality of simulated particles of the plurality of energy levels;

estimating, based on the initial phase space file, a plurality of component percentage depth-dose (PDD) curves in a phantom, wherein each of the plurality of component PDD curves corresponds to one of the plurality of energy levels;

obtaining a measured PDD curve in the phantom corresponding to the radiation of the energy spectrum traversing the phantom;

for each of the plurality of energy levels, determining, based on the plurality of component PDD curves and the measured PDD curve, a weight for the each energy level indicating a percentage of simulated particles of the each energy level, among the plurality of simulated particles, present in the radiation; and

determining a target multi-source model of the radiation source corresponding to the energy spectrum based at least in part on the initial multi-source model and the weights.

2. The system of claim 1 , wherein

the initial phase space file includes a group of initial weights each of which corresponds to one of the plurality of energy levels, and

the determining a weight for each of the plurality of energy levels includes:

determining, based on the plurality of component PDD curves, a combined PDD curve by adjusting at least one of the group of initial weights until a first difference between the combined PDD curve and the measured PDD curve is below a first threshold; and

determining, based on the adjusted group of weights, the weight for each of the plurality of energy levels.

3. The system of claim 1 , wherein

the radiation source includes a primary source and an electron applicator, and

the radiation includes primary electrons and secondary electrons, the primary electrons being generated by the primary source, the primary electrons including a first portion that exit the radiation source without being scattered and a second portion, and the secondary electrons being generated by the second portion of the primary electrons impinging on the electron applicator.

4. The system of claim 3 , wherein

a first portion of the simulated particles correspond to the first portion of primary electrons;

the radiation further includes photons, and

a second portion of the simulated particles correspond to the photons.

5. The system of claim 4 , wherein the photons are generated by the primary source.

6. The system of claim 4 , wherein the initial multi-source model of the radiation source includes a primary virtual source corresponding to the first portion of the simulated particles and the second portion of the simulated particles for simulating the primary source.

7. The system of claim 6 , wherein the initial phase space file includes positions of the first portion of the simulated particles and directions of the first portion of the simulated particles.

8. The system of claim 7 , wherein the positions of the first portion of the simulated particles or the directions of the first portion of the simulated particles are determined by a direct sampling based on a first distribution function, a particle flux distribution of the first portion of the simulated particles on a plane perpendicular to an axis of the primary virtual source conforming to the first distribution function.

9. The system of claim 6 , wherein

a third portion of the simulated particles corresponding to the secondary electrons, and

the initial multi-source model of the radiation source includes a secondary virtual source corresponding to the third portion of the simulated particles for simulating the electron applicator.

10. The system of claim 9 , wherein the initial phase space file includes positions of the third portion of the simulated particles and directions of the third portion of the simulated particles.

11. The system of claim 9 , wherein the secondary virtual source includes at least one of a second point source or a plane source,

the third portion of simulated particles includes a first sub-portion of simulated particles corresponding to the second point source and a second sub-portion of simulated particles corresponding to the plane source,

positions of the first sub-portion of simulated particles or directions of the first sub-portion of simulated particles being determined by a second direct sampling based on a second distribution function, wherein a particle flux distribution of the first sub-portion of simulated particles on a second plane perpendicular to the axis of the primary virtual source conforms to the second distribution function, and

positions of the second sub-portion of simulated particles or directions of the second sub-portion of simulated particles being determined by a third direct sampling based on a third distribution function, wherein a flux distribution of the second sub-portion of simulated particles on a third plane perpendicular to the axis of the primary virtual source conforms to the third distribution.

12. The system of claim 6 , wherein the at least one processor is further configured to cause the system to perform the operations including:

obtaining a measured off-axis ratio (OAR) curve corresponding to the radiation of the energy spectrum traversing the phantom;

determining a simulated OAR curve in the phantom based on parameters of the primary virtual source of the initial multi-source model and the initial phase space file;

adjusting the parameters of the primary virtual source until a second difference between a penumbra region of the simulated OAR curve and a penumbra region of the measured OAR curve is below a second threshold; and

determining the target multi-source model of the radiation source based further on the adjusted parameters of the primary virtual source.

13. The system of claim 12 , wherein the parameters of the primary virtual source include at least one of a size of the primary virtual source, a vertical position of the primary virtual source along the axis of the primary virtual source, or a particle flux distribution of the primary virtual source.

14. The system of claim 3 , wherein

the primary source is configured to be operably coupled to one of a plurality of second electron applicators, and

the at least one processor is further configured to cause the system to perform the operations including:

for each of the plurality of second electron applicators,

obtaining an output coefficient corresponding to the second electron applicator;

determining a simulated output coefficient corresponding to the second electron applicator based on structural parameters of the second electron applicator; and

determining a correction coefficient for the second electron applicator based on the output coefficient and the simulated output coefficient.

15. The system of claim 1 , wherein the determining a weight for each of the plurality of energy levels includes:

receiving a user input relating to the weight for the energy level; and

determining the weight for the energy level based at least in part on the user input.

16. The system of claim 4 , wherein the at least one processor is further configured to direct the system to perform operations including:

estimating an average energy level of the photons based on the energy spectrum;

determining a measured average energy level of the photons based on the measured PDD curve;

determining, based on the estimated average energy level and the measured energy level, a weight for the photons; and

determining the target multi-source model based further on the weight for the photons.

17. The system of claim 9 , wherein the secondary virtual source includes at least one of a second point source or a plane source,

the electron applicator includes at least an upper part, a middle part, and a lower part,

the second point source corresponds to the upper part and the middle part of the electron applicator, and

the plane source corresponds to the lower part of the electron applicator.

18. The system of claim 1 , wherein

the radiation source includes a collimation component, and

the initial multi-source model of the radiation source further includes a third virtual source corresponding to the collimation component.

19. A method for modeling a radiation source configured to emit radiation of an energy spectrum that includes a plurality of energy levels, the method being implemented on a computing device including at least one processor and at least one storage device, the method comprising:

obtaining an initial multi-source model of the radiation source, wherein the initial multi-source model includes an initial phase space file that includes information of a plurality of simulated particles of the plurality of energy levels;

estimating, based on the initial phase space file, a plurality of component percentage depth-dose (PDD) curves in a phantom, wherein each of the plurality of component PDD curves corresponds to one of the plurality of energy levels;

obtaining a measured PDD curve in the phantom corresponding to the radiation of the energy spectrum traversing the phantom;

for each of the plurality of energy levels, determining, based on the plurality of component PDD curves and the measured PDD curve, a weight for the each energy level indicating a percentage of simulated particles of the each energy level, among the plurality of simulated particles, present in the radiation; and

determining a target multi-source model of the radiation source corresponding to the energy spectrum based at least in part on the initial multi-source model and the weights.

20. A non-transitory computer readable medium, comprising executable instructions that, when executed by at least one processor, direct the at least one processor to perform a method for modeling a radiation source configured to emit radiation of an energy spectrum that includes a plurality of energy levels, the method comprising:

obtaining an initial multi-source model of the radiation source, wherein the initial multi-source model includes an initial phase space file that includes information of a plurality of simulated particles of the plurality of energy levels;

estimating, based on the initial phase space file, a plurality of component percentage depth-dose (PDD) curves in a phantom, wherein each of the plurality of component PDD curves corresponds to one of the plurality of energy levels;

obtaining a measured PDD curve in the phantom corresponding to the radiation of the energy spectrum traversing the phantom;

for each of the plurality of energy levels, determining, based on the plurality of component PDD curves and the measured PDD curve, a weight for the each energy level indicating a percentage of simulated particles of the each energy level, among the plurality of simulated particles, present in the radiation; and

determining a target multi-source model of the radiation source corresponding to the energy spectrum based at least in part on the initial multi-source model and the weights.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jul 22, 2023
From: SHI, ZHI
To: SHANGHAI UNITED IMAGING HEALTHCARE CO., LTD.
Reel/Frame 064347/0337 →
Continuity (2)
Continuation PCTCN2020109005 · Aug 13, 2020
Related Publication 20230149742A1 · May 18, 2023
References Cited (20)
US 20100082294A1 · Adnani · 2010 [cited by applicant]
US 20200066409A1 · Gasser et al. · 2020 [cited by applicant]
CN 1438602A · 2003 [cited by applicant]
CN 1927123A · 2007 [cited by applicant]
CN 101477202A · 2009 [cited by applicant]
CN 101477205A · 2009 [cited by applicant]
CN 102921115A · 2013 [cited by applicant]
CN 104043203A · 2014 [cited by applicant]
CN 105457170A · 2016 [cited by examiner]
CN 105468883A · 2016 [cited by applicant]
CN 105825067A · 2016 [cited by applicant]
CN 105866821A · 2016 [cited by applicant]
CN 106291650A · 2017 [cited by applicant]
CN 107072624A · 2017 [cited by applicant]
JP 2016223876A · 2016 [cited by applicant]
Townson et al. “GPU-based Monte Carlo radiotherapy dose calculation using phase-space sources”, IOP Publishing, Phys. Med. Biol. 58, 2013, p. 4341-4356. (Year: 2013). [cited by examiner]
Translation of CN105457170A (Year: 2016). [cited by examiner]
The First Office Action in Chinese Application No. 202080104244.8 mailed on Apr. 29, 2024, 41 pages. [cited by applicant]
International Search Report in PCT/CN2020/109005 mailed on May 17, 2021, 7 pages. [cited by applicant]
Written Opinion in PCT/CN2020/109005 mailed on May 17, 2021, 6 pages. [cited by applicant]
Cited By (2)
US 12,564,731 US 12,728,285