BNCT radiation dose acquisition method, radiation termination time acquisition method, terminal and medium

Abstract

The present invention provides a BNCT radiation dose acquisition method, a radiation termination time acquisition method, a terminal and a medium. The method acquires a radiation field model and a characteristic element concentration model, and performs spatiotemporal alignment on the radiation field model and the characteristic element concentration model to obtain a BNCT comprehensive model; based on the BNCT comprehensive model, the radiation field intensity and the characteristic element concentration corresponding to each voxel point in the tumor area at each time are acquired, so as to acquire the radiation dose of the voxel point at the corresponding time based on the radiation field intensity and the characteristic element concentration; the present invention can quickly and conveniently acquire the radiation dose value corresponding to each voxel point at each time, which not only improves the efficiency of obtaining the radiation dose, but also improves the precision and accuracy of obtaining the radiation dose.

Description

Technical Field

[0001] The present invention relates to the technical field of boron neutron capture, and in particular to a method for obtaining a BNCT radiation dose, a method for obtaining a radiation termination time, a terminal and a computer storage medium. Background Art

[0002] During boron neutron capture therapy (BNCT), it is often necessary to obtain the radiation dose distribution in the tumor area so as to accurately set or adjust the radiation time and radiation plan based on the radiation dose distribution, so as to minimize the radiation damage to normal tissues in the human body while achieving the treatment effect.

[0003] Existing radiation dose monitoring methods are mainly based on the monitored blood boron concentration, combined with the concentration of related elements associated with the radiation dose (such as hydrogen, carbon, oxygen, etc.) and the radiation intensity in the tumor area, to obtain the radiation dose distribution in the tumor area through calculation; however, this acquisition method usually requires multiple blood draws and blood tests, and multiple collections of radiation intensity in the tumor area, resulting in a cumbersome acquisition method, low acquisition efficiency and high cost, and unable to achieve real-time monitoring of the radiation dose in the tumor area; and this acquisition method is often based only on the blood boron concentration at a single collection moment, and does not fully consider the impact of individual differences (such as differences in boron drug absorption efficiency and migration efficiency, etc.) on the spatiotemporal distribution of blood boron concentration, that is, it is impossible to accurately obtain the blood boron concentration corresponding to each moment, nor is it possible to obtain the blood boron concentration corresponding to each sampling point, resulting in the radiation dose obtained based on the blood boron concentration often being less accurate, which in turn leads to poor results of the radiation plan formulated / adjusted based on the radiation dose. Summary of the invention

[0004] In view of the shortcomings of the prior art described above, the purpose of the present invention is to provide a BNCT radiation dose acquisition method, a radiation termination time acquisition method, a terminal and a computer storage medium, which can solve the problems of low acquisition efficiency, inability to monitor in real time and low acquisition accuracy of the existing radiation dose acquisition methods.

[0005] To achieve the above-mentioned and other related purposes, the present invention provides a BNCT radiation dose acquisition method in a first aspect, comprising:

[0006] Acquire a radiation field model, and acquire a characteristic element concentration model; wherein the radiation field model is a model for characterizing the radiation intensity distribution in the tumor region at each moment, and the characteristic element concentration model is a model for characterizing the characteristic element concentration distribution in the tumor region at each moment; perform spatiotemporal alignment on the radiation field model and the characteristic element concentration model to obtain a BNCT comprehensive model; based on the BNCT comprehensive model, acquire the radiation field intensity and characteristic element concentration corresponding to each voxel point in the tumor region at each moment, so as to acquire the radiation dose of the voxel point at the corresponding moment based on the radiation field intensity and the characteristic element concentration; wherein the characteristic elements include associated elements and boron elements, and the associated elements are basic elements of the human body associated with the BNCT radiation dose.

[0007] In some embodiments of the first aspect of the present invention, the implementation method of performing spatiotemporal alignment on the radiation field model and the characteristic element concentration model includes:

[0008] Based on a preset target resolution, a spatial grid of the tumor area is constructed; the grid points in the spatial grid are used as the voxel points; the radiation field intensity distribution and the characteristic element concentration distribution at the same time are extracted; based on the target resolution, the radiation field intensity distribution and the characteristic element concentration distribution at the same time are resampled respectively; the resampled radiation field intensity distribution and the characteristic element concentration distribution are spatially aligned with the spatial grid respectively to obtain the radiation field intensity and the characteristic element concentration corresponding to each voxel point.

[0009] In some embodiments of the first aspect of the present invention, the method for obtaining the radiation field model includes:

[0010] Determine an associated device in a BNCT working environment, wherein the associated device is a radiation device associated with the radiation field intensity; among various device parameters of the associated device, determine the device parameter related to the radiation field intensity as a target parameter; based on a preset radiation scheme, in combination with the three-dimensional model information of the associated device, obtain the spatial information and target parameter information of the associated device at each moment; based on the spatial information and the target parameter information corresponding to each moment, generate the radiation field model using twin technology; wherein the radiation scheme is a pre-set operating scheme of the associated device; the radiation field model is a model that uses a radiation field intensity calculation method to obtain the radiation intensity value corresponding to each first pixel point at each moment according to the spatial information and target parameter information of the associated device.

[0011] In some embodiments of the first aspect of the present invention, the method for acquiring the characteristic element concentration model includes:

[0012] Acquire the tumor area; acquire the associated element concentration corresponding to each second pixel point in the tumor area at each moment, and acquire the boron element concentration corresponding to each second pixel point in the tumor area at each moment; combine the associated element concentration and the boron element concentration corresponding to the same moment to obtain the characteristic element concentration corresponding to each moment; wherein, the method for acquiring the associated element concentration corresponding to each second pixel point in the tumor area at each moment includes: extracting the associated element concentration corresponding to each second pixel point at each moment according to the CT / MRI value corresponding to each second pixel point in the tumor area.

