(19)
(11)EP 3 369 458 B1

(12)EUROPEAN PATENT SPECIFICATION

(45)Mention of the grant of the patent:
10.06.2020 Bulletin 2020/24

(21)Application number: 16872229.6

(22)Date of filing:  18.10.2016
(51)Int. Cl.: 
A61N 5/10  (2006.01)
(86)International application number:
PCT/CN2016/102312
(87)International publication number:
WO 2017/097035 (15.06.2017 Gazette  2017/24)

(54)

METHOD FOR EVALUATING IRRADIATION ANGLE OF BEAM

VERFAHREN ZUR AUSWERTUNG DES BESTRAHLUNGSWINKELS EINES STRAHLS

PROCÉDÉ D'ÉVALUATION DE L'ANGLE D'IRRADIATION D'UN FAISCEAU


(84)Designated Contracting States:
AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

(30)Priority: 08.12.2015 CN 201510894752

(43)Date of publication of application:
05.09.2018 Bulletin 2018/36

(73)Proprietor: Neuboron Medtech Ltd.
Nanjing, Jiangsu 211112 (CN)

(72)Inventors:
  • LIU, Yuan-Hao
    Nanjing, Jiangsu 211112 (CN)
  • CHEN, Wei-Lin
    Nanjing, Jiangsu 211112 (CN)

(74)Representative: Sun, Yiming 
HUASUN Patent- und Rechtsanwälte Friedrichstraße 33
80801 München
80801 München (DE)


