Proton dose calculation method based on tissue segmentation of delineated regions and related equipment

Abstract

The present application relates to a proton dose calculation method for tissue segmentation based on delineated areas and related equipment, the method comprising: constructing a basic database; the basic database comprises a first database and a second database; obtaining CT information and radiotherapy delineation information; determining the density of each voxel and the tissue organ to which it belongs in combination with the basic database according to the CT information and the radiotherapy delineation information; performing proton dose calculation according to the radiotherapy plan through a preset algorithm based on the density and the tissue organ to which each voxel belongs; obtaining a dose distribution map of the area of ​​interest; comparing the target dose with the calculated dose, and performing plan evaluation and adjustment guidance based on the comparison deviation; being able to comprehensively consider tissue distribution at multiple levels, more accurately reflect the diversity and complexity of the chemical composition of different tissues, accurately simulate the interaction between protons and various tissues, improve the accuracy of dose calculation, and ensure the accuracy of the treatment plan.

Description

Technical Field

[0001] The present application relates to the field of particle radiotherapy technology, and in particular to a proton dose calculation method based on tissue segmentation of delineated areas and related equipment. Background Art

[0002] In the field of particle radiotherapy, the accuracy of dose calculation is one of the decisive factors for achieving high-quality treatment results. For proton radiotherapy, proton dose calculation methods involve two main approaches: analytical method and Monte Carlo simulation. The analytical method often converts the Hounsfield Units (HU) value obtained from CT scans into a stopping power ratio for dose calculation. Although this does not involve the specific tissue composition of the voxel, this analytical algorithm is difficult to provide accurate dose calculation when dealing with heterogeneous media (such as the lungs). Because these algorithms can only perform approximate calculations in terms of multi-elastic Coulomb scattering, elastic and inelastic nuclear reactions, and cannot achieve accurate modeling.

[0003] The Monte Carlo simulation method is hailed as the gold standard for particle dose calculation. When performing Monte Carlo simulation of proton dose calculation, the HU value is generally mapped to the density information of the voxel, and the voxels with the same HU value or within a specific HU range are classified as the same organ to determine the chemical composition of the voxel. However, this traditional method faces an important problem. When converting the HU value to density, the same HU value (density value) may correspond to different tissues and organs. Because different tissues and organs differ in elemental composition, even if the density is the same, they interact with protons differently. This results in the inability to accurately consider the differences in physical properties between tissues when performing dose calculations at the same density, which in turn leads to inaccurate dose distribution.

[0004] Specifically, although two or more regions show the same density (HU value) on CT images, they may actually be composed of different types of tissue. This difference leads to different proton interaction modes and dose absorption characteristics. This error can have an adverse effect on the formulation and implementation of treatment plans. In particle radiotherapy, it is usually hoped to accurately calculate the interaction of particle beams in the patient's body in order to accurately deliver the required dose. However, the complexity of density value conversion leads to uncertainty in dose calculation.

[0005] In order to solve the problem of this source of error, this paper proposes a proton dose calculation method based on tissue segmentation of the outlined area. This method builds a detailed human tissue and organ database based on multiple factors such as gender, human anatomical position and tissue distribution. By integrating CT information and tissue and organ outline data, the density of each CT voxel and the corresponding tissue and organ composition can be accurately determined, thereby improving the dose calculation to a higher level of accuracy and reliability. Summary of the invention

[0006] The purpose of this application is to provide a proton dose calculation method and related equipment based on tissue segmentation of delineated areas. The method is classified based on the location of human anatomical structures and tissue distribution, so it can comprehensively consider tissue distribution at multiple levels. The differences in anatomical structures in different regions are fully considered, which can more accurately reflect the diversity of chemical compositions of different tissues. The complexity of tissues is fully considered, the interaction between protons and various tissues is accurately simulated, the accuracy of dose calculation is improved, and the accuracy of treatment plans is ensured. This method allows personalized treatment planning and customizes dose distribution according to the actual anatomical structure and tissue distribution of each patient. This helps to minimize the dose of healthy tissue while ensuring that the tumor area receives sufficient treatment dose. Dose calculation based on more accurate tissue and organ information can better predict the effect of radiotherapy. By more accurately simulating dose distribution, treatment plans can be optimized, treatment effects can be improved, and treatment risks can be reduced.

[0007] The purpose of this application is achieved by the following technical solutions:

[0008] In a first aspect, the present application provides a proton dose calculation method for tissue segmentation based on delineated regions, characterized in that the method comprises:

[0009] S1. Constructing a basic database of human tissues and organs; the basic database includes a first database and a second database, both of which include names and densities of human tissues and organs, the density value in the first database is not unique, while the density value in the second database is unique, the first database has radiotherapy delineation information, while the second database does not have radiotherapy delineation information;

[0010] S2, obtaining CT information and radiotherapy delineation information; determining the density of each voxel and the tissue organ to which it belongs based on the CT information and the radiotherapy delineation information and in combination with the corresponding first database and second database in the basic database;

[0011] S3. Obtain the proton absorption dose of voxel therapy according to the density of each voxel and the tissue and organ to which it belongs through a preset algorithm;

[0012] S4. Obtain the proton absorption dose distribution of all voxels for treatment and generate a dose distribution map; compare the actual dose with the target dose, and perform plan evaluation and adjustment guidance based on the comparison deviation.

[0013] Preferably, the proton dose calculation method based on tissue segmentation of the delineated area, said S1 comprises:

[0014] Building a basic database based on reports and documents; the basic database includes a first database and a second database;

[0015] The first database includes parameters of human tissues and organs;

[0016] The second database forms multiple sub-databases according to the division criteria; the division criteria include a first criterion and a second criterion; the first criterion is gender; within each first criterion, further division is performed according to the second criterion; the second criterion is the location and tissue distribution of the human anatomical structure, and multiple regions are divided according to the location and tissue distribution of the human anatomical structure; the second database includes materials of different densities in the region.

[0017] Preferably, in the proton dose calculation method based on tissue segmentation of the delineated area, the basic database also includes α and β parameters obtained through calculation;

[0018] The α parameter is used to scale the proton track, and the β parameter is used to scale the deflection angle of the track;

[0019] Among them, the α parameter is the ratio of the stopping power of water to the stopping power of tissues and organs. Each energy corresponds to an α parameter, and the energy range is determined according to the energy range of protons:

[0020]

[0021] in, The proton energy is E i When the α parameter is The proton energy is E i The stopping power in water, The proton energy is E i The ability to stop in the organization;

[0022] The β parameter is obtained by calculating the ratio of the Coulomb scattering angle of tissues and organs to the Coulomb scattering angle of water at different energies and then averaging them. One tissue and organ corresponds to only one β parameter:

[0023]

[0024] Where n is the number of energies; When the proton energy is E i Multiple Coulomb scattering angles of tissues and organs; When the proton energy is E i The multiple Coulomb scattering angle of water.

[0025] Preferably, the proton dose calculation method based on tissue segmentation of the delineated area, said S2 comprises:

[0026] Obtain CT information and radiotherapy delineation information; determine the voxel range covered by each delineated tissue and organ through the radiotherapy delineation information;

