Dose calculation methods, devices, media, and products

Abstract

The application provides a dose calculation method, device, medium and product, and relates to the technical field of dose calculation. The method comprises the following steps: determining a dose calculation point and a cylinder string corresponding to the dose calculation point; calculating total dose contribution of the cylinder string to the dose calculation point according to a convolution kernel corresponding to the cylinder string and a line segment of the cylinder string, wherein parameters of the convolution kernel are obtained by fitting dose data calculated based on a Monte Carlo method; and calculating dose data of each dose calculation point based on the total dose contribution. In the application, the dose contribution of each cylinder string to the dose calculation point is directly calculated by using the convolution kernel and the line segment of the cylinder string, so that the introduction of errors can be effectively reduced, the calculation amount can be effectively reduced, the calculation precision can be improved, the time and the calculation resources required for dose calculation can be reduced, and the demand of clinical application can be effectively met.

Description

Technical Field

[0001] This application relates to the field of dose calculation technology, and more specifically, to a dose calculation method, device, medium, and product. Background Technology

[0002] Radiotherapy is one of the main methods for treating malignant tumors and is widely used in the treatment of cancer patients. Dosage calculation is the core of radiotherapy planning, and the speed and accuracy of dosage calculation have a significant impact on the efficiency and quality of radiotherapy plan development.

[0003] Currently, the most accurate and reliable method for dose calculation is the Monte Carlo method. However, the Monte Carlo method has a long calculation time, requires a large amount of computation and computational resources, making it difficult to effectively meet the needs of clinical applications. Summary of the Invention

[0004] This application provides a dose calculation method, device, medium, and product, which can solve the problems of long calculation time, large calculation volume, and high consumption of computing resources in dose calculation methods. To achieve this objective, this application provides the following solutions.

[0005] According to one aspect of the embodiments of this application, a dose calculation method is provided, including:

[0006] Determine the dose calculation point and the corresponding string of tubes for the dose calculation point;

[0007] The total dose contribution of the tube string to the dose calculation point is calculated based on the convolution kernel corresponding to the tube string and the line segments of the tube string. The parameters of the convolution kernel are obtained by fitting the dose data calculated based on the Monte Carlo method.

[0008] The dose data for each dose calculation point is calculated based on the total dose contribution.

[0009] In one possible implementation, determining the dose calculation point and the corresponding string of tubes includes:

[0010] A dose calculation space is determined, and the dose calculation space is discretized into a three-dimensional grid. Each voxel in the three-dimensional grid is determined as a dose calculation point.

[0011] Obtain the dose contribution direction of the dose contribution point corresponding to the dose calculation point to be calculated, and determine the tube string corresponding to the dose calculation point based on the dose contribution direction. The dose contribution point is the dose calculation point that sends energy to the dose calculation point to be calculated.

[0012] In one possible implementation, determining the convolution kernel includes:

[0013] Obtain dose data corresponding to the dose calculation space, and fit the parameters of the convolution kernel based on the particle swarm optimization algorithm and the dose data;

[0014] The analytical expression of the convolution kernel is determined based on the parameters.

[0015] In one possible implementation, the analytical expression of the convolution kernel is:

[0016] k(r)=Ae -ar

[0017] In the formula, k(r) is the convolution kernel, r represents the distance from the dose contribution point to the dose calculation point, and a and A are the parameters of the convolution kernel;

[0018] The formula for calculating the first dose contribution of the UV segment in the middle of the tube string to the dose calculation point is:

[0019]

[0020] In the formula, r u T is the effective distance from the line segment uv in the tube string to the dose calculation point Q. u T is the TERMA value of endpoint u in line segment uv. v It is the TERMA value of endpoint v in line segment uv, ρ u ρ is the density value of endpoint u in line segment uv. v It is the density value of endpoint v in line segment uv.

[0021] In one possible implementation, the method includes:

[0022] The absorbed dose at each dose calculation point in the tube string is calculated based on the first dose contribution, and the formula for calculating the absorbed dose is:

[0023] D i =E i +S i

[0024] In the formula, D i E represents the dose absorbed by the i-th voxel in the UV segment. i This represents the dose contribution of the i-th voxel to itself. L is the length of the path traversed by the ray through the i-th voxel. S i This represents the dose contribution of other voxels in the tube string to the i-th voxel. It represents the effective distance between the (i-1)th voxel and the ith voxel.

