Radiation therapy system and method for generating treatment plans therefor - Patent application

Abstract

A radiation therapy system and a treatment plan generation method thereof, the radiation therapy system including a beam irradiation device, a treatment planning module, and a control module. The beam irradiation device generates a treatment beam and irradiates the treatment beam onto a target body to form an irradiated area. The treatment planning module generates a treatment plan based on parameters of the treatment beam and medical image data of the target body. The control module retrieves a treatment plan corresponding to the target body from the treatment planning module and controls the beam irradiation device to sequentially irradiate the target body according to at least two irradiation angles and irradiation times corresponding to each irradiation angle determined by the treatment plan generation method. The radiation therapy system and treatment plan generation method according to the present invention can distribute the radiation dose in shallow parts of the irradiated area and increase the radiation dose in deep parts of the target body so as to reduce the maximum dose to normal tissue and improve the minimum dose to diseased tissue, while ensuring a uniform dose distribution within the diseased tissue.

Description

[Technical field]

[0001] One aspect of the present invention relates to a radiation therapy system, and another aspect of the present invention relates to a method for generating a treatment plan, in particular, a method for generating a treatment plan for a radiation therapy system. [Background technology]

[0002] With the development of atomic science, radiation therapy, such as cobalt-60, linear accelerator, electron beam, etc., has become one of the main means of cancer treatment. However, traditional photon or electron therapy, due to the limitations of the physical conditions of radiation itself, will not only kill tumor cells but also damage many normal tissues along the beam path. In addition, because tumor cells have different degrees of sensitivity to radiation, traditional radiation therapy has poor therapeutic effects on malignant tumors with high radiation resistance (e.g., glioblastoma multiforme, melanoma).

[0003] To reduce radiation damage to normal tissues surrounding tumors, the concept of targeted therapy in chemotherapy has been applied to radiotherapy. In addition, radiation sources with high relative biological effectiveness (RBE) are now being actively developed for tumor cells with high radiation resistance (e.g., proton therapy, heavy ion therapy, neutron capture therapy, etc.). Among these, neutron capture therapy combines the above two concepts. For example, in boron neutron capture therapy (BNCT), boron-containing drugs are specifically concentrated in tumor cells, and combined with precise beam control, it provides a better cancer treatment option compared to conventional radiation.

[0004] Boron neutron capture therapy is a method to 10 B) Taking advantage of the property that the contained drug has a large capture cross section for thermal neutrons, 10 B(n,α) 7 Li neutron capture and fission reactions 4He and 7 It produces two heavy charged particles, Li and Li. The combined range of the two particles is close to the size of a cell, so the radiation damage to the living body can be limited to the cellular level. By selectively collecting boron-containing drugs in tumor cells, combined with an appropriate neutron radiation source, it is possible to achieve the purpose of partially killing tumor cells without causing significant damage to normal tissues.

[0005] Radiotherapy uses high-energy radiation to destroy tumor cells and prevent their growth and division, with minimal side effects that normal tissues and organs can tolerate or recover from. Thus, the dose received by the tumor is limited to the radiation dose that normal tissues and organs can tolerate. At the same time, the effectiveness of boron neutron capture therapy depends on the distribution and accumulation of boron-containing drugs in the tumor and the amount of neutrons accumulated. The distribution and accumulation of boron-containing drugs are affected by tumor characteristics and the metabolic absorption capacity of the patient, and currently, patients suitable for boron neutron capture therapy are scanned and selected by positron emission tomography (PET). The neutron flux decreases with increasing depth through the patient's body, and is suppressed by interacting with the boron element in the drug, resulting in less neutron accumulation in the tumor located further in the direction of the neutrons.

[0006] 3D models are widely applied in the fields of scientific experiment analysis and simulation. For example, in the field of nuclear radiation and protection, in order to simulate the absorbed dose of the human body under certain radiation conditions to help doctors formulate treatment plans, it is always necessary to use computer technology to perform various processing on medical image data to set up the accurate lattice model required for Monte Carlo software, and then combine it with Monte Carlo software to perform simulation calculations. In the current neutron capture therapy planning system, the optimal irradiation angle of the neutron beam is selected by evaluation of the irradiation angle. Since the accumulation of neutrons deep in the tumor is low, it is necessary to increase the radiation dose, but on the other hand, due to the limit of the allowable radiation dose of normal tissues and organs, it is necessary to control the radiation dose, which greatly reduces the therapeutic effect.