[0027] The delineated voxels are assigned as the types of corresponding tissues and organs through the first database;

[0028] The region to which the voxels not outlined in the radiotherapy outline information belong is determined according to the position of the planned target area; and the tissue organ to which it belongs is determined according to the voxel density in combination with the second database.

[0029] Preferably, the proton dose calculation method based on tissue segmentation of the delineated area, said S3 comprises:

[0030] Calculate voxel parameters according to the density of each voxel and the tissue and organ it belongs to;

[0031] The proton absorbed dose is calculated according to the radiotherapy plan through voxel parameters and preset algorithms; the normalized absorbed proton dose of each voxel can be expressed as:

[0032]

[0033] Among them, D i is the absorbed dose of voxel i, E ki is the energy deposited by the kth particle in voxel i, m i is the mass of voxel i.

[0034] Preferably, the proton dose calculation method based on tissue segmentation of the delineated area, said S4 comprises:

[0035] Obtain the proton dose of all voxels and generate the dose distribution map of the region of interest, and calculate the dose index;

[0036] The target dose index is compared with the calculated dose index, and the plan is evaluated and adjusted based on the comparison deviation.

[0037] In a second aspect, the present application proposes a proton dose calculation system for tissue segmentation based on delineated regions, the system comprising:

[0038] A database construction module, used to construct a basic database of human tissues and organs; the basic database includes a first database and a second database; the first database and the second database both include names and densities of human tissues and organs, the density value in the first database is not unique, the density value in the second database is unique, the first database has radiotherapy delineation information, and the second database does not have radiotherapy delineation information;

[0039] A first determination module is used to obtain CT information and radiotherapy delineation information; according to the CT information and radiotherapy delineation information, combined with the corresponding first database and second database in the basic database, determine the density of each voxel and the tissue organ to which it belongs;

[0040] The dose acquisition module is used to obtain the proton absorption dose of the voxel therapy according to the density of each voxel and the tissue organ to which it belongs through a preset algorithm;

[0041] The adjustment module is used to obtain the proton absorption dose distribution of all voxels and generate a dose distribution map; compare the actual dose with the target dose, and provide dose adjustment guidance based on the comparison deviation.

[0042] In a third aspect, the present application provides an electronic device, comprising a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of any of the above methods when executing the computer program.

[0043] In some optional embodiments, the electronic device is also provided with a display screen.

[0044] In a fourth aspect, the present application provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of any of the above methods are implemented.

[0045] The benefits and effects of this application are: using the constructed basic database, combined with CT information and radiotherapy delineation information, the density of each voxel and the tissue organ to which it belongs are determined, and the proton dose of the voxel treatment is obtained through a preset algorithm. This method takes into account the individual differences and tissue characteristics of the human body. The method is classified based on the location of the human anatomical structure and the distribution of tissues, so it can comprehensively consider the distribution of tissues at multiple levels. The differences in the anatomical structures of different regions are fully considered, which can more accurately reflect the chemical composition diversity of different tissues; the complexity of the tissue is fully considered, and the interaction between protons and various tissues is accurately simulated, so that the proton dose distribution of each patient can be more accurately calculated, thereby improving the accuracy of the dose calculation and ensuring the accuracy of the treatment plan; by obtaining the proton dose distribution of all voxel treatments and generating a dose distribution map, the doctor can intuitively understand the dose distribution of the treatment area. Comparing the actual dose with the simulated dose and adjusting the dose according to the comparison results can ensure that the dose distribution of radiotherapy is consistent with the plan, thereby guiding precise radiotherapy, improving the treatment effect and reducing damage to normal tissues. Since this method takes into account the individual differences and tissue characteristics of the human body, a personalized treatment plan can be tailored for each patient. This method allows personalized treatment planning, customizing the dose distribution according to the actual anatomical structure and tissue distribution of each patient. This helps to minimize the dose to healthy tissue while ensuring that the tumor area receives sufficient treatment dose; by combining the patient's CT information, radiotherapy delineation information, and information in the basic database, the proton dose distribution of each patient can be determined more accurately, thereby guiding doctors to provide patients with more personalized treatment plans. By improving the accuracy of dose calculations and guiding precise radiotherapy, this method can improve treatment outcomes. Doctors can adjust the dose distribution according to the actual situation of each patient to ensure that the treatment effect is maximized. At the same time, reducing damage to normal tissue can reduce the incidence of complications and improve the quality of life of patients. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] The present application is further described below in conjunction with the accompanying drawings and embodiments.

[0047] Figure 1 It is a flowchart of a proton dose calculation method for tissue segmentation based on delineated regions provided in an embodiment of the present application;

[0048] Figure 2 It is a schematic diagram of a basic database structure provided in an embodiment of the present application;

[0049] Figure 3 It is a schematic diagram of a proton dose calculation system for tissue segmentation based on delineated regions provided in an embodiment of the present application;

[0050] Figure 4 It is a structural schematic diagram of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0051] Below, the present application is further described in conjunction with the accompanying drawings and specific implementation methods. It should be noted that, under the premise of no conflict, the various embodiments or technical features described below can be arbitrarily combined to form a new embodiment.

[0052] See also Figure 1 , an embodiment of the present application provides a proton dose calculation method for tissue segmentation based on delineated regions, the method comprising:

[0053] S1. Constructing a basic database of human tissues and organs; the basic database includes a first database and a second database; the first database and the second database both include names and densities of human tissues and organs, the density value in the first database is not unique, the density value in the second database is unique, the first database has radiotherapy delineation information, and the second database does not have radiotherapy delineation information;

[0054] S2, obtaining CT information and radiotherapy delineation information; determining the density of each voxel and the tissue organ to which it belongs based on the CT information and the radiotherapy delineation information and in combination with the corresponding first database and second database in the basic database;

[0055] S3. Obtain the proton absorption dose of voxel therapy according to the density of each voxel and the tissue and organ to which it belongs through a preset algorithm;

[0056] S4. Obtain the proton absorption dose distribution of all voxels and generate a dose distribution map; compare the actual dose with the target dose, and conduct plan evaluation and adjustment guidance based on the comparison deviation;

[0057] The outlining information is the outlining data, which includes outlining different tissues and organs of the human body on the CT image, for example, area A is the lung, area B is the heart, and area C is the breast.

[0058] The working principle of the above technical solution is: construct a basic database, which includes a first database and a second database. This information is used to determine the density of each voxel and the tissue and organ to which it belongs. Obtain CT information and radiotherapy delineation information. CT information provides the patient's body structure and density information, while radiotherapy delineation information provides doctors with delineations of tumors and other key structures. Combined with the basic database, the density and tissue and organ to which each voxel belongs can be determined. According to the density and tissue and organ to which each voxel belongs, the proton dose of voxel treatment is obtained through a preset algorithm. This step is based on the density and tissue and organ information of the voxel, and the proton dose to which each voxel is subjected is calculated through a preset algorithm. Obtain the proton dose distribution of all voxel treatments and generate a dose distribution map. This step integrates the proton doses of all voxels to generate a dose distribution map, which can intuitively display the dose distribution of the treatment area.