[0013] In some embodiments of the first aspect of the present invention, the method for extracting the concentration of the associated element corresponding to each second pixel point at each moment according to the CT / MRI value corresponding to each second pixel point in the tumor region includes:

[0014] According to the conversion relationship between the CT / MRI value and the tissue density, the tissue density corresponding to each second pixel point is obtained; using the interval classification method, combined with the tissue density corresponding to each second pixel point, each second pixel point is classified into each type of tissue; based on the tissue type corresponding to each second pixel point, the concentration of the associated element corresponding to each second pixel point is determined.

[0015] In some embodiments of the first aspect of the present invention, the method for obtaining the boron concentration corresponding to each second pixel point in the tumor area at each time includes:

[0016] The signed distance field method is used to calculate the distance function from each second pixel point in the tumor area to the tumor boundary; the spatial position corresponding to each second pixel point is determined based on the distance function; based on the spatial position information corresponding to each second pixel point and the blood boron concentration value in the human body at each moment, the diffusion equation is used to obtain the boron element concentration value corresponding to each voxel point at each moment.

[0017] In some embodiments of the first aspect of the present invention, the method for obtaining the boron concentration value corresponding to each voxel point at each moment by using the diffusion equation based on the spatial position information corresponding to each second pixel point and the blood boron concentration value in the human body at each moment includes:

[0018] ;

[0019] ;

[0020] Among them, C is the blood boron concentration, D is the diffusion coefficient, R is the reaction coefficient of neutron radiation dose to boron concentration, and k is the reaction coefficient.

[0021] In order to solve the problems in the prior art of low efficiency and low accuracy in obtaining radiation doses, which leads to low efficiency and low accuracy in obtaining the radiation duration corresponding to the radiation scheme, the present invention also provides a method for obtaining the radiation termination time in the second aspect, including: based on the current radiation scheme, constructing a radiation field model corresponding to the current radiation scheme; based on the current CT / MRI image and the current blood boron concentration value, constructing the characteristic element concentration model; based on the radiation field model and the characteristic element concentration model, using any of the BNCT radiation dose acquisition methods described above, to obtain the radiation dose corresponding to each voxel point in the tumor area at each time; according to the preset termination conditions, the radiation dose corresponding to each voxel point at each time is sequentially detected based on the termination conditions, and the time when the termination conditions are met is used as the radiation termination time.

[0022] In some embodiments of the second aspect of the present invention, the implementation method of sequentially performing detection based on the termination condition on the radiation dose corresponding to each voxel point at each moment includes: obtaining the radiation dose corresponding to each voxel point at the current moment; detecting whether the radiation dose corresponding to each voxel point satisfies the termination condition; if so, determining that the radiation dose at the current moment satisfies the termination condition; if not, determining that the radiation dose at the current moment does not satisfy the termination condition, and re-executing the detection process at the next moment.

[0023] The present invention provides a terminal in a third aspect, comprising: a processor and a memory, wherein the memory is communicatively connected to the processor; the memory is used to store a computer program, and the processor is used to execute the computer program stored in the memory, so that the terminal executes any of the BNCT radiation dose acquisition methods described above or any of the radiation termination time acquisition methods described above.

[0024] In a fourth aspect, the present invention provides a computer storage medium storing a computer program, wherein when the computer program is executed by a processor, the computer program implements any of the above-described BNCT radiation dose acquisition methods or any of the above-described radiation termination time acquisition methods.

[0025] As described above, the BNCT radiation dose acquisition method, the radiation termination time acquisition method, the terminal and the computer storage medium provided by the present invention obtain a radiation field model for characterizing the spatiotemporal distribution of the radiation field intensity, and obtain a characteristic element concentration model for characterizing the spatiotemporal distribution of the characteristic element concentration, and align the two in time and space to obtain a BNCT comprehensive model, so as to obtain the radiation dose of each voxel point at each time based on the BNCT comprehensive model, so that it is impossible to perform the measurement of blood boron concentration value and radiation field intensity multiple times. Based on the model, the radiation dose value corresponding to each voxel point at each time can be quickly and conveniently obtained, which not only improves the efficiency of obtaining the radiation dose, realizes real-time monitoring of the radiation dose, but also improves the precision and accuracy of obtaining the radiation dose, and provides effective support and guarantee for the adjustment / formulation of subsequent radiation plans. BRIEF DESCRIPTION OF THE DRAWINGS

[0026] Figure 1 Shown is a flow chart of a method for obtaining a BNCT radiation dose according to an embodiment of the present invention;

[0027] Figure 2 It is a flowchart showing an implementation method of step S200 in one embodiment of the present invention;

[0028] Figure 3 Shown is a schematic flow chart of a radiation field model acquisition process according to an embodiment of the present invention;

[0029] Figure 4 Shown is a schematic flow chart of a process for acquiring a characteristic element concentration model in one embodiment of the present invention;

[0030] Figure 5 Shown is a schematic diagram of a process of extracting the concentration of the associated element corresponding to each of the second pixel points at each time by using the interval classification method in one embodiment of the present invention;

[0031] Figure 6 Shown is a flow chart of a method for obtaining the radiation suspension time in one embodiment of the present invention;

[0032] Figure 7 It is a flowchart diagram showing a process of sequentially performing detection based on the termination condition on the radiation dose corresponding to each voxel point at each time in an embodiment of the present invention;

[0033] Figure 8 Shown is a schematic diagram of the structure of the terminal in one embodiment of the present invention. DETAILED DESCRIPTION

