Dose calculation method and device, medium and product

Abstract

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

Description

Technical Field

[0001] The present application relates to the field of dose calculation technology. Specifically, the present application relates to a dose calculation method, device, medium and product. Background Art

[0002] Radiation therapy is currently one of the main treatments for malignant tumors and is widely used in the treatment of cancer patients. Dose calculation is the core of radiotherapy planning, and the speed and accuracy of dose calculation have a significant impact on the efficiency and quality of radiotherapy planning.

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

[0004] The present invention provides a dose calculation method, device, medium, and product that can solve the problems of long calculation time, large amount of calculation, and high computing resource consumption in dose calculation methods. To achieve this goal, the present invention provides the following solutions.

[0005] According to one aspect of an embodiment of the present application, a dose calculation method is provided, comprising:

[0006] Determining a dose calculation point and a cartridge string corresponding to the dose calculation point;

[0007] 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 segment of the tube string, wherein the parameters of the convolution kernel are obtained by fitting the dose data calculated based on the Monte Carlo method;

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

[0009] In one possible implementation, determining a dose calculation point and a cartridge string corresponding to the dose calculation point includes:

[0010] determining a 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;

[0011] The dose contribution direction of the dose contribution point corresponding to the dose calculation point to be calculated is obtained, and the tube string corresponding to the dose calculation point is determined according to 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] Acquiring dose data corresponding to the dose calculation space, and fitting the parameters of the convolution kernel based on the particle swarm algorithm and the dose data;

[0014] An analytical expression of the convolution kernel is determined according to the parameters.

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

[0016] k(r)=Ae -ar

[0017] Where 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 calculation formula for the first dose contribution of the line segment uv in the tube string to the dose calculation point is:

[0019]

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

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

[0022] The absorbed dose of each dose calculation point in the cartridge string is calculated based on the first dose contribution. The calculation formula for the absorbed dose is:

[0023] D i =E i +S i

[0024] Where D i is the dose absorbed by the i-th voxel in the line segment uv, E i represents the dose contribution of the i-th voxel to itself, L is the length of the path that the ray passes through the i-th voxel, S i represents the dose contribution of other voxels in the barrel string to the i-th voxel, It represents the effective distance between the i-1th voxel and the i-th voxel.

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

[0026] A first dose contribution of each line segment is obtained according to the convolution kernel, and a total dose contribution of each cartridge string to the dose calculation point is calculated using the first dose contribution.

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

[0028] Acquiring dose data of a plurality of dose calculation points in the three-dimensional grid;

[0029] The dose distribution in the dose calculation space is calculated according to the dose data of the plurality of dose calculation points.

[0030] According to one aspect of an embodiment of the present application, an embodiment of the present application provides an electronic device, 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 above methods.

[0031] According to one aspect of an embodiment of the present application, an embodiment of the present application provides a computer-readable storage medium having a computer program stored thereon, and the computer program implements the steps of the above-described method when executed.

[0032] According to one aspect of an embodiment of the present application, an embodiment of the present application provides a computer program product, including a computer program, which implements the steps of the above method when executed by a processor.

[0033] The beneficial effects of the technical solution provided by the embodiments of the present application are:

[0034] The dose calculation method provided in the present application determines the dose calculation point and the barrel string corresponding to the dose calculation point; calculates the total dose contribution of the barrel string to the dose calculation point based on the convolution kernel corresponding to the barrel string and the line segment of the barrel string. The parameters of the convolution kernel are obtained by fitting the dose data calculated based on the Monte Carlo method, and the dose data of each dose calculation point is calculated based on the total dose contribution. The embodiment of the present application directly uses the convolution kernel and the line segment of the barrel string to calculate the dose contribution of each barrel 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. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following is a brief introduction to the drawings required for describing the embodiments of the present application.

[0036] Figure 1A flow chart of a dosage calculation method provided in an embodiment of the present application;

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