[0059] This method can help doctors more accurately calculate the proton dose distribution to tumors and other critical structures, thereby better guiding the radiation therapy process, improving treatment effectiveness and reducing damage to normal tissues.

[0060] The effect of the above technical solution is: the method uses the constructed basic database, combined with CT information and radiotherapy delineation information, to determine the density and tissue organ of each voxel, and obtains the absorbed dose of the patient's area of ​​interest according to the radiotherapy plan through a preset algorithm. This method takes into account the individual differences and tissue characteristics of the human body. The method is classified based on the location of the human anatomical structure and the distribution of tissues, so it can comprehensively consider the distribution of tissues at multiple levels. The differences in the anatomical structures of different regions are fully considered, which can more accurately reflect the chemical composition diversity of different tissues; the complexity of the tissue is fully considered, and the interaction between protons and various tissues is accurately simulated, so that the proton dose distribution of each patient can be more accurately calculated, thereby improving the accuracy of dose calculation and ensuring the accuracy of the treatment plan; by obtaining the proton dose distribution of all voxels and generating a dose distribution map, the doctor can intuitively understand the dose distribution of the treatment area. Comparing the target dose with the calculated dose and adjusting the dose according to the comparison results can ensure that the dose distribution of radiotherapy is consistent with the plan, thereby guiding precise radiotherapy, improving the treatment effect and reducing damage to normal tissues. Since this method takes into account the individual differences and tissue characteristics of the human body, a personalized treatment plan can be tailored for each patient. This method allows personalized treatment planning, customizing the dose distribution according to the actual anatomical structure and tissue distribution of each patient. This helps to minimize the dose to healthy tissue while ensuring that the tumor area receives sufficient treatment dose; by combining the patient's CT information, radiotherapy delineation information, and information in the basic database, the proton dose distribution of each patient can be determined more accurately, thereby guiding doctors to provide patients with more personalized treatment plans. By improving the accuracy of dose calculations and guiding precise radiotherapy, this method can improve treatment outcomes. At the same time, reducing damage to normal tissue can reduce the incidence of complications and improve the quality of life of patients. This method can be applied to different types of proton therapy equipment and techniques, including fixed beams, scanning beams, and pencil beams.

[0061] The proton dose calculation method for tissue segmentation based on delineated regions in this embodiment is characterized in that S1 comprises:

[0062] A basic database is constructed based on human tissues and organs in the International Commission on Radiological Protection (ICRP) report and related literature; the basic database includes a first database and a second database.

[0063] The first database includes human tissue and organ parameters that have been collected based on relevant reports and literature; for example, the density is 1.04 g / cm 3 The tissues and organs include the prostate, pancreas, gray matter and white matter of the brain, etc. The density is 1.05g / cm 3The tissues and organs include head, muscle, heart, liver, stomach, etc.; the density value (density size) of the first database is not unique, that is, the same density value may have different tissues and organs.

[0064] The second database contains multiple sub-databases; the division criteria of the sub-databases include a first criterion and a second criterion; the first criterion is gender; within each first criterion, further division is performed according to the second criterion; the second criterion is the location and tissue distribution of the human anatomical structure, and multiple regions are divided according to the location and tissue distribution of the human anatomical structure; the second database includes materials of different densities in the region; each sub-database includes materials of different densities belonging to the region, and if there are tissues and organs of the same density in the region, the most representative tissue is selected based on factors such as the volume occupied by the tissues and organs, the importance of affecting the treatment effect in radiotherapy, and the similarity of tissue components. For example, for the head, both gray matter and white matter belong to this region, and they have the same density and similar elemental composition, but the volume of white matter is about 1.1 to 1.5 times that of gray matter. Therefore, for the head and neck region, the density is 1.05 g / cm 3 Only the white matter of the brain is retained; the density value of the second database is unique, that is, the same density value will have the same tissue organ.

[0065] The working principle of the above technical solution is as follows: First, a basic database can be constructed by collecting and collating relevant reports and documents. This database includes two main parts, the first database and the second database. The first database contains various human tissues and organs, which are only identified according to their types; the second database contains multiple sub-databases, and the division criteria of these sub-databases include gender (first criterion) and the location and tissue distribution of human anatomical structures (second criterion). According to the first criterion of gender, the database can be divided into two sub-databases for males and females. Then, within each gender, further division is performed according to the second criterion of the location and tissue distribution of human anatomical structures. This ensures that each sub-database contains data of a specific gender, a specific part, and a specific tissue distribution. In the second database, each sub-database contains material information of different densities in the area. This information can help calculate the proton dose of different parts and tissues and organs. After determining the density of each voxel and the tissue and organ to which it belongs, the proton dose of each voxel can be calculated by a preset algorithm. This algorithm takes into account the density of the voxel, the characteristics of the tissue and organ, and other relevant factors in the treatment plan. Finally, the proton dose of the region of interest is obtained and a dose distribution map is generated. By comparing the deviation between the calculated dose and the target dose, the plan can be evaluated and adjusted to ensure that the dose distribution of radiotherapy is consistent with the plan. Through such a workflow, the proton dose calculation method based on tissue segmentation of the outlined area provided in this embodiment can more accurately calculate the proton dose distribution of each patient, thereby improving the treatment effect and reducing damage to normal tissues.

[0066] The effect of the above technical solution is: this method uses all human tissue and organ parameters contained in the basic database and the material information of different densities in the region, combined with CT information and radiotherapy delineation information, to more accurately determine the density of each voxel and the tissue and organ to which it belongs, and then obtain the proton dose of each voxel through a preset algorithm. This method takes into account the individual differences and tissue characteristics of the human body, and can more accurately calculate the proton dose distribution of each patient, thereby improving the accuracy and precision of dose calculation. By obtaining the dose of the region of interest and generating a dose distribution map, the doctor can intuitively understand the dose distribution of the treatment area. Comparing the calculated dose with the target dose, and evaluating and adjusting the plan according to the comparison results, it can ensure that the dose distribution of radiotherapy is consistent with the plan, thereby guiding precise radiotherapy. This method helps to improve the treatment effect, reduce the incidence of complications, and improve the quality of life of patients. Since this method takes into account the individual differences and tissue characteristics of the human body, a personalized treatment plan can be tailored for each patient. By combining the patient's CT information, radiotherapy delineation information and information in the basic database, the proton dose distribution of each patient can be more accurately determined, thereby guiding doctors to provide patients with more personalized treatment plans. This method can be applied to different types of proton therapy equipment and technologies, including fixed beam, scanning beam and pencil beam. In addition, this method can also provide accurate guidance for different types of tumor treatment, such as head and neck cancer, lung cancer, prostate cancer, etc. Therefore, this method has broad application prospects and market potential; by accurately calculating the proton dose distribution of each patient, unnecessary irradiation can be reduced, thereby improving treatment efficiency. This helps to shorten the patient's treatment cycle and reduce the consumption of medical resources. This method can promote the development of radiotherapy technology by improving the dose calculation method and establishing a more complete basic database. This method provides doctors with more accurate and personalized dose calculation tools, which helps to improve the overall level of radiotherapy.