[0007] Therefore, there is a need to provide a radiation therapy system and method for generating treatment plans therefor. Summary of the Invention

[0008] In order to overcome the deficiencies of the prior art, in one aspect of the present invention, a radiation therapy system is provided, including a beam irradiation device, a treatment planning module, and a control module. The beam irradiation device generates a treatment beam and irradiates the treatment beam to an irradiated body to form an irradiated area. The treatment planning module generates a treatment plan that determines at least two irradiation angles defined as vector directions from an irradiation point of the treatment beam to a predetermined point in the diseased tissue of the irradiated area and irradiation times corresponding to each irradiation angle based on parameters of the treatment beam generated by the beam irradiation device and medical image data of the irradiated area. The control module calls up the treatment plan corresponding to the irradiated body from the treatment planning module, and controls the beam irradiation device to sequentially irradiate the irradiated body according to the at least two irradiation angles and irradiation times corresponding to each irradiation angle determined by the treatment plan. Based on the distribution of the diseased tissue, radiation therapy is performed at multiple irradiation angles to distribute and reduce the radiation dose in the shallow part of the irradiated area, reduce the radiation dose received by normal tissue and the maximum dose of normal tissue, and further reduce the probability of side effects occurring in normal tissue after radiation therapy. At the same time, the total radiation dose can be appropriately increased to increase the dose in the diseased tissue, especially the radiation dose deep in the diseased tissue, improve the minimum dose in the diseased tissue, and the multiple incidence directions can further make the dose in the diseased tissue more uniform.

[0009] Furthermore, the treatment planning module uses a Monte Carlo simulation program to simulate a radiation dose distribution when the irradiated area is irradiated with the treatment beam, and generates the treatment plan in combination with a mathematical algorithm. Furthermore, the treatment planning module sets an objective function of a region of interest based on the simulated radiation dose distribution, optimizes the objective function, and calculates the at least two irradiation angles and the irradiation times corresponding to each irradiation angle.

[0010] Preferably, the treatment planning module sets a three-dimensional voxel prosthetic tissue model based on the medical image data of the irradiated area, inputs the parameters of the treatment beam and the three-dimensional voxel prosthetic tissue model into the Monte Carlo simulation program to simulate sampling at different irradiation angles, and calculates the radiation dose D received by each voxel unit i per unit time at the sampled irradiation angle k. ki Calculate.

[0011] Preferably, the 3D voxel prosthetic tissue model has information of tissue type and tissue density, which can provide more accurate tissue type, elemental composition and density, and the geometric model set up can more closely match the actual situation reflected by medical image data. Furthermore, the above radiation therapy system is a boron neutron capture therapy system, and the above 3D voxel prosthetic tissue model further has information of tissue boron concentration, which can clearly understand the concentration of boron-containing drug in each tissue, and can more realistically reflect the actual situation when performing irradiation simulation of boron neutron capture therapy.

[0012] Preferably, when sampling different irradiation angles, or after sampling and calculating, the sampled Beam angle may be further selected.

[0013] Furthermore, the above objective function is expressed by Equation 1:

number

number

number

number

number

number

[0014] Furthermore, the treatment planning module uses an optimization algorithm to optimize the objective function according to Equation 4;

number

number

[0015] The treatment planning module further sets constraints on the optimization solution of the objective function, the constraints further including selecting one or more normal organs or tissues M and determining a total dose d for all voxels i in each normal organ or tissue M. i Sum of d M satisfies Equation 6.

number

[0016] Preferably, the units of the dose, total dose and prescribed dose in the above formulas 1 to 6 are eq-Gy, and the unit of the irradiation time is s.

[0017] More preferably, the treatment planning module evaluates or selects results of an optimization solution of the objective function with a dosimetric review.