[0038] Figure 3 This is a structural diagram of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0039] The following describes the embodiments of the present application in conjunction with the accompanying drawings. It should be understood that the embodiments described below in conjunction with the accompanying drawings are exemplary descriptions for explaining the technical solutions of the embodiments of the present application and do not constitute a limitation on the technical solutions of the embodiments of the present application.

[0040] Those skilled in the art will understand that, unless otherwise stated, the singular forms "a," "an," "said," and "the" used herein may also include plural forms. It should be further understood that the terms "including" and "comprising" used in the embodiments of the present application mean that the corresponding features can be implemented as the presented features, information, data, steps, operations, elements, and / or components, but do not exclude implementation as other features, information, data, steps, operations, elements, components, and / or combinations thereof supported by the present technical field. It should be understood that when we refer to an element as being "connected" or "coupled" to another element, the element can be directly connected or coupled to the other element, or it can refer to the element and the other element establishing a connection relationship through an intermediate element. In addition, the "connection" or "coupling" used herein may include wireless connection or wireless coupling. The term "and / or" 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] In order to make the objectives, technical solutions and advantages of this application clearer, the implementation methods of this application will be further described in detail below with reference to the accompanying drawings.

[0042] The following describes several exemplary embodiments to illustrate the technical solutions of the embodiments of the present application and the technical effects produced by the technical solutions of the present application. It should be noted that the following embodiments can refer to, draw on, or combine with each other, and the same terms, similar features, and similar implementation steps in different embodiments will not be repeated.

[0043] The dose calculation method, device, 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 executes the dose calculation method of the present application may be a mobile phone, a tablet computer, a server, or other smart terminal capable of performing dose calculation.

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

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

[0047] Optionally, determining the dose calculation point and the tube string corresponding to the dose calculation point 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 according to the dose contribution direction, the dose contribution point being the dose calculation point that sends energy to the dose calculation point.

[0048] Optionally, range information of a region for which dose distribution is to be calculated may be acquired, and a space corresponding to the range information may be determined as a dose calculation space.

[0049] Optionally, grid points in a three-dimensional grid can be identified as dose calculation points. Each dose calculation point receives dose contributions from other dose calculation points in the three-dimensional grid. The direction of the dose calculation points contributing to it can be determined based on the position of the dose calculation point. The corresponding barrel string can then be derived based on this direction, and the barrel string can be used to represent the dose contribution received in a particular direction. The dose calculation points in each barrel string are located in the same direction as the dose calculation point to be calculated.

[0050] Optionally, the dose calculation points to be calculated may 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 may be determined based on 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 the dose data calculated based on the Monte Carlo method. The convolution kernel is used to represent the probability distribution of the dose contribution of a voxel to other voxels in the three-dimensional grid.

[0053] Optionally, dose data corresponding to the dose calculation space is obtained, and parameters of the convolution kernel are fitted based on the particle swarm algorithm and the dose data; and an analytical expression of the convolution kernel is determined according to the parameters.

[0054] In one embodiment, Figure 2As shown, there are voxels P and Q, and an incident light ray (which may be a therapeutic ray) passes through voxel P. There is a deviation angle between the direction of the dose contribution from voxel P to voxel Q and the direction of the light ray. After determining the dose calculation point to be calculated, the Monte Carlo method can be used to calculate the dose data for multiple voxels corresponding to the dose calculation point. Based on this dose data, the particle swarm algorithm is used to fit the parameters at the current deviation angle, and the analytical expression of the convolution kernel is derived based on these parameters.

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

[0056] k(r)=Ae -ar

[0057] Where 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 of the convolution kernel can also be k(r)=Ae -ar +Be -br , B and b are also the parameters of the convolution kernel. These parameters can also be fitted using the particle swarm 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 scattered 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 tube string corresponding to the dose calculation point is determined as an energy transmission line, and any point on the line will be the dose calculation point to be calculated. Any point on the line will have a dose contribution to the dose calculation point, so the integral of all dose contributions is the total dose contribution of the tube string to the dose calculation point. When calculating the total dose contribution, the energy transmission line segment can be used. For a certain line segment uv, its first dose contribution D to the dose calculation point Q uv for