[0067] In summary, the proton dose calculation method based on tissue segmentation of the outlined area in this embodiment can improve the accuracy and precision of dose calculation, guide precise radiotherapy, and provide personalized treatment plans. It has the benefits and effects of wide applicability, improving treatment efficiency, and promoting the development of radiotherapy technology.

[0068] This embodiment is based on the proton dose calculation method for tissue segmentation of the delineated area, and the basic database includes the name, density and calculated α and β parameters of human tissues and organs.

[0069] The α parameter is used to scale the proton's track, and the β parameter is used to scale the deflection angle of the track.

[0070] The α parameter is the ratio of the stopping power of water to the stopping power of tissues and organs. Each energy corresponds to an α parameter, and the energy range is determined by the energy range of protons:

[0071]

[0072] in, The proton energy is E i When the α parameter is The proton energy is E i The stopping power in water,

[0073] The proton energy is E i The ability to stop in an organization; and Obtained through FLUKA Monte Carlo simulation;

[0074] The β parameter is obtained by calculating the ratio of the Coulomb scattering angle of tissues and organs to the Coulomb scattering angle of water at different energies and then averaging them. One tissue and organ corresponds to only one β parameter:

[0075]

[0076] Where n is the number of energies; and When the proton energy is E i Multiple Coulomb scattering angles of tissues and water.

[0077] The working principle of the above technical solution is as follows: First, by collecting and collating relevant reports and documents, a basic database can be constructed. This database includes the names, densities, and calculated α and β parameters of human tissues and organs. These parameters are very important for accurately calculating radiation doses. The density values ​​of the first database are not unique, while the density values ​​of the second database are unique.

[0078] In the basic database, the α parameter is used to scale the proton track, while the β parameter is used to scale the deflection angle of the track.

[0079]

[0080] in, The proton energy is E i When the α parameter is The proton energy is E i The stopping power in water,

[0081] The proton energy is E i The ability to stop in an organization; and Obtained through FLUKA Monte Carlo simulation.

[0082] The β parameter is obtained by calculating the ratio of the Coulomb scattering angle of tissues and organs to the Coulomb scattering angle of water at different energies and then averaging them. One tissue and organ corresponds to only one β parameter:

[0083]

[0084] Where n is the number of energies; and When the proton energy is E i Multiple Coulomb scattering angles of tissues and water.

[0085] Among them, the proton energy range is 1-230MeV, and the energy interval is 1MeV, that is, one material corresponds to 230 α values ​​and one β value.

[0086] After determining the density and tissue of each voxel, the proton absorbed dose of each voxel can be calculated using the α and β parameters in the basic database through a preset algorithm. This algorithm takes into account the density of the voxel, the characteristics of the tissue and organ, and other relevant factors in the treatment plan.

[0087] Finally, the dose of the region of interest is obtained and a dose distribution map is generated. By comparing the deviation between the calculated dose and the target dose, the plan can be evaluated and adjusted to ensure that the dose distribution of radiotherapy is consistent with the plan.

[0088] The effect of the above technical solution is: by constructing a basic database, which contains the parameters of human tissues and organs, accurate tissue segmentation results can be provided. Each tissue and organ has a specific density and α, β parameters, which can be used to scale the length and deflection angle of the proton track, thereby achieving accurate proton dose calculation. FLUKA Monte Carlo simulation is used to obtain the stopping power and Coulomb scattering angle of protons in water and tissues. FLUKA is a simulation tool widely used in particle transport and radiation physics research, which can provide highly reliable results. This method can be applied to proton therapy planning and dose assessment in the medical field. By accurately calculating the dose distribution of protons in different tissues, doctors can be better guided in treatment planning and minimize side effects for patients. In short, this proton dose calculation method based on tissue segmentation based on delineation area can provide accurate and reliable proton dose calculation results by constructing a basic database, considering the influence of different regions and materials, and using FLUKA Monte Carlo simulation. It is widely used in the medical field and helps to improve the effect and safety of proton therapy.

[0089] In this embodiment, the proton dose calculation method for tissue segmentation based on the delineated area, S2 includes:

[0090] Obtain CT information and radiotherapy delineation information; determine the voxel range covered by each delineated tissue and organ through the radiotherapy delineation information;

[0091] The delineated voxels are assigned as the types of corresponding tissues and organs through the first database;

[0092] The region to which the voxels not outlined in the radiotherapy outline information belong is determined according to the position of the planned target area; and the tissue organ to which it belongs is determined according to the voxel density in combination with the second database.

[0093] The working principle of the above technical solution is as follows: Radiotherapy technology is widely used in the treatment of diseases. In the practice of radiotherapy, accurately determining the material properties of the patient's CT geometry is a key step to ensure the accuracy of the treatment plan. Only by deeply understanding the physical properties of the tissue, such as density and atomic composition, can we more accurately calculate the distribution of radiation doses and thus formulate a more accurate treatment plan.

[0094] First, the patient's CT scan image and radiotherapy delineation information are obtained. The CT image provides the anatomical structure and density information of the tissues and organs, while the radiotherapy delineation information determines the specific area of ​​interest in the radiotherapy treatment plan. According to the radiotherapy delineation information, the voxels of the specific organs are delineated on the CT image. These voxels represent the position and shape of the organs in three-dimensional space. The delineated voxels are matched to the type of the corresponding tissues and organs using the human tissue and organ parameters contained in the first database. The first database includes the name, density, and calculated α and β parameters of human tissues and organs. According to the location of the planned target area, the region to which the voxels not delineated in the radiotherapy delineation belong is determined. Usually, the tissues and organs around the planned target area may not be delineated, and it is necessary to determine the region to which they belong. According to the density of the undelined voxels and the regional and density material information contained in the second database, the specific tissue and organ to which they belong is determined. The second database divides multiple regions according to the location of the human anatomical structure and tissue distribution, and includes materials of different densities in the region. By comparing the density and regional information of the voxels, the tissue and organ to which they belong can be determined.

[0095] Through the above steps, the proton dose calculation method based on tissue segmentation of the outlined area can more accurately determine the tissue and organ type of each voxel in the CT image. This provides accurate basic data for subsequent proton dose calculation to guide the design of radiotherapy treatment plans and dose assessment.

[0096] The effect of the above technical solution is: by acquiring CT information and radiotherapy delineation information, the voxels of specific organs are delineated, and the delineated voxels are designated as the type of corresponding tissue organs according to the first database. At the same time, the area to which the voxels not delineated in the radiotherapy delineation belong is determined according to the position of the planned target area, and then the tissue organs to which they belong are determined according to the voxel density combined with the second database. By combining the delineation information and the area of ​​the tumor target area, the material component information of the voxels can be determined more accurately, so that the proton dose can be calculated more accurately to achieve better treatment effects. This information can guide the design of radiotherapy treatment plans, making the treatment plans more personalized, refined, and effective. By accurately calculating the proton dose, excessive or insufficient proton radiotherapy can be avoided, thereby improving the safety of treatment.

[0097] In this embodiment, the proton dose calculation method for tissue segmentation based on the delineated area, S3 includes:

[0098] According to the density of each voxel and the tissue organ to which it belongs, the voxel parameters are calculated, and the voxel parameters include α i and β i ;

[0099]

[0100] Among them, α r and β r are the reference α and β values ​​corresponding to the tissues in the database, ρ r and ρ i are the density of the tissue in the database and the density of the current voxel, α i and β i They are the α and β parameters of the current voxel obtained after scaling the parameters of the database according to the density; it can be seen that α scales inversely with the density, while β scales directly with the density.

[0101] The proton dose is calculated according to the radiotherapy plan through the voxel parameters and the preset algorithm. The normalized absorbed dose of each voxel can be expressed as;

[0102]

[0103] Among them, D i is the absorbed dose of voxel i, E ki is the energy deposited by the kth particle in voxel i, m i is the mass of voxel i.

[0104] The working principle of the above technical solution is: first, the voxel parameter α is calculated according to the density of each voxel and the tissue organ to which it belongs. i and β i Among them, α irepresents the alpha parameter of the current voxel, β i Represents the beta parameter of the current voxel. These parameters are calculated based on the reference alpha and beta values ​​and density information corresponding to the tissue in the database. The voxel parameters are used to calculate the proton dose through the FDC (Fastdose calculation) algorithm. Based on the calculated voxel dose, the proton dose distribution of the region of interest is compared with the target dose distribution, and the proton therapy plan can be evaluated and adjusted to ensure adequate dose coverage within the tumor tissue while minimizing radiation damage to surrounding normal tissues.

[0105] In general, this method achieves accurate calculation of proton dose based on tissue segmentation by calculating voxel parameters and applying preset algorithms, thereby better guiding the design of radiotherapy treatment plans and improving the accuracy and safety of treatment.

[0106] The effect of the above technical solution is: by considering the density of each voxel and the tissue organ to which it belongs, the voxel parameter α is calculated. i and β i , making the proton dose calculation more accurate. Since proton therapy is very sensitive to changes in tissue density, a more refined division of tissue components can more realistically reflect the absorption of protons by different tissue types, thereby improving the accuracy of dose calculation. This method performs dose calculation based on the patient's actual tissue segmentation, and can provide patients with more personalized treatment plans. By considering the density of different tissue types and their absorption of protons, the treatment plan can be better evaluated and adjusted to minimize radiation damage to surrounding normal tissues while ensuring that tumor tissues receive sufficient doses. Accurate proton dose calculation can help doctors better evaluate the safety of treatment plans and avoid excessive or insufficient treatment. This helps to maximize the safety and effectiveness of treatment and reduce complications and adverse reactions during treatment. Overall, the proton dose calculation method based on tissue segmentation of the outlined area can improve the personalization and accuracy of proton therapy, thereby guiding doctors to provide patients with safer and more effective treatment plans.

[0107] In this embodiment, the proton dose calculation method for tissue segmentation based on the delineated area, the S4 includes:

[0108] Obtain doses for all voxels and generate dose-volume histograms for planning targets and organs at risk; calculate dose indices

[0109] The target dose index is compared with the calculated dose index, and the plan is evaluated and adjusted according to the comparison deviation.

[0110] The working principle of the above technical solution is: according to the treatment plan and simulation parameters, the dose algorithm is used to perform simulation calculations to obtain the dose of each voxel during the treatment process. The obtained dose value is mapped to the corresponding voxel position to generate a three-dimensional dose distribution data set. This data set can be used to visualize the dose distribution and help doctors analyze and evaluate. This can optimize the treatment plan and improve the accuracy and effectiveness of proton therapy. This helps to improve the accuracy and personalization of proton therapy and provide patients with safer and more effective treatment options.

[0111] The effect of the above technical solution is: by obtaining and analyzing the proton dose distribution or dose assessment index of the region of interest, the distribution of the dose in the patient's body can be understood more accurately. This helps doctors better evaluate the treatment effect and make adjustments as needed to ensure maximum destruction of tumor cells while minimizing damage to healthy tissues. This method uses CT information to map the dose value to the corresponding voxel position to generate a three-dimensional dose distribution data set. This provides doctors with detailed dose distribution information, enabling them to develop personalized treatment plans and adjust the dose according to the characteristics and requirements of different tissues in specific cases.

[0112] In summary, the proton dose calculation method based on tissue segmentation of the outlined area can improve the accuracy and personalization of proton therapy, optimize the dose distribution, and provide patients with safer and more effective treatment plans, while improving the treatment effect and minimizing side effects.

[0113] See also Figure 3 The present application also provides a proton dose calculation system for tissue segmentation based on delineated regions, the system comprising:

[0114] A database construction module, used to construct a basic database of human tissues and organs; the basic database includes a first database and a second database; the first database and the second database both include names and densities of human tissues and organs, the density value in the first database is not unique, the density value in the second database is unique, the first database has radiotherapy delineation information, and the second database does not have radiotherapy delineation information;

[0115] A first determination module is used to obtain CT information and radiotherapy delineation information; according to the CT information and radiotherapy delineation information, combined with the corresponding first database and second database in the basic database, determine the density of each voxel and the tissue organ to which it belongs;

[0116] The dose acquisition module is used to obtain the absorbed dose of each voxel according to the density of each voxel and the tissue and organ to which it belongs, through a preset algorithm and according to the treatment plan;

[0117] The adjustment module is used to obtain the proton dose distribution of all voxels for treatment and generate dose distribution maps; compare the actual dose with the target dose, and perform plan evaluation and adjustment guidance based on the comparison deviation.

[0118] The specific execution steps and contents of the data construction module include:

[0119] Building a basic database based on reports and documents; the basic database includes a first database and a second database;

[0120] The first database includes parameters of human tissues and organs;

[0121] The second database forms multiple sub-databases according to the division criteria; the division criteria include a first criterion and a second criterion; the first criterion is gender; within each first criterion, further division is performed according to the second criterion; the second criterion is the location and tissue distribution of the human anatomical structure, and multiple regions are divided according to the location and tissue distribution of the human anatomical structure; the second database includes materials of different densities in the region.

[0122] The basic database also includes the calculated α and β parameters;

[0123] The α parameter is used to scale the proton track, and the β parameter is used to scale the deflection angle of the track;

[0124] Among them, the α parameter is the ratio of the stopping power of water to the stopping power of tissues and organs. Each energy corresponds to an α parameter, and the energy range is determined according to the energy range of protons:

[0125]

[0126] in, The proton energy is E i When the α parameter is The proton energy is E i The stopping power in water, The proton energy is E i The ability to stop in the organization;

[0127] The β parameter is obtained by calculating the ratio of the Coulomb scattering angle of tissues and organs to the Coulomb scattering angle of water at different energies and then averaging them. One tissue and organ corresponds to only one β parameter:

[0128]

[0129] Where n is the number of energies; When the proton energy is E i Multiple Coulomb scattering angles of tissues and organs; When the proton energy is Ei The multiple Coulomb scattering angle of water.

[0130] The specific execution steps and contents of the first determination module include:

[0131] Obtain CT information and radiotherapy delineation information; determine the voxel range covered by each delineated tissue and organ through the radiotherapy delineation information;