[0034] The following describes the embodiments of the present application through specific examples, and those skilled in the art can easily understand other advantages and effects of the present application from the contents disclosed in this specification. The present application can also be implemented or applied through other different specific embodiments, and the details in this specification can also be modified or changed in various ways based on different viewpoints and applications without departing from the spirit of the present application. It should be noted that the following embodiments and features in the embodiments can be combined with each other without conflict.

[0035] It should be noted that in the following description, reference is made to the accompanying drawings, which describe several embodiments of the present application. It should be understood that other embodiments may be used, and mechanical, structural, electrical, and operational changes may be made without departing from the spirit and scope of the present application.

[0036] Also, as used herein, the singular forms "a", "an", and "the" are intended to include the plural forms as well, unless the context indicates otherwise. It should be further understood that the terms "comprises", "includes", "includes" indicate the presence of the described features, operations, elements, components, items, kinds, and / or groups, but do not exclude the presence, occurrence, or addition of one or more other features, operations, elements, components, items, kinds, and / or groups. The terms "or" and "and / or" used herein are interpreted as inclusive, or mean any one or any combination.

[0037] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the technical solution in the embodiments of the present invention is further described in detail through the following embodiments and in combination with the accompanying drawings. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the invention.

[0038] In order to solve the technical problems existing in the prior art, the present application first provides a BNCT radiation dose acquisition method, which is used for the radiation dose corresponding to each voxel point at each moment; wherein the radiation dose is a simulation value.

[0039] The voxel points are spatial pixel points in the tumor area.

[0040] See also Figure 1 , showing a flow chart of the BNCT radiation dose acquisition method provided in an embodiment of the present application; Figure 1 As shown, the method comprises the following steps:

[0041] S100, obtaining a radiation field model and a characteristic element concentration model;

[0042] Among them, the radiation field model is used to characterize the neutron source radiation field intensity distribution (hereinafter referred to as "radiation field intensity distribution") in the tumor area at each time; the radiation intensity distribution is the radiation intensity corresponding to each first pixel point in the tumor area; the first pixel point is each sampling point determined based on the first sampling resolution.

[0043] The characteristic element concentration model is a model used to characterize the characteristic element concentration distribution in the tumor area at each moment; the characteristic element concentration distribution is the characteristic element concentration corresponding to each second pixel point in the tumor area; the second pixel point is each sampling point determined based on the first sampling resolution.

[0044] The characteristic elements include associated elements and boron elements. The associated elements are basic elements of the human body associated with the BNCT radiation dose, including C elements, N elements, H elements, and the like.

[0045] Specifically, the radiation field model includes device data corresponding to each moment, and the device data is radiation device data related to the radiation field intensity, including but not limited to the radiation source power, radiation source position and beam angle of the radiation device. In the radiation field model, based on the device data corresponding to each moment, the radiation intensity value corresponding to each of the first pixel points at the corresponding moment is obtained using a radiation field intensity calculation method.

[0046] Among them, the radiation field intensity calculation method includes but is not limited to Monte Carlo simulation method, finite element analysis method, finite difference analysis method or other existing radiation field intensity calculation methods; exemplarily, according to the radiation source position, beam angle and radiation source power at the i-th moment, the Monte Carlo simulation method is used to obtain the radiation intensity value corresponding to each of the first pixel points in the tumor area at the i-th moment.

[0047] The characteristic element concentration model includes biological tissue data corresponding to each moment in the tumor region, and the biological tissue data is data related to the associated element, including but not limited to CT value, MRI value or other existing biological tissue value. In the characteristic element concentration model, based on the biological tissue data corresponding to each moment, the associated element value corresponding to each second pixel point in the tumor region at the corresponding moment is calculated;

[0048] Furthermore, the characteristic element concentration model also includes the blood boron concentration value of the human body at each time in the tumor area. Based on the blood boron concentration value corresponding to each time, the diffusion equation is used to obtain the boron element concentration value corresponding to each second pixel point in the tumor area.

[0049] It should be noted that, in the embodiment of the present application, the first sampling resolution and the second sampling resolution may be the same or different.

[0050] S200, performing spatiotemporal alignment on the radiation field model and the characteristic element concentration model to obtain a BNCT comprehensive model;

[0051] Wherein, the BNCT comprehensive model is a model that simultaneously characterizes the spatiotemporal distribution characteristics of the radiation field intensity and the characteristic element concentration;

[0052] Specifically, the radiation field model data and the characteristic element concentration model data are aligned in the time dimension and the space dimension respectively to ensure the consistency of each data in the time and space dimensions, so as to obtain the radiation field intensity signal distribution and the characteristic element concentration distribution in the tumor area at each time as the BNCT comprehensive model.

[0053] In a specific embodiment, the implementation of step S200 is as follows: Figure 2 As shown, including:

[0054] S201, constructing a spatial grid of the tumor area based on a preset target resolution;

[0055] The target resolution is a sampling resolution that characterizes the radiation dose in the tumor region, that is, the spatial distance between two adjacent radiation dose sampling points.

[0056] Specifically, three-dimensional spatial data of the patient's tumor area is obtained; based on a preset target resolution, grid division is performed on the three-dimensional spatial data to construct a spatial grid of the patient's tumor area, and each grid point in the spatial grid is set as a voxel point.

[0057] S202, extracting the radiation field intensity distribution and the characteristic element concentration distribution at the same time to achieve time alignment of each data;

[0058] S203, based on the target resolution, resample the radiation field intensity distribution and the characteristic element concentration distribution at the same time respectively; spatially align the resampled radiation field intensity distribution and the characteristic element concentration distribution with each of the spatial grids respectively to achieve spatial alignment of each data.