[0060]

[0061] Where x is the coordinate of the energy transfer line, the first dose contribution is calculated by linearly varying the TERMA (total energy released per unit mass of material) within the line segment and the density is a constant. -ar Substituting into the calculation formula of the first dose contribution, we can obtain the calculation formula of the first dose contribution of the line segment uv in the tube string to the dose calculation point:

[0062]

[0063] Where r u is the effective distance from the line segment uv in the cartridge string to the dose calculation point Q. This effective distance can be the equivalent distance from point Q to the endpoint u of the line segment uv. This distance can be equivalent to the distance from point Q to the endpoint u under water density. T u is the TERMA value of endpoint u in line segment uv, T v is the TERMA value of endpoint v in line segment uv, ρ u is the density value of endpoint u in line segment uv, ρ v 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 may be calculated by a recursive method.

[0065] Optionally, the ray advances along the energy transmission line. When passing through each voxel on the energy transmission line, the ray accumulates the dose contributions of all voxels along the path and gives the dose contributions of these dose contributions to the current voxel, while attenuating each time it advances. Therefore, the dose calculation method of the present application also includes: calculating the absorbed dose of each dose calculation point in the barrel string based on the first dose contribution (i.e., the dose contribution of other voxels on the energy transmission line to the dose calculation point). The calculation formula for the absorbed dose is:

[0066] D i =E i +S i

[0067] Where D i is the dose absorbed by the i-th voxel in the line segment uv, E i represents the dose contribution of the i-th voxel to itself, L is the length of the path that the ray passes through the i-th voxel (the length of the path that the ray corresponding to the line segment uv passes through along the energy transmission line to reach the i-th voxel), E i The first dose contribution of voxel s on the energy transfer line is equal to e -aru The ratio (i.e. D 1uv represents the first dose contribution of line segment uv to voxel s), and the effective distance from line segment uv to voxel s is S i represents the dose contribution of other voxels in the barrel string to the i-th voxel, It represents the effective distance between the i-1th voxel and the i-th voxel.

[0068] Optionally, the total dose contribution of the barrel string to the dose calculation point is calculated based on the convolution kernel corresponding to the barrel string and the line segments of the barrel string, including: obtaining the first dose contribution of each line segment based on the convolution kernel, calculating the total dose contribution of each barrel string to the dose calculation point using the first dose contribution; and determining the dose data for the dose calculation point based on the sum of the total dose contributions. Summing the first contributions of all line segments yields the dose contribution of the current barrel string to the dose calculation point, eliminating the need for discretization of the convolution kernel. This avoids the problem of a large gradient in the convolution kernel when r→0, which can easily lead to an increased difference between the calculated value and the actual value, thereby reducing calculation accuracy.

[0069] In one embodiment, the energy transmission line corresponding to the cartridge string is split into multiple line segments, and the dose contribution of each line segment uv to the current point is D uv , so using the analytical expression given in the mathematical model, then Among them D i (u, v) represents the dose contribution of a line segment on the current energy transmission line to the i-th voxel. Specifically, the length of the line segment can be set to 3 mm based on computational accuracy and speed requirements. Alternatively, the length can be set to 4 mm, 2.5 mm, or other values based on practical requirements. This is not a limitation.

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

[0071] Optionally, after obtaining the total dose contribution of each cartridge 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 Among them D i is the dose contribution of the i-th cartridge string to the point.

[0072] Optionally, the method further includes: acquiring dose data of multiple dose calculation points in a three-dimensional grid; and calculating the dose distribution in the dose calculation space according to the dose data of the multiple dose calculation points.

[0073] Optionally, the multiple dose calculation points may be determined as a group of dose results, and interpolation processing may be performed based on the dose results to obtain dose data at different positions in the dose calculation space.