[0025] In one possible implementation, calculating the total dose contribution of the tube string to the dose calculation point based on the convolution kernel corresponding to the tube string and the line segments of the tube string includes:

[0026] The first dose contribution of each line segment is obtained based on the convolution kernel, and the total dose contribution of each tube string to the dose calculation point is calculated using the first dose contribution.

[0027] In one possible implementation, the method includes:

[0028] Obtain dose data from multiple dose calculation points in the three-dimensional grid;

[0029] The dose distribution of the dose calculation space is calculated based on the dose data from the multiple dose calculation points.

[0030] According to one aspect of the embodiments of this application, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of any of the methods described above.

[0031] According to one aspect of the embodiments of this application, the embodiments of this application provide a computer-readable storage medium having a computer program stored thereon, which, when executed, implements the steps of the method described above.

[0032] According to one aspect of the embodiments of this application, the embodiments of this application provide a computer program product, including a computer program that, when executed by a processor, implements the steps of the method described above.

[0033] The beneficial effects of the technical solutions provided in this application are:

[0034] The dose calculation method provided in this application determines the dose calculation point and the corresponding tube string. Based on the convolution kernel corresponding to the tube string and the line segment of the tube string, the total dose contribution of the tube string to the dose calculation point is calculated. The parameters of the convolution kernel are obtained by fitting the dose data calculated by the Monte Carlo method. The dose data of each dose calculation point is calculated based on the total dose contribution. The embodiments of this application directly use the convolution kernel and the line segment of the tube string to calculate the dose contribution of each tube string to the dose calculation point, thereby effectively reducing the amount of calculation and the introduction of errors, improving the calculation accuracy, reducing the time and computing resources required for dose calculation, and effectively meeting the needs of clinical applications. Attached Figure Description

[0035] To more clearly illustrate the technical solutions in the embodiments of this application, the accompanying drawings used in the description of the embodiments of this application will be briefly introduced below.

[0036] Figure 1A flowchart of the dose calculation method provided in the embodiments of this application;

[0037] Figure 2 A schematic diagram illustrating convolution kernel calculation provided in an embodiment of this application;

[0038] Figure 3 This is a structural diagram of an electronic device provided in an embodiment of this application. Detailed Implementation

[0039] The embodiments of this application are described below with reference to the accompanying drawings. It should be understood that the embodiments described below with reference to the accompanying drawings are exemplary descriptions for explaining the technical solutions of the embodiments of this application, and do not constitute a limitation on the technical solutions of the embodiments of this application.

[0040] Those skilled in the art will understand that, unless specifically stated otherwise, the singular forms “a,” “an,” “the,” and “the” used herein may also include the plural forms. It should be further understood that the terms “comprising” and “including” as used in embodiments of this application mean that the corresponding feature can be implemented as the presented feature, information, data, step, operation, element, and / or component, but do not exclude implementation as other features, information, data, step, operation, element, component, and / or combinations thereof supported by the art. It should be understood that when we say that an element is “connected” or “coupled” to another element, the one element can be directly connected or coupled to the other element, or it can mean that the one element and the other element establish a connection relationship through an intermediate element. Furthermore, “connected” or “coupled” as used herein can include wireless connection or wireless coupling. The term “and / or” as used herein indicates at least one of the items defined by the term; for example, “A and / or B” indicates implementation as “A,” or implementation as “A,” or implementation as “A and B.”

[0041] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0042] The technical solutions of this application and their effects are described below through several exemplary embodiments. It should be noted that the following embodiments can be referenced, borrowed from, or combined with each other. Identical terms, similar features, and similar implementation steps in different embodiments will not be repeated.

[0043] The dosage calculation method, equipment, medium, and product provided in this application are intended to solve at least one technical problem existing in the prior art.

[0044] Optionally, the device that performs the dose calculation method of this application can be a mobile phone, tablet computer, server, or other smart terminal capable of dose calculation.

[0045] Optionally, such as Figure 1 , Figure 2 As shown, the dosage calculation method of this application includes:

[0046] S101: Determine the dose calculation point and the corresponding tube string.

