CN115025403A - Dose prediction method and device based on radiotherapy - Google Patents
Dose prediction method and device based on radiotherapy Download PDFInfo
- Publication number
- CN115025403A CN115025403A CN202210496970.XA CN202210496970A CN115025403A CN 115025403 A CN115025403 A CN 115025403A CN 202210496970 A CN202210496970 A CN 202210496970A CN 115025403 A CN115025403 A CN 115025403A
- Authority
- CN
- China
- Prior art keywords
- dose
- volume histogram
- value
- predicted
- distribution
- 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
Images
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/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1071—Monitoring, verifying, controlling systems and methods for verifying the dose delivered by the treatment plan
-
- 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
Abstract
The invention discloses a dose prediction method and a device based on radiotherapy, comprising the following steps: acquiring an original integral dose volume histogram of a structure to be predicted, and acquiring reference dose distribution and actual dose distribution according to a preset measurement plane; acquiring the projection weight distribution of the structure to be predicted on the preset measuring plane according to the preset measuring plane and the control point parameters; calculating and obtaining a reference dose volume histogram and an actual dose volume histogram based on projection according to the projection weight distribution, the original integral dose volume histogram, the reference dose distribution and the actual dose distribution; and calculating a dose difference value of the reference dose volume histogram and the actual dose volume histogram based on the projection, and calculating and obtaining a prediction integral dose volume histogram according to the dose difference value, thereby completing the prediction of the dose. The invention solves the technical problem that complex modeling and model testing processes are required in the prior art.
Description
Technical Field
The invention relates to the technical field of medical software, in particular to a dose prediction method and device based on radiotherapy.
Background
Due to the complexity of radiation treatment, pre-treatment dose verification is required before the treatment plan is applied to the actual treatment of the patient in order to ensure medical safety and to ensure that the plan used for treatment meets clinical requirements, to ensure consistency of the treatment plan being performed with the plan design. The dose verification is to compare the dose distribution measured during the actual execution of the plan with the dose distribution calculated by the plan, and evaluate the effect of the actual execution of the plan by evaluating the difference between the actual measurement value and the calculated value, and generally using an ionization chamber matrix, a semiconductor matrix, and an Electronic Portal Imaging Device (EPID) as a dose measuring tool for the pre-treatment intensity-modulated treatment plan.
Common Dose distribution comparison and analysis methods include a Dose Difference method (Dose Difference Test), a Distance-To-agent method (Distance-To-agent), a Gamma analysis method and a three-dimensional Dose reconstruction comparison method, and the first three methods aim at the comparison of direct measurement results and require that a treatment plan is transplanted To a measurement device for calculation so as To obtain the Dose distribution of a measurement plane. The general measuring equipment is in a uniform cubic, cylindrical or polygonal shape, but not in a shape of a human body or a shape of an actually treated human body, because the measuring equipment is different from the actual human body, and the obtained difference result is inconsistent with the actual human body dose difference result, the direct comparison has the biggest defect that clinically required data, such as information of integral dose volume histograms, average doses and the like of different anatomical structures, cannot be intuitively reflected; the three-dimensional dose reconstruction comparison method can obtain the result closest to the treatment, but on one hand, the accurate reconstruction can be completed only by high requirements on measuring equipment, for example, high resolution is required, and the complete capture of accurate flux distribution is ensured; a faster data acquisition speed is required to ensure that the acquisition process can acquire complete plan execution information in real time. On the other hand, complex modeling and model testing processes are needed, and the quality of a modeling result directly influences a final result.
Therefore, there is a need for a dose prediction method and apparatus that avoids the modeling and model testing processes.
Disclosure of Invention
The invention provides a dose prediction method and a dose prediction device based on radiotherapy, and aims to solve the technical problem that complex modeling and model testing processes are required in the prior art.
In order to solve the above technical problem, an embodiment of the present invention provides a dose prediction method based on radiation therapy, including:
acquiring an original integral dose volume histogram of a structure to be predicted, and acquiring reference dose distribution and actual dose distribution according to a preset measurement plane;
acquiring the projection weight distribution of the structure to be predicted on the preset measuring plane according to the preset measuring plane and the control point parameters; the control point parameters are acquired through a medical linear accelerator;
calculating and obtaining a reference dose volume histogram and an actual dose volume histogram based on projection according to the projection weight distribution, the original integral dose volume histogram, the reference dose distribution and the actual dose distribution;
and calculating a dose difference value of the reference dose volume histogram and the actual dose volume histogram based on the projection, and calculating and obtaining a prediction integral dose volume histogram according to the dose difference value, thereby completing the prediction of the dose.
Compared with the prior art, the method can acquire the original integral dose volume histogram of the structure to be predicted by the existing dose distribution comparative analysis method, acquire the reference dose distribution and the actual dose distribution according to the measurement plane, acquire the projection weight distribution of the structure to be predicted on the preset measurement plane according to the control point parameters, further ensure that the reference dose volume histogram and the actual dose volume histogram which are acquired by calculation can be based on the projection weight, improve the accuracy of the dose volume histogram, and accurately calculate and obtain the predicted integral dose volume histogram by the dose difference value between the reference dose volume histogram and the actual dose volume histogram, thereby avoiding the existing process of needing complex modeling and model testing and improving the efficiency of dose prediction.
Preferably, the control point parameters include: the total number of the control points of the medical linear accelerator, the boom angle of each control point and the beam-out hop ratio of each control point.
It can be understood that the control point parameters are determined by the data such as the boom angle of the execution control point of the medical linear accelerator, so that the projection weight distribution can be accurately calculated, and the complex data required in modeling is not required.
As a preferred scheme, the obtaining of the projection weight distribution of the structure to be predicted on the preset measurement plane according to the preset measurement plane and the control point parameters specifically includes:
dividing the preset measuring plane into a plurality of measuring points, respectively obtaining the track length of the straight line corresponding to each measuring point passing through the structure to be predicted according to the straight line determined by each measuring point and the beam outlet source of the linear accelerator, and respectively calculating and obtaining the projection weight distribution of each measuring point according to each measuring point, the track length corresponding to each measuring point, the total number of control points and the beam outlet hop ratio of each control point.
It can be understood that the track length corresponding to each measuring point is obtained by dividing the preset measuring plane according to the obtained measuring points and the straight line determined by the beam-emitting source of the linear accelerator, and the structural projection weight distribution of the structure to be predicted on the measuring plane can be accurately obtained according to the execution control point parameters of the linear accelerator, so that the process of prediction calculation is simplified.
Preferably, after the obtaining the original integral dose volume histogram of the structure to be predicted, the method further includes:
and respectively carrying out normalization operation on the volume value and the dose value of the original integral dose volume histogram according to the preset volume value and the preset dose value of the structure to be predicted, thereby obtaining the normalized original integral dose volume histogram.
It can be understood that the normalization operation is respectively performed on the volume value and the dose value of the integral dose volume histogram through the preset volume value and the preset dose value of the structure to be predicted, so that the normalized original integral dose volume histogram can avoid errors in the calculation process caused by data of different magnitudes, and the calculation accuracy is improved.
As a preferred scheme, the calculating and obtaining a projection-based reference dose volume histogram and an actual dose volume histogram according to the projection weight distribution, the original integral dose volume histogram, the reference dose distribution, and the actual dose distribution specifically includes:
acquiring and adding projection weights of the measuring points which are more than 0 in the projection weight distribution to obtain a total weight of the measuring points;
respectively carrying out normalization operation on the reference dose distribution and the actual dose distribution according to a preset dose value;
according to the dose value interval in the original integral dose volume histogram, adding the projection weights of the measurement points which are larger than the dose value interval in the normalized reference dose distribution and the normalized actual dose distribution respectively to obtain a first weight sum and a second weight sum of each dose value interval respectively; wherein the first weighted sum is a weighted sum of a reference dose distribution and the second weighted sum is a weighted sum of an actual dose distribution;
and respectively carrying out normalization processing on the first weight sum and the second weight sum according to the total weight of the measuring points, and respectively obtaining a reference dose volume histogram and an actual dose reference histogram based on projection according to the normalized first weight sum and the normalized second weight sum.
It can be understood that, the total weight of the measurement points whose projection weights are greater than zero are obtained and added through statistics of the projection weight distribution, the projection weights of the measurement points corresponding to the interval greater than the dose value are added through the normalized reference dose distribution and the normalized actual dose distribution, so as to obtain a first weight sum and a second weight sum, and the first weight sum and the second weight sum are normalized through the total weight of the measurement points, so that the reference dose volume histogram and the actual dose reference histogram based on projection are accurately obtained.
As a preferred scheme, the calculating of the dose difference value is performed on the reference dose volume histogram and the actual dose volume histogram based on projection, and the predicted integral dose volume histogram is calculated and obtained according to the dose difference value, specifically:
calculating the difference of corresponding dose values when the first weight sum and the second weight sum are the same as each other according to the first weight sum and the second weight sum, and taking the difference as a dose difference value;
and calculating and obtaining a predicted dose value according to the dose difference value and the original integral dose volume histogram, and calculating and obtaining a predicted integral dose volume histogram by taking the weight sum with the same first weight sum and second weight sum as a predicted volume value.
It can be understood that, by calculating the difference between the corresponding dose values when the first weighted sum and the second weighted sum are the same as the dose difference value, the dose difference value between the reference dose volume histogram and the actual dose reference histogram can be accurately calculated, the predicted dose value is calculated by the dose difference value, and the weighted sum when the first weighted sum and the second weighted sum are the same as the predicted volume value, the predicted integral dose volume histogram is accurately calculated and obtained, so that the calculation efficiency of the overall predicted integral dose volume histogram is improved.
Preferably, after the calculating and acquiring the predicted integrated dose volume histogram, the method further includes:
and respectively carrying out reduction operation on the predicted volume value and the predicted dosage value according to the preset volume value and the preset dosage value of the structure to be predicted, thereby obtaining a predicted integral dosage volume histogram reduced to an absolute value.
It can be understood that the predicted volume value and the predicted dose value are respectively reduced through the preset volume value and the preset dose value of the structure to be predicted, so that the predicted integral dose volume histogram reduced to an absolute value is accurately obtained, the intermediate calculation complexity is simplified, the modeling process is avoided, and the calculation efficiency and the accuracy are improved.
Accordingly, the present invention also provides a radiation therapy-based dose prediction device comprising: the device comprises an acquisition module, a projection weight module, a histogram module and a prediction module;
the acquisition module is used for acquiring an original integral dose volume histogram of a structure to be predicted and acquiring reference dose distribution and actual dose distribution according to a preset measurement plane;
the projection weight module is used for acquiring the projection weight distribution of the structure to be predicted on the preset measuring plane according to the preset measuring plane and the control point parameters; the control point parameters are obtained through a medical linear accelerator;
the histogram module is used for calculating and acquiring a reference dose volume histogram and an actual dose volume histogram according to the projection weight distribution, the original integral dose volume histogram, the reference dose distribution and the actual dose distribution;
the prediction module is used for calculating a dose difference value of the reference dose volume histogram and the actual dose volume histogram based on projection, and calculating and obtaining a prediction integral dose volume histogram according to the dose difference value, so that the prediction of the dose is completed.
Preferably, the method further comprises the following steps: a normalization module;
and the normalization module is used for respectively carrying out normalization operation on the volume value and the dose value of the original integral dose volume histogram according to the preset volume value and the preset dose value of the structure to be predicted, so as to obtain the normalized original integral dose volume histogram.
Preferably, the method further comprises the following steps: a reduction module;
and the reduction module is used for respectively carrying out reduction operation on the predicted volume value and the predicted dosage value according to the preset volume value and the preset dosage value of the structure to be predicted, so as to obtain a predicted integral dosage volume histogram reduced to an absolute value.
Drawings
FIG. 1: a flowchart of the steps of a method for radiation therapy based dose prediction is provided in accordance with an embodiment of the present invention;
FIG. 2 is a schematic diagram: the embodiment of the invention provides a dose prediction method based on radiotherapy, which is a schematic diagram of an outgoing beam source of a linear accelerator for Chinese medicine;
FIG. 3: the embodiment of the invention provides a flow chart of calculation steps for obtaining a projection-based dose volume histogram in a radiation therapy-based dose prediction method;
FIG. 4: the embodiment of the invention provides a flow chart of the calculation steps of the prediction integral dose volume histogram in the dose prediction method based on radiation treatment;
FIG. 5: the embodiment of the invention provides a structural schematic diagram of a dose prediction device based on radiation therapy.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example one
Referring to fig. 1, a dose prediction method based on radiation therapy provided by an embodiment of the present invention includes the following steps S101-S104:
s101: and acquiring an original integral dose volume histogram of the structure to be predicted, and acquiring reference dose distribution and actual dose distribution according to a preset measuring plane.
Preferably, after the obtaining the original integral dose volume histogram of the structure to be predicted, the method further includes:
and respectively carrying out normalization operation on the volume value and the dose value of the original integral dose volume histogram according to the preset volume value and the preset dose value of the structure to be predicted, thereby obtaining the normalized original integral dose volume histogram.
In this embodiment, a normalized integrated dose volume histogram DV of the anatomical structure S to be predicted to be analyzed in the original Treatment plan is derived, preferably by a radiation Treatment Planning System (TPS) or third party independent dose verification softwareH plan The anatomical structures to be predicted include, but are not limited to, the target area for irradiation, the brain stem, the spinal cord, and the like. In this embodiment, the normalized original integrated dose volume histogram is calculated by comparing the volume value V of the original integrated dose volume histogram plan Divided by the volume value V of the anatomical structure S to be predicted plan，max And the dose value D of the original integrated dose volume histogram plan Divided by the maximum dose value D in the original treatment plan plan，max To obtain a normalized original integrated dose volume histogram.
It can be understood that the normalization operation is respectively carried out on the volume value and the dose value of the integral dose volume histogram through the preset volume value and the preset dose value of the structure to be predicted, so that the normalized original integral dose volume histogram can avoid errors in the calculation process caused by data of different magnitudes, and the calculation accuracy is improved.
Further, in this embodiment, the reference dose distribution and the actual dose distribution are calculated by transplanting the target treatment plan to a phantom of an electronic portal imaging system (EPID) through a radiation treatment planning system or a third-party independent dose verification software, and calculating a reference dose distribution D of an EPID measurement plane under an ideal condition ref And executing a target treatment plan, obtaining an actual dose distribution D of the medical linear accelerator executing the plan through EPID measurement eva 。
S102: acquiring the projection weight distribution of the structure to be predicted on the preset measuring plane according to the preset measuring plane and the control point parameters; and acquiring the control point parameters through a medical linear accelerator.
It should be noted that, the medical linear accelerator uses the microwave electromagnetic field to accelerate the electrons and has the accelerating device of the linear motion orbit, a medical apparatus used for the radiotherapy of the tumor or other focuses of the patient, it can produce high-energy X-ray and electron beam, it has high dose rate, the time of irradiation is short, the radiation field is large, the uniformity and stability of the dose are good, and the half shadow area is small, etc. The control point parameters are obtained from a plan file (generally, a standardized Dicom file) of a radiation treatment plan in the medical linear accelerator, and may also be obtained from a log file during the execution of the medical linear accelerator, and the content of the control point parameters is specific accelerator device parameters describing the execution of the treatment plan by the medical linear accelerator.
Preferably, the control point parameters include: the total number of control points of the medical linear accelerator, the boom (Gantry) angle of each control point and the beam-out hop ratio of each control point.
It can be understood that the control point parameters are determined by the data such as the boom angle of the execution control point of the medical linear accelerator, so that the projection weight distribution can be accurately calculated, and the complex data required in modeling is not required. Meanwhile, the control point parameters are the total number of control points of the medical linear accelerator for executing a treatment plan, the boom (Gantry) angle of each control point (the boom angle is used for determining the position of an accelerator source, and the medical linear accelerator rotates around an isocenter when executing the treatment plan), and the beam-out hop ratio of each control point.
Specifically, the preset measurement plane is divided into a plurality of measurement points, the length of the track of the straight line corresponding to each measurement point, which passes through the structure to be predicted, is respectively obtained according to the straight line determined by each measurement point and the beam outlet source of the linear accelerator, and the projection weight distribution of each measurement point is respectively calculated and obtained according to each measurement point, the track length corresponding to each measurement point, the total number of control points and the beam outlet hop ratio of each control point. The specific calculation formula is as follows:
It can be understood that the track length corresponding to each measuring point is obtained by dividing the preset measuring plane according to the obtained measuring points and the straight line determined by the beam-emitting source of the linear accelerator, and the structural projection weight distribution of the structure to be predicted on the measuring plane can be accurately obtained according to the execution control point parameters of the linear accelerator, so that the process of prediction calculation is simplified.
S103: and calculating and obtaining a reference dose volume histogram and an actual dose volume histogram based on projection according to the projection weight distribution, the original integral dose volume histogram, the reference dose distribution and the actual dose distribution.
As a preferred solution, referring to fig. 3, step S103 specifically includes the following steps S201 to S204:
s201: and acquiring and adding the projection weights of the measuring points which are more than 0 in the projection weight distribution to obtain the total weight of the measuring points. The specific calculation formula is as follows:
N＝∑ ij W ij wherein W is ij >And 0 and N are the total weight of the measuring points.
S202: and respectively carrying out normalization operation on the reference dose distribution and the actual dose distribution according to a preset dose value.
In this embodiment, a reference dose distribution D is preferably chosen ref The maximum dose value of (A) as a dose normalization value D max . Taking the reference dose distribution as an example, the reference dose distribution D ref Dose value D in ij Is normalized to D ij ’＝D ij /D max To obtain a normalized reference dose distribution D ref '. Similarly, for the actual dose distribution D eva Normalized operation of (2), actual dose distribution D eva Normalizing the dose value in (A) to D ij ’＝D ij /D max To obtain a normalized actual dose distribution D eva ’。
S203: according to the dose value interval in the original integral dose volume histogram, adding the projection weights of the measurement points which are larger than the dose value interval in the normalized reference dose distribution and the normalized actual dose distribution respectively to obtain a first weight sum and a second weight sum of each dose value interval respectively; wherein the first weighted sum is a weighted sum of a reference dose distribution, and the second weighted sum is a weighted sum of an actual dose distribution.
In this embodiment, the original integrated dose volume histogram DVH is preferably used plan Dose interval D of (1) plan For statistical dose value intervals, the normalized reference dose distribution D is counted ref ' the middle dose is larger than D of the k dose interval plan，k Is measured at each measurement point, and a first weight sum N corresponding to each measurement point k1 And statistically normalized actual dose distribution D eva ' the middle dose is larger than D of the k dose interval plan，k Each measuring point of (2) corresponds to a second weighted sum N k2 。
S204: and respectively carrying out normalization processing on the first weight sum and the second weight sum according to the total weight of the measuring points, and respectively obtaining a reference dose volume histogram and an actual dose reference histogram based on projection according to the normalized first weight sum and the normalized second weight sum.
In this embodiment, the total weight of the measurement points is calculated by N k ’＝N k and/N respectively carrying out normalization processing on the first weight sum and the second weight sum so as to obtain a projection-based reference dose volume histogram and an actual dose reference histogram.
It can be understood that, the total weight of the measurement points whose projection weights are greater than zero are obtained and added through the statistics of the projection weight distribution, the projection weights of the measurement points corresponding to the interval greater than the dose value are added through the normalized reference dose distribution and the normalized actual dose distribution, so as to obtain a first weight sum and a second weight sum, and the first weight sum and the second weight sum are normalized through the total weight of the measurement points, so as to accurately obtain the reference dose volume histogram DVH based on projection ref And actual dose reference histogram DVH eva 。
S104: and calculating a dose difference value of the reference dose volume histogram and the actual dose volume histogram based on projection, and calculating and obtaining a prediction integral dose volume histogram according to the dose difference value so as to complete the prediction of the dose.
As a preferred scheme, referring to fig. 4, step S104 specifically includes the following steps S301 to S302:
s301: and calculating the difference of corresponding dose values when the first weight sum and the second weight sum are the same according to the first weight sum and the second weight sum, and using the difference as a dose difference value.
In this embodiment, the DVH is calculated by fitting a projection-based reference dose volume histogram ref And the actual dose volume histogram DVH eva Comparing to calculate the same N k ' time DVH ref And DVH eva The difference between the dose values as the dose difference value PD k 。
S302: and calculating and obtaining a predicted dose value according to the dose difference value and the original integral dose volume histogram, and calculating and obtaining a predicted integral dose volume histogram by taking the weight sum with the same first weight sum and second weight sum as a predicted volume value.
In this embodiment, the dose interval D in the original integrated dose volume histogram is used as a basis plan，k And dose differential value PD k Through D pre，k ＝D plan，k +PD k Calculating an acquired predicted dose value D pre，k And in DVH ref And DVH eva Respectively correspond to the same N k ' As a predictive volume value V pre，k According to the predicted dose value D pre，k And a predicted volume value V pre，k To obtain a predicted integrated dose volume histogram.
It can be understood that by calculating the difference between the corresponding dose values when the first weighted sum and the second weighted sum are the same as the dose difference value, the dose difference value between the reference dose volume histogram and the actual dose reference histogram can be accurately calculated, the predicted dose value is calculated by the dose difference value, and the predicted integral dose volume histogram is accurately calculated and obtained by using the weighted sum when the first weighted sum and the second weighted sum are the same as the predicted volume value, so that the calculation efficiency of the overall predicted integral dose volume histogram is improved.
Preferably, after the calculating and obtaining the predicted integrated dose volume histogram, the method further includes:
and respectively carrying out reduction operation on the predicted volume value and the predicted dosage value according to the preset volume value and the preset dosage value of the structure to be predicted, thereby obtaining a predicted integral dosage volume histogram reduced to an absolute value.
In the present embodiment, the volume value V of the anatomical structure S to be predicted in step S101 is calculated plan，max And the maximum dose value D in the original treatment plan plan，max To predict the volume value V pre And predicted dose value D pre And carrying out reduction operation. In particular by predicting the volume value V pre Multiplied by the volume value V of the anatomical structure S to be predicted plan，max And by predicting the dose value D pre Multiplied by the maximum dose value D in the original treatment plan plan，max Thereby obtaining a predicted integrated dose volume histogram that is reduced to absolute values.
It can be understood that the predicted volume value and the predicted dose value are respectively reduced through the preset volume value and the preset dose value of the structure to be predicted, so that a predicted integral dose volume histogram reduced to an absolute value is accurately obtained, the intermediate calculation complexity is simplified, the modeling process is avoided, and the calculation efficiency and accuracy are improved.
As a preferable scheme of this embodiment, the above steps are repeated, and the calculation and acquisition of the predicted integrated dose volume histogram may be performed on other anatomical structures to be analyzed, so as to obtain dose prediction of other anatomical structures.
The embodiment has the following effects:
compared with the prior art, the embodiment of the invention can obtain the original integral dose volume histogram of the structure to be predicted under the original treatment plan through the existing dose distribution comparative analysis method, obtain the reference dose distribution and the actual dose distribution according to the measurement plane, and obtain the projection weight distribution of the structure to be predicted on the preset measurement plane according to the control point parameters, thereby ensuring that the reference dose volume histogram and the actual dose volume histogram which are obtained through calculation can be based on the projection weight, improving the accuracy of the dose volume histogram, and accurately calculating to obtain the predicted integral dose volume histogram through the dose difference value between the reference dose volume histogram and the actual dose volume histogram, avoiding the existing process of needing complex modeling and model testing, and improving the efficiency of dose prediction.
Example two
Accordingly, referring to fig. 5, the present invention also provides a dose prediction device based on radiation therapy, comprising: an acquisition module 401, a projection weight module 402, a histogram module 403, and a prediction module 404.
The obtaining module 401 is configured to obtain an original integral dose volume histogram of a structure to be predicted, and obtain a reference dose distribution and an actual dose distribution according to a preset measurement plane.
As a preferable scheme, the embodiment of the present invention further includes: a normalization module 4011; the normalization module 4011 is configured to, according to the preset volume value and the preset dose value of the structure to be predicted, perform normalization operation on the volume value and the dose value of the original integral dose volume histogram, respectively, so as to obtain a normalized original integral dose volume histogram.
The projection weight module 402 is configured to obtain a projection weight distribution of the structure to be predicted on the preset measurement plane according to the preset measurement plane and the control point parameter; and acquiring the control point parameters through a medical linear accelerator.
Preferably, the control point parameters include: the total number of control points of the medical linear accelerator, the arm support angle of each control point and the beam-out hop ratio of each control point.
As a preferred scheme, the projection weight module 402 is specifically configured to: dividing the preset measuring plane into a plurality of measuring points, respectively obtaining the track length of the straight line corresponding to each measuring point penetrating through the structure to be predicted according to the straight line determined by each measuring point and the beam outlet source of the linear accelerator, and respectively calculating and obtaining the projection weight distribution of each measuring point according to each measuring point, the track length corresponding to each measuring point, the total number of control points and the beam outlet hop ratio of each control point.
The histogram module 403 is configured to calculate and obtain a reference dose volume histogram and an actual dose volume histogram according to the projection weight distribution, the original integral dose volume histogram, the reference dose distribution, and the actual dose distribution.
Preferably, the histogram module 403 is specifically configured to: acquiring and adding projection weights of the measuring points which are more than 0 in the projection weight distribution to obtain a total weight of the measuring points; respectively carrying out normalization operation on the reference dose distribution and the actual dose distribution according to a preset dose value; according to the dose value interval in the original integral dose volume histogram, adding the projection weights of the measurement points which are larger than the dose value interval in the normalized reference dose distribution and the normalized actual dose distribution respectively to obtain a first weight sum and a second weight sum of each dose value interval respectively; wherein the first weighted sum is a weighted sum of a reference dose distribution and the second weighted sum is a weighted sum of an actual dose distribution; and respectively carrying out normalization processing on the first weight sum and the second weight sum according to the total weight of the measuring points, and respectively obtaining a reference dose volume histogram and an actual dose reference histogram based on projection according to the normalized first weight sum and the normalized second weight sum.
The prediction module 404 is configured to perform dose difference value calculation on the reference dose volume histogram and the actual dose volume histogram based on the projection, and calculate and obtain a predicted integral dose volume histogram according to the dose difference value, thereby completing prediction of dose.
Preferably, the prediction module 404 is specifically configured to: calculating the difference of corresponding dose values when the first weight sum and the second weight sum are the same as each other according to the first weight sum and the second weight sum, and taking the difference as a dose difference value; and calculating and obtaining a predicted dose value according to the dose difference value and the original integral dose volume histogram, and calculating and obtaining a predicted integral dose volume histogram by taking the weight sum with the same first weight sum and second weight sum as a predicted volume value.
Preferably, the method further comprises the following steps: a reduction module 4041; the reducing module 4041 is configured to perform reducing operations on the predicted volume value and the predicted dose value according to the preset volume value and the preset dose value of the structure to be predicted, so as to obtain a predicted integral dose volume histogram reduced to an absolute value.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working process of the apparatus described above may refer to the corresponding process in the foregoing method embodiment, and is not described herein again.
The above embodiment is implemented, and has the following effects:
compared with the prior art, the embodiment of the invention can obtain the original integral dose volume histogram of the structure to be predicted by the existing dose distribution comparative analysis method, obtain the reference dose distribution and the actual dose distribution according to the measuring plane, obtain the projection weight distribution of the structure to be predicted on the preset measuring plane according to the control point parameters, further ensure that the reference dose volume histogram and the actual dose volume histogram which are obtained by calculation can be based on the projection weight, improve the accuracy of the dose volume histogram, and accurately calculate and obtain the predicted integral dose volume histogram by the dose difference value between the reference dose volume histogram and the actual dose volume histogram, thereby avoiding the existing process of needing complex modeling and model testing and improving the efficiency of dose prediction.
The above-mentioned embodiments are provided to further explain the objects, technical solutions and advantages of the present invention in detail, and it should be understood that the above-mentioned embodiments are only examples of the present invention and are not intended to limit the scope of the present invention. It should be understood that any modifications, equivalents, improvements and the like, which come within the spirit and principle of the invention, may occur to those skilled in the art and are intended to be included within the scope of the invention.
Claims (10)
1. A radiation therapy based dose prediction method, comprising:
acquiring an original integral dose volume histogram of a structure to be predicted, and acquiring reference dose distribution and actual dose distribution according to a preset measurement plane;
acquiring the projection weight distribution of the structure to be predicted on the preset measuring plane according to the preset measuring plane and the control point parameters; the control point parameters are acquired through a medical linear accelerator;
calculating and obtaining a reference dose volume histogram and an actual dose volume histogram based on projection according to the projection weight distribution, the original integral dose volume histogram, the reference dose distribution and the actual dose distribution;
and calculating a dose difference value of the reference dose volume histogram and the actual dose volume histogram based on the projection, and calculating and obtaining a prediction integral dose volume histogram according to the dose difference value, thereby completing the prediction of the dose.
2. A radiation therapy based dose prediction method as claimed in claim 1, wherein said control point parameters comprise: the total number of control points of the medical linear accelerator, the boom angle of each control point and the beam-out hop ratio of each control point.
3. The radiation therapy-based dose prediction method according to claim 2, wherein the obtaining of the projection weight distribution of the structure to be predicted on the preset measurement plane according to the preset measurement plane and the control point parameters comprises:
dividing the preset measuring plane into a plurality of measuring points, respectively obtaining the track length of the straight line corresponding to each measuring point passing through the structure to be predicted according to the straight line determined by each measuring point and the beam outlet source of the linear accelerator, and respectively calculating and obtaining the projection weight distribution of each measuring point according to each measuring point, the track length corresponding to each measuring point, the total number of control points and the beam outlet hop ratio of each control point.
4. A radiation therapy based dose prediction method as claimed in claim 3, further comprising, after said obtaining a raw integrated dose volume histogram of the structure to be predicted:
and respectively carrying out normalization operation on the volume value and the dose value of the original integral dose volume histogram according to the preset volume value and the preset dose value of the structure to be predicted, thereby obtaining the normalized original integral dose volume histogram.
5. The radiation therapy based dose prediction method as claimed in claim 4, wherein said calculating to obtain a projection based reference dose volume histogram and an actual dose volume histogram based on said projection weight distribution, said original integrated dose volume histogram, said reference dose distribution and said actual dose distribution comprises:
acquiring and adding projection weights of the measuring points which are more than 0 in the projection weight distribution to obtain a total weight of the measuring points;
respectively carrying out normalization operation on the reference dose distribution and the actual dose distribution according to a preset dose value;
according to the dose value interval in the original integral dose volume histogram, adding the projection weights of the measurement points which are larger than the dose value interval in the normalized reference dose distribution and the normalized actual dose distribution respectively to obtain a first weight sum and a second weight sum of each dose value interval respectively; wherein the first weighted sum is a weighted sum of a reference dose distribution and the second weighted sum is a weighted sum of an actual dose distribution;
and respectively carrying out normalization processing on the first weight sum and the second weight sum according to the total weight of the measuring points, and respectively obtaining a reference dose volume histogram and an actual dose reference histogram based on projection according to the normalized first weight sum and the normalized second weight sum.
6. The radiation therapy based dose prediction method as claimed in claim 5, wherein said calculating of the dose difference values for said projection-based reference dose volume histogram and actual dose volume histogram and calculating the predicted integral dose volume histogram according to said dose difference values comprises:
calculating the difference of corresponding dose values when the first weight sum and the second weight sum are the same as each other according to the first weight sum and the second weight sum, and taking the difference as a dose difference value;
and calculating and obtaining a predicted dose value according to the dose difference value and the original integral dose volume histogram, and calculating and obtaining a predicted integral dose volume histogram by taking the weight sum with the same first weight sum and second weight sum as a predicted volume value.
7. A radiation therapy based dose prediction method as claimed in claim 6, wherein after said calculating obtains a predicted integrated dose volume histogram, further comprising:
and respectively carrying out reduction operation on the predicted volume value and the predicted dosage value according to the preset volume value and the preset dosage value of the structure to be predicted, thereby obtaining a predicted integral dosage volume histogram reduced to an absolute value.
8. A radiation therapy based dose prediction device, comprising: the device comprises an acquisition module, a projection weight module, a histogram module and a prediction module;
the acquisition module is used for acquiring an original integral dose volume histogram of a structure to be predicted and acquiring reference dose distribution and actual dose distribution according to a preset measurement plane;
the projection weight module is used for acquiring the projection weight distribution of the structure to be predicted on the preset measuring plane according to the preset measuring plane and the control point parameters; the control point parameters are acquired through a medical linear accelerator;
the histogram module is used for calculating and acquiring a reference dose volume histogram and an actual dose volume histogram according to the projection weight distribution, the original integral dose volume histogram, the reference dose distribution and the actual dose distribution;
the prediction module is used for calculating a dose difference value of the reference dose volume histogram and the actual dose volume histogram based on projection, and calculating and obtaining a prediction integral dose volume histogram according to the dose difference value, so that the prediction of the dose is completed.
9. A radiation therapy based dose prediction device as defined in claim 8, further comprising: a normalization module;
and the normalization module is used for respectively carrying out normalization operation on the volume value and the dose value of the original integral dose volume histogram according to the preset volume value and the preset dose value of the structure to be predicted, so as to obtain the normalized original integral dose volume histogram.
10. A radiation therapy based dose prediction device as defined in claim 8, further comprising: a reduction module;
and the reduction module is used for obtaining a prediction integral dose volume histogram reduced to an absolute value according to the preset volume value and the preset dose value of the structure to be predicted.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210496970.XA CN115025403A (en) | 2022-05-07 | 2022-05-07 | Dose prediction method and device based on radiotherapy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210496970.XA CN115025403A (en) | 2022-05-07 | 2022-05-07 | Dose prediction method and device based on radiotherapy |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115025403A true CN115025403A (en) | 2022-09-09 |
Family
ID=83119313
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210496970.XA Pending CN115025403A (en) | 2022-05-07 | 2022-05-07 | Dose prediction method and device based on radiotherapy |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115025403A (en) |
-
2022
- 2022-05-07 CN CN202210496970.XA patent/CN115025403A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11266858B2 (en) | Systems, devices, and methods for quality assurance of radiation therapy | |
| US10369383B2 (en) | Time-resolved pre-treatment portal dosimetry systems, devices, and methods | |
| US11020615B2 (en) | Computing radiotherapy dose distribution | |
| CN109464756B (en) | Method and device for verifying radiation therapy dosage and radiation therapy equipment | |
| JP5567058B2 (en) | System for performing dosimetry in the body (VIVO) | |
| Nijsten et al. | Routine individualised patient dosimetry using electronic portal imaging devices | |
| CN104548372A (en) | Radiotherapy planning method and device, radiotherapy dose determining method and device and radiotherapy quality guaranteeing method and device | |
| CN109621228A (en) | The calculating unit and storage medium of radiological dose | |
| US20100213394A1 (en) | Optimization of control parameters for a particle irradiation system | |
| CN110237445B (en) | EPID-based in-vivo three-dimensional dose monitoring and verifying method | |
| CN111001097B (en) | Radiotherapy dose evaluation system, device and storage medium | |
| US9364686B2 (en) | Planning a treatment beam aimed at one or more target regions | |
| US10881877B2 (en) | Systems and methods for evaluating motion tracking for radiation therapy | |
| CN112133440A (en) | Dose calculation system, dose quality assurance system and storage medium | |
| CN115025403A (en) | Dose prediction method and device based on radiotherapy | |
| CN114420252A (en) | Method, device and medium for determining intensity modulated radiotherapy plan evaluation parameter matrix | |
| US20220249867A1 (en) | Independent stereotactic radiotherapy dose calculation and treatment plan verification | |
| Giordanengo et al. | RIDOS: A new system for online computation of the delivered dose distributions in scanning ion beam therapy | |
| Li et al. | Monte Carlo as a QA tool for advanced radiation therapy | |
| CN114146329B (en) | Radiation treatment plan optimizing system introducing gamma passing rate optimizing target | |
| CN113426030B (en) | Proton dosage calculation method and device | |
| CN115154926A (en) | Data processing method, device and equipment | |
| CN117727423A (en) | Particle beam range real-time verification method and system | |
| CN116832344A (en) | Coefficient acquisition method for dose calculation, radiation therapy system and related device | |
| Karlsson et al. | Symposia/Proffered papers Monday, 17 September 2001 $25 |
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 |