[0132] The delineated voxels are assigned as the types of corresponding tissues and organs through the first database;

[0133] The region to which the voxels not outlined in the radiotherapy outline information belong is determined according to the position of the planned target area; and the tissue organ to which it belongs is determined according to the voxel density in combination with the second database.

[0134] The specific execution steps and contents of the dose acquisition module include:

[0135] Calculate voxel parameters according to the density of each voxel and the tissue and organ it belongs to;

[0136] The proton absorbed dose is calculated according to the radiotherapy plan through voxel parameters and preset algorithms; the normalized absorbed proton dose of each voxel can be expressed as:

[0137]

[0138] Among them, D i is the absorbed dose of voxel i, E ki is the energy deposited by the kth particle in voxel i, m i is the mass of voxel i.

[0139] The specific execution steps and contents of the adjustment module include:

[0140] Obtain the proton dose of all voxels and generate the dose distribution map of the region of interest, and calculate the dose index;

[0141] The target dose index is compared with the calculated dose index, and the plan is evaluated and adjusted based on the comparison deviation.

[0142] The working principle of the above technical solution is: construct a basic database, which includes a first database and a second database. This information is used to determine the density of each voxel and the tissue and organ to which it belongs. Obtain CT information and radiotherapy delineation information. CT information provides the patient's body structure and density information, while radiotherapy delineation information provides doctors with delineations of tumors and other key structures. Combined with the basic database, the density and tissue and organ to which each voxel belongs can be determined. According to the density and tissue and organ to which each voxel belongs, the proton dose of voxel treatment is obtained through a preset algorithm. This step is based on the density and tissue and organ information of the voxel, and the proton dose to which each voxel is subjected is calculated through a preset algorithm. Obtain the proton dose distribution of all voxel treatments and generate a dose distribution map. This step integrates the proton doses of all voxels to generate a dose distribution map, which can intuitively display the dose distribution of the treatment area.

[0143] This method can help doctors more accurately calculate the proton dose distribution to tumors and other critical structures, thereby better guiding the radiation therapy process, improving treatment effectiveness and reducing damage to normal tissues.

[0144] The effect of the above technical solution is: the method uses the constructed basic database, combined with CT information and radiotherapy delineation information, to determine the density and tissue organ of each voxel, and obtains the absorbed dose of the patient's area of ​​interest according to the radiotherapy plan through a preset algorithm. This method takes into account the individual differences and tissue characteristics of the human body. The method is classified based on the location of the human anatomical structure and the distribution of tissues, so it can comprehensively consider the distribution of tissues at multiple levels. The differences in the anatomical structures of different regions are fully considered, which can more accurately reflect the chemical composition diversity of different tissues; the complexity of the tissue is fully considered, and the interaction between protons and various tissues is accurately simulated, so that the proton dose distribution of each patient can be more accurately calculated, thereby improving the accuracy of dose calculation and ensuring the accuracy of the treatment plan; by obtaining the proton dose distribution of all voxels and generating a dose distribution map, the doctor can intuitively understand the dose distribution of the treatment area. Comparing the target dose with the calculated dose and adjusting the dose according to the comparison results can ensure that the dose distribution of radiotherapy is consistent with the plan, thereby guiding precise radiotherapy, improving the treatment effect and reducing damage to normal tissues. Since this method takes into account the individual differences and tissue characteristics of the human body, a personalized treatment plan can be tailored for each patient. This method allows personalized treatment planning, customizing the dose distribution according to the actual anatomical structure and tissue distribution of each patient. This helps to minimize the dose to healthy tissue while ensuring that the tumor area receives sufficient treatment dose; by combining the patient's CT information, radiotherapy delineation information, and information in the basic database, the proton dose distribution of each patient can be determined more accurately, thereby guiding doctors to provide patients with more personalized treatment plans. By improving the accuracy of dose calculations and guiding precise radiotherapy, this method can improve treatment outcomes. At the same time, reducing damage to normal tissue can reduce the incidence of complications and improve the quality of life of patients. This method can be applied to different types of proton therapy equipment and techniques, including fixed beams, scanning beams, and pencil beams.

[0145] See also Figure 4 The embodiment of the present application further provides an electronic device 200, which includes at least one memory 210, at least one processor 220, and a bus 230 connecting different platform systems. The at least one processor is configured to construct a basic database in the following manner when executing the computer program:

[0146] Building a basic database based on reports and documents; the basic database includes a first database and a second database;

[0147] The first database includes parameters of human tissues and organs;

[0148] The second database forms multiple sub-databases according to the division criteria; the division criteria include a first criterion and a second criterion; the first criterion is gender; within each first criterion, further division is performed according to the second criterion; the second criterion is the location and tissue distribution of the human anatomical structure, and multiple regions are divided according to the location and tissue distribution of the human anatomical structure; the second database includes materials of different densities in the region.

[0149] The basic database also includes the calculated α and β parameters;

[0150] The α parameter is used to scale the proton track, and the β parameter is used to scale the deflection angle of the track;

[0151] Among them, the α parameter is the ratio of the stopping power of water to the stopping power of tissues and organs. Each energy corresponds to an α parameter, and the energy range is determined according to the energy range of protons:

[0152]

[0153] in, The proton energy is E i When the α parameter is The proton energy is E i The stopping power in water, The proton energy is E i The ability to stop in the organization;

[0154] The β parameter is obtained by calculating the ratio of the Coulomb scattering angle of tissues and organs to the Coulomb scattering angle of water at different energies and then averaging them. One tissue and organ corresponds to only one β parameter:

[0155]

[0156] Where n is the number of energies; When the proton energy is E i Multiple Coulomb scattering angles of tissues and organs; When the proton energy is E i The multiple Coulomb scattering angle of water.

[0157] The at least one processor is configured to determine the density and the tissue organ to which each voxel belongs in the following manner when executing the computer program:

[0158] Obtain CT information and radiotherapy delineation information; determine the voxel range covered by each delineated tissue and organ through the radiotherapy delineation information;

[0159] The delineated voxels are assigned as the types of corresponding tissues and organs through the first database;

[0160] The region to which the voxels not outlined in the radiotherapy outline information belong is determined according to the position of the planned target area; and the tissue organ to which it belongs is determined according to the voxel density in combination with the second database.

[0161] The at least one processor is configured to acquire the proton dose of the voxel therapy in the following manner when executing the computer program:

[0162] Calculate voxel parameters according to the density of each voxel and the tissue and organ it belongs to;

[0163] The proton dose is calculated according to the radiotherapy plan using voxel parameters and a preset algorithm; the normalized absorbed proton dose of each voxel can be expressed as:

[0164]

[0165] Among them, D i is the absorbed dose of voxel i, E ki is the energy deposited by the kth particle in voxel i, m i is the mass of voxel i.

[0166] The at least one processor is configured to perform plan evaluation and adjustment guidance according to the comparison deviation in the following manner when executing the computer program:

[0167] Obtain treatment proton doses for all voxels and generate dose distribution maps for regions of interest; calculate dose indicators;

[0168] The target dose index is compared with the calculated dose index, and the plan is evaluated and adjusted based on the comparison deviation.