[0059] Specifically, after resampling, the new first sampling resolution of the radiation field intensity distribution and the new second sampling resolution of the characteristic element concentration distribution are the same as the target resolution; the first pixel point and the second pixel point representing the same sampling space are respectively aligned with the voxel point corresponding to the sampling space, thereby achieving spatial alignment of each data.

[0060] It should be noted that, in other specific embodiments, multimodal data fusion technology can also be used to fuse the radiation field model data with the characteristic element concentration model data in the time dimension and the space dimension; wherein the multimodal data fusion technology includes but is not limited to principal component analysis PCA or independent component analysis ICA.

[0061] S300, based on the BNCT comprehensive model, according to the radiation field intensity and characteristic element concentration corresponding to each voxel point in the tumor area at each moment, the radiation dose corresponding to each voxel point at the corresponding moment is obtained.

[0062] Specifically, according to the BNCT comprehensive model, the radiation field intensity value and the characteristic element concentration value corresponding to each voxel point at each treatment moment are obtained; based on the characteristic element concentration value and the radiation field intensity value, the relative biological effectiveness (RBE) calculation method is used to obtain the radiation dose corresponding to each voxel point.

[0063] Among them, the implementation method of obtaining the radiation dose corresponding to each voxel point using the RBE calculation method based on the characteristic element concentration value combined with the radiation field intensity value is a prior art and will not be expanded here.

[0064] The BNCT radiation dose acquisition method provided in the embodiment of the present application can quickly and conveniently acquire the radiation intensity and characteristic element concentration of the same voxel point in the tumor area by separately constructing a radiation field model and a characteristic element concentration model, and by performing spatiotemporal alignment of the radiation field model and the characteristic element concentration model, that is, the radiation intensity and characteristic element concentration at the same spatial position can be quickly acquired, so that the radiation dose at each spatial position can be accurately acquired, which greatly improves the accuracy and precision of radiation dose acquisition and provides convenience for the subsequent analysis and adjustment of the BNCT radiation plan.

[0065] In order to quickly and conveniently obtain the radiation field model, in some optional embodiments, the radiation field model is obtained by a modeling process using twin technology, such as Figure 3 As shown, including:

[0066] S101, determining associated devices in the BNCT working environment;

[0067] The associated device is a radiation device associated with the radiation field intensity distribution in the BNCT working environment, such as a collimator or a radiation platform.

[0068] S102, determining, among the device parameters of the associated devices, a device parameter related to the radiation field intensity distribution as a target parameter;

[0069] Specifically, after determining the associated device, various device parameters of the associated device are obtained, including configuration parameters, working parameters, and state parameters, etc.; among the device parameters, the device parameters associated with the radiation field intensity distribution are extracted, and the extracted device parameters are set as the target parameters, such as radiation energy, radiation angle, etc.

[0070] S103, based on the preset radiation plan and combined with the three-dimensional model information of the associated equipment, obtain the spatial information and target parameter information of the associated equipment at each moment; based on the spatial information and the target parameter information corresponding to each moment, use the twin technology to generate a twin model, and use the twin model as the radiation field model.

[0071] Among them, the radiation scheme is a pre-set operation scheme of the associated equipment, such as the operation position, state and parameters of the associated equipment at each time; illustratively, the radiation scheme includes the collimator having a first beam angle from time T1 to time T2, a second beam angle from time T2 to time T3, a third beam angle from time T3 to time T4, and so on.

[0072] Specifically, using deep time camera technology or three-dimensional scanning technology, the associated equipment in the BNCT working environment is scanned and imaged to obtain three-dimensional model data corresponding to the associated equipment; based on the operation plan at each moment in the radiation plan, combined with the three-dimensional model data of the associated equipment, the position information, size information, posture information, etc. of the associated equipment are obtained as the spatial information of the associated equipment;

[0073] And, according to the operating parameters or operating status of the associated device at each moment in the radiation plan, the target parameter information of the associated device is obtained.

[0074] After extracting the spatial information and target parameter information of the associated equipment at each moment in the radiation plan, the spatial information and the target parameter information corresponding to each moment are respectively imported into a pre-built data processing platform; wherein, the data processing platform is a software platform built based on twin technology for realizing BNCT simulation.

[0075] In the data processing platform, a twin model for performing radiation field intensity calculation is generated, and the twin model is used as the radiation field model corresponding to the radiation scheme; the radiation field model is a model that obtains the radiation intensity value corresponding to each of the first pixel points at each moment based on the spatial information and target parameter information of the associated equipment and using the radiation field intensity calculation method.

[0076] More specifically, the spatial information includes the radiation distance of the associated device; the target parameter information includes the radiation angle and the radiation energy; the radiation distance, radiation angle and radiation energy corresponding to each moment are respectively imported into a pre-built data processing platform to generate the corresponding radiation field model; based on this, the radiation field model is based on the radiation distance and radiation angle corresponding to the same moment, combined with the radiation energy, and using the radiation field intensity calculation method to obtain the radiation intensity value corresponding to each of the first pixel points at each moment.

[0077] In some optional embodiments, since part of the target parameter information is associated with spatial information such as the position information or posture information of the associated device; for example, the radiation distance (target parameter) is associated with the position and emission angle of the collimator (associated device), then when the spatial information and the target parameter information corresponding to each moment are respectively imported into a preset data processing platform, it also includes:

[0078] Based on the spatial information of the associated device, performing a numerical check on the target parameter associated with the spatial information;

[0079] Specifically, based on the spatial information of the associated device, parameters related to the spatial information are extracted from each of the target parameters as related target parameters; according to the spatial information of the associated device, parameter values ​​of the related target parameters are obtained, and the target parameter values ​​are verified using the parameter values ​​to further improve the accuracy of each of the target parameter values ​​input into the data processing platform.