[0018] In another aspect of the present invention, a radiation therapy system including a beam irradiation device, a treatment planning module, and a control module is provided. The beam irradiation device generates a treatment beam and irradiates the treatment beam to an irradiated body to form an irradiated area. The treatment planning module generates a treatment plan for determining a plurality of irradiation angles defined as vector directions from an irradiation point of the treatment beam to a predetermined point in the diseased tissue of the irradiated area and a planned irradiation dose corresponding to each of the irradiation angles based on parameters of the treatment beam generated by the beam irradiation device and medical image data of the irradiated area. The control module calls up the treatment plan corresponding to the irradiated body from the treatment planning module, and controls the beam irradiation device to irradiate in sequence according to the plurality of irradiation angles and the planned irradiation dose corresponding to each of the irradiation angles in a process of performing one radiation therapy on the irradiated body according to the treatment plan. Based on the distribution of the diseased tissue, radiation therapy is performed at a plurality of irradiation angles to distribute and reduce the radiation dose in the shallow part of the irradiated area, reduce the radiation dose received by normal tissue and the maximum dose of normal tissue, and further reduce the probability of side effects occurring in normal tissue after radiation therapy. At the same time, the total radiation dose can be appropriately increased to increase the dose in the diseased tissue, especially the radiation dose deep in the diseased tissue, improve the minimum dose in the diseased tissue, and the multiple incidence directions can further make the dose in the diseased tissue more uniform.

[0019] In yet another aspect of the present invention, a method for generating a treatment plan is provided, comprising the steps of: establishing a three-dimensional voxel prosthesis tissue model based on medical image data; defining beam parameters in a Monte Carlo simulation program; simulating sampling at different irradiation angles; and calculating a radiation dose D received by each voxel unit i per unit time at the sampled irradiation angle k. kiand setting an objective function of the region of interest, optimizing the objective function, and calculating at least two irradiation angles and irradiation times corresponding to each irradiation angle, the irradiation angle being defined as a vector direction from the irradiation point of the beam to a predetermined point in the diseased tissue of the three-dimensional voxel prosthetic tissue model. By performing radiation therapy at multiple irradiation angles, the radiation dose in the shallow part of the irradiated area can be distributed and reduced, the radiation dose received by the normal tissue and the maximum dose of the normal tissue can be reduced, and the probability of side effects occurring in the normal tissue after radiation therapy can be reduced, and at the same time, the total radiation dose can be appropriately increased so as to increase the dose of the diseased tissue, especially the radiation dose in the deep part of the diseased tissue, and the minimum dose of the diseased tissue can be improved, and the dose in the diseased tissue can be more uniform through multiple incidence directions.

[0020] Preferably, the objective function is expressed using Equation 1 as follows:

number

number

number

number

number

number

[0021] Furthermore, an optimization algorithm is used to optimize the objective function according to Equation 4,

number

number

[0022] Furthermore, the method for generating a treatment plan further includes a step of setting a constraint on the optimization solution of the objective function. The constraint further includes selecting one or more normal organs or tissues M, and determining a total dose d for all voxels i in each normal organ or tissue M. i Sum of d M satisfies Equation 6.

number

[0023] Preferably, the units of the dose, total dose and prescribed dose in the above formulas 1 to 6 are eq-Gy, and the unit of the irradiation time is s.

[0024] Preferably, the method for generating a treatment plan further comprises a step of dose verification in which results of an optimization solution of the objective function are evaluated or screened by dose verification.

[0025] Preferably, the step of setting the three-dimensional voxel prosthesis tissue model based on the medical image data further includes the steps of reading the medical image data, setting the three-dimensional medical image voxel model, defining or reading the boundary of the region of interest, defining the tissue type (element composition) and tissue density of each voxel unit, and setting the three-dimensional voxel prosthesis tissue model. The three-dimensional voxel prosthesis tissue model is set based on the conversion relationship between the medical image data and the tissue type and tissue density, which provides the tissue type (element composition) and tissue density more accurately, and the set geometric model more closely matches the actual situation reflected by the medical image data. Furthermore, the above treatment plan generation method is applied to boron neutron capture therapy, and the step of setting the three-dimensional voxel prosthesis tissue model based on the medical image data further includes the step of defining the tissue boron concentration of each voxel unit, which can clearly understand the concentration of the boron-containing drug in each tissue, and can more realistically reflect the actual situation when performing the irradiation simulation of boron neutron capture therapy.