[0169] The memory 210 may include a readable medium in the form of a volatile memory, such as a random access memory (RAM) 211 and / or a cache memory 212 , and may further include a read-only memory (ROM) 213 .

[0170] Among them, the memory 210 also stores a computer program, which can be executed by the processor 220, so that the processor 220 executes the steps of the proton dose calculation method for tissue segmentation based on the outlined area in the embodiment of the present application. Its specific implementation method is consistent with the implementation method and the technical effect achieved in the embodiment of the proton dose calculation method for tissue segmentation based on the outlined area, and some contents will not be repeated here.

[0171] The memory 210 may also include a utility 214 having at least one program module 215, such program module 215 including but not limited to: an operating system, one or more application programs, other program modules, and program data, each of which or some combination may include an implementation of a network environment.

[0172] Accordingly, the processor 220 may execute the above-mentioned computer program and may execute the utility 214 .

[0173] The bus 230 may be a bus representing one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, a processor, or a local bus using any of a variety of bus architectures.

[0174] The electronic device 200 may also communicate with one or more external devices 240 such as keyboards, pointing devices, Bluetooth devices, etc., and may also communicate with one or more devices that can interact with the electronic device 200, and / or communicate with any device (such as a router, a modem, etc.) that enables the electronic device 200 to communicate with one or more other computing devices. This communication may be performed through an input / output interface 250. In addition, the electronic device 200 may also communicate with one or more networks (such as a local area network (LAN), a wide area network (WAN) and / or a public network, such as the Internet) through a network adapter 260. The network adapter 260 may communicate with other modules of the xxxx device 200 through a bus 230. It should be understood that, although not shown in the figure, other hardware and / or software modules may be used in conjunction with the electronic device 200, including but not limited to: microcode, device drivers, redundant processors, external disk drive arrays, RAID systems, tape drives, and data backup storage platforms, etc.

[0175] The embodiment of the present application also provides a computer-readable storage medium, which is used to store a computer program. When the computer program is executed, the steps of the proton dose calculation method for tissue segmentation based on delineated areas in the embodiment of the present application are implemented. The specific implementation method is consistent with the implementation method and the technical effect recorded in the embodiment of the proton dose calculation method for tissue segmentation based on delineated areas, and some contents are not repeated here.

[0176] An embodiment of the present application also provides a computer-readable storage medium, which is used to store a computer program. When the computer program is executed, the steps of a charger position deviation visualization method in the embodiment of the present application are implemented. The specific implementation method is consistent with the implementation method and the technical effect achieved recorded in the embodiment of the above-mentioned charger position deviation visualization method, and some contents are not repeated here.

[0177] In the present application, a readable storage medium can be any tangible medium containing or storing a program, which can be used by an instruction execution system, device or device or used in combination with it. The program product 300 can adopt any combination of one or more readable media. The readable medium can be a readable signal medium or a readable storage medium. The readable storage medium can be, for example, but not limited to, a system, device or device of electricity, magnetism, light, electromagnetic, infrared, or semiconductor, or any combination of the above. More specific examples (non-exhaustive list) of readable storage media include: an electrical connection with one or more wires, a portable disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above.

[0178] Computer-readable storage media may include data signals propagated in baseband or as part of a carrier wave, wherein readable program codes are carried. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. The readable storage medium may also be any readable medium, which may send, propagate, or transmit a program for use by an instruction execution system, an apparatus, or a device or for use in combination with it. The program code contained on the readable storage medium may be transmitted with any appropriate medium, including but not limited to wireless, wired, optical cable, RF, etc., or any suitable combination of the above. The program code for performing the operation of the present invention may be written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Java, C++, etc., and also conventional procedural programming languages ​​such as C language or similar programming languages. The program code may be executed entirely on a user computing device, partially on an associated device, as an independent software package, partially on a user computing device, partially on a remote computing device, or entirely on a remote computing device or server. Where a remote computing device is involved, the remote computing device may be connected to the user computing device through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computing device (e.g., through the Internet using an Internet service provider).

[0179] This application is explained from the perspectives of purpose of use, effectiveness, progress and novelty, and has met the functional enhancement and usage requirements emphasized by the Patent Law. The above description and drawings of this application are only the preferred embodiments of this application, and are not intended to limit this application. Therefore, all structures, devices, features, etc. that are similar or identical to this application, that is, all equivalent replacements or modifications made in accordance with the scope of the patent application of this application, should fall within the scope of protection of the patent application of this application.

[0180] Although the embodiments of the present invention have been shown and described above, it is to be understood that the above embodiments are exemplary and are not to be construed as limitations on the present invention. A person skilled in the art may change, modify, substitute and modify the above embodiments within the scope of the invention without departing from the principles and purpose of the present invention. All such changes should fall within the scope of protection of the claims of the present invention.

Claims

1. A proton dose calculation method for tissue segmentation based on delineated regions, characterized in that: The method comprises: S1. Constructing a basic database of human tissues and organs; the basic database includes a first database and a second database, both of which include names and densities of human tissues and organs, the density value in the first database is not unique, while the density value in the second database is unique, the first database has radiotherapy delineation information, while the second database does not have radiotherapy delineation information; S2, obtaining CT information and radiotherapy delineation information; determining the density of each voxel and the tissue organ to which it belongs based on the CT information and the radiotherapy delineation information and in combination with the corresponding first database and second database in the basic database; S3. Obtain the proton absorption dose of voxel therapy according to the density of each voxel and the tissue and organ to which it belongs through a preset algorithm; S4. Obtain the proton absorption dose distribution of all voxels and generate a dose distribution map; compare the actual dose with the target dose, and conduct plan evaluation and adjustment guidance based on the comparison deviation; Said S1 comprises: said first database comprises parameters of human tissues and organs; The second database forms a plurality of sub-databases according to the division criteria; the division criteria include a first criterion and a second criterion; the first criterion is gender; within each first criterion, further division is performed according to the second criterion; the second criterion is the location and tissue distribution of the human anatomical structure, and a plurality of regions are divided according to the location and tissue distribution of the human anatomical structure; the second database includes materials of different densities in the regions; The basic database also includes α parameters and β parameters; The α parameter is the ratio of the stopping power of water to the stopping power of tissues and organs, and is used to scale the proton track; The β parameter is used to scale the deflection angle of the proton track; The β parameter is obtained by averaging the ratio of the Coulomb scattering angle of the tissue organ to the Coulomb scattering angle of water at different energies; One tissue or organ corresponds to only one β parameter; The S3 includes: Calculate voxel parameters according to the density of each voxel and the tissue and organ it belongs to; The proton absorbed dose is calculated according to the radiotherapy plan using voxel parameters and a preset algorithm; the normalized proton absorbed dose of each voxel is: Among them, D i is the absorbed dose of voxel i, E ki is the energy deposited by the kth particle in voxel i, m i is the mass of voxel i.

2. The proton dose calculation method for tissue segmentation based on delineated regions according to claim 1, characterized in that: The S1 further comprises: Build a basic database based on reports and literature.