[0080] The BNCT radiation dose acquisition method provided in the embodiment of the present application is based on the three-dimensional model information of the radiation plan and the associated equipment, and uses twin technology to construct a radiation field twin model. It can not only realize the rapid integration and spatiotemporal synchronous display of radiation field related data, but also can realize the rapid acquisition of the radiation intensity information corresponding to each first pixel point at each moment based on the radiation field twin model, thereby improving the efficiency of BNCT radiation dose acquisition and providing support for real-time monitoring of BNCT radiation dose.

[0081] In order to quickly and accurately obtain the characteristic element concentration model, in an optional embodiment, the characteristic element concentration model is obtained in a manner such as Figure 4 As shown, including:

[0082] S201, acquiring the tumor region;

[0083] Specifically, CT / MRI images covering the tumor area are collected; and the tumor area is extracted based on the CT / MRI images, which is a three-dimensional spatial range.

[0084] S202, obtaining the concentration of the associated element corresponding to each of the second pixel points in the tumor area at each time;

[0085] Specifically, according to the CT / MRI value corresponding to each second pixel point in the tumor region, the concentration of the associated element corresponding to each second pixel point at each time is extracted using an interval classification method.

[0086] S203, obtaining the boron concentration corresponding to each of the second pixel points in the tumor area at each time;

[0087] Specifically, the blood boron concentration value in the human body at each moment is obtained; based on the blood boron concentration value, the boron concentration value corresponding to each of the second pixel points in the tumor area is obtained using the diffusion formula of the boron element.

[0088] It is understandable that those skilled in the art may also adopt, for example, an empirical formula for boron element distribution to obtain the boron element concentration value corresponding to each second pixel point, which will not be elaborated here.

[0089] S204, combining the associated element concentration and the boron element concentration corresponding to the same moment to obtain the characteristic element concentration corresponding to each moment.

[0090] It should be noted that the step S202 and the step S203 may be performed simultaneously or sequentially, which is not specifically limited here.

[0091] When the radiation field model is constructed based on the twin technology, in order to further improve the execution efficiency and effect of the fusion of the radiation field model and the characteristic element concentration model, in some optional embodiments, the acquisition process of the characteristic element concentration model is a modeling process using the twin technology, that is, obtaining the CT / MRI images at each moment, extracting the CT / MRI values ​​corresponding to each second pixel point in the CT / MRI images; importing the CT / MRI values ​​corresponding to each second pixel point at each moment into a preset data processing platform, so as to calculate the associated element concentration corresponding to each second pixel point at each moment in the data processing platform according to the CT / MRI values ​​corresponding to each second pixel point, using the interval classification method, thereby generating a twin associated element concentration model.

[0092] In a specific implementation, the interval classification method is used to extract the concentration of the associated element corresponding to each of the second pixel points at each time, such as Figure 5 As shown, including:

[0093] S202A, acquiring the tissue density corresponding to each second pixel point according to the conversion relationship between the CT / MRI value and the tissue density;

[0094] S202B, using an interval classification method and combining the tissue density corresponding to each second pixel point, classifying each second pixel point into each type of tissue;

[0095] S202C, based on the tissue type corresponding to each second pixel point, determine the concentration of the associated element corresponding to each second pixel point.

[0096] There is a preset relationship between the CT / MRI value and the tissue density, that is, the tissue density value corresponding to each second pixel point can be determined according to the CT / MRI value corresponding to the second pixel point;

[0097] The tissue types include air, soft tissue, bone, muscle, etc., and different tissue types have elemental compositions and elemental concentrations of specific elements; the specific elements are used as associated elements corresponding to the tissue types, and the elemental concentrations of the specific elements are used as associated element concentrations.

[0098] Exemplarily, the associated elements include H, C, O, Ca, etc.

[0099] In some preferred embodiments, the step S203, when specifically performed, includes:

[0100] Based on the blood boron concentration value at each moment, the signed distance field (SDF) method is used to obtain the boron concentration of each second pixel point in the tumor area.

[0101] Specifically, using the SDF method, a distance function from each second pixel point in the tumor region to the tumor boundary is calculated; based on the distance function, a spatial position corresponding to each second pixel point is determined, that is, coordinate information of each second pixel point in the tumor region;

[0102] Based on the spatial position information corresponding to each second pixel point and based on the blood boron concentration value in the human body at each moment, the boron element concentration value corresponding to each second pixel point at each moment is obtained by constructing a diffusion equation.

[0103] More specifically, the implementation method of obtaining the boron element concentration value corresponding to each second pixel at each time by constructing a diffusion equation based on the spatial position information corresponding to each second pixel and the blood boron concentration value of the patient at each time is:

[0104] ;

[0105] ;

[0106] in, for C The rate at which a value changes; C is the blood boron concentration value at each time;D is the diffusion coefficient; is the value corresponding to each of the second pixels at each moment R value; R is the response coefficient of neutron radiation dose to boron concentration, which is used to characterize the response mechanism of boron in the tumor area; k is the reaction coefficient; The corresponding I value; I is the neutron radiation dose distribution; ( x, y, z ) is the spatial position information of each of the second pixel points.

[0107] It should be noted that when the tumor area is small, the radiation area range of BNCT is relatively fixed, that is, the radiation area range of BNCT does not change during the radiation process; accordingly, the concentration information of each associated element in the radiation area is relatively fixed; when the step S202 is executed, it can also be:

[0108] Before the radiation scheme is started, a CT / MRI image corresponding to the tumor area at a certain moment is obtained; based on the CT / MRI image, the CT / MRI value corresponding to each second pixel point in the tumor area is obtained; using the interval classification method, the associated element concentration corresponding to each second pixel point at that moment is extracted; and the associated element concentration corresponding to that moment is used as the associated element concentration corresponding to each moment after the radiation scheme is started.