[0026] Preferably, when sampling different irradiation angles, or after sampling and calculating, the sampled Beam angle may be further selected. Effect of the Invention

[0027] The radiation therapy system and the method for generating a treatment plan according to the present invention can distribute the radiation dose in the shallow part of the irradiated area and increase the radiation dose in the deep part of the diseased tissue, so as to reduce the maximum dose in normal tissue and improve the minimum dose in diseased tissue, while ensuring uniform distribution of the dose within the diseased tissue. [Brief description of the drawings]

[0028] [Figure 1] FIG. 1 is a schematic diagram of a boron neutron capture reaction. [Diagram 2] This is the neutron capture nuclear reaction formula for 10B(n,α)7Li. [Diagram 3] FIG. 1 is a block diagram of a neutron capture therapy system in an embodiment of the present invention. [Figure 4] 1 is a flow chart of a method for generating a treatment plan by a treatment planning module in an embodiment of the present invention. [Diagram 5] 1 is a flow chart of a method for establishing a three-dimensional voxel prosthetic tissue model in an embodiment of the present invention. [Figure 6] 4 is a flowchart of a method for setting an objective function of a region of interest and calculating it through an optimization solution in an embodiment of the present invention. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0029] The embodiments of the present invention will be described in more detail below with reference to the drawings, so that those skilled in the art can practice the invention by referring to the specification.

[0030] Preferably, a neutron capture therapy system and a method for generating a treatment plan therefor are embodiments of the present invention.Neutron capture therapy, and in particular, boron neutron capture therapy, will now be briefly described.

[0031] Neutron capture therapy has been increasingly applied in recent years as an effective means of cancer treatment, among which boron neutron capture therapy has become the most common. Neutrons used in boron neutron capture therapy can be supplied by a nuclear reactor or an accelerator. In the embodiment of the present invention, accelerator boron neutron capture therapy is taken as an example. The basic module of accelerator boron neutron capture therapy generally includes an accelerator used to accelerate charged particles (e.g., protons, deuterium nuclei, etc.), a target, a heat removal system, and a beam shaping assembly. Neutrons are generated by the action of the accelerated charged particles and a metal target, and an appropriate nuclear reaction is selected depending on the required neutron yield and energy, the energy and current magnitude of the accelerated charged particles that can be provided, the physical and chemical properties of the metal target, etc. Nuclear reactions that have been well studied are: 7 Li(p,n) 7 Be and 9Be(p,n) 9 B, and these two reactions are both endothermic. These two nuclear reactions have energy thresholds of 1.881 MeV and 2.055 MeV, respectively. The ideal neutron source for boron neutron capture therapy is epithermal neutrons at the keV energy level, so theoretically, bombardment of a metallic lithium target with protons whose energy is slightly higher than the threshold would produce neutrons with relatively low energy, making clinical application possible without requiring too much moderation. However, the two targets, lithium metal (Li) and beryllium metal (Be), do not have a large cross section for interacting with protons of the threshold energy, so nuclear reactions are generally induced with protons with relatively high energy to generate sufficient neutron flux.

[0032] Boron Neutron Capture Therapy (BNCT) utilizes the property that boron (10B)-containing drugs have a large capture cross section for thermal neutrons. 10 B(n,α) 7 Li neutron capture and fission reactions 4 He and 7 The boron neutron capture reaction produces two types of heavy charged particles, Li and Li. Figures 1 and 2 show a schematic diagram of the boron neutron capture reaction, respectively. 10 B(n,α) 7 The nuclear reaction formula for Li neutron capture is shown below. The two types of charged particles have an average energy of about 2.33 MeV, and are characterized by a high linear energy transfer (LET) and a short range. The linear energy transfer and range of the α particle are 150 keV / μm and 8 μm, respectively. 7 For Li heavy charged particles, the range is 175keV / μm and 5μm, respectively. The combined range of the two types of particles is close to the size of a cell, so that the radiation damage to the living body can be limited to the cellular level. By selectively collecting boron-containing drugs in tumor cells, combined with an appropriate neutron radiation source, it is possible to achieve the purpose of partially killing tumor cells without causing significant damage to normal tissues.