(56)References cited: : 
WO-A1-2011/042819
CN-A- 102 068 763
US-A- 5 976 066
CN-A- 101 014 383
US-A- 5 341 292
US-A1- 2007 034 812
  
      
    Note: Within nine months from the publication of the mention of the grant of the European patent, any person may give notice to the European Patent Office of opposition to the European patent granted. Notice of opposition shall be filed in a written reasoned statement. It shall not be deemed to have been filed until the opposition fee has been paid. (Art. 99(1) European Patent Convention).


    Description

    FIELD OF THE DISCLOSURE



    [0001] The present disclosure relates to a method for evaluating irradiating angles, and, more particularly, to a method for evaluating irradiating angles of a beam.

    [0002] In particular the disclosure relates to a method for determining an irradiation angle for treating a tumor in a body including evaluating sample irradiation angles of a beam defined as being the vector direction from an irradiation point of the beam to a pre-set point of the tumor, the method comprising: a step of sampling the sample irradiation angles of the beam by selecting the irradiation points; a step of calculating the track of the beam passing through the body along an irradiation track to an effective treatment depth of each sample irradiation angle; and a step of determining the irradiation angle from the sample irradiation angles.

    BACKGROUND OF THE DISCLOSURE



    [0003] US 2007/0034812 A1 discloses a method of the generic type defined in the preceding chapter. Specifically US 2007/0034812 A1 discloses a method of optimizing a laser-accelerated proton radiation dose to a targeted region are disclosed. Disclosed methods include providing a plurality of modulated poly-energetic proton beamlets and irradiating the targeted region with the plurality of modulated beamlets. An intensity or weight of each beamlet may be adjusted accordingly to achieve a uniform and conform dose distribution.

    [0004] CN 106474634 A discloses a geometric model building method based on medical image data. The method comprises steps of reading medical image data, defining a tissue type through the conversion relationship between the medical image data and tissue types, deciding a tissue clustering number, defining the tissue density through the conversion relationship between the medical image data and densities, building a 3D coding matrix with the tissue and the density information, generating a geometric model and judging whether the medical image voxel is in an ROI boundary.

    [0005] As atomics moves ahead, such radiotherapy as Cobalt-60, linear accelerators and electron beams has been one of major means to cancer therapy. However, conventional photon or electron therapy has been undergone physical restrictions of radioactive rays; for example, many normal tissues on a beam path will be damaged as tumor cells are destroyed. On the other hand, sensitivity of tumor cells to the radioactive rays differs greatly, so in most cases, conventional radiotherapy falls short of treatment effectiveness on radioresistant malignant tumors (such as glioblastoma multiforme and melanoma).

    [0006] For the purpose of reducing radiation damage to the normal tissue surrounding a tumor site, target therapy in chemotherapy has been employed in the radiotherapy. While for high-radioresistant tumor cells, radiation sources with high RBE (relative biological effectiveness) including such as proton, heavy particle and neutron capture therapy have also developed. Among them, the neutron capture therapy combines the target therapy with the RBE, such as the boron neutron capture therapy (BNCT). By virtue of specific grouping of boronated pharmaceuticals in the tumor cells and precise neutron beam regulation, BNCT is provided as a better cancer therapy choice than conventional radiotherapy.

    [0007] BNCT takes advantage that the boron (10B)-containing pharmaceuticals have high neutron capture cross section and produces 4He and 7Li heavy charged particles through 10B(n,α)7Li neutron capture and nuclear fission reaction. As illustrated in FIGS. 1 and 2, a schematic drawing of BNCT and a nuclear reaction formula of 10B (n, α) 7Li neutron capture are shown, the two charged particles, with average energy at about 2.33MeV, are of linear energy transfer (LET) and short-range characteristics. LET and range of the alpha particle are 150keV/micrometer and 8 micrometers respectively while those of the heavy charged particle 7Li are 175keV/micrometer and 5 micrometers respectively, and the total range of the two particles approximately amounts to a cell size. Therefore, radiation damage to living organisms may be restricted at the cells' level. When the boronated pharmaceuticals are gathered in the tumor cells selectively, only the tumor cells will be destroyed locally with a proper neutron source on the premise of having no major normal tissue damage.

    [0008] In the existing neutron capture treatment planning system, the geometric angle of irradiation is determined and defined artificially based on experience. Due to the structure of the human body is quite complex and the sensitivity of various tissues or organs to radiation is also different. Therefore, a better angle of irradiation is likely to be ignored by relying on human judgment only, resulting in poor therapeutic effect. In order to achieve the optimization of the therapeutic effect, the irradiation angle of the beam is an essential part to be considered.

    [0009] Therefore, it is necessary to propose a method for evaluating the irradiation angle of a beam.

    SUMMARY



    [0010] The invention is defined in claim 1. Further aspects and preferred embodiments are defined in the dependent claims. Aspects, embodiments and examples of the present disclosure which do not fall under the scope of the appended claims do not form part of the claimed invention and are merely provided for illustrative purposes.

    [0011] In order to overcome defects of prior art, the inventor of the present disclosure developed an executable method for optimizing the irradiation angle, which can be used as a favorable reference, and in combination with doctor's experience, for finding out the best irradiation angle as far as possible. Implementation of the routing optimization is carried out by assigning multiple lines of tracks from the tumor to the body surface, and calculating the evaluation coefficient of individual track taking into account the proportion of organs in the track and their radiation sensitive coefficients, by a forward calculation method (from the body surface to the tumor) or a reverse calculation method (from the tumor to the body surface) with selecting points sequentially or randomly. Through the above method, the evaluation score for the incidence into the tumor from each point on the surface can be calculated, which can be fused with 2D or 3D images to help the user find the best point of incidence.

    [0012] The present invention provides a method for determining an irradiation angle for treating a tumor in a body including evaluating sample irradiation angles of a beam defined as being the vector direction from an irradiation point of the beam to a pre-set point of the tumor, the method comprising: a step of sampling the sample irradiation angles of the beam by selecting the irradiation points; a step of calculating the track of the beam passing through the body along an irradiation track to an effective treatment depth for each sample irradiation angle; and a step of determining the irradiation angle from the sample irradiation angles. In addition the beam is a neutron beam, and it is determined for each sample irradiation angle whether the tumor is fully covered within the effective treatment depth: if yes, entering steps of calculating an evaluation coefficient of the sampled irradiation angle, recording irradiation conditions and calculating results, and returning to the step of sampling the irradiation angle of the beam; and if no, entering a step of setting a worst evaluation coefficient of the irradiation angle, and returning to the step of sampling the irradiation angle of the beam.

    [0013] According to the method for evaluating the irradiation angle of the beam, it can be clearly recognized whether the performance of the beam irradiating at a certain position and at a certain angle is good or bad, thus providing strong data support for the doctor or the physicists to decide on the irradiation mode. The "pre-set point of the tumor" can be set as the centroid of the tumor or the deepest point in the tumor, and the location of the specific pre-set point of the tumor can be adjusted according to the user's needs. As preferred, the irradiation angle of the beam is defined as the vector direction from the irradiation point of the beam to the centroid of the tumor or the deepest point in the tumor. Certainly, it is well known to those skilled in the art that the irradiation angle of the beam can also be customized according to user's needs.

    [0014] Certainly, it is well known to those skilled in the art that the above vector direction includes the direction of the positive vector from the irradiation point of the beam to the pre-set point of the tumor, and the direction of the negative vector.

    [0015] In order to more accurately calculate the evaluation coefficient, calculating the evaluation coefficient based on the beam characteristics, the radiation sensitivity coefficient, and the boron concentration of an organ.

    [0016] in the sampled irradiation angle and irradiation track, a weighting factor (W(i)) of the organ (i) is calculated using the formula I:

    Wherein, I(i), S(i) and C(i) are the beam intensity, the radiation sensitive coefficient of the organ (i) and the boron concentration of the organ (i), respectively.

    [0017] Further, the I(i) is calculated using Formula II that integrates the depth intensity or dose curve of the human body according to the beam used:

    Wherein i(x) is the depth intensity or dose curve function of the therapeutic beam in an approximate human body and x0-x is the depth range of the organ (i) in the beam track.

    [0018] The evaluation coefficient is calculated using Formula III:

    Wherein Q(x, y, z, Φ, θ) as the evaluation coefficient is the sum of the weighting factors of each organ in the organ track.

    [0019] Certainly, another way of presenting the evaluation coefficient, that is, the ratio of organ evaluation coefficient to tumor evaluation coefficient, can also be used.

    [0020] The ratio (QR(x, y, z, Φ, θ)) of the evaluation coefficient to the tumor evaluation coefficient is calculated using Formula IV:

    Wherein, W(tumor) is the weighting factor of the tumor.

    [0021] The "medical image data" may be Magnetic Resonance Imaging (MRI) or Computed Tomography (CT) or Positron Emission computed Tomography (PET-CT). However, it is well known to those skilled in the art that other medical image data can also be used as long as the medical image data can define the materials and densities of organs, tissues and tumors, it can be applied to the method for evaluating the irradiation angle of the beam disclosed in the present disclosure. As preferred, the step of reading medical image data is a step of reading at least one kind of medical image data selected from the group consisting of CT image data, MRI image data, or PET-CT image data.

    [0022] "Organ track" refers to a beam track through skin, bone, tissue and tumor when irradiated a beam at a certain position and at a certain angle. For example, the organ track is a track of a beam sequentially passing through the skin, bone, tissue, tumor, tissue, bone, and skin. In some calculations, for example, when the weighting factors of the organs are added to a total, the weighting factor of the tumor may not be included. As preferred, the organ track is a track of a beam sequentially through the skin, skull, brain tissue, tumor, brain tissue, skull, and skin. This clearly shows the track of the beam through the brain. Certainly, it is well known to those skilled in the art that the organ track may be a beam sequentially through other part of the human body other than the brain, for example, the liver.

    [0023] As preferred, the method for evaluating the irradiation angle of the beam further includes a step of displaying each evaluation coefficient in a 3D image. Certainly, those skilled in the art can also display each evaluation coefficient in other ways, as long as the doctor or the physicists can identify each evaluation coefficient displayed.

    [0024] The "beam" may be one or more radioactive beams. Preferably, the beam is a mixed beam of neutron beam and gamma beam, or may be an individual neutron beam, an individual proton beam or an individual heavy particle beam.

    [0025] The method for evaluating the irradiation angle of the beam further includes a step of reading the medical image data; a step of defining or reading the profile of an organ, tissue and tumor; and a step of defining the material and the density of the organ, tissue and tumor.

    [0026] Other advantages, objects, and features of the disclosure will be apparent to those skilled in the art from the following description, taken in part by the research and practice of the disclosure.

    BRIEF DESCRIPTION OF THE DRAWINGS



    [0027] 

    FIG. 1 is a boron neutron capture reaction schematic.

    FIG. 2 is a 10B(n,α) 7Li neutron capture nuclear reaction formula.

    FIG. 3 is a logic block diagram of the method for evaluating the irradiation angle of a beam in an embodiment of the present disclosure.

    FIG. 4 is a schematic diagram of an organ track when a beam is irradiated according to an embodiment of the present disclosure.


    DETAILED DESCRIPTION



    [0028] The following description of the preferred embodiments is merely exemplary in nature and is in no way intended to limit the invention, its application, or uses.

    [0029] The embodiments of the present disclosure are further described in detail below with reference to the accompanying drawings, so that those skilled in the art can implement the technical solutions according to the description.

    [0030] As preferred, a method for evaluating the irradiation angle of a beam for neutron capture therapy is taken as an embodiment of the present disclosure. The following will briefly introduce neutron capture therapy, especially boron neutron capture therapy.

    [0031] Neutron capture therapy (NCT) has been increasingly practiced as an effective cancer curing means in recent years, and BNCT is the most common. Neutrons for NCT may be supplied by nuclear reactors or accelerators. Take AB-BNCT for example, its principal components include, in general, an accelerator for accelerating charged particles (such as protons and deuterons), a target, a heat removal system and a beam shaping assembly. The accelerated charged particles interact with the metal target to produce the neutrons, and suitable nuclear reactions are always determined according to such characteristics as desired neutron yield and energy, available accelerated charged particle energy and current and materialization of the metal target, among which the most discussed two are 7Li (p, n) 7Be and 9Be (p, n) 9B and both are endothermic reaction. Their energy thresholds are 1.881MeV and 2.055MeV respectively. Epithermal neutrons at a keV energy level are considered ideal neutron sources for BNCT. Theoretically, bombardment with lithium target using protons with energy slightly higher than the thresholds may produce neutrons relatively low in energy, so the neutrons may be used clinically without many moderations. However, Li (lithium) and Be (beryllium) and protons of threshold energy exhibit not high action cross section. In order to produce sufficient neutron fluxes, high-energy protons are usually selected to trigger the nuclear reactions.

    [0032] BNCT takes advantage that the boron (10B)-containing pharmaceuticals have high neutron capture cross section and produces 4He and 7Li heavy charged particles through 10B(n,α)7Li neutron capture and nuclear fission reaction. As illustrated in FIGS. 1 and 2, a schematic drawing of BNCT and a nuclear reaction formula of 10B (n, α) 7Li neutron capture are shown, the two charged particles, with average energy at about 2.33MeV, are of linear energy transfer (LET) and short-range characteristics. LET and range of the alpha particle are 150keV/micrometer and 8 micrometers respectively while those of the heavy charged particle 7Li are 175keV/micrometer and 5 micrometers respectively, and the total range of the two particles approximately amounts to a cell size. Therefore, radiation damage to living organisms may be restricted at the cells' level. When the boronated pharmaceuticals are gathered in the tumor cells selectively, only the tumor cells will be destroyed locally with a proper neutron source on the premise of having no major normal tissue damage.

    [0033] No matter BNCT neutron sources are from the nuclear reactor or the nuclear reactions between the accelerator charged particles and the target, only mixed radiation fields are produced, that is, beams include neutrons and photons having energies from low to high. As for BNCT in the depth of tumors, except the epithermal neutrons, the more the residual quantity of radiation ray is, the higher the proportion of nonselective dose deposition in the normal tissue is. Therefore, radiation causing unnecessary dose should be lowered down as much as possible. Besides air beam quality factors, dose is calculated using a human head tissue prosthesis in order to understand dose distribution of the neutrons in the human body. The prosthesis beam quality factors are later used as design reference to the neutron beams, which is elaborated hereinafter.

    [0034] The International Atomic Energy Agency (IAEA) has given five suggestions on the air beam quality factors for the clinical BNCT neutron sources. The suggestions may be used for differentiating the neutron sources and as reference for selecting neutron production pathways and designing the beam shaping assembly, and are shown as follows:

    Epithermal neutron flux > 1 x 109 n/cm2s

    Fast neutron contamination < 2 x 10-13 Gy-cm2/n

    Photon contamination < 2 x 10-13 Gy-cm2/n

    Thermal to epithermal neutron flux ratio < 0.05

    Epithermal neutron current to flux ratio > 0.7



    [0035] Note: the epithermal neutron energy range is between 0.5eV and 40keV, the thermal neutron energy range is lower than 0.5eV, and the fast neutron energy range is higher than 40keV.

    1. Epithermal neutron flux



    [0036] The epithermal neutron flux and the concentration of the boronated pharmaceuticals at the tumor site codetermine clinical therapy time. If the boronated pharmaceuticals at the tumor site are high enough in concentration, the epithermal neutron flux may be reduced. On the contrary, if the concentration of the boronated pharmaceuticals in the tumors is at a low level, it is required that the epithermal neutrons in the high epithermal neutron flux should provide enough doses to the tumors. The given standard on the epithermal neutron flux from IAEA is more than 109 epithermal neutrons per square centimeter per second. In this flux of neutron beams, therapy time may be approximately controlled shorter than an hour with the boronated pharmaceuticals. Thus, except that patients are well positioned and feel more comfortable in shorter therapy time, and limited residence time of the boronated pharmaceuticals in the tumors may be effectively utilized.

    2. Fast neutron contamination



    [0037] Unnecessary dose on the normal tissue produced by fast neutrons are considered as contamination. The dose exhibit positive correlation to neutron energy, hence, the quantity of the fast neutrons in the neutron beams should be reduced to the greatest extent. Dose of the fast neutrons per unit epithermal neutron flux is defined as the fast neutron contamination, and according to IAEA, it is supposed to be less than 210-13Gy-cm2/n.

    3. Photon contamination (gamma-ray contamination)



    [0038] Gamma-ray long-range penetration radiation will selectively result in dose deposit of all tissues in beam paths, so that lowering the quantity of gamma-ray is also the exclusive requirement in neutron beam design. Gamma-ray dose accompanied per unit epithermal neutron flux is defined as gamma-ray contamination which is suggested being less than 210-13Gy-cm2/n according to IAEA.

    4. Thermal to epithermal neutron flux ratio



    [0039] The thermal neutrons are so fast in rate of decay and poor in penetration that they leave most of energy in skin tissue after entering the body. Except for skin tumors like melanocytoma, the thermal neutrons serve as neutron sources of BNCT, in other cases like brain tumors, the quantity of the thermal neutrons has to be lowered. The thermal to epithermal neutron flux ratio is recommended at lower than 0.05 in accordance with IAEA.

    5. Epithermal neutron current to flux ratio



    [0040] The epithermal neutron current to flux ratio stands for beam direction, the higher the ratio is, the better the forward direction of the neutron beams is, and the neutron beams in the better forward direction may reduce dose surrounding the normal tissue resulted from neutron scattering. In addition, treatable depth as well as positioning posture is improved. The epithermal neutron current to flux ratio is better of larger than 0.7 according to IAEA.

    [0041] The prosthesis beam quality factors are deduced by virtue of the dose distribution in the tissue obtained by the prosthesis according to a dose-depth curve of the normal tissue and the tumors. The three parameters as follows may be used for comparing different neutron beam therapy effects.

    1. Advantage depth



    [0042] Tumor dose is equal to the depth of the maximum dose of the normal tissue. Dose of the tumor cells at a position behind the depth is less than the maximum dose of the normal tissue, that is, boron neutron capture loses its advantages. The advantage depth indicates penetrability of neutron beams. Calculated in cm, the larger the advantage depth is, the larger the treatable tumor depth is.

    2. Advantage depth dose rate



    [0043] The advantage depth dose rate is the tumor dose rate of the advantage depth and also equal to the maximum dose rate of the normal tissue. It may have effects on length of the therapy time as the total dose on the normal tissue is a factor capable of influencing the total dose given to the tumors. The higher it is, the shorter the irradiation time for giving a certain dose on the tumors is, calculated by Gy/mA-min.

    3. Advantage ratio



    [0044] The average dose ratio received by the tumors and the normal tissue from the brain surface to the advantage depth is called as advantage ratio. The average ratio may be calculated using dose-depth curvilinear integral. The higher the advantage ratio is, the better the therapy effect of the neutron beams is.

    [0045] To provide comparison reference to design of the beam shaping assembly, we also provide the following parameters for evaluating expression advantages and disadvantages of the neutron beams in the embodiments of the present disclosure except the air beam quality factors of IAEA and the abovementioned parameters.
    1. 1. Irradiation time <=30min (proton current for accelerator is 10mA)
    2. 2. 30.0RBE-Gy treatable depth >=7cm
    3. 3. The maximum tumor dose>=60.0RBE-Gy
    4. 4. The maximum dose of normal brain tissue<=12.5RBE-Gy
    5. 5. The maximum skin dose<=11.0RBE-Gy


    [0046] Note: RBE stands for relative biological effectiveness. Since photons and neutrons express different biological effectiveness, the dose above should be multiplied with RBE of different tissues to obtain equivalent dose.

    [0047] Please refer to FIG. 3 and FIG. 4, in order to provide reference basis for the irradiation angle of the beam for neutron capture therapy, in the following, calculate each possible irradiation angle by random or one by one, and calculate the evaluation coefficients by Formula I and Formula II, and display each evaluation coefficient as a 3D image, which is convenient for the doctor or physician to judge the irradiation angle of the treatment.

    [0048] The establishment step of the 3D evaluation coefficient image will be described in detail below, which is a preferred implementation step, and it is well known to those skilled in the art that it is not necessary to follow the preferred implementation step. Specifically, the step is: read the images of patients with definite human anatomy, such as CT/MRI/PET-CT, to define the profile of each organ, tissue and tumor one by one, and to assign the type and density of the material. After the completion of the definition of the geometric material and the density, it is decided to calculate the starting position and the angle of the beam. The determination of the starting position and angle in the calculation can be a forward algorithm or a reverse algorithm. In the forward algorithm, the starting position is determined in an in-vitro position and can be sampled sequentially at a fixed angle or distance interval, or sampled randomly; regarding the beam angle, it can be set as the vector direction from the irradiation point of the beam to the centroid of the tumor or the deepest point in the tumor, and the location of the specific endpoint of the tumor can be adjusted according to user's needs. In the reverse algorithm, the starting position is determined in the range of a tumor, and the starting position can be the tumor centroid, the deepest or a random point within the tumor, and the beam angle can be sampled at random intervals or at specified intervals. After the position and angle of the beam are decided, calculate the organ track of the beam passing through the organs, that is, calculate the type and thickness of the organ which the beam passes after entering the human body. After obtaining the track information of the beam axis passing through the human body, judge whether the tumor range falls within the maximum treatable depth range. If yes, calculate the evaluation coefficient of the track based on the track information, the boron concentration in the organ, the radiation sensitive coefficient of the organ and the beam characteristic information set by the user; If not, give the worst evaluation coefficient, and re-sample the beam position and irradiation angle. After the calculation of the evaluation coefficient is completed, record of the irradiation position, angle and evaluation coefficient. After a certain number of calculations are repeated, outputting a report and provided to the doctor or physician as a reference to determine the irradiation pattern. The data can also be expressed in the form of 3D images to more conveniently judge the pros and cons based on the evaluation coefficients around the beam axis.

    [0049] The evaluation coefficient is calculated based on the beam characteristics, the organ radiation sensitive coefficient and the boron concentration in the organ. Corresponding to a position and an irradiation track, the weighting factor (W(i)) of organ i is calculated as shown in Formula I, wherein, I(i), S(i) and C(i) are the beam intensity, the radiation sensitive coefficient of the organ i and the boron concentration of the organ i, respectively.



    [0050] Wherein, the I(i) is obtained by integrating the depth intensity or dose curve of the human body according to the beam used, as shown in Formula II. In the formula, i(x) is the depth intensity or dose curve function of the therapeutic beam in an approximate human body and x0-x is the depth range of the organ i in the beam track.



    [0051] Through the above algorithm, calculating the weighting factors of each organ in the organ track sequentially calculated and then summing up to obtain the evaluation coefficient corresponding to the beam, as shown in Formula 3. In this calculation, the weighting factor of the tumor should not be included in the calculation.



    [0052] According to the level of the above evaluation coefficient, the degree of damage to the normal tissue during the treatment can be more intuitively determined. In addition to the evaluation coefficient, the evaluation ratio coefficient can also be used to evaluate the irradiation position and angle, which is defined as a ratio of the evaluation coefficient to the weighting factor of the tumor, as shown in Formula IV. In this way, the expected therapeutic effect of the irradiation position and angle can be more fully revealed.



    [0053] The steps involved in the above embodiments - to read the images of patients with definite human anatomy, such as CT/MRI/PET-CT, to define the profile of each organ, tissue and tumor one by one, and to assign the type and density of the material - may refer to the patent application published as CN 106474634 A, filed by the applicant on November17, 2015 to the State Intellectual Property Office with the application No. 201510790248.7, and entitled "Geometric Model Establishment Method Based on Medical Image Data".

    [0054] It is well known to persons skilled in the art that some simple transformations in the above Formula I to Formula IV are still within the claimed scope of the present disclosure, for example, the multiplication of I(i), S(i) and C(i) changes to addition; I(i), S(i) and C(i) are respectively multiplied by the n-th power, n may be an integer multiple of 1 or other multiplicities, depending on the situation; i(x) may be the average or the medium of x0-x times (x0-x), or any calculation method that can achieve matching to the result of intensity integral calculation.

    [0055] While illustrative embodiments of the disclosure have been described above in order to provide a person of ordinary skill in the art with an understanding of the disclosure, it should be clear that the disclosure is not to be limited in scope by the specific embodiments.


    Claims

    1. A computer-implemented method for determining an irradiation angle in a treatment planning system for treating a tumor in a body including evaluating sample irradiation angles of a beam defined as being each a vector direction from an irradiation point of the beam to a pre-set point of the tumor, the method comprising:

    a step of sampling the sample irradiation angles of the beam by selecting the irradiation points;

    a step of calculating the track of the beam passing through the body along an irradiation track to an effective treatment depth for each sample irradiation angle; and

    a step of determining the irradiation angle from the sample irradiation angles; wherein

    for each sample irradiation angle it is determined whether the tumor is fully covered within the effective treatment depth:

    if yes, entering steps of calculating an evaluation coefficient of the sample irradiation angle, recording irradiation conditions and calculating results, and returning to the step of sampling the irradiation angle of the beam; and

    if no, entering a step of setting a worst evaluation coefficient of the irradiation angle, and returning to the step of sampling the irradiation angle of the beam;

    characterized in that the beam is a neutron beam.


     
    2. The method according to claim 1, wherein the evaluation coefficient is calculated based on beam characteristics, radiation sensitivity coefficients of tissue within the body and irradiated by the beam, and a boron concentration within the body.
     
    3. The method according to claim 2 wherein in the sampled irradiation angle and irradiation
    track, a weighting factor (W(i)) of an irradiated organ (i) is calculated using the formula I:

    Wherein, I(i), S(i) and C(i) are the beam intensity, the radiation sensitive coefficient of the organ (i) and the boron concentration of the organ (i), respectively.
     
    4. The method according to claim 3, wherein the I(i) is calculated using Formula II that integrates a depth intensity or dose curve of the body according to the beam used:

    Wherein i(x) is the depth intensity or dose curve function of the beam in an approximate body and x0-x is the depth range of the organ (i) in the irradiation track.
     
    5. The method according to claim 3 or 4, wherein the evaluation coefficient is calculated using Formula III:

    Wherein Q(x, y, z, Φ, θ) as the evaluation coefficient is the sum of the weighting factors of each organ in the irradiation track.
     
    6. The method according to claim 5, wherein the ratio (QR(x, ,z, Φ, θ)) of the evaluation coefficient to a tumor evaluation coefficient of the tumor is calculated using Formula IV:

    Wherein, W(tumor) is a weighting factor of the tumor.
     
    7. The method according to any one of claims 1 to 4, wherein the method further comprises a step of reading medical image data selected from the group consisting of CT image data, MRI image data, or PET-CT image data.
     
    8. The method according to any one of claims 1 to 4, wherein the method further comprises a step of displaying each evaluation coefficient in a 3D image.
     
    9. The method according to claim 7, wherein the method further comprises a step of defining or reading the profile of an organ, tissue and tumor in the body; and a step of defining the material and the density of the organ, tissue and tumor.
     


    Ansprüche

    1. Computerrealisiertes Verfahren zum Bestimmen eines Bestrahlungswinkels in einem Behandlungsplanungssystem zum Behandeln eines Tumors in einem Körper, der das Beurteilen von Probebestrahlungswinkeln eines Strahls umfasst, die jeweils als eine Vektorrichtung von einer Bestrahlungsstelle des Strahls zu einer voreingestellten Stelle des Tumors definiert sind, wobei das Verfahren umfasst:

    einen Schritt des Ausprobierens der Probebestrahlungswinkel des Strahls durch Auswählen der Bestrahlungsstellen;

    einen Schritt des Berechnens der Bahn des Strahls, der durch den Körper einer Bestrahlungsbahn entlang bis zu einer wirksamen Behandlungstiefe für jeden Probebestrahlungswinkels hindurchgeht; und

    einen Schritt des Bestimmens des Bestrahlungswinkels von den Probebestrahlungswinkeln

    wobei für jeden Probebestrahlungswinkel bestimmt wird, ob der Tumor innerhalb der wirksamen Behandlungstiefe vollständig bedeckt ist:

    falls ja, Eingeben von Schritten zum Berechnen der Berechnungsschritte eines Bewertungskoeffizienten des Probebestrahlungswinkels, Aufzeichnen der Bestrahlungsbedingungen und Berechnen der Ergebnisse und Zurückkehren zum Schritt des Ausprobierens des Bestrahlungswinkels des Strahls; und

    falls nein, Eingeben eines Schritts des Einstellens eines schlechtesten Bewertungskoeffizienten des Bestrahlungswinkels und Zurückkehren zum Schritt des Ausprobierens des Bestrahlungswinkels des Strahls;

    dadurch gekennzeichnet, dass der Strahl ein Neutronenstrahl ist.


     
    2. Verfahren nach Anspruch 1, wobei der Bewertungskoeffizient auf der Basis von Strahlcharakteristiken, Strahlungsempfindlichkeitskoeffizienten von Gewebe innerhalb des Körpers und von dem Strahl bestrahlt und einer Borkonzentration innerhalb des Körpers berechnet wird.
     
    3. Verfahren nach Anspruch 2, wobei für den ausprobierten Bestrahlungswinkel und die Bestrahlungsbahn ein Gewichtungsfaktor (W(i)) eines bestrahlten Organs (i) unter Anwendung der Formel I berechnet wird:

    wobei I(i), S(i) und C(i) die Strahlintensität, der Strahlungsempfindlichkeitskoeffizient des Organs (i) bzw. die Borkonzentration des Organs (i) sind.
     
    4. Verfahren nach Anspruch 3, wobei I(i) unter Anwendung der Formel II berechnet wird, die eine Tiefenintensität oder Dosiskurve des Körpers dem verwendeten Strahl entsprechend integriert:

    wobei i(x) die Tiefenintensität oder Dosiskurvenfunktion des Strahls in einem angenäherten Körper ist und xo-x der Tiefenbereich des Organs (i) in der Bestrahlungsbahn ist.
     
    5. Verfahren nach Anspruch 3 oder 4, wobei der Bewertungskoeffizient unter Anwendung der Formel III berechnet wird:

    wobei Q(x, y, z, Φ, θ) als Bewertungskoeffizient die Summe der Gewichtungsfaktoren jedes Organs in der Bestrahlung Bahn ist.
     
    6. sich Verfahren nach Anspruch 5, wobei das Verhältnis von (QR(x, z, Φ, θ)) des Bewertungskoeffizienten zu einem Tumorbewertungskoeffizienten des Tumors unter Anwendung der Formel IV berechnet wird:

    wobei W(Tumor) ein Gewichtungsfaktor des Tumors ist.
     
    7. Verfahren nach einem der Ansprüche 1 bis 4, wobei das Verfahren ferner einen Schritt des Lesens medizinischer Bilddaten ausgewählt aus der Gruppe bestehend aus CT-Bilddaten, MRI-Bilddaten oder PET-CT Bilddaten umfasst.
     
    8. Verfahren nach einem der Ansprüche 1 bis 4, wobei das Verfahren ferner einen Schritt des Anzeigens jedes Bewertungskoeffizienten in einem 3D-Bild umfasst.
     
    9. Verfahren nach Anspruch 7, wobei das Verfahren ferner einen Schritt des Definierens oder Lesens des Profils eines Organs, Gewebes oder Tumors in dem Körper und einen Schritt des Definierens des Materials und der Dichte des Organs, Gewebes oder Tumors umfasst.
     


    Revendications

    1. Procédé mis en œuvre par un ordinateur permettant de déterminer un angle de rayonnement dans un système de planification de traitement destiné au traitement d'une tumeur dans un corps comprenant des angles de rayonnement d'échantillon d'évaluation d'un faisceau défini comme étant chacun une direction de vecteur à partir d'un point de rayonnement du faisceau vers un point prédéfini de la tumeur, le procédé comprenant :

    une étape d'échantillonnage des angles de rayonnement d'échantillon du faisceau en choisissant les points de rayonnement ;

    une étape de calcul de la piste du faisceau passant à travers le corps le long d'une piste de rayonnement à une profondeur de traitement efficace pour chaque angle de rayonnement d'échantillon ; et

    une étape de détermination de l'angle de rayonnement à partir des angles de rayonnement d'échantillon ;

    dans lequel pour chaque angle de rayonnement d'échantillon, il est déterminé si la tumeur est entièrement recouverte à l'intérieur de la profondeur de traitement efficace ;

    si oui, l'entrée d'étapes de calcul d'un coefficient d'évaluation de l'angle de rayonnement d'échantillon, l'enregistrement des conditions de rayonnement et le calcul des résultats,

    et le retour à l'étape d'échantillonnage de l'angle de rayonnement du faisceau ; et

    si non, l'entrée d'une étape de réglage d'un pire coefficient d'évaluation de l'angle de rayonnement, et le retour à l'étape d'échantillonnage de l'angle de rayonnement du faisceau ;

    caractérisé en ce que le faisceau est un faisceau de neutrons.


     
    2. Procédé selon la revendication 1, dans lequel le coefficient d'évaluation est calculé en fonction des caractéristiques de faisceau, des coefficients de sensibilité de rayonnement de tissu à l'intérieur du corps et irradié par le faisceau, et de la concentration en bore à l'intérieur du corps.
     
    3. Procédé selon la revendication 2, dans lequel dans l'angle de rayonnement échantillonné et la piste de rayonnement, un facteur de pondération (W(i)) d'un organe irradié (i) est calculé en utilisant la formule I :

    dans laquelle, I(i), S(i) et C(i) représentent l'intensité de faisceau, le coefficient sensible aux rayonnements de l'organe (i) et la concentration en bore de l'organe (i), respectivement.
     
    4. Procédé selon la revendication 3, dans lequel le I(i) est calculée en utilisant la formule Il qui intègre une intensité de profondeur ou une courbe de dosage du corps conformément au faisceau utilisé :

    dans laquelle i(x) représente la fonction d'intensité de profondeur ou la fonction de courbe de dosage du faisceau dans d'un corps approximatif et x0 - x représente la plage de profondeur de l'organe (i) dans la piste de rayonnement.
     
    5. Procédé selon la revendication 3 ou la revendication 4, dans lequel le coefficient d'évaluation est calculé en utilisant la Formule III :

    dans laquelle, Q(x, y, z, Φ, θ) en tant que coefficient d'évaluation représente la somme des facteurs de pondération de chaque organe dans la piste de rayonnement.
     
    6. Procédé selon la revendication 5, dans lequel le rapport (QR(x, z, Φ, θ)) du coefficient d'évaluation à un coefficient d'évaluation de tumeur de la tumeur est calculé en utilisant la formule IV :


    dans laquelle, W(tumeur) est un facteur de pondération de la tumeur.
     
    7. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel le procédé comprend en outre une étape de lecture de données d'image médicale choisie dans le groupe constitué de données d'image de tomodensitométrie, de données d'image d'I.R.M. ou de données d'image TEP-TDM.
     
    8. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel le procédé comprend en outre une étape d'affichage de chaque coefficient d'évaluation dans une image tridimensionnelle.
     
    9. Procédé selon la revendication 7, dans lequel le procédé comprend en outre une étape de définition ou de lecture du profil, d'un organe, d'un tissu ou d'une tumeur dans le corps ; et une étape de définition du matériau et de la densité de l'organe, du tissu et de la tumeur.
     




    Drawing












    REFERENCES CITED IN THE DESCRIPTION



    This list of references cited by the applicant is for the reader's convenience only. It does not form part of the European patent document. Even though great care has been taken in compiling the references, errors or omissions cannot be excluded and the EPO disclaims all liability in this regard.

    Patent documents cited in the description