[0109] In addition, it should be noted that, in some embodiments, since the blood boron concentration value in the human body is usually obtained by blood drawing and testing; therefore, when step S203 is executed, the blood boron concentration value in the human body at each moment is obtained based on the blood boron concentration value actually measured at a certain moment, deduced using the boron concentration diffusion model.

[0110] Based on the same inventive concept, the present application also provides a method for obtaining the radiation termination time in the second aspect, which is used to obtain the radiation termination time corresponding to the current radiation plan according to the current radiation plan, the current CT / MRI image and the current blood boron concentration value, and combined with the radiation termination conditions, so as to assist in setting the radiation duration corresponding to the current radiation plan based on the radiation termination time and realize dynamic adjustment of the current radiation plan.

[0111] The current radiation scheme is a pre-built scheme for characterizing the operating state of the radiation equipment corresponding to each moment after the current moment; for example, for T i The current radiation plan corresponding to the time includes T i+1 , T i+2 ,……T i+n The operating status of the radiation equipment corresponding to each moment.

[0112] The current CT / MRI image is the latest CT / MRI image acquired before the current radiation scheme is initiated; and the current blood boron concentration value is the latest blood boron concentration value in the human body acquired before the current radiation scheme is initiated.

[0113] The radiation termination time is the time when the corresponding radiation dose in the tumor area satisfies the radiation termination condition after the radiation process is started according to the current radiation plan.

[0114] See also Figure 6 , showing a flow chart of a method for obtaining the radiation termination time provided in an embodiment of the present application; Figure 6 As shown, the method comprises the following steps:

[0115] S10, based on the current radiation scheme, constructing a radiation field model corresponding to the current radiation scheme;

[0116] Specifically, after determining the current radiation plan, the spatial information and target parameter information of the associated equipment at each moment in the radiation plan are extracted, and the spatial information and the target parameter information corresponding to each moment are respectively imported into a preset data processing platform, so that in the data processing platform, according to the spatial information and the target parameter information of the associated equipment, a radiation field intensity calculation method is used to obtain the radiation intensity value corresponding to each of the first pixel points at each moment, thereby generating a twin radiation field model corresponding to the radiation plan as the current radiation field model.

[0117] Among them, the data processing platform is a software platform built based on twin technology for realizing BNCT simulation.

[0118] S20, constructing the characteristic element concentration model based on the current CT / MRI image and the current blood boron concentration value;

[0119] Specifically, based on the current CT / MRI image, a tumor region of the patient is obtained; and according to the CT / MRI value corresponding to each second pixel point in the tumor region, an interval classification method is used to extract the concentration of the associated element corresponding to each second pixel point at each time;

[0120] And, based on the current blood boron concentration value, obtaining the boron element concentration corresponding to each of the second pixel points in the tumor area at each time;

[0121] The associated element concentration and the boron element concentration corresponding to the same moment are combined to obtain the characteristic element concentration corresponding to each moment, thereby obtaining the current characteristic element concentration model.

[0122] It should be noted that the step S10 and the step S20 may be performed simultaneously or sequentially, which is not specifically limited here.

[0123] S30, based on the radiation field model and the current characteristic element concentration model, using the BNCT radiation dose acquisition method to acquire the radiation dose corresponding to each voxel point in the tumor area at each time;

[0124] Specifically, in this embodiment, the BNCT radiation dose acquisition method is the same as the BNCT radiation dose acquisition method described in any one of the above embodiments, and will not be repeated here.

[0125] S40, according to a preset termination condition, sequentially performing detection based on the termination condition on the radiation dose corresponding to each voxel point at each moment, and taking the moment when the termination condition is satisfied as the radiation termination moment.

[0126] The termination time is the time when the corresponding radiation dose meets the termination condition.

[0127] In the present application, the treatment termination condition is a determination condition related to the radiation dose, which is used to determine whether the radiation dose meets the preset radiation level requirement.

[0128] In a specific implementation, the radiation dose corresponding to each voxel point at each moment is sequentially detected based on the termination condition, such as Figure 7 As shown, including:

[0129] S41, obtaining the radiation dose corresponding to each voxel point at the current moment;

[0130] S42, detecting whether the radiation dose corresponding to each voxel point satisfies the termination condition;

[0131] S43a, if yes, determining that the radiation dose at the current moment satisfies the termination condition;

[0132] S43b, if not, it is determined that the radiation dose at the current moment does not meet the termination condition, and the above steps are re-executed at the next moment.

[0133] It can be understood that the termination condition may include a variety of radiation dose determination conditions; for example, the termination condition includes that the radiation dose is less than the maximum radiation dose; for example, the termination condition is that the total radiation dose corresponding to each voxel point is not greater than the design value of the radiation dose, which will not be elaborated here.

[0134] The method for acquiring the radiation termination time provided in the embodiment of the present application constructs a radiation field model corresponding to the current radiation plan and a characteristic element concentration model corresponding to the current CT / MRI image and the current blood boron concentration value, respectively, and obtains the radiation dose distribution of the tumor area at each time based on the radiation field model and the characteristic element concentration model using the BNCT radiation dose acquisition method, and extracts the termination time that meets the radiation termination conditions based on the radiation dose distribution corresponding to each time, so that the radiation termination time can be quickly and accurately acquired to assist in adjusting / formulating the radiation plan.

[0135] Based on the same technical concept, the radiation dose acquisition method based on BNCT provided in the above embodiment of the present invention can be implemented on the terminal side or the server side.