3. The proton dose calculation method for tissue segmentation based on delineated regions according to claim 1, characterized in that: Each energy corresponds to an α parameter, and the energy range is determined by the energy range of the proton: in, The proton energy is E i When the α parameter is The proton energy is E i The stopping power in water, The proton energy is E i The ability to stop in the organization; Where n is the number of energies; When the proton energy is E i Multiple Coulomb scattering angles of tissues and organs; When the proton energy is E i The multiple Coulomb scattering angle of water.

4. The proton dose calculation method for tissue segmentation based on delineated regions according to claim 1, characterized in that: The S2 includes: Obtain CT information and radiotherapy delineation information; determine the voxel range covered by each delineated tissue and organ through the radiotherapy delineation information; The delineated voxels are assigned as the types of corresponding tissues and organs through the first database; The region to which the voxels not outlined in the radiotherapy outline information belong is determined according to the position of the planned target area; and the tissue organ to which it belongs is determined according to the voxel density in combination with the second database.

5. The proton dose calculation method for tissue segmentation based on delineated regions according to claim 1, characterized in that: The S4 includes: Obtain the proton dose of all voxels and generate the dose distribution map of the region of interest, and calculate the dose index; The target dose index is compared with the calculated dose index, and the plan is evaluated and adjusted based on the comparison deviation.

6. A proton dose calculation system for tissue segmentation based on delineated regions, characterized in that: The system comprises: A database construction module, used to construct a basic database of human tissues and organs; the basic database includes a first database and a second database; the first database and the second database both include names and densities of human tissues and organs, the density value in the first database is not unique, the density value in the second database is unique, the first database has radiotherapy delineation information, and the second database does not have radiotherapy delineation information; A first determination module is used to obtain CT information and radiotherapy delineation information; according to the CT information and radiotherapy delineation information, combined with the corresponding first database and second database in the basic database, determine the density of each voxel and the tissue organ to which it belongs; The dose acquisition module is used to obtain the proton absorption dose of the voxel therapy according to the density of each voxel and the tissue organ to which it belongs through a preset algorithm; The adjustment module is used to obtain the proton absorption dose distribution of all voxels and generate a dose distribution map; compare the actual dose with the target dose, and provide dose adjustment guidance based on the comparison deviation; The first database includes parameters of human tissues and organs; The second database forms a plurality of sub-databases according to the division criteria; the division criteria include a first criterion and a second criterion; the first criterion is gender; within each first criterion, further division is performed according to the second criterion; the second criterion is the location and tissue distribution of the human anatomical structure, and a plurality of regions are divided according to the location and tissue distribution of the human anatomical structure; the second database includes materials of different densities in the regions; The basic database also includes α parameters and β parameters; The α parameter is the ratio of the stopping power of water to the stopping power of tissues and organs, and is used to scale the proton track; The β parameter is used to scale the deflection angle of the proton track; The β parameter is obtained by averaging the ratio of the Coulomb scattering angle of the tissue organ to the Coulomb scattering angle of water at different energies; One tissue or organ corresponds to only one β parameter; The dose acquisition module comprises: Calculate voxel parameters according to the density of each voxel and the tissue and organ it belongs to; The proton absorbed dose is calculated according to the radiotherapy plan using voxel parameters and a preset algorithm; the normalized proton absorbed dose of each voxel is: Among them, D i is the absorbed dose of voxel i, E ki is the energy deposited by the kth particle in voxel i, m i is the mass of voxel i.

7. An electronic device, characterized in that: The electronic device comprises a memory and at least one processor, the memory stores a computer program, and the at least one processor is configured to implement the following steps when executing the computer program: S1. Constructing a basic database of human tissues and organs; the basic database includes a first database and a second database; the first database and the second database both include names and densities of human tissues and organs, the density value in the first database is not unique, the density value in the second database is unique, the first database has radiotherapy delineation information, and the second database does not have radiotherapy delineation information; S2, obtaining CT information and radiotherapy delineation information; determining the density of each voxel and the tissue and organ to which it belongs based on the CT information and radiotherapy delineation information and in combination with the basic database; S3. Obtain the proton absorption dose of voxel therapy according to the density of each voxel and the tissue and organ to which it belongs through a preset algorithm; S4. Obtain the proton absorption dose distribution of all voxels and generate a dose distribution map; compare the actual dose with the simulated dose, and provide dose adjustment guidance based on the comparison deviation; The first database includes parameters of human tissues and organs; The second database forms a plurality of sub-databases according to the division criteria; the division criteria include a first criterion and a second criterion; the first criterion is gender; within each first criterion, further division is performed according to the second criterion; the second criterion is the location and tissue distribution of the human anatomical structure, and a plurality of regions are divided according to the location and tissue distribution of the human anatomical structure; the second database includes materials of different densities in the regions; The basic database also includes α parameters and β parameters; The α parameter is the ratio of the stopping power of water to the stopping power of tissues and organs, and is used to scale the proton track; The β parameter is used to scale the deflection angle of the proton track; The β parameter is obtained by averaging the ratio of the Coulomb scattering angle of the tissue organ to the Coulomb scattering angle of water at different energies; One tissue or organ corresponds to only one β parameter; The S3 includes: Calculate voxel parameters according to the density of each voxel and the tissue and organ it belongs to; The proton absorbed dose is calculated according to the radiotherapy plan using voxel parameters and a preset algorithm; the normalized proton absorbed dose of each voxel is: Among them, D i is the absorbed dose of voxel i, E ki is the energy deposited by the kth particle in voxel i, m i is the mass of voxel i.

8. The electronic device according to claim 7, characterized in that: The at least one processor is configured to construct a basic database in the following manner when executing the computer program: Build a basic database based on reports and literature.

9. The electronic device according to claim 7, characterized in that: Each energy corresponds to an α parameter, and the energy range is determined by the energy range of the proton: in, The proton energy is E i When the α parameter is The proton energy is E i The stopping power in water, The proton energy is E i The ability to stop in the organization; Where n is the number of energies; When the proton energy is E i Multiple Coulomb scattering angles of tissues and organs; When the proton energy is E i The multiple Coulomb scattering angle of water.

10. The electronic device according to claim 7, characterized in that: The at least one processor is configured to determine the density and the tissue organ to which each voxel belongs in the following manner when executing the computer program: Obtain CT information and radiotherapy delineation information; determine the voxel range covered by each delineated tissue and organ through the radiotherapy delineation information; The delineated voxels are assigned as the types of corresponding tissues and organs through the first database; The region to which the voxels not outlined in the radiotherapy outline information belong is determined according to the position of the planned target area; and the tissue organ to which it belongs is determined according to the voxel density in combination with the second database.

11. The electronic device according to claim 7, characterized in that: The at least one processor is configured to perform plan evaluation and adjustment guidance according to the comparison deviation in the following manner when executing the computer program: Obtain treatment proton doses for all voxels and generate dose distribution maps for regions of interest; calculate dose indicators; The target dose index is compared with the calculated dose index, and the plan is evaluated and adjusted based on the comparison deviation.

12. A computer-readable storage medium storing computer instructions, wherein when a computer reads the computer instructions, the computer executes the method according to any one of claims 1 to 5.

Images