[0033] 3, the radiation therapy system in this embodiment is preferably a neutron capture therapy system 100, and includes a neutron beam irradiation device 10, a treatment planning module 20, and a control module 30. The neutron beam irradiation device 10 includes: Neutron beam generator 11 and treatment table 12. Neutron beam generator 11 The treatment planning module 20 generates a therapeutic neutron beam and irradiates the patient on the treatment table 12 to form an irradiated area. In neutron capture therapy, in order to simulate the absorbed dose of the living body under specific radiation conditions and help the doctor to formulate a treatment plan, it is always necessary to use computer technology to perform various processing on medical images to set an accurate lattice model required for the Monte Carlo software, and perform simulation calculations in combination with the Monte Carlo software. The treatment planning module 20 Neutron beam generator Based on the neutron beam parameters generated by the neutron beam irradiation module 11 and the medical image data of the patient's irradiated area, a Monte Carlo simulation program is used to simulate the radiation dose distribution when the patient undergoes radiation therapy, and a treatment plan is generated in combination with a mathematical algorithm. In one embodiment, the treatment plan module 20 sets an objective function of the region of interest based on the simulated radiation dose distribution, optimizes the objective function, and calculates at least two irradiation angles and the corresponding irradiation times for each irradiation angle. As can be understood, the at least two irradiation angles and the corresponding irradiation times for each irradiation angle may also be calculated by other methods. The control module 30 retrieves a treatment plan corresponding to the current patient from the treatment plan module 20, and controls the irradiation of the neutron beam irradiation device 10 based on the treatment plan, for example, to generate a neutron beam and sequentially irradiate the patient on the treatment couch 12 according to the at least two irradiation angles and the corresponding irradiation times determined by the treatment plan. Neutron beam generator 11. As can be understood, the irradiation time corresponding to each irradiation angle may be the planned irradiation dose corresponding to each irradiation angle, and may be obtained by conversion through simulation calculation.

[0034] Referring to FIG. 4, the method for generating a treatment plan by the treatment plan module 20 according to this embodiment specifically includes steps S410 to S440.

[0035] In S410, a 3D voxel prosthesis tissue model is set based on the medical image data. In S420, the parameters of the beam are defined in a Monte Carlo simulation program (e.g., Monte Carlo N Particle Transport Code, MCNP), and the sampling at different irradiation angles k is simulated to obtain the dose distribution D at different irradiation angles k in unit time of voxel i in the 3D voxel prosthesis tissue model. ki Calculate In S430, an objective function of the region of interest is set, and at least two irradiation angles and irradiation times corresponding to each irradiation angle are calculated by optimizing the solution. The so-called region of interest may be an important organ such as an eye or a liver, an important tissue such as a bone tissue or a brain tissue, or a tumor cell. In S440, a dose review is performed and the results of the optimization solution of the objective function are evaluated or selected.

[0036] Referring to FIG. 5, in one embodiment, the step S410 of establishing a three-dimensional voxel prosthetic tissue model based on medical image data further includes steps S510 to S550.

[0037] In S510, the medical image data is read, In the S520, a 3D medical image voxel model is set up, In S530, a boundary of the region of interest is defined or read; In S540, the tissue type (element composition) and tissue density of each voxel unit are defined, and may be automatically defined based on the conversion relationship between the CT image data and the tissue type and tissue density, or may be manually defined by a user, for example, by giving a specific tissue type and tissue density to a voxel unit within the boundary of each region of interest; In S550, a three-dimensional voxel prosthetic tissue model is established.