[0136] See also Figure 8 , is an optional hardware structure diagram of an electronic terminal 700 provided in an embodiment of the present invention. The electronic terminal 700 may be a live broadcast machine, a camera, a mobile phone, a computer device, a tablet device, a personal digital processing device, a factory background processing device, etc. that integrates photo / video functions. The electronic terminal 700 includes: at least one processor 701, a memory 702, at least one network interface 704 and a user interface 706. The various components in the device are coupled together through a bus system 705. It can be understood that the bus system 705 is used to realize the connection and communication between these components. In addition to the data bus, the bus system 705 also includes a power bus, a control bus and a status signal bus.

[0137] The user interface 706 may include a display, a keyboard, a mouse, a trackball, a click gun, keys, buttons, a touch pad or a touch screen.

[0138] It is understood that the memory 702 can be a volatile memory or a non-volatile memory, and can also include both volatile and non-volatile memories. Among them, the non-volatile memory can be a read-only memory (ROM), a programmable read-only memory (PROM), which is used as an external cache. By way of exemplary but not limiting explanation, many forms of RAM are available, such as static random access memory (SRAM), synchronous static random access memory (SSRAM). The memory described in the embodiments of the present invention is intended to include but is not limited to these and any other suitable categories of memory.

[0139] The memory 702 in the embodiment of the present invention is used to store various types of data to support the operation of the electronic terminal 700. Examples of these data include: any executable program used to operate on the electronic terminal 700, such as an operating system 7021 and an application 7022; the operating system 7021 includes various system programs, such as a framework layer, a core library layer, a driver layer, etc., for implementing various basic services and processing hardware-based tasks. The application 7022 can include various applications, such as a media player (MediaPlayer), a browser (Browser), etc., for implementing various application services. The radiation dose acquisition method based on BNCT that implements the embodiment of the present invention can be included in the application 7022.

[0140] The method disclosed in the above embodiment of the present invention can be applied to the processor 701, or implemented by the processor 701. The processor 701 may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by the hardware integrated logic circuit in the processor 701 or the instruction in the form of software. The above processor 701 may be a general processor, a digital signal processor (DSP, Digital Signal Processor), or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. The processor 701 can implement or execute the various methods, steps and logic block diagrams disclosed in the embodiment of the present invention. The general processor 701 may be a microprocessor or any conventional processor, etc. In combination with the steps of the accessory optimization method provided in the embodiment of the present invention, it can be directly embodied as a hardware decoding processor to execute, or it can be executed by a combination of hardware and software modules in the decoding processor. The software module can be located in a storage medium, which is located in a memory, and the processor reads the information in the memory and completes the steps of the above method in combination with its hardware.

[0141] In an exemplary embodiment, the electronic terminal 700 may be implemented by one or more application specific integrated circuits (ASIC), DSP, programmable logic device (PLD), complex programmable logic device (CPLD) to execute the aforementioned method.

[0142] An embodiment of the present invention further provides a computer-readable storage medium on which a computer program is stored. When the program is called by a processor, the method for acquiring radiation dose based on BNCT provided by the present invention is implemented.

[0143] Among them, the computer-readable storage medium can be a tangible device that can hold and store instructions used by the instruction execution device. The computer-readable storage medium can be, for example, (but not limited to) an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the above. More specific examples of computer-readable storage media (a non-exhaustive list) include: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a static random access memory (SRAM), a portable compact disk read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, and a mechanical encoding device.

[0144] The computer-readable program described herein can be downloaded from a computer-readable storage medium to each computing / processing device, or downloaded to an external computer or external storage device via a network, such as the Internet, a local area network, a wide area network, and / or a wireless network. A network adapter card or a network interface in each computing / processing device receives computer-readable program instructions from the network and forwards the computer-readable program instructions for storage in a computer-readable storage medium in each computing / processing device.

[0145] It should be noted that, in the various embodiments of the present application, the serial numbers of the above steps do not indicate the order of execution. The execution order of the steps should be determined by their functions and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0146] The above embodiments are merely illustrative of the principles and effects of the present invention, and are not intended to limit the present invention. Anyone familiar with the art may modify or alter the above embodiments without departing from the spirit and scope of the present invention. Therefore, all equivalent modifications or alterations made by a person of ordinary skill in the art without departing from the spirit and technical concept disclosed by the present invention shall still be covered by the claims of the present invention.

Claims

1. A BNCT radiation dose acquisition method, characterized in that: include: Acquire a radiation field model and a characteristic element concentration model; wherein the radiation field model is a model for characterizing the radiation intensity corresponding to each first pixel point in the tumor region at each moment, and the characteristic element concentration model is a model for characterizing the associated element concentration and the boron element concentration corresponding to each second pixel point in the tumor region at each moment; the method for acquiring the boron element concentration corresponding to each second pixel point in the tumor region at each moment includes: acquiring the boron element concentration based on the blood boron concentration value in the human body at each moment using the diffusion formula of the boron element; performing spatiotemporal alignment on the radiation field model and the characteristic element concentration model to obtain a BNCT comprehensive model; Based on the BNCT comprehensive model, the radiation field intensity and characteristic element concentration corresponding to each voxel point in the tumor area at each time are obtained, so as to obtain the radiation dose of the voxel point at the corresponding time based on the radiation field intensity and the characteristic element concentration; The associated element is a basic element of the human body associated with the BNCT radiation dose, and the first pixel point and the second pixel point correspond to a first sampling resolution and a second sampling resolution, respectively.