[0047] Optionally, determining the dose calculation point and the corresponding tube string includes: determining the dose calculation space, discretizing the dose calculation space into a three-dimensional grid, and determining each voxel in the three-dimensional grid as a dose calculation point; obtaining the dose contribution direction of the dose contribution point corresponding to the dose calculation point to be calculated, and determining the tube string corresponding to the dose calculation point based on the dose contribution direction, wherein the dose contribution point is the dose calculation point that sends energy to the dose calculation point.

[0048] Optionally, the range information of the region where the dose distribution to be calculated can be obtained, and the space corresponding to the range information can be determined as the dose calculation space.

[0049] Optionally, grid points in the 3D mesh can be defined as dose calculation points, each of which is affected by dose contributions from other dose calculation points in the same mesh. The direction of the dose calculation points contributing to the dose calculation point can be determined based on its position, and a corresponding string of points can be obtained based on this direction. This string represents the dose contribution received in a certain direction. The dose calculation points in each string are located in the same direction as the dose calculation point to be calculated.

[0050] Optionally, the dose calculation points to be calculated can be a set of scattered points in the dose calculation space, or multiple points on a three-dimensional grid. The dose calculation points to be calculated can be determined according to the calculation requirements of the dose distribution and the calculation method used.

[0051] S102: Calculate the total dose contribution of the tube string to the dose calculation point based on the convolution kernel corresponding to the tube string and the line segment of the tube string.

[0052] Optionally, the parameters of the convolution kernel are obtained by fitting dose data calculated using the Monte Carlo method. This convolution kernel is used to represent the probability distribution of the dose contribution of a voxel to other voxels in a 3D mesh.

[0053] Optionally, obtain the dose data corresponding to the dose calculation space, fit the parameters of the convolution kernel based on the particle swarm optimization algorithm and the dose data, and determine the analytical expression of the convolution kernel based on the parameters.

[0054] In one embodiment, such as Figure 2As shown, there are voxels P and Q, and the incident light ray (which could be a therapeutic ray) passes through voxel P. The direction of the dose contribution of voxel P to voxel Q has an angle with the direction of the light ray. After determining the dose calculation point, the dose data of multiple voxels corresponding to the dose calculation point can be calculated using the Monte Carlo method. Based on this dose data, the parameters at the current angle are obtained using a particle swarm optimization algorithm, and the analytical expression of the convolution kernel is derived based on these parameters.

[0055] Alternatively, the analytical expression of the convolution kernel is:

[0056] k(r)=Ae -ar

[0057] In the formula, k(r) is the convolution kernel, r represents the distance from the dose contribution point to the dose calculation point, a and A are the parameters of the convolution kernel, A can be the weight or intensity of energy diffusion, and a can be the attenuation coefficient.

[0058] Alternatively, the analytical expression for the convolution kernel can also be k(r) = Ae -ar +Be -br B and b are also parameters of the convolution kernel. These parameters can also be obtained by fitting using the particle swarm optimization algorithm. In this case, A can be the weight or intensity corresponding to the low-angle energy diffusion (close to the vertical incident direction), reflecting the concentration of energy in the main propagation direction; B can be the weight or intensity corresponding to the high-angle energy diffusion (deviating from the vertical direction), reflecting the distribution of energy in the lateral or scattering direction; a can be the attenuation coefficient in the low-angle direction; and b can be the attenuation coefficient in the high-angle direction.

[0059] Optionally, the string of tubes corresponding to the dose calculation point can be defined as a straight energy transfer line. Any point on this line represents the dose calculation point. Since every point on the line contributes to the dose calculation point, the integral of all dose contributions represents the total dose contribution of this string of tubes to the dose calculation point. When calculating the total dose contribution, for a specific segment of the energy transfer line (uv), its first dose contribution D to the dose calculation point Q can be calculated. uv for

[0060]

[0061] In the formula, x is the coordinate of the energy transfer line, and the first dose contribution is calculated by linearly varying the TERMA (total energy released by the ray in a unit mass of matter) within the line segment, with density being a constant. k(r) = Ae -ar Substituting these values ​​into the formula for calculating the first dose contribution, we can obtain the formula for calculating the first dose contribution of the line segment uv in the tube string to the dose calculation point:

[0062]