[0074] The dose calculation method provided by this application does not make constant or linear assumptions for the convolution kernel within the three-dimensional grid, and directly uses analytical expressions to derive the dose contribution of the dose calculation point (the dose is directly obtained using the analytical expression of the convolution kernel without making assumptions such as linearity). This avoids the problem of excessive gradient of the convolution kernel near zero, thereby providing higher calculation accuracy. In addition, a set of dose data is obtained by Monte Carlo simulation, and the analytical expression of the convolution kernel is obtained by fitting this set of data. This set of analytical expressions is then used to calculate the dose. Since this application maintains the original analytical expression without making further assumptions, no new errors are introduced in the dose calculation step, effectively improving the calculation accuracy.

[0075] Based on the derived dose calculation formula, the present invention can directly calculate the dose of scattered points in space without the help of a three-dimensional grid, thereby improving the calculation speed of the scattered point dose calculation task.

[0076] Based on the same inventive concept, the embodiment of the present application provides an electronic device, such as Figure 3 As shown, Figure 3 The electronic device 2000 shown 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 can be a CPU (Central Processing Unit), a general-purpose processor, a DSP (Digital Signal Processor), an ASIC (Application Specific Integrated Circuit), an 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 can also be a combination that implements computing functions, such as a combination of one or more microprocessors, a combination of a DSP and a microprocessor, etc.

[0078] Bus 2002 may include a path 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, for example. Bus 2002 may be divided into an address bus, a data bus, a control bus, and so on. For ease of illustration, the figure shows only one thick line, but this does not indicate that there is only one bus or only one type of bus.

[0079] The memory 2003 may be a ROM (Read-Only Memory) or other types of static storage devices that can store static information and instructions, a RAM (random access memory) or other types of dynamic storage devices that can store information and instructions, or an EEPROM (Electrically Erasable Programmable Read Only Memory), a CD-ROM (Compact Disc Read-Only Memory) or other optical disk storage, optical disk storage (including compact discs, laser discs, optical discs, digital versatile discs, Blu-ray discs, etc.), a magnetic disk storage medium or other magnetic storage device, or any other medium that can be used to carry or store desired program code in the form of instructions or data structures and can be accessed by a computer, but is not limited to these.

[0080] Optionally, the electronic device 2000 may further include a communication unit 2004. The communication unit 2004 may be used to receive and transmit signals. The communication unit 2004 may allow the electronic device 2000 to communicate with other devices wirelessly or by wire to exchange data. It should be noted that in actual applications, the number of communication units 2004 is not limited to one.

[0081] Optionally, the electronic device 2000 may further include an input unit 2005. The input unit 2005 may be configured to receive input digital, character, image, and / or sound information, or to generate key signal input 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 a touch screen, a physical keyboard, function keys (such as a volume control key, a power key, etc.), a trackball, a mouse, a joystick, a camera, a microphone, and the like.

[0082] Optionally, the electronic device 2000 may further include an output unit 2006. The output unit 2006 may 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, and the like.

[0083] Although the electronic device 2000 is shown with various devices, it should be understood that it is not required to implement or possess all the devices shown, and more or fewer devices may be implemented or possessed instead.

[0084] Optionally, the memory 2003 is used to store a computer program for executing the solution of the present application, and the 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 the present application.

[0085] Based on the same inventive concept, an embodiment of the present application provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by an electronic device / processor, it implements the steps of any method provided in the present application / implements the steps of various optional implementation methods of the method provided in the present application.

[0086] Based on the same inventive concept, an embodiment of the present application provides a computer program product, which includes a computer program, which, when executed by an electronic device / processor, implements the steps of any method provided in the present application / implements the steps of various optional implementation methods of the method provided in the present application.

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

[0088] In the description of the present application, the directions or positional relationships indicated by words such as "center", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", and "outside" are based on the exemplary directions or positional relationships shown in the accompanying drawings. They are for the convenience of describing or simplifying the description of the embodiments of the present application, and do not 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 understood as limitations on the present application.