2. The BNCT radiation dose acquisition method according to claim 1, characterized in that: The implementation method of performing spatiotemporal alignment on the radiation field model and the characteristic element concentration model includes: Based on a preset target resolution, construct a spatial grid of the tumor region; and use grid points in the spatial grid as voxel points; Extract the radiation field intensity distribution and characteristic element concentration distribution at the same time; Based on the target resolution, the radiation field intensity distribution and the characteristic element concentration distribution at the same time are resampled respectively; the resampled radiation field intensity distribution and the characteristic element concentration distribution are spatially aligned with the spatial grid respectively to obtain the radiation field intensity and the characteristic element concentration corresponding to each voxel point.

3. The BNCT radiation dose acquisition method according to claim 1, characterized in that: The method for obtaining the radiation field model comprises: Determining associated equipment in a BNCT working environment, wherein the associated equipment is a radiation device associated with the radiation field intensity; Among the device parameters of the associated devices, determining a device parameter related to the radiation field intensity as a target parameter; Based on the preset radiation scheme, combined with the three-dimensional model information of the associated device, the spatial information and target parameter information of the associated device at each moment are obtained; Based on the spatial information and the target parameter information corresponding to each moment, the radiation field model is generated by using the twin technology; Among them, the radiation plan is a pre-set operation plan of the associated device; the radiation field model is a model that uses a radiation field intensity calculation method to obtain the radiation intensity value corresponding to each first pixel point at each moment based on the spatial information and target parameter information of the associated device.

4. The BNCT radiation dose acquisition method according to claim 1, characterized in that: The method for obtaining the characteristic element concentration model includes: acquiring the tumor region; Obtain the associated element concentration corresponding to each second pixel point in the tumor area at each time, and obtain the boron element concentration corresponding to each second pixel point in the tumor area at each time; combine the associated element concentration and the boron element concentration corresponding to the same time to obtain the characteristic element concentration corresponding to each time; wherein, The method for obtaining the concentration of the associated element corresponding to each of the second pixel points in the tumor area at each time includes: According to the CT / MRI value corresponding to each second pixel point in the tumor area, the concentration of the associated element corresponding to each second pixel point at each time is extracted.

5. The BNCT radiation dose acquisition method according to claim 4, characterized in that: The method of extracting the concentration of the associated element corresponding to each second pixel point at each time according to the CT / MRI value corresponding to each second pixel point in the tumor area includes: According to the conversion relationship between the CT / MRI value and the tissue density, the tissue density corresponding to each second pixel point is obtained; Using an interval classification method, combined with the tissue density corresponding to each second pixel point, each second pixel point is classified into each type of tissue; Based on the tissue type corresponding to each second pixel point, the concentration of the associated element corresponding to each second pixel point is determined.

6. The BNCT radiation dose acquisition method according to claim 4, characterized in that: The method for obtaining the boron element concentration corresponding to each second pixel point in the tumor area at each time includes: Calculating the distance function from each second pixel point in the tumor area to the tumor boundary by using a signed distance field method; Determine the spatial position corresponding to each second pixel point based on the distance function; Based on the spatial position information corresponding to each second pixel point and based on the blood boron concentration value in the human body at each moment, the boron element concentration value corresponding to each voxel point at each moment is obtained using the diffusion equation.

7. The BNCT radiation dose acquisition method according to claim 6, characterized in that: The method of obtaining the boron concentration value corresponding to each voxel point at each moment by using the diffusion equation based on the spatial position information corresponding to each second pixel point and the blood boron concentration value in the human body at each moment includes: R(x,y,z,t)=-k·I(x,y,z)·C in, is the rate of change of C value; C is the blood boron concentration at each moment; D is the diffusion coefficient; R(x, y, z, t) is the R value corresponding to each second pixel point at each moment; R is the reaction coefficient of neutron radiation dose to boron concentration; k is the reaction coefficient; I(x, y, z) is the I value corresponding to each second pixel point; I is the neutron radiation dose distribution; (x, y, z) is the spatial position information of each second pixel point.

8. A method for obtaining the radiation termination time, characterized in that: include: Based on the current radiation plan, construct a radiation field model corresponding to the current radiation plan; Based on the current CT / MRI image and the current blood boron concentration value, construct the characteristic element concentration model; Based on the radiation field model and the characteristic element concentration model, the radiation dose corresponding to each voxel point in the tumor area at each time is obtained using the BNCT radiation dose acquisition method according to any one of claims 1 to 7; the radiation dose is a simulated value; According to a preset termination condition, the radiation dose corresponding to each voxel point at each moment is sequentially detected based on the termination condition, and the moment when the termination condition is satisfied is used as the radiation termination moment; Among them, the current CT / MRI image is the latest CT / MRI image obtained before the current radiation plan is started; and the current blood boron concentration value is the latest blood boron concentration value in the human body obtained before the current radiation plan is started.

9. The method for obtaining the radiation termination time according to claim 8, characterized in that: The implementation method of sequentially performing detection based on the termination condition on the radiation dose corresponding to each voxel point at each moment includes: Obtain the radiation dose corresponding to each voxel point at the current moment; Detecting whether the radiation dose corresponding to each voxel point satisfies the termination condition; If yes, it is determined that the radiation dose at the current moment meets the termination condition; if no, it is determined that the radiation dose at the current moment does not meet the termination condition, and the detection process is re-executed at the next moment.

10. A terminal, characterized in that: include: A processor and a memory, wherein the memory is communicatively connected to the processor; The memory is used to store a computer program, and the processor is used to execute the computer program stored in the memory, so that the terminal executes the BNCT radiation dose acquisition method according to any one of claims 1 to 7 or executes the radiation suspension time acquisition method according to claim 8 or 9.

11. A computer storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the computer program implements the BNCT radiation dose acquisition method according to any one of claims 1 to 7 or executes the radiation suspension time acquisition method according to claim 8 or 9.

Images