[0063] In the formula, r u This is the effective distance from the line segment uv in the tube string to the dose calculation point Q. This effective distance can be considered as the equivalent distance from point Q to the endpoint u of the line segment uv, which can be equivalent to the distance from point Q to endpoint u at water density. u T is the TERMA value of endpoint u in line segment uv. v It is the TERMA value of endpoint v in line segment uv, ρ u ρ is the density value of endpoint u in line segment uv. v It is the density value of endpoint v in line segment uv.

[0064] Optionally, after obtaining the dose contribution of the line segment to the dose calculation point to be calculated, the dose contribution of the line segment to other dose calculation points can also be calculated by recursion.

[0065] Optionally, the ray travels along the energy transport line. As it passes each voxel along the energy transport line, the ray accumulates the dose contribution of all voxels along its path while simultaneously providing the dose contribution of these dose contributions to the current voxel, attenuating with each forward movement. Therefore, the dose calculation method of this application further includes: calculating the absorbed dose at each dose calculation point in the tube string (i.e., the dose contribution of other voxels on the energy transport line to that dose calculation point) based on the first dose contribution. The formula for calculating the absorbed dose is:

[0066] D i =E i +S i

[0067] In the formula, D i E represents the dose absorbed by the i-th voxel in the UV segment. i This represents the dose contribution of the i-th voxel to itself. L is the length of the path traversed by the ray through the i-th voxel (the length of the path traversed by the ray corresponding to line segment uv along the energy transfer line to reach the i-th voxel). E i The first dose contribution of voxel s to e along the energy transfer line. -aru The ratio (i.e.) D 1uv (This represents the first dose contribution of line segment uv to voxel s), where the effective distance from line segment uv to voxel s is... S i This represents the dose contribution of other voxels in the tube string to the i-th voxel. It represents the effective distance between the (i-1)th voxel and the ith voxel.

[0068] Optionally, the total dose contribution of the tube string to the dose calculation point is calculated based on the convolution kernel corresponding to the tube string and the line segments of the tube string. This includes: obtaining the first dose contribution of each line segment based on the convolution kernel, calculating the total dose contribution of each tube string to the dose calculation point using the first dose contribution, and determining the dose data of the dose calculation point based on the sum of the total dose contributions. By summing the first contributions of all line segments, the dose contribution of the current tube string to the dose calculation point can be obtained. This eliminates the need for discretization of the convolution kernel, avoiding the problem that when r→0, the convolution kernel has a large gradient, which can easily lead to an increased difference between the calculated value and the actual value, reducing the calculation accuracy.

[0069] In one embodiment, the energy transfer line corresponding to the cylinder string is divided into multiple line segments, and the dose contribution of each line segment uv to the current point is D. uv Therefore, using the analytical expression given in the mathematical model, then Where D i (u,v) represents the dose contribution of a certain segment on the current energy transfer line to the i-th voxel. Specifically, the length of the segment can be determined to be 3mm according to the requirements of calculation accuracy and speed, or it can be determined to be 4mm, 2.5mm, or other values ​​according to actual requirements; no limitation is made here.

[0070] S103: Calculate dose data for each dose calculation point based on the total dose contribution.

[0071] Optionally, after obtaining the total dose contribution of each tube string to the dose calculation point, the sum of the total dose contributions is determined as the dose data of the dose calculation point. Specifically, for any given point in space, its dose... Where D i It is the dose contribution of the i-th tube string to that point.

[0072] Optionally, it also includes: acquiring dose data from multiple dose calculation points in a three-dimensional mesh; and calculating the dose distribution in the dose calculation space based on the dose data from the multiple dose calculation points.

[0073] Optionally, the multiple dose calculation points can be defined as a set of dose results, and interpolation can be performed on these dose results to obtain dose data at different locations in the dose calculation space.

[0074] The dose calculation method provided in this application does not make constant or linear assumptions about the convolution kernel within the 3D mesh. It directly derives the dose contribution at the dose calculation point using analytical expressions (the dose is directly obtained using the analytical expressions of the convolution kernel without making assumptions such as linearity), thus avoiding the problem of excessive gradients of the convolution kernel near zero and providing higher computational accuracy. Furthermore, the Monte Carlo simulation yields a set of dose data, and the analytical expressions of the convolution kernel are obtained by fitting this set of data. The dose is then calculated using these analytical expressions. Since this application maintains the original analytical expressions without making additional assumptions, no new errors are introduced in the dose calculation step, effectively improving computational accuracy.

[0075] Based on the derived analytical formula for dose calculation, this invention can directly calculate the dose of scattered points in space without the need for a three-dimensional mesh, thereby improving the calculation speed of scattered dose calculation tasks.

[0076] Based on the same inventive concept, embodiments of this application provide an electronic device, such as... Figure 3 As shown, Figure 3 The illustrated electronic device 2000 includes a processor 2001 and a memory 2003. The processor 2001 and the memory 2003 are communicatively connected, for example, via a bus 2002.

[0077] Processor 2001 may be a CPU (Central Processing Unit), a general-purpose processor, a DSP (Digital Signal Processor), an ASIC (Application Specific Integrated Circuit), a FPGA (Field-Programmable Gate Array), or other programmable logic devices, transistor logic devices, hardware components, or any combination thereof. It can implement or execute the various exemplary logic blocks, modules, and circuits described in conjunction with the disclosure of this application. Processor 2001 may also be a combination that implements computing functions, such as including one or more microprocessor combinations, a combination of a DSP and a microprocessor, etc.

[0078] Bus 2002 may include a pathway for transmitting information between the aforementioned components. Bus 2002 may be a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus, etc. Bus 2002 can be divided into address bus, data bus, control bus, etc. For ease of illustration, only one thick line is used to represent it in the figure, but this does not indicate that there is only one bus or one type of bus.

[0079] The memory 2003 may be ROM (Read-Only Memory) or other types of static storage devices capable of storing static information and instructions, RAM (Random Access Memory) or other types of dynamic storage devices capable of storing information and instructions, or EEPROM (Electrically Erasable Programmable Read Only Memory), CD-ROM (Compact Disc Read-Only Memory) or other optical disc storage, optical disc storage (including compressed optical discs, laser discs, optical discs, digital universal optical discs, Blu-ray discs, etc.), magnetic disk storage media or other magnetic storage devices, or any other medium capable of carrying or storing desired program code in the form of instructions or data structures and accessible by a computer, but not limited thereto.

[0080] Optionally, the electronic device 2000 may also include a communication unit 2004. The communication unit 2004 can be used for receiving and transmitting signals. The communication unit 2004 allows the electronic device 2000 to communicate wirelessly or wiredly with other devices to exchange data. It should be noted that in practical applications, the communication unit 2004 is not limited to one.

[0081] Optionally, the electronic device 2000 may further include an input unit 2005. The input unit 2005 can be used to receive input numbers, characters, images, and / or sound information, or to generate key signal inputs related to user settings and function control of the electronic device 2000. The input unit 2005 may include, but is not limited to, one or more of the following: a touchscreen, a physical keyboard, function keys (such as volume control buttons, power buttons, etc.), a trackball, a mouse, a joystick, a camera, a microphone, etc.

[0082] Optionally, the electronic device 2000 may also include an output unit 2006. The output unit 2006 can be used to output or display information processed by the processor 2001. The output unit 2006 may include, but is not limited to, one or more of a display device, a speaker, a vibration device, etc.

[0083] Although the electronic device 2000 with various devices is shown in the figure, it should be understood that it is not required to implement or have all of the devices shown. More or fewer devices may be implemented or have alternatively.

[0084] Optionally, the memory 2003 is used to store a computer program for executing the scheme of this application, and its execution is controlled by the processor 2001. The processor 2001 is used to execute the computer program stored in the memory 2003 to implement the steps of any method provided in the embodiments of this application.

[0085] Based on the same inventive concept, embodiments of this application provide a computer-readable storage medium storing a computer program that, when executed by an electronic device / processor, implements the steps of any method provided in this application or the steps of various optional implementations of the method provided in this application.

[0086] Based on the same inventive concept, embodiments of this application provide a computer program product, which includes a computer program that, when executed by an electronic device / processor, implements the steps of any method provided in this application or the steps of various optional implementations of the method provided in this application.

[0087] Those skilled in the art will understand that the steps, measures, and solutions in the various operations, methods, and processes discussed in this application can be alternated, modified, combined, or deleted. Furthermore, other steps, measures, and solutions in the various operations, methods, and processes discussed in this application can also be alternated, modified, rearranged, decomposed, combined, or deleted. Furthermore, steps, measures, and solutions in related technologies that are similar to those disclosed in this application can also be alternated, modified, rearranged, decomposed, combined, or deleted.