[0038] The three-dimensional voxel prosthetic tissue model is set based on the conversion relationship between the medical image data and the tissue type and tissue density, which provides more accurate tissue type (element composition) and tissue density, and the set geometric model more closely matches the actual situation reflected by the medical image data. If the radiation therapy system is a boron neutron capture therapy system, the step S410 of setting the three-dimensional voxel prosthetic tissue model based on the medical image data may further include S560, after S540, defining the tissue boron concentration of each voxel unit. As can be understood, S560 may be performed before S540. The geometric model labeled with the information of the tissue boron concentration can clearly understand the concentration of the boron-containing drug in each tissue, which will more realistically reflect the actual situation when the next neutron irradiation simulation is performed.

[0039] For detailed process of setting a 3D voxel prosthetic tissue model based on medical image data, reference can be made to the patent application disclosed on March 8, 2017, publication number CN106474634A, and title "Method for setting a geometric model based on medical image data", the entire contents of which are incorporated herein by reference.

[0040] The Monte Carlo method is currently a tool that can accurately simulate the collision trajectory and energy distribution of nuclear particles in a three-dimensional space inside a radiation irradiation target, and the combination of a human body model and a Monte Carlo simulation program can accurately calculate and evaluate the absorbed dose of the human body under a radiation environment. In step S420, the parameters of the beam (e.g., energy, intensity, radius, etc. of the beam) are defined in the Monte Carlo simulation program, and the dose distribution at different irradiation angles of the three-dimensional voxel prosthetic tissue model is simulated and calculated by sampling at different irradiation angles, i.e., the radiation dose D received by each voxel unit i in a unit time by the defined beam irradiation at each sampled irradiation angle k is calculated. kiis simulated and calculated.

[0041] If sampling, the starting position of the beam and Beam angle The starting position and angle in the calculation can be determined by a forward algorithm or a backward algorithm. In the forward algorithm, the starting position is determined as an extracorporeal position, and the calculation can be performed by sampling in sequence at fixed angles or distance intervals. It can also be performed by random sampling. Beam angle For the part, it can be set as a vector direction from the irradiation point to the tumor center of gravity or the deepest point of the tumor, and the specific tumor end point position can be adjusted according to the user's needs. In the backward algorithm, the starting position is determined within the tumor range, and the starting position can be the tumor center of gravity, the deepest point or a random point within the tumor range; Beam angle may be determined in a random sampling or interval sampling manner.

[0042] When sampling, Beam angle For example, Beam angle is evaluated and used in subsequent calculations based on the evaluation result. Beam angle Alternatively, we can select Beam angle Further selection may be made, for example, based on the results of radiation dose distribution or Beam angle The selection may be made based on the evaluation results. Beam angle The evaluation method will not be described in detail in this specification, and reference may be made to the patent application disclosed on June 16, 2017, with publication number CN106853272A, and the invention title "Method for Evaluating Beam Irradiation Angle", the entire contents of which are incorporated herein by reference.

[0043] Referring to FIG. 6, step S430 of setting an objective function for an area of ​​interest and calculating at least two irradiation angles and irradiation times corresponding to each irradiation angle through an optimization solution will be described in detail below. In one embodiment, step S430 further includes steps S610 to S630.

[0044] In S610, for a region of interest N (which is a tumor cell in this embodiment), an objective function for the region of interest N is set. In one embodiment, in order to uniformly distribute the dose received by all voxels in the region of interest N, the objective function is the square of the difference between the desired dose (prescription dose) and the calculated dose, and as can be appreciated, other objective functions may also be used. The objective function in this embodiment is expressed using Equation 1 as follows:

number

number

[0045] For a voxel i, the total dose received is d i may be calculated using Equation 2:

number

[0046] Prescribed dose for voxel i

number

number

number

number

[0047] In S620, an optimization algorithm is used to solve and calculate an objective function according to Equation 4 so as to reduce the distribution difference between the desired dose and the calculated dose as much as possible;

number

number

[0048] By using an appropriate optimization algorithm such as support vector machine, response surface methodology, least squares vector regression, etc., the optimal solution for the design variable X can be obtained, and based on the obtained optimal solution, multiple irradiation angles k and irradiation times w corresponding to different irradiation angles k can be calculated. k As can be appreciated, the optimization solution for the objective function may be performed in other ways.

