CN117224864B - Method and device for determining radiotherapy dosage, storage medium and electronic equipment - Google Patents
Method and device for determining radiotherapy dosage, storage medium and electronic equipment Download PDFInfo
- Publication number
- CN117224864B CN117224864B CN202311503258.9A CN202311503258A CN117224864B CN 117224864 B CN117224864 B CN 117224864B CN 202311503258 A CN202311503258 A CN 202311503258A CN 117224864 B CN117224864 B CN 117224864B
- Authority
- CN
- China
- Prior art keywords
- dose
- reference point
- determining
- target
- distribution map
- 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.)
- Active
Links
- 238000001959 radiotherapy Methods 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000008569 process Effects 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 5
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 15
- 230000000694 effects Effects 0.000 description 5
- 239000007787 solid Substances 0.000 description 4
- 230000003902 lesion Effects 0.000 description 3
- 238000002729 3-dimensional conformal radiation therapy Methods 0.000 description 2
- 238000011455 3D conformal radiation therapy Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Abstract
The invention discloses a method and a device for determining radiotherapy dosage, a storage medium and electronic equipment. The method comprises the following steps: for each of a plurality of beams generated by the radiotherapy apparatus, determining a target plane corresponding to the beam; determining a reference point matched with the beam from the target plane according to the dose of the beam on the target plane, wherein the ratio of the dose at the reference point to the maximum dose of the beam on the target plane is greater than or equal to a preset value; a total dose profile is determined from the dose profile of each beam, the dose at the reference point to which each beam is matched, and the prescribed dose, wherein the total dose profile is used to characterize the desired dose distribution during radiation therapy. The invention solves the technical problem that the accuracy of the determined radiotherapy dosage is low because the relative error at the fixed point is difficult to control when the radiotherapy dosage is determined according to the dosage of the beam at the fixed point in the related technology.
Description
Technical Field
The present invention relates to the field of data information processing, and in particular, to a method and apparatus for determining a radiotherapy dose, a storage medium, and an electronic device.
Background
Three-dimensional conformal radiotherapy (3D-Conformal Radiation Therapy, 3D-CRT) is a high-precision radiotherapy technology, which adopts a virtual three-dimensional reconstruction technology to project a lesion part of a patient and normal tissues around the lesion part into a three-dimensional image, and sets a plurality of radiation fields in different directions, and adjusts multi-leaf gratings according to the shapes of target areas in the three-dimensional image under different visual angles to change the shapes of the radiation fields, so that the shapes of dose distribution of a high-dose area are consistent with the shapes of the target areas in a three-dimensional direction, thereby being convenient for better protecting the normal tissues and improving the radiotherapy effect.
Since the dose distribution in three-dimensional conformal radiation therapy directly affects the effect of radiation therapy, how to determine an accurate dose distribution is a matter of great concern. Currently, in the related art, a uniform fixed point is generally selected as a reference point after the beam is set, and then a final required dose distribution during radiotherapy is determined according to the dose of the beam at the reference point. However, since the reference point is fixed, for some beams the relative error at the reference point is large, resulting in problems with low accuracy of the finally determined radiotherapy dose.
In view of the above problems, no effective solution has been proposed at present.
Disclosure of Invention
The embodiment of the invention provides a method, a device, a storage medium and electronic equipment for determining radiotherapy dosage, which at least solve the technical problem that the accuracy of the determined radiotherapy dosage is low because the relative error at a fixed point is difficult to control when the radiotherapy dosage is determined according to the dosage of a beam at the fixed point in the related technology.
According to an aspect of an embodiment of the present invention, there is provided a method for determining a radiotherapy dose, including: for each of a plurality of beams generated by the radiotherapy apparatus, determining a target plane corresponding to the beam; determining a reference point matched with the beam from the target plane according to the dose of the beam on the target plane, wherein the ratio of the dose at the reference point to the maximum dose of the beam on the target plane is greater than or equal to a preset value; a total dose profile is determined from the dose profile of each beam, the dose at the reference point to which each beam is matched, and the prescribed dose, wherein the total dose profile is used to characterize the desired dose distribution during radiation therapy.
Further, the target plane is a plane passing through the isocenter of the radiotherapy apparatus and perpendicular to the isocenter of the beam.
Further, the method for determining the radiotherapy dosage further comprises the following steps: determining a dose of the beam at the isocenter, resulting in a first dose; determining a maximum dose of the beam on the target plane to obtain a second dose; a reference point for beam matching is determined from the target plane based on a proportional relationship between the first dose and the second dose.
Further, the method for determining the radiotherapy dosage further comprises the following steps: determining the isocenter as a reference point in case the ratio between the first dose and the second dose is greater than or equal to a preset value; in case the ratio between the first dose and the second dose is smaller than a preset value, a reference point is determined from points on the target plane other than the isocenter.
Further, the method for determining the radiotherapy dosage further comprises the following steps: for each beam, adjusting the dose distribution map of the beam according to the dose and the prescription dose at the reference point matched with the beam to obtain an adjusted dose distribution map; a total dose distribution map is determined from the adjusted dose distribution maps of the plurality of beams.
Further, the method for determining the radiotherapy dosage further comprises the following steps: determining an initial total dose distribution map based on the weights of the individual beams and the adjusted dose distribution maps of the plurality of beams; acquiring a target reference point, wherein the target reference point is in an initial total dose distribution map; and adjusting the initial total dose distribution map according to the dose at the target reference point and the prescription dose to obtain a total dose distribution map.
Further, the method for determining the radiotherapy dosage further comprises the following steps: acquiring a preset reference point, and determining the preset reference point as a target reference point; according to the dose of each point in the initial total dose distribution map, sequencing each point in the initial total dose distribution map to obtain a target sequence, and determining the Nth point in the target sequence as a target reference point, wherein N is a positive integer.
According to another aspect of the embodiment of the present invention, there is also provided a radiotherapy dose determining apparatus, including: a first determining module, configured to determine, for each of a plurality of beams generated by the radiotherapy apparatus, a target plane corresponding to the beam; a second determining module, configured to determine a reference point for beam matching from the target plane according to a dose of the beam on the target plane, where a ratio between the dose at the reference point and a maximum dose of the beam on the target plane is greater than or equal to a preset value; and a third determining module, configured to determine a total dose distribution map according to the dose distribution map of each beam, the dose at the reference point matched by each beam, and the prescribed dose, where the total dose distribution map is used to characterize the dose distribution required in the radiotherapy process.
According to another aspect of the embodiments of the present invention, there is also provided a computer readable storage medium having a computer program stored therein, wherein the computer program is arranged to perform the above-described method of determining a radiation dose when run.
According to another aspect of an embodiment of the present invention, there is also provided an electronic device including one or more processors; and a memory for storing one or more programs which, when executed by the one or more processors, cause the one or more processors to implement a method for running the program, wherein the program is configured to perform the method of determining a radiation therapy dose as described above when run.
In the embodiment of the invention, a reference point is determined according to the dose of the beam on a target plane, the required dose distribution in the radiotherapy process is determined according to the reference point of each beam, the target plane corresponding to the beam is determined for each beam in a plurality of beams generated by the radiotherapy equipment, and then the reference point matched with the beam is determined from the target plane according to the dose of the beam on the target plane, so that a total dose distribution map is determined according to the dose distribution map of each beam, the dose at the reference point matched with each beam and the prescription dose, wherein the ratio between the dose at the reference point and the maximum dose of the beam on the target plane is larger than or equal to a preset value, and the total dose distribution map is used for representing the required dose distribution in the radiotherapy process.
In the above process, for each beam, by determining the beam matching reference point from the target plane according to the ratio between the dose of the point on the target plane and the maximum dose, on the one hand, the independent determination of the reference point of each beam is realized, and on the other hand, the relatively large dose at the reference point of each beam is ensured, so that the relative error at the reference point of each beam can be controlled within a certain range, and the accuracy of the determined radiotherapy dose can be effectively improved when the total dose distribution map is determined according to the dose distribution map of each beam, the dose at the reference point of each beam matching and the prescription dose.
Therefore, the scheme provided by the application achieves the purposes of determining the reference point according to the dose of the beam on the target plane and determining the required dose distribution in the radiotherapy process according to the reference points of the beams, thereby realizing the technical effect of improving the accuracy of the determined radiotherapy dose, and further solving the technical problem that the accuracy of the determined radiotherapy dose is low due to the fact that the relative error at the fixed point is difficult to control when the radiotherapy dose is determined according to the dose of the beam at the fixed point in the related technology.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiments of the invention and together with the description serve to explain the invention and do not constitute a limitation on the invention. In the drawings:
FIG. 1 is a schematic illustration of an alternative method of determining a radiation therapy dose according to an embodiment of the invention;
FIG. 2 is a schematic illustration of an alternative beam in accordance with an embodiment of the present invention;
FIG. 3 is a schematic illustration of the relationship between an alternative isocenter and a dose profile of a beam according to an embodiment of the present invention;
FIG. 4 is a schematic illustration of another relationship between an alternative isocenter and a dose profile of a beam according to an embodiment of the present invention;
fig. 5 is a schematic view of an alternative radiotherapy dose determination device according to an embodiment of the present invention;
fig. 6 is a schematic diagram of an alternative electronic device according to an embodiment of the invention.
Detailed Description
In order that those skilled in the art will better understand the present invention, a technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in which it is apparent that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present invention without making any inventive effort, shall fall within the scope of the present invention.
It should be noted that the terms "first," "second," and the like in the description and the claims of the present invention and the above figures are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used may be interchanged where appropriate such that the embodiments of the invention described herein may be implemented in sequences other than those illustrated or otherwise described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
Example 1
In accordance with an embodiment of the present invention, there is provided an embodiment of a method of determining a radiation therapy dose, it being noted that the steps illustrated in the flowchart of the figures may be performed in a computer system, such as a set of computer executable instructions, and that although a logical sequence is illustrated in the flowchart, in some cases the steps illustrated or described may be performed in a different order than what is illustrated herein.
Fig. 1 is a schematic diagram of an alternative method of determining a radiation therapy dose according to an embodiment of the invention, as shown in fig. 1, the method comprising the steps of:
step S101, for each of a plurality of beams generated by the radiotherapy apparatus, determining a target plane corresponding to the beam.
Alternatively, an electronic device, an application system, a server, or the like may be used as an execution subject of the present application, and in this embodiment, a target processing system is used as an execution subject to determine the target plane corresponding to each beam.
Optionally, in this embodiment, in order to determine the dose distribution in the three-dimensional conformal radiotherapy, the target processing system may first acquire a target image, where the target image is a three-dimensional image corresponding to a lesion part of a patient and a surrounding normal tissue, and the target image may be an image such as CT (Computed Tomography, electronic computed tomography) \mri (Magnetic resonance imaging ) \cbct (cone beam CT) \ Positron Emission Tomography, positron emission tomography). Further, the target processing system may add a plurality of beams for irradiating the region in the target image in response to the beam setting operation by the user, and adjust the beams according to the parameter information such as the beam angle, the multileaf grating shape, and the like in the beam setting operation.
Alternatively, after each beam is adjusted, the target processing system may determine the target plane to which each beam corresponds. Wherein one beam corresponds to one target plane, which for each beam is a plane through the beam, and further the target plane may be a plane through the beam and through an isocenter of the radiotherapy apparatus, and further the target plane may be a plane through the isocenter of the radiotherapy apparatus and perpendicular to the isocenter of the beam, i.e. the target plane is an isocenter plane. In the radiotherapy apparatus, the reference axes of the various movements move around a common center point, and the beam passes through a minimum sphere centered on this point, which is the isocenter point, and the isocenter axis is the line between the center of the light source emitting the beam and the isocenter point. For example, FIG. 2 is a schematic illustration of an alternative beam according to an embodiment of the invention, shown in FIG. 2Go out to have B 1 、B 2 B, B 3 The three beams, as shown in fig. 2, are unique in isocenter for the beams emitted by the same radiotherapy apparatus, and the isocenter corresponds to the beams one by one.
Step S102, determining a reference point matched with the beam from the target plane according to the dose of the beam on the target plane, wherein the ratio of the dose at the reference point to the maximum dose of the beam on the target plane is greater than or equal to a preset value.
Currently, a uniform fixed point is generally selected as a reference point in the related art, however, this reference point is not necessarily located in all fields, and thus the relative error at the reference point is large for some beams. Alternatively, the above principle is described taking the reference point as an isocenter as an example. According to the dose distribution curve, the dose distribution of the beam on the isocenter plane can be divided into a highlight area, a penumbra area and a shadow area, wherein the highlight area is defined in an area of 80% -100% of the dose distribution curve on the isocenter plane, the penumbra area is defined in an area of 20% -80% of the dose distribution curve on the isocenter plane, and the shadow area is defined in an area of 0% -20% of the dose distribution curve on the isocenter plane. Fig. 3 is a schematic diagram of a relationship between an optional isocenter and a dose distribution diagram of a beam according to an embodiment of the present invention, as shown in fig. 3, a dashed curve in fig. 3 represents a dose size, a dashed square frame represents a beam, three solid circles in fig. 3 represent dose distribution diagrams, wherein different solid circles represent different dose distribution intensities, and a dose distribution intensity in an innermost region is higher than a dose distribution intensity in a region between a middle circle and an innermost circle, a dose distribution intensity in a region between a middle circle and an innermost circle is higher than a dose distribution intensity in a region between an outermost circle and a middle circle, when the isocenter is in a field, an isocenter (i.e., an intersection point of the isocenter and the isocenter plane in fig. 3) is in a highlight region or is adjacent to the highlight region, a dose value at the isocenter is relatively high, and a relative error= (|measured value-true value/true value) is 100% at this time, for example, and a relative error at the isocenter is 1.25 is small, for example, and a relative error at the isocenter is 0.25. Fig. 4 is a schematic diagram of another relationship between an optional isocenter and a dose distribution diagram of a beam according to an embodiment of the present invention, as shown in fig. 4, the dashed curve in fig. 4 represents a dose size, the dashed square frame represents the beam, three solid circles in fig. 4 represent the dose distribution diagram, wherein different solid circles represent different dose distribution intensities, and the dose distribution intensity in the innermost region is higher than that in the region between the middle circle and the innermost circle, the dose distribution intensity in the region between the middle circle and the innermost circle is higher than that in the region between the outermost circle and the middle circle, and when the isocenter is outside the shot, the isocenter (i.e. the intersection of the isocenter and the isocenter plane in fig. 4) is in or adjacent to the shadow region, the dose value at the isocenter is relatively low, and the relative error of the dose value at the isocenter is relatively large, for example, the measured value is 0.02, and the true value is 0.01, and the relative error is 1.
Further, since the dose distribution finally required in the radiotherapy process is determined according to the dose at the reference point, if the relative error of the dose value at the reference point is large, the deviation between the finally determined dose distribution and the actual dose distribution is large, which results in the problem of low accuracy of the finally determined radiotherapy dose.
Thus, in this embodiment, a reference point of relatively large dose may be selected from the target plane to which the beam corresponds to control that the relative error at the reference point is small. Alternatively, the target processing system may employ a dose calculation model such as a Monte Carlo or pencil beam, determine a dose distribution map of the beam based on the parameter information of the beam, and then determine the dose of the beam on the target plane from the dose distribution map based on the coordinate information of the target plane. Further, the target processing system may determine a maximum dose of the beam on the target plane, thereby determining the reference point based on a ratio between the dose and the maximum dose at each point on the target plane. The preset values may be freely set according to different actual requirements, which is not limited in this embodiment, and in addition, in the case that there are multiple points on the target plane that satisfy the reference point condition, the target processing system may select any point as the reference point.
Step S103, determining a total dose distribution map according to the dose distribution map of each beam, the dose at the reference point matched by each beam and the prescription dose, wherein the total dose distribution map is used for representing the dose distribution required in the radiotherapy process.
Alternatively, for each beam, after determining the reference point at which the beam matches, the target processing system may determine the dose of the beam at the reference point and then adjust the dose profile of the beam based on the prescribed dose and the dose of the beam at the reference point. For example, the ratio between the prescribed dose and the dose of the beam at the reference point is calculated and then directly multiplied by the individual dose values in the dose profile to obtain an adjusted dose profile. For another example, the contrast value is adjusted upward or downward, and then the adjusted ratio is multiplied by each dose value in the dose profile to obtain the adjusted dose profile. The prescription dose is preset by a user, and the prescription dose is a unique value.
Further, after determining the adjusted dose distribution map corresponding to each beam, the target processing system may directly add the dose values at the same point in each adjusted dose distribution map, and then adjust the dose according to the prescribed dose, so as to obtain a total dose distribution map, for example, add the dose values at the same point in each adjusted dose distribution map to obtain a total dose distribution map to be adjusted, and then make the maximum dose in the total dose distribution map to be adjusted be the prescribed dose, so as to obtain the total dose distribution map. Optionally, the target processing system may also perform weighted summation on the dose values at the same point in each adjusted dose distribution map to obtain a total dose distribution map. Optionally, the target processing system may further adjust the adjusted dose distribution map according to the new reference point after performing weighted summation on the adjusted dose distribution maps, so as to obtain a total dose distribution map.
Based on the above-defined schemes of step S101 to step S103, it may be known that, in the embodiment of the present invention, the reference point is determined according to the dose of the beam on the target plane, the target plane corresponding to the beam is determined for each of the plurality of beams generated by the radiotherapy apparatus by determining the dose distribution required in the radiotherapy process according to the reference point of each beam, and then the reference point matched with the beam is determined from the target plane according to the dose of the beam on the target plane, so as to determine the total dose distribution according to the dose distribution of each beam, the dose at the reference point matched with each beam, and the prescribed dose, wherein the ratio between the dose at the reference point and the maximum dose of the beam on the target plane is greater than or equal to a preset value, and the total dose distribution is used to characterize the dose distribution required in the radiotherapy process.
In the above process, for each beam, by determining the beam matching reference point from the target plane according to the ratio between the dose of the point on the target plane and the maximum dose, on the one hand, the independent determination of the reference point of each beam is realized, and on the other hand, the relatively large dose at the reference point of each beam is ensured, so that the relative error at the reference point of each beam can be controlled within a certain range, and the accuracy of the determined radiotherapy dose can be effectively improved when the total dose distribution map is determined according to the dose distribution map of each beam, the dose at the reference point of each beam matching and the prescription dose.
Therefore, the scheme provided by the application achieves the purposes of determining the reference point according to the dose of the beam on the target plane and determining the required dose distribution in the radiotherapy process according to the reference points of the beams, thereby realizing the technical effect of improving the accuracy of the determined radiotherapy dose, and further solving the technical problem that the accuracy of the determined radiotherapy dose is low due to the fact that the relative error at the fixed point is difficult to control when the radiotherapy dose is determined according to the dose of the beam at the fixed point in the related technology.
In an alternative embodiment, the target plane is a plane passing through the isocenter of the radiotherapy apparatus and perpendicular to the isocenter of the beam.
Wherein in this case the target plane corresponds to the isocenter plane to which the beam corresponds.
It should be noted that, by setting the target plane as the isocenter plane of beam matching, it is ensured that the maximum dose on the target plane can be relatively high at all times, that is, the dose at the determined reference point can be relatively high, and the relative error is small, so that when any plane passing through the beam is selected as the target plane, the situation that the maximum dose on the target plane is relatively low is avoided, thereby affecting the accuracy of the determined radiotherapy dose.
In an alternative embodiment, in determining the beam matching reference point from the target plane based on the dose of the beam on the target plane, the target processing system may determine the dose of the beam at the isocenter, obtain a first dose, and then determine the maximum dose of the beam on the target plane, obtain a second dose, thereby determining the beam matching reference point from the target plane based on the proportional relationship between the first dose and the second dose.
The target processing system can determine the dose of the beam at the isocenter from the target plane according to the coordinate information of the isocenter, and the first dose is obtained.
Alternatively, after determining the first dose and the second dose, the target processing system may calculate a ratio between the first dose and the second dose, thereby determining a beam matching reference point from the target plane based on the magnitude of the ratio. For example, the ratio is compared with the above-described preset value, then in the case where the ratio is greater than or equal to the preset value, the isocenter is taken as the reference point, in the case where the ratio is less than the preset value, the point other than the isocenter is taken as the reference point, and for example, the ratio is compared with the target value, then in the case where the ratio is greater than or equal to the target value, the isocenter is taken as the reference point, and in the case where the ratio is less than the target value, the point other than the isocenter is taken as the reference point. The target value may be set smaller than the preset value, and the difference between the target value and the preset value is smaller than the preset difference, so that the reference points corresponding to the beams can be unified as much as possible under the condition that the dose value of the control reference point is in a relatively high interval, thereby reducing the processing time consumption of the target processing system and improving the efficiency of determining the total dose distribution map.
It should be noted that, by determining the reference point from the target plane according to the proportional relationship between the dose of the beam at the isocenter and the maximum dose, there is a case that the isocenter can be directly used as the reference point, so that the efficiency of determining the reference point can be effectively improved, and the problem of resource waste caused by comparing the dose values of all the points with the maximum dose is avoided.
In an alternative embodiment, in determining the beam matching reference point from the target plane according to the proportional relationship between the first dose and the second dose, the target processing system may determine the isocenter as the reference point if the ratio between the first dose and the second dose is greater than or equal to a preset value, and determine the reference point from points other than the isocenter on the target plane if the ratio between the first dose and the second dose is less than the preset value.
Optionally, the ratio between the first dose and the second dose is the first dose divided by the second dose. In this embodiment, the preset value may be 0.5, and when the preset value is 0.5, if the first dose is greater than or equal to 1/2 of the second dose, the isocenter may be determined as the reference point, otherwise, if the first dose is less than 1/2 of the second dose, the reference point is determined from points other than the isocenter on the target plane.
Wherein in determining the reference point from points other than the isocenter on the target plane, the target processing system may determine the point corresponding to the second dose from the target plane to obtain the target point, thereby determining the target point as the reference point, and thereby controlling the relative error at the reference point as small as possible.
It should be noted that, through the above procedure, it is ensured that the dose of the determined reference point is always relatively large, so that the accuracy of the finally determined dose distribution can be effectively improved.
In an alternative embodiment, in determining the total dose distribution map from the dose distribution map of each beam, the dose at the reference point to which each beam is matched, and the prescribed dose, the target processing system may adjust the dose distribution map of the beam for each beam from the dose at the reference point to which the beam is matched and the prescribed dose, resulting in an adjusted dose distribution map, thereby determining the total dose distribution map from the adjusted dose distribution maps of the plurality of beams.
Optionally, in determining the adjusted dose distribution map, the target processing system may calculate a ratio between the prescribed dose and the dose at the reference point to obtain a first target ratio, thereby determining the adjusted dose distribution map based on a product of the dose distribution map of the beam and the first target ratio. Wherein the first target ratio is obtained by dividing the prescribed dose by the dose at the reference point.
For example, after the first target ratio is obtained, the dose distribution map multiplied by the first target ratio is directly used as an adjusted dose distribution map, that is, the dose value at the reference point is equal to the prescription dose, for example, the dose distribution map multiplied by the first target ratio is used as an adjusted initial dose distribution map, and then the adjusted initial dose distribution map is adjusted according to a preset adjustment strategy to obtain an adjusted dose distribution map. The foregoing adjustment strategy may be to round up the dose value in the adjusted initial dose distribution map, uniformly up-adjust or uniformly down-adjust the dose value in the adjusted initial dose distribution map, and the like.
Further, after obtaining the adjusted dose distribution diagrams of the plurality of beams, the target processing system may directly add the dose values at the same point in each adjusted dose distribution diagram, and then adjust the dose according to the prescribed dose, so as to obtain a total dose distribution diagram, and optionally, the target processing system may also perform weighted summation on the dose values at the same point in each adjusted dose distribution diagram, so as to obtain the total dose distribution diagram. Optionally, the target processing system may further adjust the adjusted dose distribution map according to the new reference point after performing weighted summation on the adjusted dose distribution maps, so as to obtain a total dose distribution map.
It should be noted that, the dose distribution map is adjusted according to the prescribed dose and the dose at the reference point, so that the total dose distribution map is determined according to the adjusted dose distribution map, so that the dose distribution in the total dose distribution map can meet the prescription requirement as much as possible, and the accuracy of the total dose distribution map is improved.
In an alternative embodiment, in determining the total dose distribution map from the adjusted dose distribution maps of the plurality of beams, the target processing system may determine an initial total dose distribution map from the weights of the respective beams and the adjusted dose distribution maps of the plurality of beams, and then obtain a target reference point, thereby adjusting the initial total dose distribution map to obtain the total dose distribution map based on the dose and the prescribed dose at the target reference point, wherein the target reference point is within the initial total dose distribution map.
Alternatively, the weights of the individual beams may be user preset and the target processing system may determine the initial total dose profile according to the following formula:
；
wherein D is all Represents the initial total dose distribution map, n represents the number of beams, i represents the ith beam, w i Representing the weight of the ith beam, D i Representing the adjusted dose distribution map of the ith beam.
Further, the target processing system may adjust the initial total dose profile based on the dose at the target reference point and the prescribed dose to obtain a total dose profile. Alternatively, the target processing system may calculate the ratio between the prescribed dose and the dose at the target reference point, i.e. dividing the prescribed dose by the dose at the target reference point, to obtain a second target ratio, after which the target processing system may determine the total dose profile from the product between the initial total dose profile and the second target ratio.
For example, after the second target ratio is obtained, the initial total dose distribution map multiplied by the second target ratio is directly taken as a total dose distribution map, that is, the dose value at the target reference point is equal to the prescription dose, for example, the initial total dose distribution map multiplied by the second target ratio is taken as a new initial dose distribution map, and then the new initial dose distribution map is adjusted according to a preset adjustment strategy, so as to obtain the total dose distribution map. The foregoing adjustment strategy may be to round up the dose value in the new initial dose distribution map, uniformly up-regulate the dose value in the new initial dose distribution map, uniformly down-regulate the dose value, and the like.
It should be noted that, by combining the new reference point to perform the secondary adjustment in the process of determining the total dose distribution map according to the adjusted dose distribution map, the satisfaction degree of the total dose distribution map on the prescription requirement can be further improved, so that the accuracy of the determined total dose distribution map is improved.
In an alternative embodiment, the target processing system may obtain the target reference point by one of the following: acquiring a preset reference point, and determining the preset reference point as a target reference point; according to the dose of each point in the initial total dose distribution map, sequencing each point in the initial total dose distribution map to obtain a target sequence, and determining the Nth point in the target sequence as a target reference point, wherein N is a positive integer.
Optionally, in the first manner, the preset reference point is determined by the user according to the actual determination, and the target processing system may acquire the preset reference point, so as to use the preset reference point as the target reference point.
Alternatively, in the second manner, the target processing system may determine the target reference point according to preset requirement information. For example, the foregoing demand information may be that the dose value to which 95% of the volume of the target volume is planned is required to be above the prescribed dose, in which case the target processing system may order the points in the initial total dose profile in order of the dose from high to low, resulting in a target sequence, whereby the point at the first 5% of the target sequence is determined to be the target reference point. As another example, the aforementioned demand information may be that the dose value to which 90% of the volume of the target volume is planned is required to be above the prescribed dose, in which case the target processing system may determine the point at the first 10% position in the target sequence as the target reference point, thereby enabling the final determined total dose profile to meet the content of the demand information.
By the aid of the method, the target reference point is flexibly determined, so that various practical application requirements are met conveniently, and applicability of the method and accuracy of a total dose distribution map are improved.
Therefore, the scheme provided by the application achieves the purposes of determining the reference point according to the dose of the beam on the target plane and determining the required dose distribution in the radiotherapy process according to the reference points of the beams, thereby realizing the technical effect of improving the accuracy of the determined radiotherapy dose, and further solving the technical problem that the accuracy of the determined radiotherapy dose is low due to the fact that the relative error at the fixed point is difficult to control when the radiotherapy dose is determined according to the dose of the beam at the fixed point in the related technology.
Example 2
According to an embodiment of the present invention, there is provided an embodiment of a radiotherapy dose determination device, wherein fig. 5 is a schematic diagram of an alternative radiotherapy dose determination device according to an embodiment of the present invention, as shown in fig. 5, the device comprises:
a first determining module 501, configured to determine, for each of a plurality of beams generated by a radiotherapy apparatus, a target plane corresponding to the beam;
A second determining module 502, configured to determine a reference point for beam matching from the target plane according to a dose of the beam on the target plane, where a ratio between the dose at the reference point and a maximum dose of the beam on the target plane is greater than or equal to a preset value;
a third determining module 503 is configured to determine a total dose distribution map according to the dose distribution map of each beam, the dose at the reference point matched by each beam, and the prescribed dose, where the total dose distribution map is used to characterize the dose distribution required in the radiotherapy process.
It should be noted that the first determining module 501, the second determining module 502, and the third determining module 503 correspond to steps S101 to S103 in the above embodiment, and the three modules are the same as examples and application scenarios implemented by the corresponding steps, but are not limited to those disclosed in the above embodiment 1.
Optionally, the target plane is a plane passing through an isocenter of the radiotherapy apparatus and perpendicular to an isocenter of the beam.
Optionally, the second determining module 502 further includes: a first determining sub-module for determining a dose of the beam at the isocenter, resulting in a first dose; a second determining sub-module for determining a maximum dose of the beam on the target plane, resulting in a second dose; a third determination sub-module for determining a reference point for beam matching from the target plane according to a proportional relationship between the first dose and the second dose.
Optionally, the third determining sub-module further comprises: a first determining unit configured to determine an isocenter as a reference point in a case where a ratio between the first dose and the second dose is greater than or equal to a preset value; and a second determining unit for determining a reference point from points other than the isocenter on the target plane in the case where the ratio between the first dose and the second dose is smaller than a preset value.
Optionally, the third determining module 503 further includes: the adjusting submodule is used for adjusting the dose distribution map of each beam according to the dose and the prescription dose at the reference point matched with the beam to obtain an adjusted dose distribution map; a fourth determination sub-module for determining a total dose profile from the adjusted dose profiles of the plurality of beams.
Optionally, the fourth determining submodule further includes: a third determining unit for determining an initial total dose distribution map based on the weights of the respective beams and the adjusted dose distribution maps of the plurality of beams; the acquisition unit is used for acquiring a target reference point, wherein the target reference point is positioned in the initial total dose distribution map; and the adjusting unit is used for adjusting the initial total dose distribution map according to the dose at the target reference point and the prescription dose to obtain a total dose distribution map.
Optionally, the radiotherapy dose determining device further comprises: the acquisition module is used for acquiring a preset reference point and determining the preset reference point as a target reference point; and a fourth determining module, configured to sort the points in the initial total dose distribution map according to the doses of the points in the initial total dose distribution map, obtain a target sequence, and determine an nth point in the target sequence as a target reference point, where N is a positive integer.
Example 3
According to another aspect of embodiments of the present invention, there is also provided a computer readable storage medium having a computer program stored therein, wherein the computer program is arranged to perform the above-described method of determining a radiation dose when run.
Example 4
According to another aspect of an embodiment of the present invention, there is also provided an electronic device, wherein fig. 6 is a schematic diagram of an alternative electronic device according to an embodiment of the present invention, as shown in fig. 6, the electronic device including one or more processors; and a memory for storing one or more programs which, when executed by the one or more processors, cause the one or more processors to implement a method for running the program, wherein the program is configured to perform the method of determining a radiation therapy dose as described above when run.
The foregoing embodiment numbers of the present invention are merely for the purpose of description, and do not represent the advantages or disadvantages of the embodiments.
In the foregoing embodiments of the present invention, the descriptions of the embodiments are emphasized, and for a portion of this disclosure that is not described in detail in this embodiment, reference is made to the related descriptions of other embodiments.
In the several embodiments provided in the present application, it should be understood that the disclosed technology content may be implemented in other manners. The above-described embodiments of the apparatus are merely exemplary, and the division of units may be a logic function division, and there may be another division manner in actual implementation, for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be through some interfaces, units or modules, or may be in electrical or other forms.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed over a plurality of units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.
The integrated units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in essence or a part contributing to the prior art or all or part of the technical solution in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server or a network device, etc.) to perform all or part of the steps of the method of the various embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a removable hard disk, a magnetic disk, or an optical disk, or other various media capable of storing program codes.
The foregoing is merely a preferred embodiment of the present invention and it should be noted that modifications and adaptations to those skilled in the art may be made without departing from the principles of the present invention, which are intended to be comprehended within the scope of the present invention.
Claims (9)
1. A method of determining a radiation dose, comprising:
for each of a plurality of beams generated by a radiotherapy apparatus, determining a target plane corresponding to the beam;
determining a reference point for matching the beam from the target plane according to the dose of the beam on the target plane, wherein the ratio between the dose at the reference point and the maximum dose of the beam on the target plane is greater than or equal to a preset value;
determining a total dose distribution map according to the dose distribution map of each beam, the dose at the reference point matched by each beam and the prescription dose, wherein the total dose distribution map is used for representing the dose distribution required in the radiotherapy process;
wherein determining the total dose distribution map from the dose distribution map for each beam, the dose at the reference point for each beam match, and the prescribed dose comprises:
Dividing the prescribed dose by the dose at the beam-matched reference point for each beam to obtain a first target ratio;
determining an adjusted dose distribution map from the product between the dose distribution map of the beam and the first target ratio;
determining the total dose profile from the adjusted dose profiles of the plurality of beams.
2. The method of claim 1, wherein the target plane is a plane passing through an isocenter of the radiotherapy apparatus and perpendicular to an isocenter of the beam.
3. The method of claim 2, wherein determining the reference point for beam matching from the target plane based on the dose of the beam on the target plane comprises:
determining a dose of the beam at the isocenter, resulting in a first dose;
determining a maximum dose of the beam on the target plane, resulting in a second dose;
a reference point for the beam matching is determined from the target plane according to a proportional relationship between the first dose and the second dose.
4. A method according to claim 3, wherein determining the beam matching reference point from the target plane based on the proportional relationship between the first dose and the second dose comprises:
Determining the isocenter as the reference point in case the ratio between the first dose and the second dose is greater than or equal to the preset value;
in the case where the ratio between the first dose and the second dose is smaller than the preset value, the reference point is determined from points other than the isocenter on the target plane.
5. The method of claim 1, wherein determining the total dose profile from the adjusted dose profiles of the plurality of beams comprises:
determining an initial total dose distribution map based on the weights of the individual beams and the adjusted dose distribution maps of the plurality of beams;
obtaining a target reference point, wherein the target reference point is within the initial total dose distribution map;
and adjusting the initial total dose distribution map according to the dose at the target reference point and the prescription dose to obtain the total dose distribution map.
6. The method of claim 5, wherein the target reference point is obtained by one of:
acquiring a preset reference point, and determining the preset reference point as the target reference point;
And sequencing all points in the initial total dose distribution map according to the dose of all points in the initial total dose distribution map to obtain a target sequence, and determining the Nth point in the target sequence as the target reference point, wherein N is a positive integer.
7. A radiotherapy dosage determination device, comprising:
a first determining module, configured to determine, for each beam of a plurality of beams generated by a radiotherapy apparatus, a target plane corresponding to the beam;
a second determining module, configured to determine a reference point for matching the beam from the target plane according to a dose of the beam on the target plane, where a ratio between a dose at the reference point and a maximum dose of the beam on the target plane is greater than or equal to a preset value;
a third determining module, configured to determine a total dose distribution map according to the dose distribution map of each beam, the dose at the reference point matched by each beam, and the prescribed dose, where the total dose distribution map is used to characterize the dose distribution required in the radiotherapy process;
wherein the radiotherapy dosage determining device is further used for:
dividing the prescribed dose by the dose at the beam-matched reference point for each beam to obtain a first target ratio;
Determining an adjusted dose distribution map from the product between the dose distribution map of the beam and the first target ratio;
determining the total dose profile from the adjusted dose profiles of the plurality of beams.
8. A computer readable storage medium, characterized in that a computer program is stored in the computer readable storage medium, wherein the computer program is arranged to perform the method of determining a radiation therapy dose according to any one of claims 1 to 6 when run.
9. An electronic device, the electronic device comprising one or more processors; a memory for storing one or more programs which, when executed by the one or more processors, cause the one or more processors to implement a method for running a program, wherein the program is arranged to perform the method of determining a radiation therapy dose as claimed in any one of claims 1 to 6 when run.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311503258.9A CN117224864B (en) | 2023-11-13 | 2023-11-13 | Method and device for determining radiotherapy dosage, storage medium and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311503258.9A CN117224864B (en) | 2023-11-13 | 2023-11-13 | Method and device for determining radiotherapy dosage, storage medium and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117224864A CN117224864A (en) | 2023-12-15 |
| CN117224864B true CN117224864B (en) | 2024-02-27 |
Family
ID=89086490
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311503258.9A Active CN117224864B (en) | 2023-11-13 | 2023-11-13 | Method and device for determining radiotherapy dosage, storage medium and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117224864B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101011617A (en) * | 2006-12-29 | 2007-08-08 | 成都川大奇林科技有限责任公司 | Method for determining radiating field output dose accurately in conformalradiotherapy |
| CN104189993A (en) * | 2014-08-14 | 2014-12-10 | 广州医科大学附属肿瘤医院 | Radiotherapy plan designing simulation and three-dimensional dose distribution testing device |
| CN104645500A (en) * | 2015-02-12 | 2015-05-27 | 上海联影医疗科技有限公司 | Method and system for optimizing intensity modulated radiation therapy |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11065471B2 (en) * | 2015-11-20 | 2021-07-20 | Duke University | Systems and methods for automatic, customized radiation treatment plan generation for cancer |
-
2023
- 2023-11-13 CN CN202311503258.9A patent/CN117224864B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101011617A (en) * | 2006-12-29 | 2007-08-08 | 成都川大奇林科技有限责任公司 | Method for determining radiating field output dose accurately in conformalradiotherapy |
| CN104189993A (en) * | 2014-08-14 | 2014-12-10 | 广州医科大学附属肿瘤医院 | Radiotherapy plan designing simulation and three-dimensional dose distribution testing device |
| CN104645500A (en) * | 2015-02-12 | 2015-05-27 | 上海联影医疗科技有限公司 | Method and system for optimizing intensity modulated radiation therapy |
Non-Patent Citations (6)
| Title |
|---|
| 宫颈癌腔内后装治疗剂量学研究进展;张书旭, 林生趣, 徐海荣;国外医学.肿瘤学分册(第11期);全文 * |
| 放射性粒子定向植入治疗计划系统在口腔颌面部肿瘤中的初步应用;宋铁砾, 张建国, 于世平, 周刚, 申文江;现代口腔医学杂志(第02期);全文 * |
| 放射治疗中常规剂量的测算(之二)――临床处方剂量的计算;张绍刚;;医疗装备(第05期);全文 * |
| 放射治疗技术和仪器综述;付礼霞;孙仲轩;冯念伦;;中国医学装备(第05期);全文 * |
| 鼻咽癌动态调强放疗计划的剂量学验证;肖锋;孙朝阳;石梅;张学成;陆宙;陆军;刘晓莉;;现代肿瘤医学(第02期);全文 * |
| 鼻咽癌近距离放疗剂量参考点的设置;李智华, 李瑞英, 赵春艳;中国肿瘤临床(第07期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117224864A (en) | 2023-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11413474B2 (en) | System and method for modelling of dose calculation in radiotherapy treatment planning | |
| US11383102B2 (en) | Three-dimensional radiotherapy dose distribution prediction | |
| Engwall et al. | Effectiveness of different rescanning techniques for scanned proton radiotherapy in lung cancer patients | |
| Bortfeld et al. | Effects of motion on the total dose distribution | |
| Shepard et al. | Direct aperture optimization: a turnkey solution for step‐and‐shoot IMRT | |
| US9956428B2 (en) | Radiosurgical planning and treatment | |
| Dias et al. | Simulated annealing applied to IMRT beam angle optimization: a computational study | |
| EP2994195A2 (en) | An interactive dose gradient based optimization technique to control imrt delivery complexity | |
| Yuan et al. | Lung IMRT planning with automatic determination of beam angle configurations | |
| Fernandez et al. | Interplay effects in highly modulated stereotactic body radiation therapy lung cases treated with volumetric modulated arc therapy | |
| Lee et al. | Dosimetric effects of roll rotational setup errors on lung stereotactic ablative radiotherapy using volumetric modulated arc therapy | |
| CN117224864B (en) | Method and device for determining radiotherapy dosage, storage medium and electronic equipment | |
| Leal et al. | MLC leaf width impact on the clinical dose distribution: a Monte Carlo approach | |
| Huesa-Berral et al. | A method using 4D dose accumulation to quantify the interplay effect in lung stereotactic body radiation therapy | |
| Mutaf et al. | Optimization of internal margin to account for dosimetric effects of respiratory motion | |
| Shang et al. | Impact of multiple beams on plan quality, linear energy transfer distribution, and plan robustness of intensity modulated proton therapy for lung cancer | |
| van den Dobbelsteen et al. | Dosimetric evaluation of off-axis fields and angular transmission for the 1.5 T MR-linac | |
| CN115797262A (en) | Method and device for recommending portal information, processor and electronic equipment | |
| Lee et al. | Intra‐arc binary collimation with dynamic axes trajectory optimization for the SRS treatment of multiple metastases with multiple prescriptions | |
| Arjunan et al. | Dosimetric impact of random spot positioning errors in intensity modulated proton therapy plans of small and large volume tumors | |
| Tang et al. | Stochastic versus deterministic kernel-based superposition approaches for dose calculation of intensity-modulated arcs | |
| Giles et al. | Radiosurgery treatment planning using conformal arc informed volumetric modulated arc therapy | |
| Vrančić et al. | Experimental evaluation of a robust optimization method for IMRT of moving targets | |
| Wang et al. | A dose calculation framework for dynamic electron arc radiotherapy (DEAR) using VirtuaLinac Monte Carlo simulation tool | |
| Rojas-López et al. | Assessment of PTV margin adequacy for single isocenter multiple brain metastases using genetic algorithms |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |