EP3941576A1 - Dosimetric projection of collimator geometry - Google Patents
Dosimetric projection of collimator geometryInfo
- Publication number
- EP3941576A1 EP3941576A1 EP20709334.5A EP20709334A EP3941576A1 EP 3941576 A1 EP3941576 A1 EP 3941576A1 EP 20709334 A EP20709334 A EP 20709334A EP 3941576 A1 EP3941576 A1 EP 3941576A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- dosimetric
- calculating
- collimator
- collimating device
- geometry
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000005855 radiation Effects 0.000 claims abstract description 95
- 238000000034 method Methods 0.000 claims abstract description 66
- 238000004364 calculation method Methods 0.000 claims description 36
- 238000002719 stereotactic radiosurgery Methods 0.000 claims description 18
- 230000008021 deposition Effects 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 8
- 238000005094 computer simulation Methods 0.000 claims description 4
- 230000032258 transport Effects 0.000 description 25
- 238000004422 calculation algorithm Methods 0.000 description 23
- 230000015654 memory Effects 0.000 description 17
- 238000001959 radiotherapy Methods 0.000 description 13
- 230000006870 function Effects 0.000 description 10
- 238000009826 distribution Methods 0.000 description 6
- 230000005291 magnetic effect Effects 0.000 description 6
- 238000003860 storage Methods 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000011455 3D conformal radiation therapy Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000002728 volumetric modulated arc therapy Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 1
- 229910052721 tungsten Inorganic materials 0.000 description 1
- 239000010937 tungsten Substances 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
- A61N5/1031—Treatment planning systems using a specific method of dose optimization
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1042—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head
- A61N5/1045—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head using a multi-leaf collimator, e.g. for intensity modulated radiation therapy or IMRT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1001—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy using radiation sources introduced into or applied onto the body; brachytherapy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1077—Beam delivery systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
Definitions
- This application relates generally to radiation methods and apparatuses.
- various embodiments of methods of dosimetric projection of collimator geometry and calculation of dose deposition for treatment planning are described.
- External beam radiotherapy requires treatment planning to determine a desired radiation dose to the target and a maximal dose which can be safely absorbed by healthy tissue or organs within the treatment volume near the target.
- dose deposited within a treatment volume is calculated or predicted using a suitable dose calculation algorithm, which may involve tracking of the radiation beam from the source to the treatment volume.
- collimating devices are used in a treatment machine to size, shape and/or modulate the intensity of the beam.
- MLCs multileaf collimators
- An MLC includes a plurality of pairs of beam-blocking leaves which can be independently moved in and out of a radiation beam. In use, a number of selected beam-blocking leaves are positioned in the radiation beam forming an aperture through which the unblocked beam passes.
- the shape and size of the aperture in the MLC define a shaped field at the target in the isocenter plane.
- the aperture in an MLC is configured sufficiently small for small-field radiotherapy or stereotactic radiosurgery (SRS).
- SRS cones may be used to perform stereotactic radiosurgery.
- An SRS cone is a collimating device typically made from tungsten having a small conical hole to allow a radiation beam passing through for focused beam delivery.
- a radiation beam from a source passes through an MLC, SRS cone, and/or other collimating devices, the characteristics of the beam, including the shape or fluence of the beam, are modified and thus need be determined to provide input for a dose calculation algorithm to compute or predict the dose deposited in the treatment volume.
- One conventional approach of calculating the shape or fluence of a beam exiting a collimating device such as an MLC is to project the MLC geometry along the beamline onto a plane, essentially flattening the three-dimensional (3D) MLC geometry into a two-dimensional (2D) opening, where the structural details of the MLC leaf tips are not taken into account in the calculation.
- Embodiments of this disclosure provide for a method of tracking or calculating beam fluence passing through a collimating device.
- the three-dimensional (3D) shape of the collimating device such as an MLC
- the three-dimensional (3D) shape of the collimating device is taken into account by a local dosimetric projection of the geometry.
- the dimension of an MLC part along the beam is projected into a dosimetric opacity value or dosimetric thickness.
- the dosimetric projection is based on the realistic 3D geometry of the MLC, but an algorithm can be made faster than a full 3D model by making the dosimetric projection local (e.g., around the MLC leaf tip and edges) and by defining bounding points or surfaces.
- the dosimetric opacity can be calculated at the geometrical boundaries by e.g. ray-tracing or analytical calculation, and by interpolating between them or using analytical approximations.
- the dependence on variables such as collimator shape and position can be taken into account by, e.g., parameterizing the dosimetric projection as a function of the variable or by solving the dependence by analytical approximation.
- the dosimetrically projected opacity values or thicknesses can be then used as input in a dose calculation algorithm within a fluence calculation module.
- the disclosed method requires a computational effort comparable to the current 2D approach, but improves the accuracy of dose calculation significantly to a level comparable to a full 3D model.
- a method comprises the steps of defining a 3D geometry of a collimating device which defines an aperture configured to allow a radiation beam passing through, projecting the collimating device along the radiation beam into a 2D geometry in a plane, calculating dosimetric opacity values of the collimating device at locations adjacent to the aperture based on the 3D geometry of the collimating device, and calculating transport of the radiation beam through the collimating device based on the 2D geometry projected in the plane and using the dosimetric opacity values of the collimating device at the locations adjacent to the aperture.
- a computer product comprises a non-transitory computer-readable medium storing instructions executable by a computer system.
- the instructions comprises defining a 3D geometry of a collimating device which defines an aperture configured to allow a radiation beam passing through, projecting the collimating device along the radiation beam into a 2D geometry in a plane, calculating dosimetric opacity values of the collimating device at locations adjacent to the aperture based on the 3D geometry of the collimating device, and calculating transport of the radiation beam through the collimating device based on the 2D geometry projected in the plane and using the dosimetric opacity values of the collimating device at the locations adjacent to the aperture.
- a computer system comprises a processor and memory storing instructions executable by the processor.
- the instructions comprises defining a 3D geometry of a collimating device which defines an aperture configured to allow a radiation beam passing through, projecting the collimating device along the radiation beam into a 2D geometry in a plane, calculating dosimetric opacity values of the collimating device at locations adjacent to the aperture based on the 3D geometry of the collimating device, and calculating transport of the radiation beam through the collimating device based on the 2D geometry projected in the plane and using the dosimetric opacity values of the collimating device at the locations adjacent to the aperture.
- FIG. 1 is a simplified illustration of a radiation system according to embodiments of this disclosure.
- FIG. 1A is an enlarged perspective view of the multileaf collimator (MLC) shown in FIG. 1.
- MLC multileaf collimator
- FIG. 2 illustrates a three-dimensional (3D) model for calculating beam transport through an MLC.
- FIG. 3 illustrates a two-dimensional (2D) model for calculating beam fluence through an MLC.
- FIGS. 4 and 5 illustrate a dosimetric projection model for calculating beam fluence through an MLC according to embodiments of the disclosure.
- FIG. 6 is a flowchart illustrating a method of calculating dose deposition using dosimetric projection of a collimator geometry according to embodiments of the disclosure.
- FIG. 7 illustrates an example method of calculating dosimetric opacity values for an MLC leaf according to embodiments of the disclosure.
- FIG. 8 is a flowchart illustrating example steps of dosimetric projection of a collimator geometry according to embodiments of the disclosure.
- FIG. 9 is an image showing a result of dosimetric projection of an example collimator according to embodiments of the disclosure.
- FIG. 10 shows plots comparing the accuracies of dose calculation based on various models.
- FIG. 11 is a simplified illustration of a computing system according to embodiments of the disclosure.
- FIG. 1 is a simplified illustration of a radiation system 100 in which the methods of the disclosure can be implemented.
- the radiation system 100 may include a radiation source 102 producing or emitting a beam 104 of radiation such as photons, electrons, protons, or other types of radiation.
- the radiation source 102 may include a metallic target configured to produce a beam of x-rays or photons upon impingement of electrons.
- the principle of dosimetric projection of a collimator geometry and dose calculation method described in this disclosure can also be practiced in other types of radiation system such as ion, electron, carbon, or proton systems.
- the radiation system 100 may include various collimating devices or components configured to modify the size, shape, fluence, and other characteristics of the beam.
- a primary collimator 106 adjacent to the source 102 generally limits the extent of the beam 104 as it travels away from the source 102 toward a subject such as a patient 108.
- motorized secondary collimators or collimation jaws 109 may be included to define the field size.
- a multileaf collimator (MLC) 110 may be disposed between the source 102 and the patient 108 to shape the beam, as indicated by the shaped field 112 shown in FIG. 1A.
- the MLC 110 may be rotated about the central beamline or axis 114, placing the MLC 110 in various orientations.
- the radiation system 100 may optionally include a flattening filter 116 to provide a uniform beam profile.
- the radiation system does not include any flattening filter or is flattening-filter-free (FFF) to enhance dose rates for treatment.
- Ion chamber 118 monitors the parameters of the beam from the source 102.
- the source 102, primary collimator 106, secondary collimators 109, MLC 110, and other devices or components may be enclosed in the treatment head 120, which can be rotated by a gantry (not shown) about an axis such as a horizontal axis.
- the system 100 can deliver radiation to a target in the subject 108 from various beam angles.
- the shape, size, and/or intensity of the beam 104 can be adjusted or dynamically adjusted by the MLC 110 as the beam angle is stepped or swept around the target.
- the operation of the source 102, MLC 110, and other devices of the radiation system 100 can be controlled by a control system 122.
- the control system 122 may comprise a computer processor which receives and executes a treatment plan generated by a treatment planning system 1100 to be described in greater detail below.
- the MLC 110 may be a single level MLC as shown in FIGS. 1 and 1A or a multi level MLC.
- the term“MLC” or“multileaf collimator” refers to a collection of a plurality of beam-blocking leaves 111 each of which can be longitudinally moved in and out of a beam to modify one or more parameters of the beam such as the beam shape, size, energy, or fluence etc.
- Each beam-blocking leaf 111 may be driven by a motor with a lead screw or other suitable means.
- the beam-blocking leaves 111 may be arranged in pairs. The beam blocking leaves of each pair may be brought in contact or retracted from each other to close or open a path for a radiation beam to pass through the MLC.
- the beam-blocking leaves may be arranged in opposing banks and supported by a frame, box, carriage or other support structure, which has features allowing the individual beam-blocking leaves to extend into and retract from the beam.
- the frame, box, carriage or other support structure can be further moved or translated in addition to the individual leaf travel.
- the beam-blocking leaves 111 may have various leaf tip profiles or end portion configurations, and leaf edge profiles or leaf side configurations.
- the beam-blocking leaves 111 of the MLC 110 may have a flat front-end surface and a flat side surface.
- a beam-blocking leaf with a flat front- end surface may be shown to have a straight line orthogonal to the leaf longitudinal moving direction and two right angles at each side of the straight line.
- the term“top view” refers to a view observed from the source or in a direction along the central beamline.
- the beam-blocking leaves of the MLC 110 may have a curved front- end surface.
- a beam-blocking leaf with a curved front-end surface may be shown to have a curved line with a radius and two parallel lines at either side of the curved line.
- a leaf with a curved front-end surface may be shown to have a straight line orthogonal to the leaf longitudinal moving direction and two right angles on each side of the straight line.
- the term“square shape” may be used to describe a leaf tip profile which, in a top view, has a straight line orthogonal to the leaf longitudinal moving direction and two right angles at each side of the straight line.
- the term “non-square shape” may be used to describe any leaf tip profile which, in a top view, does not have a square shape.
- the non-square shape in a top view may include a curved or elliptic shape or a shape of a chamfer leaf as will be described in greater detail below.
- the MLC leaf tip may have a circular or non-circular profile.
- the MLC leaf side may have a tongue-and-groove feature, a tilted surface, or some other type of profile.
- the beam-blocking leaves 111 of the MLC 110 may have a tip profile that includes a combination of a curved surface portion and a beveled or flat surface portion at each side of the curved surface portion.
- the term“chamfer leaf” may be used herein to refer to a leaf which includes a combination of a curved end surface portion and beveled or flat end surface portions at either side of the curved end surface portion.
- the end portion or tip of a chamfer leaf may be shown to have a straight middle line section orthogonal to the leaf longitudinal moving direction and beveled or angled lines sections at each side of the straight middle line section.
- the MLC 110 may be a multi-level MLC.
- the MLC 110 may include a first MLC in a first level distal to the source 102 and a second MLC in a second level proximal to the source 102.
- the the first and second MLCs may be arranged such that the moving direction of individual beam-blocking leaves of the first and second MLCs are generally in parallel.
- the first and second MLCs may be arranged such that the moving direction of the beam-blocking leaves of the first MLC is non parallel e.g. perpendicular to the moving direction of the beam-blocking leaves of the second MLC.
- the first and second MLCs may be arranged such that the leaves of the second MLC may laterally offset the leaves of the first MLC in a top view or as viewed in a direction from the source.
- the MLC may be absent or combined with other collimating or beam shaping devices, including, but not limited to, collimation jaws, cones, wedges, and filters.
- the MLC 110 may be configured to form an aperture defining a shaped field 112 approximating the target geometry at the isocenter plane.
- the MLC 110 may be configured to define differently shaped fields at different MLC orientations and/or beam angles, and the doses of multiple fields may be summed to build up a desired dose distribution in the target.
- Radiation may be delivered intermittently or statically wherein the MLC leaves are in positions while radiation is being delivered. Radiation may also be delivered dynamically wherein the MLC leaves are moving or the MLC is rotating while radiation is being delivered.
- the aperture of the MLC is formed substantially small for small-field radiotherapy such as stereotactic radiosurgery (SRS).
- SRS stereotactic radiosurgery
- dose deposited within the target or treatment volume is calculated or predicted using a suitable dose calculation algorithm.
- various factors are considered including the composition of tissue and characteristics of the radiation beam from the source.
- the characteristics of the beam exiting the MLC need be determined and provided to the dose calculation algorithm for computing or predicting dose deposited in the target or treatment volume.
- FIG. 2 illustrates a method or model for calculating radiation transport through a collimating device such as an MLC.
- a 3D collimator geometry model is constructed and computer software based on the 3D model created.
- the 3D model for calculating radiation transport can be very accurate because the structural details of the 3D collimator geometry, including the fine features in the region near the aperture, are taken into account in calculation.
- constructing a 3D model and computer software based on the 3D model in full detail are extremely time-consuming. Calculation of radiation transport based on a fully 3D model is computationally too slow to be used in a treatment planning system in a clinical environment, especially in treatment plan optimization.
- FIG. 3 illustrates another method or model for calculating radiation transport through a collimating device such as an MLC.
- a 3D collimator geometry is flattened or projected along the beamline into a 2D geometry and computer software based on the 2D geometry is created.
- the 2D model for calculating radiation transport is fast and suitable for many clinical cases.
- conventional 2D models for calculating radiation transport are less accurate, especially around the edges of the aperture because the fine features of the aperture edges are not modeled accurately, which is especially relevant when the aperture is small for e.g. stereotactic radiosurgery.
- an actual MLC leaf end portion having a curve profile in a side view may be shown as having a square shape when the MLC leaf is flattened or projected along the beamline into a 2D geometry.
- the fine features of the actual MLC leaf tip profile which would affect attenuation or transmission of the radiation beam, are not taken into account in conventional 2D modeling.
- the fine features of chamfer leaves described above are not taken into account in conventional 2D modeling for calculating radiation transport.
- the field size may be comparable to the projected dimension of the MLC leaf tip profile. Therefore, the fine features of the aperture edges have to be taken into account to accurately calculate radiation transport.
- the geometry of the collimator 410 is flattened or projected along the beamline from 3D to 2D on a plane, and a calculation algorithm is created based on the 2D geometry. While the collimator geometry is flattened from 3D to 2D, dosimetric opacity values or dosimetric thicknesses of the collimator parts are calculated based on the actual 3D geometry of the collimator. The dosimetric opacity values or thicknesses can be encoded in the calculation algorithm based on the 2D geometry, as schematically shown in a range of shades around the edges of the aperture in FIG. 4.
- the calculation algorithm model is in 2D, the computation is fast. Further, because the dosimetric opacity values or thicknesses are obtained based on the actual 3D shape of the collimator, the calculation of radiation transport, including the calculation at the edges of the aperture, is as accurate or almost as accurate as the 3D model.
- FIG. 5 shows that the dosimetric opacity values or thicknesses depend on both the 3D shape of the collimator parts (e.g. MLC leaf tip profiles) and the positions of the collimator parts within the radiation field.
- a beam-blocking leaf 412 of the MLC 410 at a position shown in FIG. 5 exhibits a broader range of shades than does the same leaf 412 at a different position shown in FIG. 4. This can be explained as that the dosimetric opacity value or thickness is a function of the length or dimension of the collimator part along the beam path.
- the term’’dosimetric projection may be used hereafter to refer to a method of obtaining a dosimetric opacity value or thickness of a collimator based on the 3D shape or geometry of the collimator.
- the fluence and energy of the radiation produced by the source and direction at all points in a reference plane can be calculated.
- the characteristics of the radiation source can be provided as an initial spatial photon spectrum describing the fluence, energy and directions in a beam reference plane.
- the 3D geometry of a collimator may be predetermined by the design of the collimator such as in the case of a custom-molded collimation block, which may in turn be dictated by the geometry of the target to be treated.
- the 3D geometry of the MLC may be determined by the positioning of a set of beam-blocking leaves, which may readily change depending on the treatment plan.
- the beam-blocking leaves of an MLC may have various leaf tip profiles e.g. curved or rounded leaf end surfaces, whose beam’s eye view or beam-blocking effect may change depending on the positions of the beam-blocking leaves.
- a same geometry of the MLC may have different collimating or attenuation effect especially at the edges of the aperture when the MLC is in different orientations because the tips of the beam-blocking leaves defining the aperture change the angles with respect to the source.
- the 3D geometry of the MLC may be defined by specifying the spatial positions of the beam-blocking leaves and their angles with respect to the central beamline, the physical dimensions, and the leaf tip profiles of the beam-blocking leaves, and other variables.
- dosimetric projection of the collimator is calculated.
- Dosimetric projection of a collimator involves in projection of the collimator geometry from 3D to 2D along the beamline wherein the dimension of the collimator along the beamline is represented or encoded with a dosimetric opacity value.
- dosimetric opacity refers to a quality of the collimator in attenuating radiation. The value of dosimetric opacity is obtained by applying the dosimetric projection method on the given collimator.
- dosimetric opacity depends on the material the collimator is made of, the physical dimension of the collimator or the“thickness” of the material along which the radiation passes, and the position of the collimator in relation to the beam centerline. Therefore, term “dosimetric opacity” may be used interchangeably with the term“dosimetric thickness.”
- the dosimetric opacity value of a given collimator part can be calculated using a suitable method based on the realistic 3D geometry of the collimator.
- local dosimetric projection of the collimator geometry is performed.
- dosimetric opacity values at locations adjacent to the aperture defined by the collimator are calculated.
- a plurality of points or surfaces on an MLC leaf end portion which collectively may represent a boundary of the MLC leaf end or indicate the leaf tip profile, are specified.
- the dosimetric opacity values at these bounding points or surfaces may be calculated by for example ray-tracing or any other suitable analytical methods. Interpolating of dosimetric opacity values between the bounding points and surfaces or analytical approximation may be performed.
- FIG. 7 illustrates an example collimator part (an MLC beam-blocking leaf) and a method of calculating dosimetric opacity values according to embodiments of the disclosure.
- the horizontal axis indicates leaf projection positions and the vertical axis represents normalized leaf thicknesses.
- Curve 702 depicts the physical geometry or real profile of the collimator part at a position relative to the source.
- the example MLC leaf has a curved leaf tip profile.
- the dots on the curve 702 may represent specified bounding points or surfaces for calculating dosimetric opacity values.
- Curve 704 shows a calculated thickness of the collimator part.
- curve 704 substantially superimposes over curve 702, indicating that the thickness of the collimator part can be accurately calculated using an analytical or numerical method such as ray-tracing and interpolating known in the art.
- the dosimetrically projected thickness can be defined, for example, as a function approximating the ray-traced curve 704.
- the dosimetric projection can be defined in such a way that the approximation remains valid and takes into account the dependence of the dosimetric thickness on the position of the MLC leaf, as indicated by the illustrations in FIG. 4 and FIG. 5.
- the dosimetrically projected shape or thickness is indicated by curve 706.
- FIG. 8 is a flowchart illustrating example steps for calculating dosimetric projection of a collimator.
- the collimator thickness along a ray from pointA at the source to point B on the projection plane is calculated.
- the calculation of collimator thickness can be conducted numerically e.g. by ray-tracing or analytically.
- Analytical calculation can be conducted wherein the collimator is represented by a set of suitable geometric primitives, e.g., circles and lines in two dimensions or cylinders and planes in three dimensions.
- the collimator position can be fixed.
- the calculation of 802 is repeated for different points B’ on the projection plane.
- the different points B’ can be those located adjacent to the aperture defined by the collimator.
- a thickness function from the results of 806 is constructed. This can be done by e.g. interpolating between the points calculated at 806 or by an analytical approximation.
- the function based on interpolation can take the form of a pre-calculated look- up table coupled with an interpolation algorithm such as multilinear interpolation or other suitable method known in the art.
- the function based on analytical approximation can be obtained, e.g., by numerically fitting an analytical function such as a multidimensional polynomial to the points calculated at 806.
- the analytical approximation can also be obtained from the exact analytical solution, if such is available, by approximating the solution for greater numerical efficiency by using numerical computation methods known in the art.
- the thickness function is converted into dosimetric opacity.
- the radiation source characteristics such as the primarily particle type and energy spectrum and radiation transport physics can be considered in obtaining the opacity.
- This step may be implementation dependent - it can be considered to be part of the radiation transport calculation (608 in FIG. 6) if the algorithm for radiation transport uses geometrical data (thickness) instead of transmission/opacity for the collimator components.
- FIG. 9 is an image showing a result of dosimetric projection, the dosimetric opacity using a single MLC leaf as an example.
- the vertical axis indicates dosimetric opacity value
- the horizontal axes indicate the leaf tip poistion and the distance from the leaf tip respectively.
- the units of the axes in FIG. 9 are arbitrary for illustration purpose.
- FIG. 9 shows that for an exemplary MLC leaf, the dosimetric opacity can be represented by a function of two variables, the leaf position and distance from the leaf tip.
- the method may proceed to calculate radiation transport using the dosimetric opacity values or thicknesses.
- Various algorithms have been developed to calculate beam transport through collimating devices. These methods are generally known in the art and therefore their detailed description is omitted to avoid obscuring the description of various embodiments of this disclosure.
- an algorithm for calculating radiation transport takes into account factors including primary source, extra focal source, electron contamination, photon scattering etc.
- the dosimetric opacity values obtained by local dosimetric projection described above are used as input for an algorithm based on a 2D geometry model to calculate radiation transport.
- the beam exiting the collimator can be provided as a spatial fluence spectrum in a reference plane.
- the method may proceed to calculate dose deposition in the target within the treatment volume using a treatment planning software including a dose calculation algorithm.
- a treatment planning software including a dose calculation algorithm.
- a method and algorithm based on pencil beam convolution-superposition may be used.
- Other suitable methods and algorithms based on e.g. Boltzmann transport equation (BTE) can also be used.
- BTE Boltzmann transport equation
- the present disclosure is not limited to a particular dose calculation algorithm. Examples of commercially available dose calculation software include Acuros ® XB dose calculation algorithm and the Anisotropic Analytical Algorithm (AAA), both available from Varian Medical Systems, Inc. of Palo Alto, CA.
- AAAA Anisotropic Analytical Algorithm
- the spatial fluence spectrum that has been calculated in block 608 using the dosimetric opacity values or thicknesses at the boundary of the collimator defining the aperture are provided to the dose calculation algorithm.
- Embodiments of this disclosure provide for an improved method of calculating radiation transport or transmission through a collimator using a 2D or quasi-2D model to achieve a calculation accuracy that approximates a 3D model.
- the disclosed method does not require calculation of radiation transmission in 3D geometry every time the result is needed or over the whole spatial extent of the collimator.
- the disclosed method may use dosimetric projection described above, to pre-calculate radiation transmission just once in 3D or in an approximate 3D geometry, pre-calculate only locally around the aperture defined by the collimator, and pre-calculate for different collimator positions (field sizes). Using the pre-calculated values, a 2D or quasi-2D geometry or computational model, algorithm, function, etc.
- the 2D or quasi-2D model of the dosimetric opacity can be implemented into the fluence calculation module of the treatment planning software in calculation of dose desposition.
- FIG. 10 provides a comparison of accuracies of dose calculation using different methods.
- the plots shown in FIG. 10 depict dose distribution for a small field (e.g. 5 x 5 millimeters).
- Plots 1002 and 1004 show dose distributions calculated with Eclipse ® treatment planning software and Acuros ® XB dose deposition algorithm, with plot 1002 being calculated using a conventional 2D MLC modeling and plot 1004 being calculated using the dosimetric projection method according to embodiments of the disclosure.
- plot 1006 shows a dose distribution calculated with Geant4 Monte Carlo using a full 3D MLC modeling.
- FIGS. 1-10 Various embodiments of dosimetric projection of a collimator geometry and dose calculation methods are described in connection with FIGS. 1-10. It will be appreciated that more or fewer steps, actions, or processes may be incorporated into the methods without departing from the scope of the disclosure. No particular order is implied by the arrangement of blocks shown and described herein. It further will be appreciated that the methods described in conjunction with FIGS. 1-10 may be embodied in a computer-based treatment planning system which performs dosimetric projection and dose calculation methods according to this disclosure. The invention may be in the form of a computer product comprising a computer- readable medium storing or carrying instructions which, when executed by a computer processor, cause the computer processor to perform the methods described in the disclosure.
- the instructions may be implemented as software code to be executed by a processor using any suitable computer language such as, Java, C++ or Perl using, e.g. conventional or object- oriented techniques.
- the computer-readable medium may include any suitable medium that is capable of storing or encoding a sequence of instructions for execution by the computer processor and that causes the computer processor to perform any one of the methodologies of the present invention.
- the computer-readable medium shall accordingly include, but not be limited to, solid-state memories, optical and magnetic disks. Examples of computer-readable medium include volatile and nonvolatile, removable and non-removable media for storage of computer-readable instructions.
- a computer-readable medium encoded with the program code may be packaged with a compatible device or provided separately from other devices e.g. via Internet download. Further, any such computer-readable medium may reside on or within a computer product e.g. a hard drive, a CD, or an entire computer system.
- FIG. 11 schematically illustrates a treatment planning system 1100 upon which embodiments of the disclosure can be performed.
- the treatment planning system 1100 includes a computer system 1110, which generally comprises a processor 1112, memory 1114, a user interface 1116, and a network interface 1118, each of which is coupled to a system bus 1120.
- the processor 1112 may include a central processing unit (CPU) that is generally known in the art, such as an INTEL ® processor or an AMD ® processor, or a graphical processing unit (GPU), such as an NVIDIA ® GPU, or other type of processing unit.
- CPU central processing unit
- GPU graphical processing unit
- the processor 1112 may retrieve and execute computer-executable instructions from the memory 1114, which may cause the processor 1112 to perform any of the methods and/or steps according to the embodiments of this disclosure described above.
- the memory 1114 may include any one of or a combination of volatile memory elements and nonvolatile memory elements.
- the memory 1114 may include a random-access memory (RAM) or other dynamic storage device for storing information and instructions to be executed by the processor 1112, and for storing temporary variables or other intermediate information during execution of instructions by the processor 1112.
- RAM random-access memory
- the memory 1114 may also include read-only memory (ROM) or other static storage device for storing static information and instructions for the processor 1112.
- ROM read-only memory
- the memory 1114 may further include a data storage device such as a magnetic disk or optical disk, for storing information and instructions.
- the memory 1114 e.g.
- a non-transitory computer-readable medium may comprise programs (logic) for operating the computer system and for performing applications including dosimetric projection and dose calculation as described above, or other treatment planning applications.
- the memory 1114 may include a database storing any information that can be selected by a user, such as a radiation oncologist or radiation therapist.
- the network interface 1118 allows the computer system to communicate with other devices or systems over a communication network 1122 such as the Internet or an intranet (e.g., a local area network).
- the network interface 1118 may include a Wi-Fi interface, Ethernet interface, Bluetooth interface, or other wireless or wired interfaces.
- the network interface 1118 allows the computer system 1110 to receive and send electrical, electromagnetic or optical signals that carry data streams representing various types of information. For instance, the network interface 1118 may allow the computer system 1110 to receive data streams representing software programs for treatment planning over the communication network 1122.
Abstract
Description
Claims
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US16/359,562 US10940330B2 (en) | 2019-03-20 | 2019-03-20 | Dosimetric projection of collimator geometry |
PCT/IB2020/051444 WO2020188372A1 (en) | 2019-03-20 | 2020-02-20 | Dosimetric projection of collimator geometry |
Publications (1)
Publication Number | Publication Date |
---|---|
EP3941576A1 true EP3941576A1 (en) | 2022-01-26 |
Family
ID=69743627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP20709334.5A Pending EP3941576A1 (en) | 2019-03-20 | 2020-02-20 | Dosimetric projection of collimator geometry |
Country Status (3)
Country | Link |
---|---|
US (1) | US10940330B2 (en) |
EP (1) | EP3941576A1 (en) |
WO (1) | WO2020188372A1 (en) |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9931521B2 (en) | 2014-03-25 | 2018-04-03 | Varian Medical Systems, Inc. | Multi level multileaf collimator leaf tip shape effects and penumbra optimization |
-
2019
- 2019-03-20 US US16/359,562 patent/US10940330B2/en active Active
-
2020
- 2020-02-20 WO PCT/IB2020/051444 patent/WO2020188372A1/en unknown
- 2020-02-20 EP EP20709334.5A patent/EP3941576A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
WO2020188372A1 (en) | 2020-09-24 |
US20200298018A1 (en) | 2020-09-24 |
US10940330B2 (en) | 2021-03-09 |
CN113597326A (en) | 2021-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3515557B1 (en) | Controlling and shaping the dose distribution outside treatment targets in external-beam radiation treatments | |
EP3515556B1 (en) | Generating time-efficient treatment field trajectories for external-beam radiation treatments | |
EP3459595B1 (en) | Vmat treatment planning using multicriteria optimization and a progressive optimization scheme | |
US8913717B2 (en) | Dose calculation method for multiple fields | |
US7801270B2 (en) | Treatment plan optimization method for radiation therapy | |
US8085899B2 (en) | Treatment planning system and method for radiotherapy | |
US20100012859A1 (en) | Method For Treating A Target Volume With A Particle Beam And Device Implementing Same | |
US20210038914A1 (en) | System and method for radiotherapy treatment planning and delivery | |
US10500417B2 (en) | Compensating for leakage radiation in MLC modulated treatment | |
Otto et al. | Enhancement of IMRT delivery through MLC rotation | |
US8280003B2 (en) | Method for calculating head scatter phase space for radiation treatment using a multi-leaf collimator with dynamic jaws | |
EP3946584B1 (en) | Using isodose surfaces for optimizing dose distribution in radiation treatment planning | |
US20210138267A1 (en) | Expansion of search space for multicriteria optimization in radiation treatment planning | |
US10940330B2 (en) | Dosimetric projection of collimator geometry | |
CN113597326B (en) | Dosimetry projection of collimator geometry | |
JP7357059B2 (en) | Systems and methods for passive ion radiotherapy planning | |
US20240139543A1 (en) | Pre-calculated collimator model for dose calculation source modeling | |
MacDonald | Development and Implementation of Trajectory Optimization Technologies for Cranial Stereotactic Radiation Therapy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: UNKNOWN |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE |
|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
|
17P | Request for examination filed |
Effective date: 20210916 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
DAV | Request for validation of the european patent (deleted) | ||
DAX | Request for extension of the european patent (deleted) | ||
RAP3 | Party data changed (applicant data changed or rights of an application transferred) |
Owner name: SIEMENS HEALTHINEERS INTERNATIONAL AG |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
17Q | First examination report despatched |
Effective date: 20240111 |