[0088] In the description of this application, the terms "center," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," and "outer," etc., indicate directions or positional relationships based on the exemplary directions or positional relationships shown in the accompanying drawings. They are used to facilitate the description or simplification of the embodiments of this application and are not intended to indicate or imply that the device or component referred to must have a specific orientation or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on this application.

[0089] The terms "first," "second," "third," "fourth," "1," "2," etc. (if present) in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of this application described herein can be implemented in a sequence other than that shown in the figures or text.

[0090] It should be understood that although arrows indicate various operation steps in the flowcharts of this application's embodiments, the order in which these steps are implemented is not limited to the order indicated by the arrows. Unless explicitly stated herein, in some implementation scenarios of this application's embodiments, the implementation steps in each flowchart can be executed in other orders as required. Furthermore, some or all steps in each flowchart, based on the actual implementation scenario, may include multiple sub-steps or multiple stages. Some or all of these sub-steps or stages can be executed at the same time, and each sub-step or stage can also be executed at different times. In scenarios where execution times differ, the execution order of these sub-steps or stages can be flexibly configured according to requirements, and this application's embodiments do not limit this.

[0091] The above description is only an optional implementation method for some implementation scenarios of this application. It should be noted that for those skilled in the art, other similar implementation methods based on the technical concept of this application without departing from the technical concept of this application also fall within the protection scope of the embodiments of this application.

Claims

1. A dosage calculation method, characterized in that, include: Determining the dose calculation point and the corresponding cylinder string includes: determining the dose calculation space, discretizing the dose calculation space into a three-dimensional grid, and determining each voxel in the three-dimensional grid as a dose calculation point; Obtain the dose contribution direction of the dose contribution point corresponding to the dose calculation point to be calculated, and determine the tube string corresponding to the dose calculation point based on the dose contribution direction. The dose contribution point is the dose calculation point that sends energy to the dose calculation point to be calculated. The total dose contribution of the tube string to the dose calculation point is calculated based on the convolution kernel corresponding to the tube string and the line segments of the tube string. The parameters of the convolution kernel are obtained by fitting the dose data calculated based on the Monte Carlo method. Calculate dose data for each dose calculation point based on the total dose contribution; The determination of the convolution kernel includes: Obtain dose data corresponding to the dose calculation space, and fit the parameters of the convolution kernel based on the particle swarm optimization algorithm and the dose data; The analytical expression of the convolution kernel is determined based on the parameters, and the analytical expression of the convolution kernel is: In the formula, Here, r represents the convolution kernel, and r represents the distance from the dose contribution point to the dose calculation point. A and A are the parameters of the convolution kernel; The formula for calculating the first dose contribution of the UV segment in the middle of the tube string to the dose calculation point is: In the formula, It is the middle segment of the tube string. To the dose calculation point Effective distance, It is a line segment The TERMA value of the midpoint u, It is a line segment The TERMA value of the midpoint v. For line segments Density value at the midpoint u, It is a line segment The density value at the midpoint v; The method includes: The absorbed dose at each dose calculation point in the tube string is calculated based on the first dose contribution, and the formula for calculating the absorbed dose is: In the formula, Let be the dose absorbed by the i-th voxel in the UV segment. This represents the dose contribution of the i-th voxel to itself. L is the length of the path traversed by the ray through the i-th voxel. , This represents the dose contribution of other voxels in the tube string to the i-th voxel. , It represents the effective distance between the (i-1)th voxel and the ith voxel.

2. The dosage calculation method according to claim 1, characterized in that, The step of calculating the total dose contribution of the tube string to the dose calculation point based on the convolution kernel corresponding to the tube string and the line segments of the tube string includes: The first dose contribution of each line segment is obtained based on the convolution kernel, and the total dose contribution of each tube string to the dose calculation point is calculated using the first dose contribution.

3. The dosage calculation method according to claim 2, characterized in that, The method includes: Obtain dose data from multiple dose calculation points in the three-dimensional grid; The dose distribution of the dose calculation space is calculated based on the dose data from the multiple dose calculation points.

4. An electronic device comprising a memory, a processor, and a computer program stored in the memory, characterized in that, The processor executes the computer program to implement the steps of the method as described in any one of claims 1-3.

5. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed, it implements the steps of the method as described in any one of claims 1-3.

6. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1-3.

Images