[0049] In S630, a constraint is set so that the optimization solution of the objective function can better meet the treatment needs. In this embodiment, the dose limit value of a normal organ or tissue is set as a constraint, one or more normal organs or tissues M are selected, and the total dose d i Sum of d M satisfies Equation 6.

number

[0050] As can be appreciated, no constraints may be set, or different constraints may be set to seek different optimum solutions and generate different treatment planning tools for selection by an operator, such as a physician.

[0051] After the objective function is optimized, the results of the optimization of the objective function are evaluated or selected in step S440 of the dose inspection. For example, a dose volume histogram (DVH) is used to evaluate or select the multiple irradiation angles k determined by the optimal solution of the design variable X and the irradiation times w corresponding to the different irradiation angles k. k The overlapping dose distribution obtained by simulating the three-dimensional voxel prosthetic tissue model may be evaluated, and the above-mentioned irradiation angle evaluation may be performed. Different optimal solutions obtained under different constraint conditions may be simultaneously evaluated so that an operator such as a doctor can select a treatment planning method that best meets his / her needs. As can be understood, it is not necessary to perform a dose inspection.

[0052] According to the distribution of the tumor, radiation therapy can be performed at multiple irradiation angles to distribute and reduce the neutron dose in the shallow part of the patient's irradiated area, reduce the radiation dose received by normal tissue and the maximum dose of normal tissue, and further reduce the probability of side effects occurring in normal tissue after radiation therapy; at the same time, the total radiation dose can be appropriately increased so as to increase the tumor dose, especially the neutron dose deep in the tumor, improve the minimum dose of the tumor, and the multiple incidence directions can make the dose in the tumor more uniform.

[0053] The units of dose, total dose, and prescribed dose in the above formulas 1 to 6 are eq-Gy, and the unit of irradiation time is s. It can be understood that some simple conversions, simple conversions of dose and time units in the above formulas 1 to 6 are still within the scope of the claims of the present invention. The total number of irradiation angles k is at least two, and the specific number can be manually set or automatically obtained by an algorithm or continuously adjusted in an arc shape. The sampling of irradiation angles k may be on the same side or opposite side of the patient.

[0054] After the treatment planning module 20 determines the treatment planning means through calculation and manual selection by the operator, the control module 30 calls up the treatment plan based on an instruction, and controls the neutron beam irradiation device 10 to sequentially irradiate the patient according to the multiple irradiation angles and corresponding irradiation times determined by the treatment plan. As can be understood, the first irradiation angle and corresponding irradiation time for irradiating the patient may be the irradiation that can give the maximum dose to the tumor, followed by irradiation of other supplementary dose irradiation angles, and adjusting according to the next irradiation angle after the irradiation of the current irradiation angle is completed. Adjustment of the irradiation angle may be achieved by the control module 30 controlling the direction of the beam exit of the neutron beam generating device 11 (e.g., using a rotatable frame), or by controlling the set-up for the patient, which may be achieved by the control module 30 directly controlling the movement of the treatment couch 12 based on a treatment plan, or an operator such as a doctor may set up the patient in a simulation positioning room (not shown) based on the treatment plan, and then manually or automatically adjust the position of the treatment couch 12 and the patient based on the patient set-up determined by simulation positioning in the irradiation room (not shown).

[0055] As can be understood, the present invention can also be applied to other radiation therapy fields that can be simulated by Monte Carlo software well known to those skilled in the art, such as proton, heavy ion, X-ray or gamma ray therapy, in which the neutron beam irradiation device is another radiation beam irradiation device. The present invention can also be applied to other diseases that can be treated by radiation irradiation, such as Alzheimer's disease, rheumatoid arthritis, in which the tumor cells are other diseased tissues. The patient can also be other irradiated subjects.

[0056] The above description of exemplary specific embodiments of the present invention will facilitate understanding of the present invention for those skilled in the art. However, it is obvious that the present invention is not limited to the scope of the specific embodiments. Various changes are obvious to those skilled in the art, and are all within the scope of the claims of the present invention, provided that they are within the spirit and scope of the present invention as defined and determined by the appended claims.