[0089] The terms "first," "second," "third," "fourth," "1," "2," and the like (if any) in the specification and claims of this application and the accompanying drawings are used to distinguish similar objects and are not necessarily used to describe a particular order or sequential sequence. It should be understood that the terms used in this manner are interchangeable where appropriate, so that the embodiments of the application described herein can be implemented in an order other than that shown or described in the drawings.

[0090] It should be understood that, although each operation step is indicated by arrows in the flowchart of the embodiment of the present application, the order of implementation of these steps is not limited to the order indicated by the arrows. Unless otherwise clearly stated herein, in some implementation scenarios of the embodiment of the present application, the implementation steps in each flowchart can be performed in other orders according to demand. In addition, some or all of the steps in each flowchart can include multiple sub-steps or multiple stages based on actual implementation scenarios. Some or all of these sub-steps or stages can be executed at the same time, and each sub-step or stage in these sub-steps or stages can also be executed at different times respectively. Under different scenarios at the execution time, the execution order of these sub-steps or stages can be flexibly configured according to demand, and the embodiment of the present application does not limit this.

[0091] The above is only an optional implementation method for some implementation scenarios of this application. It should be pointed out that for ordinary technicians in this technical field, without departing from the technical concept of the solution of this application, the use of other similar implementation methods based on the technical ideas of this application also falls within the protection scope of the embodiments of this application.

Claims

1. A dosage calculation method, characterized in that: include: Determining a dose calculation point and a cartridge string corresponding to the dose calculation point; 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 segment of the tube string, wherein the parameters of the convolution kernel are obtained by fitting the dose data calculated based on the Monte Carlo method; Dose data for each dose calculation point is calculated based on the total dose contribution.

2. The dose calculation method according to claim 1, characterized in that: The step of determining a dose calculation point and a cartridge string corresponding to the dose calculation point includes: determining a 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; The dose contribution direction of the dose contribution point corresponding to the dose calculation point to be calculated is obtained, and the tube string corresponding to the dose calculation point is determined according to the dose contribution direction. The dose contribution point is the dose calculation point that sends energy to the dose calculation point to be calculated.

3. The dose calculation method according to claim 2, characterized in that: Determining the convolution kernel includes: Acquiring dose data corresponding to the dose calculation space, and fitting the parameters of the convolution kernel based on the particle swarm algorithm and the dose data; An analytical expression of the convolution kernel is determined according to the parameters.

4. The dose calculation method according to claim 3, characterized in that: The analytical formula of the convolution kernel is: k(r)=Ae -ar Where 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; The calculation formula for the first dose contribution of the line segment uv in the tube string to the dose calculation point is: Where r u is the effective distance from the line segment uv in the tube string to the dose calculation point Q, T u is the TERMA value of endpoint u in line segment uv, T v is the TERMA value of endpoint v in line segment uv, ρ u is the density value of endpoint u in line segment uv, ρ v is the density value of endpoint v in line segment uv.

5. The dose calculation method according to claim 4, characterized in that: The method comprises: The absorbed dose of each dose calculation point in the cartridge string is calculated based on the first dose contribution. The calculation formula for the absorbed dose is: D i =E i +S i Where D i is the dose absorbed by the i-th voxel in the line segment uv, E i represents the dose contribution of the i-th voxel to itself, L is the length of the path that the ray passes through the i-th voxel, S i represents the dose contribution of other voxels in the barrel string to the i-th voxel, r(x i-1 ,x i ) represents the effective distance between the i-1th voxel and the i-th voxel.

6. The dose calculation method according to claim 4, characterized in that: The calculating the total dose contribution of the tube string to the dose calculation point according to the convolution kernel corresponding to the tube string and the line segment of the tube string includes: A first dose contribution of each line segment is obtained according to the convolution kernel, and a total dose contribution of each cartridge string to the dose calculation point is calculated using the first dose contribution.

7. The dose calculation method according to claim 6, characterized in that: The method comprises: Acquiring dose data of a plurality of dose calculation points in the three-dimensional grid; The dose distribution in the dose calculation space is calculated according to the dose data of the plurality of dose calculation points.

8. 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 according to any one of claims 1 to 7.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images