Claims

1. a beam irradiation device that generates a therapeutic beam and irradiates the therapeutic beam onto an irradiated body to form an irradiated area; a treatment planning module that generates a treatment plan based on the parameters of the treatment beam generated by the beam irradiation device and medical image data of the irradiated site, the treatment plan determining at least two irradiation angles defined as vector directions from an irradiation point of the treatment beam to a predetermined point in a diseased tissue of the irradiated site, and an irradiation time corresponding to each irradiation angle; a control module that calls up the treatment plan corresponding to the irradiated object from the treatment plan module, and controls the beam irradiation device so as to sequentially irradiate the irradiated object according to at least two irradiation angles and irradiation times corresponding to each irradiation angle determined by the treatment plan; A neutron capture therapy system, characterized in that a first irradiation angle for irradiating a patient and an irradiation time corresponding to the first irradiation angle are irradiation at a maximum dose applied to a tumor within a range of radiation tolerance for normal tissues and organs, and further, irradiation is performed at other irradiation angles as a supplementary dose.

2. The neutron capture therapy system of claim 1, wherein the treatment planning module uses a Monte Carlo simulation program to simulate a radiation dose distribution when the irradiated area is irradiated with the treatment beam, and generates the treatment plan in combination with a mathematical algorithm.

3. 3. The neutron capture therapy system of claim 2, wherein the treatment planning module sets an objective function for a region of interest based on a simulated radiation dose distribution, optimizes the objective function, and calculates the at least two irradiation angles and irradiation times corresponding to each irradiation angle.

4. The treatment planning module sets a three-dimensional voxel prosthetic tissue model based on the medical image data of the irradiated area, inputs the parameters of the treatment beam and the three-dimensional voxel prosthetic tissue model into the Monte Carlo simulation program to simulate sampling at different irradiation angles, and calculates the radiation dose D received by each voxel unit i per unit time at the sampled irradiation angle k. ｋｉ 4. The neutron capture therapy system according to claim 3, wherein the neutron capture therapy system calculates:

5. The neutron capture therapy system according to claim 4, wherein the neutron capture therapy system is a boron neutron capture therapy system, and the three-dimensional voxel prosthetic tissue model includes information on tissue type, tissue density, and tissue boron concentration.

6. The neutron capture therapy system according to claim 4 , wherein the treatment planning module selects the sampled irradiation angles when sampling different irradiation angles or after sampling and calculating.

7. The objective function is expressed by using Equation 1 as follows: [0010] In the formula, d ｉ is the total dose in voxel i, [0025] is the prescribed dose for voxel i, For a voxel i, the total dose received d ｉ can be calculated using Equation 2: [0030] During the ceremony, w ｋ is the irradiation time at different irradiation angles, and D ｋｉ is the dose at irradiation angle k of voxel i in unit time, and d ｉ is the total dose in voxel i, Prescribed dose for voxel i [0045] can be calculated using Equation 3: [0050] During the ceremony, [006] is the prescribed dose in the region of interest N, and C Ｎ The neutron capture therapy system of claim 4, wherein: is the number of voxels in the region of interest.

8. The treatment planning module uses an optimization algorithm to optimize the objective function according to Equation 4; [0070] The design variable of Equation 4 is defined as X. The design variable X is shown in Equation 5. [0080] Based on the optimal solution of the design variable X, the at least two irradiation angles k and the irradiation times w corresponding to the different irradiation angles k are determined. ｋ 8. The neutron capture therapy system of claim 7, wherein:

9. The treatment planning module sets constraints on the optimization solution of the objective function, the constraints including selecting one or more normal organs or tissues M and calculating a total dose d for all voxels i in each normal organ or tissue M. ｉ Sum of d Ｍ The neutron capture therapy system according to claim 8, characterized in that: [0090]

10. The neutron capture therapy system of claim 3 , wherein the treatment planning module evaluates or screens results of the optimization solution of the objective function by dosimetry.

Images