Beam spot processing method, system, computer device and storage medium

By constructing a beam direction-equivalent water depth grid and performing beam spot truncation, the problems of excessive memory usage and computational load in proton intensity-modulated radiotherapy are solved, and efficient beam spot processing is achieved.

CN117831715BActive Publication Date: 2026-06-26SHANGHAI UNITED IMAGING HEALTHCARE

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
SHANGHAI UNITED IMAGING HEALTHCARE
Filing Date
2022-09-29
Publication Date
2026-06-26

AI Technical Summary

Technical Problem

In intensity-modulated proton therapy, traditional beam spot truncation algorithms cannot achieve a balance between memory usage and truncation accuracy, resulting in excessive consumption of computational resources.

Method used

By constructing a beam direction-equivalent water depth grid, the cutoff region of the beam spot in the mapping grid set is determined according to the preset cutoff configuration information, and the target region is obtained according to the reference grid for beam spot dose optimization processing.

Benefits of technology

While maintaining truncation accuracy, memory usage was reduced, the computational load of the proton intensity modulation program was decreased, and beam spot processing efficiency was improved.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN117831715B_ABST
    Figure CN117831715B_ABST
Patent Text Reader

Abstract

The application relates to a beam spot processing method, a system, computer equipment and a storage medium. The method comprises the following steps: acquiring a spatial grid set generated based on a medical image of a target object, and acquiring a beam direction of at least one beam for the target object; taking a spatial grid in the beam direction of the at least one beam as a target grid; constructing a mapping grid corresponding to each target grid according to a projection plane coordinate and an equivalent water depth value corresponding to each target grid, to obtain a mapping grid set; determining a truncation region of a beam spot of the at least one beam in the mapping grid set according to preset truncation configuration information, and taking a mapping grid in the truncation region as a reference grid; and obtaining a target region of the beam spot of the at least one beam in the spatial grid set according to a target grid having a mapping relationship with the reference grid. The method can reduce memory occupation and reduce calculation amount on the basis of ensuring beam spot truncation accuracy.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This application relates to the field of computer technology, and in particular to a method, system, computer device, storage medium, and computer program product for processing spot patterns. Background Technology

[0002] With the development of radiotherapy technology, intensity-modulated proton therapy (IMRT) has emerged. Pencil beam IMRT can deliver thousands to tens of thousands of weighted proton beams within the target area. However, for such a large number of proton beams, the amount of grid dose information required to be stored in memory is substantial, typically exceeding 10GB in non-robust scenarios and exceeding 100GB in robust scenarios. This enormous memory footprint and the resulting massive computational demands mean that pencil beam IMRT planning requires significant computer resources. Traditional beam truncation algorithms cannot achieve a balance between memory usage and truncation accuracy. Summary of the Invention

[0003] Therefore, it is necessary to provide a beam spot processing method, system, computer device, storage medium, and computer program product that can solve the above-mentioned technical problems.

[0004] In a first aspect, this application provides a method for processing beam spots, the method comprising:

[0005] Obtain a set of spatial grids generated from medical images of a target object, and obtain the beam direction of at least one beam directed toward the target object. Use the spatial grids located in the beam direction of the at least one beam as the target grid.

[0006] Based on the projection plane coordinates and equivalent water depth values ​​corresponding to each target grid, a mapping grid corresponding to each target grid is constructed to obtain a mapping grid set; the mapping grid set is used to characterize the material density fluctuation in the beam direction of the at least one beam;

[0007] Based on the preset truncation configuration information, the truncation region of the beam spot of the at least one beam in the mapping grid set is determined, and the mapping grid in the truncation region is used as the reference grid.

[0008] Based on the target grid that has a mapping relationship with the reference grid, the target region of the beam spot of the at least one beam in the spatial grid set is obtained; the target region is used to participate in the beam spot dose optimization process.

[0009] In one embodiment, obtaining the beam direction of at least one beam targeting the target object, and using a spatial grid along the beam direction of the at least one beam as the target grid, includes:

[0010] Based on the beam direction of the at least one beam, determine the spatial grid through which each ray in the beam direction passes;

[0011] Based on the spatial grid through which each of the rays passes, a spatial grid in the beam direction of the at least one beam is obtained, which serves as the target grid.

[0012] In one embodiment, after the step of using a spatial grid located in the beam direction of the at least one beam as the target grid, the method further includes:

[0013] For each target mesh, determine the coordinates of the target mesh projected onto the isocenter plane, and use them as the projection plane coordinates corresponding to the target mesh;

[0014] The equivalent water depth value corresponding to the target grid is obtained based on the material density of the target grid and the ray path length passing through the target grid.

[0015] In one embodiment, determining the truncated region of the beam spot of the at least one beam in the mapping grid set according to preset truncation configuration information includes:

[0016] Obtain the energy and beam projection coordinates corresponding to the beam spot of the at least one beam;

[0017] By combining the energy and beam projection coordinates corresponding to the beam spot of the at least one beam, and the truncation configuration information, the truncation region of the beam spot of the at least one beam in the mapping grid set is obtained.

[0018] In one embodiment, prior to the step of determining the truncated region of the beam spot of the at least one beam in the mapping grid set according to preset truncation configuration information, the method further includes:

[0019] The energy of the beam spot of the at least one beam is obtained, and the overall depth dose information of the beam spot of the at least one beam is obtained;

[0020] Based on the energy and the overall depth dose information, a preset cutoff configuration for the beam spot of the at least one beam is determined.

[0021] In one embodiment, the truncation configuration information includes configuration parameters for multiple different truncation sizes. The step of determining the truncation region of the beam spot of the at least one beam in the mapping grid set according to the preset truncation configuration information, and using the mapping grid within the truncation region as a reference grid, includes:

[0022] Based on the configuration parameters for each cutoff size, the cutoff regions of different sizes for the beam spots of the at least one beam in the mapping grid set are determined respectively, and the reference grid corresponding to each size cutoff region is obtained.

[0023] In one embodiment, after the step of obtaining the target region of the beam spot of the at least one beam in the spatial grid set based on the target grid having a mapping relationship with the reference grid, the method further includes:

[0024] Based on the beam spot of the at least one beam, different optimization strategies are used to optimize the beam spot dose in target regions of different sizes in the spatial grid set.

[0025] Secondly, this application also provides a beam spot processing system, the system comprising:

[0026] The target mesh acquisition module is used to acquire a set of spatial meshes generated based on medical images of a target object, and to acquire the beam direction of at least one beam for the target object, and to take the spatial meshes located in the beam direction of the at least one beam as the target meshes;

[0027] The mapping mesh acquisition module is used to construct a mapping mesh corresponding to each target mesh based on the projection plane coordinates and equivalent water depth values ​​corresponding to each target mesh, thereby obtaining a mapping mesh set; the mapping mesh set is used to characterize the material density fluctuation in the beam direction of the at least one beam;

[0028] The truncated region determination module is used to determine the truncated region of the beam spot of the at least one beam in the mapping grid set according to the preset truncated configuration information, and to use the mapping grid in the truncated region as the reference grid.

[0029] The target region acquisition module is used to obtain the target region of the beam spot of the at least one beam in the spatial grid set based on the target grid that has a mapping relationship with the reference grid; the target region is used to participate in the beam spot dose optimization process.

[0030] Thirdly, this application also provides a computer device. The computer device includes a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement the steps of the beam spot processing method described above.

[0031] Fourthly, this application also provides a computer-readable storage medium. The computer-readable storage medium stores a computer program thereon, which, when executed by a processor, implements the steps of the beam spot processing method described above.

[0032] Fifthly, this application also provides a computer program product. The computer program product includes a computer program that, when executed by a processor, implements the steps of the beam spot processing method described above.

[0033] The aforementioned beam spot processing method, system, computer device, storage medium, and computer program product acquire a spatial grid set generated from medical images of a target object, and acquire the beam direction of at least one beam targeting the target object. The spatial grid located in the beam direction of the at least one beam is used as the target grid. Based on the projection plane coordinates and equivalent water depth values ​​corresponding to each target grid, a mapping grid corresponding to each target grid is constructed, resulting in a mapping grid set. This mapping grid set is used to characterize the material density fluctuation in the beam direction of at least one beam. Then, based on preset truncation configuration information, the truncation region of the beam spot of at least one beam in the mapping grid set is determined. The mapping grid located in the truncation region is used as the reference grid. Furthermore, based on the target grids that have a mapping relationship with the reference grids, the target region of the beam spot of at least one beam in the spatial grid set is obtained. The target region is used to participate in beam spot dose optimization processing, realizing beam spot truncation processing optimization. By introducing the equivalent water depth of the beam spot to construct the corresponding mapping grid, and setting the truncation in this mapping grid, memory usage can be reduced while ensuring truncation accuracy, thus reducing the computational load of proton intensity modulation program creation and improving beam spot processing efficiency. Attached Figure Description

[0034] Figure 1 This is a schematic flowchart of a spot treatment method in one embodiment;

[0035] Figure 2 This is a schematic diagram of a beam spot truncation process in one embodiment;

[0036] Figure 3a This is a schematic diagram of a spatial grid in one embodiment;

[0037] Figure 3b This is a schematic diagram of a cutoff region in one embodiment;

[0038] Figure 4 This is a schematic diagram of another beam spot truncation process in one embodiment;

[0039] Figure 5 This is a flowchart illustrating another beam spot processing method in one embodiment;

[0040] Figure 6 This is a structural block diagram of a beam spot processing system in one embodiment;

[0041] Figure 7This is an internal structural diagram of a computer device in one embodiment. Detailed Implementation

[0042] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0043] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for display, data used for analysis, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties; correspondingly, this application also provides a corresponding user authorization entry point for users to choose to authorize or refuse.

[0044] In one embodiment, such as Figure 1 As shown, a speckle processing method is provided. This embodiment illustrates the application of this method to a terminal. It is understood that this method can also be applied to a server, and further to a system including both a terminal and a server, and implemented through interaction between the terminal and the server. In this embodiment, the method includes the following steps:

[0045] Step 101: Obtain a set of spatial grids generated from medical images of the target object, and obtain the beam direction of at least one beam for the target object, and take the spatial grid in the beam direction of the at least one beam as the target grid.

[0046] The target object can be the organ involved in the radiotherapy. For example, medical images of the organ to be radiotreated can be obtained and used as the target object.

[0047] As an example, the beam can be a proton beam, or it can include an electron beam, a photon beam, a heavy ion beam, etc. In this embodiment, no specific limitation is made, that is, the beam spot processing method can be applied to proton beam spots, electron beam spots, photon beam spots, heavy ion beam spots, etc.

[0048] In practical applications, a set of spatial grids generated from medical images of the target object can be obtained, as well as the beam direction of at least one beam targeting the target object. Then, based on the beam direction of at least one beam, the spatial grid through which each ray passes in that beam direction can be determined. Furthermore, based on the spatial grid through which each ray passes, the spatial grid in the beam direction of at least one beam can be obtained as the target grid.

[0049] Specifically, such as Figure 2As shown, a three-dimensional spatial dose grid (i.e., a spatial grid in a set of spatial grids) can be generated based on imported patient medical images using a preset grid generation algorithm. This grid can contain all organs involved in radiotherapy. The accuracy of the spatial grid can be specified by the user in advance. The voxel corresponding to a single grid is the smallest unit for dose calculation in radiotherapy. For each beam, a series of rays can be emitted from the beam source point. The grid through which each ray passes can be used as the target grid. That is, based on the beam direction of at least one beam, the spatial grid through which each ray passes in that beam direction is determined.

[0050] Step 102: Based on the projection plane coordinates and equivalent water depth values ​​corresponding to each target grid, construct a mapping grid corresponding to each target grid to obtain a mapping grid set; the mapping grid set is used to characterize the material density fluctuation in the beam direction of the at least one beam.

[0051] As an example, the projection plane coordinates can be two-dimensional BEV (Beam Eyes View) coordinates, and the equivalent water depth value can be the WEPL (water-equivalent path length) value.

[0052] In practical implementation, for each target mesh, the coordinates of the target mesh projected onto the isocentric plane can be determined as the corresponding projection plane coordinates of the target mesh, such as two-dimensional BEV coordinates. The equivalent water depth value corresponding to the target mesh, such as the WEPL value, can be obtained based on the material density of the target mesh and the length of the ray path passing through the target mesh. Then, based on the projection plane coordinates and equivalent water depth values ​​of each target mesh, a mapping mesh corresponding to each target mesh can be constructed to obtain a mapping mesh set. There is a mapping relationship between the spatial mesh set and the mapping mesh set.

[0053] For example, for each spatial dose grid traversed by a ray, the two-dimensional coordinates of the isocentric plane of that spatial dose grid (i.e., the projection plane coordinates corresponding to the target grid) can be obtained by projecting any spatial dose grid (i.e., each target grid) onto the isocentric plane, such as the two-dimensional BEV coordinates. Furthermore, the WEPL value (i.e., equivalent water depth value) corresponding to that spatial dose grid can be obtained by summing the product of the mass density of that spatial dose grid and the path length traversing the spatial dose grid. Thus, each spatial dose grid can correspond to a set of two-dimensional BEV coordinates and WEPL values, such as... Figure 2 As shown, a beam direction-equivalent water depth grid (i.e., a mapped grid) can be constructed, such as the BEV-WEPL grid.

[0054] In one example, since the three-dimensional morphology of the beam spot in the patient's body fluctuates with the material density along the range direction, the material density fluctuation of the beam spot along the range direction can be considered in the beam spot mesh truncation. By introducing the equivalent water depth of the beam spot to describe its three-dimensional morphology, a beam direction-equivalent water depth mesh is constructed to further set the truncation in this beam direction-equivalent water depth mesh. Since the distribution of the proton beam spot remains relatively stable on this mesh, the geometric truncation can be set better, thereby reducing memory usage and computational load for proton intensity modulation program optimization while ensuring truncation accuracy.

[0055] Step 103: Based on the preset truncation configuration information, determine the truncation region of the beam spot of the at least one beam in the mapping grid set, and use the mapping grid in the truncation region as the reference grid.

[0056] In this context, truncation can refer to obtaining the effective spatial dose grid for each beam spot in the beam. For example, the spatial dose grid within the truncation can be obtained based on the mapping relationship of the BEV-WEPL grid, that is, the mapping grid in the truncation region can be used as the reference grid.

[0057] In practical applications, the energy of at least one beam spot and the overall depth-dose information of at least one beam spot are obtained. Then, based on the energy and the overall depth-dose information, a preset truncation configuration for the beam spot of at least one beam is determined. Furthermore, by combining the proton energy and beam spot projection coordinates corresponding to the beam spot of at least one beam, and the preset truncation configuration information, the truncation region of the beam spot of at least one beam in the mapping grid set can be obtained.

[0058] Specifically, taking a proton beam spot as an example, such as Figure 2 As shown, for a certain beam spot, the cutoff range of the beam spot in the BEV-WEPL grid (i.e., the cutoff region in the mapped grid set) can be obtained based on the proton energy of the beam spot and the isocenter plane coordinates of the beam spot, according to the geometric cutoff settings (i.e., cutoff configuration information) of the BEV-WEPL grid.

[0059] In one optional embodiment, the equivalent water depth (WEPL) is an important parameter of the Monte Carlo dosing engine. It can be used in dose calculation. After the Monte Carlo dosing engine completes the dose calculation, the beam direction-equivalent water depth (BEV-WEPL) can be set as a cutoff, and the grid within the BEV-WEPL cutoff can be extracted as the effective dose grid, which can then participate in subsequent dose optimization processing. Compared with the traditional geometric cutoff method, the cutoff evaluation based on the BEV-WEPL grid is more efficient in obtaining the effective grid dose, thereby reducing memory usage and reducing the amount of dose optimization computation.

[0060] Step 104: Based on the target grid that has a mapping relationship with the reference grid, obtain the target region of the beam spot of the at least one beam in the spatial grid set; the target region is used to participate in the beam spot dose optimization process.

[0061] After obtaining the reference grid, the target grid can be determined from multiple spatial grids based on the mapping relationship between the reference grid and the spatial grids in the spatial grid set. Then, based on the target grid, the target region of at least one beam spot in the spatial grid set can be determined, so that the spatial grid in the target region can be used to participate in the subsequent beam spot dose optimization process.

[0062] In one example, such as Figure 2 As shown, for each beam spot, the spatial location of the beam spot can be determined in the spatial grid set, and the grid dose of the spatial grid corresponding to the beam spot can be calculated. Then, for each beam spot, after obtaining the cutoff range of the beam spot in the BEV-WEPL grid (the cutoff region in the mapping grid set), the spatial dose grid (i.e., the target grid) corresponding to the cutoff range can be obtained based on the mapping relationship between the spatial grid set and the mapping grid set, and can be kept in memory to participate in subsequent beam spot weight optimization processing.

[0063] Compared to traditional methods that rely on the spatial location of the proton beam spot for geometric truncation (e.g., truncation within a spatial dose grid), which fail to consider the fluctuations in the three-dimensional morphology of the proton beam spot within the patient's body due to variations in material density along the range direction, making it difficult to accurately predict the distribution area of ​​the spatial dose grid, resulting in low truncation accuracy and high memory consumption, this embodiment's technical solution constructs a beam direction-equivalent water depth (BEV-WEPL) grid based on the spatial dose grid and the patient's medical images, establishing a mapping relationship between it and the spatial dose grid. This allows for geometric truncation within the BEV-WEPL grid, determining the spatial dose grid within the truncation range based on the mapping relationship, and participating in subsequent weight optimization. This improves the truncation accuracy of the spatial dose grid, saves memory resources, and more effectively preserves the effective dose (e.g., reducing the memory consumption of grids with dose values ​​of 0 or ~0), further reducing the computational load for dose optimization.

[0064] In the above-described beam spot processing method, a spatial grid set generated from the medical image of the target object is obtained, and the beam direction of at least one beam targeting the target object is obtained. The spatial grid located in the beam direction of at least one beam is used as the target grid. Based on the projection plane coordinates and equivalent water depth values ​​corresponding to each target grid, a mapping grid corresponding to each target grid is constructed to obtain a mapping grid set. Then, according to the preset truncation configuration information, the truncation region of the beam spot of at least one beam in the mapping grid set is determined. The mapping grid located in the truncation region is used as the reference grid. Then, based on the target grid with the mapping relationship with the reference grid, the target region of the beam spot of at least one beam in the spatial grid set is obtained. This realizes the optimization of beam spot truncation processing. By introducing the equivalent water depth of the beam spot to construct the corresponding mapping grid and setting the truncation in the mapping grid, the memory occupation can be reduced while ensuring the truncation accuracy, reducing the computational load of proton intensity modulation program production and improving beam spot processing efficiency.

[0065] In one embodiment, acquiring a set of spatial grids generated from medical images of a target object, and acquiring the beam direction of at least one beam directed towards the target object, and using the spatial grid located in the beam direction of the at least one beam as the target grid, may include the following steps:

[0066] A set of spatial grids generated from medical images of a target object is obtained, and the beam direction of at least one beam directed at the target object is obtained; based on the beam direction of the at least one beam, the spatial grid through which each ray passes in the beam direction is determined; based on the spatial grid through which each ray passes, a spatial grid in the beam direction of the at least one beam is obtained as the target grid.

[0067] In practical applications, a three-dimensional spatial dose grid (i.e., a spatial grid in a set of spatial grids) can be generated based on imported patient medical images (i.e., medical images of the target object). This spatial grid can contain all organs involved in radiotherapy. The accuracy of the spatial grid can be specified by the user in advance. The voxel corresponding to a single grid is the smallest unit for dose calculation in radiotherapy.

[0068] For example, for each beam, a series of rays can be emitted from the beam source point, and the spatial dose grid through which each ray passes can be used as the target grid. That is, based on the beam direction of at least one beam, the spatial grid through which each ray passes in that beam direction is determined.

[0069] In this embodiment, a set of spatial grids generated from medical images of the target object is obtained, and the beam direction of at least one beam targeting the target object is obtained. Then, based on the beam direction of at least one beam, the spatial grid through which each ray passes in the beam direction is determined. Based on the spatial grid through which each ray passes, the spatial grid in the beam direction of at least one beam is obtained as the target grid, providing data support for the subsequent construction of the beam direction-equivalent water depth grid.

[0070] In one embodiment, after the step of using the spatial grid located in the beam direction of the at least one beam as the target grid, the following step may also be included:

[0071] For each target grid, the coordinates of the target grid projected onto the isocentric plane are determined as the projection plane coordinates corresponding to the target grid; based on the material density of the target grid and the length of the ray path passing through the target grid, the equivalent water depth value corresponding to the target grid is obtained.

[0072] In one example, for each spatial dose grid traversed by a ray, the two-dimensional coordinates of the isocenter plane of the spatial dose grid (i.e., the projection plane coordinates), such as two-dimensional BEV coordinates, can be obtained by projecting any spatial dose grid (i.e., each target grid) onto the isocenter plane. The WEPL value (i.e., the equivalent water depth value) corresponding to the spatial dose grid can be obtained by accumulating the product of the material density of the spatial dose grid and the path length through the spatial dose grid.

[0073] For example, such as Figure 3a As shown, the spatial dose grid is related to CT. Different spatial dose grids correspond to different material densities (i.e., the material density of the target grid), which can be determined by the HU value of the CT for that spatial dose grid. The path length of a ray passing through a certain spatial dose grid can be obtained by accumulating the following:

[0074] ∑length*ρ

[0075] Where length is the distance the ray travels through the spatial dose grid, and ρ is the material density corresponding to that spatial dose grid.

[0076] In this embodiment, for each target grid, the coordinates of the target grid projected onto the isocenter plane are determined as the projection plane coordinates corresponding to the target grid. Then, based on the material density of the target grid and the length of the ray path passing through the target grid, the equivalent water depth value corresponding to the target grid is obtained, which provides data support for constructing the beam direction-equivalent water depth grid and establishing its mapping relationship with the spatial dose grid.

[0077] In one embodiment, determining the truncated region of the beam spot of the at least one beam in the mapping grid set according to preset truncation configuration information may include the following steps:

[0078] Obtain the energy and beam projection coordinates corresponding to the beam spot of the at least one beam; the beam projection coordinates are the coordinates of the beam spot projection of the at least one beam onto the isocenter plane; combine the energy and beam projection coordinates corresponding to the beam spot of the at least one beam, and the truncation configuration information, to obtain the truncation region of the beam spot of the at least one beam in the mapping grid set.

[0079] In one example, taking a proton beam spot as an example, for a certain beam spot, the cutoff range of the beam spot in the BEV-WEPL grid can be obtained based on the proton energy of the beam spot (i.e., the energy corresponding to the beam spot) and the isocenter plane coordinates of the beam spot (i.e., the beam spot projection coordinates) according to the geometric cutoff settings of the BEV-WEPL grid (i.e., the cutoff configuration information), that is, the cutoff region in the mapped grid set.

[0080] For example, such as Figure 3b As shown, the left diagram illustrates a spatial dose grid cutoff with two densities in three-dimensional space. The beam (such as a pencil beam) appears as a cone in three-dimensional space. A spatial dose grid can be established using medical images with two uniform densities. A beam is perpendicularly injected into the spatial dose grid, and a spot in the beam passes through the boundary between the two densities of the spatial dose grid.

[0081] For example, such as Figure 3b As shown in the diagram on the right, the beam direction-equivalent water depth (BEV-WEPL) grid is truncated from the beam direction BEV perspective. The beam (such as a pencil beam) appears as a cylinder in three-dimensional space. By calculating the WEPL value corresponding to the beam passing through the spatial dose grid and establishing the beam direction-equivalent water depth (BEV-WEPL) grid, a geometric truncation can be set for the beam spot in the beam direction-equivalent water depth (BEV-WEPL) grid, as shown in the bottom part of the cylinder in the diagram on the right. This yields the truncated region in the beam direction-equivalent water depth (BEV-WEPL) grid, which is the truncated region in the mapping grid set. Furthermore, based on the mapping relationship between the spatial dose grid and the BEV-WEPL grid (i.e., the spatial grid set and the mapping grid set), the effective grid region in the spatial dose grid (i.e., the target region in the spatial grid set) can be obtained, as shown in the bottom part of the cone in the diagram on the left.

[0082] In this embodiment, by obtaining the energy and beam projection coordinates corresponding to the beam spot of at least one beam, and then combining the energy and beam projection coordinates corresponding to the beam spot of at least one beam, as well as the truncation configuration information, the truncation region of the beam spot of at least one beam in the mapping grid set is obtained. Since the distribution of the proton beam spot remains relatively stable in the beam direction-equivalent water depth, the geometric truncation can be set better, which can reduce memory usage while ensuring truncation accuracy.

[0083] In one embodiment, before the step of determining the truncated region of the beam spot of the at least one beam in the mapping grid set according to preset truncation configuration information, the following step may be included:

[0084] The energy of the beam spot of the at least one beam is obtained, and the overall depth dose information of the beam spot of the at least one beam is obtained; based on the energy and the overall depth dose information, a preset cutoff configuration information for the beam spot of the at least one beam is determined.

[0085] In practical applications, taking the proton beam spot as an example, a preset cutoff configuration can be obtained by setting a geometric cutoff in the beam direction-equivalent water depth (BEV-WEPL) grid. This geometric cutoff setting can be determined based on the energy of different proton beam spots in different beams and the overall depth dose (IDD) peak width (i.e., overall depth dose information) of the proton beam spot. The beam direction BEV can be determined by the beam spot radius of the proton beam spot.

[0086] In this embodiment, by acquiring the energy of the beam spot of at least one beam and the overall depth dose information of the beam spot of at least one beam, and then determining the preset cutoff configuration information for the beam spot of at least one beam based on the energy and the overall depth dose information, it is helpful to set the cutoff in the beam direction-equivalent water depth grid.

[0087] In one embodiment, the truncation configuration information may include configuration parameters for multiple different truncation sizes. The step of determining the truncation region of the beam spot of the at least one beam in the mapping grid set according to the preset truncation configuration information, and using the mapping grid within the truncation region as a reference grid, may include the following steps:

[0088] Based on the configuration parameters for each cutoff size, the cutoff regions of different sizes of the beam spot of the at least one beam in the mapping grid set are determined respectively, and the reference grid corresponding to each size of the cutoff region is obtained.

[0089] The step of obtaining the target region of the beam spot of the at least one beam in the spatial grid set based on the target grid that has a mapping relationship with the reference grid includes:

[0090] Based on the truncated region of each size, and according to the target grid that has a mapping relationship with the reference grid in the truncated region, the beam spot of the at least one beam is obtained in the target region of different sizes in the spatial grid set.

[0091] In one example, such as Figure 4 As shown, the geometric cutoff in the beam direction-equivalent water depth (BEV-WEPL) grid can be set as multiple geometric cutoffs to adapt to the multiple iteration algorithms in subsequent weight optimization. For example, based on multiple geometric cutoffs, multiple different geometric cutoff sizes (i.e., configuration parameters of the cutoff size) can be set, and for different geometric cutoff sizes, the spatial dose grid corresponding to each geometric cutoff size can be recorded, i.e., the spatial grid in the target region of different sizes.

[0092] In this embodiment, by configuring parameters based on each cutoff size, the cutoff regions of at least one beam spot in the mapping grid set of different sizes are determined, and the reference grid corresponding to each cutoff region is obtained. Based on each cutoff region of different sizes, the target regions of at least one beam spot in the spatial grid set of different sizes are obtained according to the target grid that has a mapping relationship with the reference grid in the cutoff region. This can adapt to the multiple iterative algorithms in the subsequent weight optimization based on multiple geometric cutoff settings.

[0093] In one embodiment, after the step of obtaining the target region of the beam spot of the at least one beam in the spatial grid set based on the target grid having a mapping relationship with the reference grid, the following step may be further included:

[0094] Based on the beam spot of the at least one beam, different optimization strategies are used to optimize the beam spot dose in target regions of different sizes in the spatial grid set.

[0095] In practical implementation, the cutoff settings based on the beam direction-equivalent water depth (BEV-WEPL) grid can be configured in the following ways to set different geometric cutoff sizes:

[0096] For the dual-threshold case, a smaller cutoff size can be set as the dose core region, in which the spatial dose grid recorded can be a high-dose grid; a larger cutoff size can be set and the dose core region subtracted to form the dose non-core region, in which the spatial dose grid recorded can be a low-dose grid. Thus, different strategies can be used for the two regions during dose optimization, which can reduce the computational load of dose optimization.

[0097] In this embodiment, by using different optimization strategies to optimize the beam spot dose based on the beam spot of at least one beam in target regions of different sizes in a spatial grid set, the computational load of dose optimization can be effectively reduced.

[0098] In one embodiment, such as Figure 5 The diagram shows a flowchart of another beam spot treatment method. In this embodiment, the method includes the following steps:

[0099] In step 501, a spatial mesh set generated from the medical image of the target object is obtained, as well as the beam direction of at least one beam targeting the target object. The spatial mesh located in the beam direction of at least one beam is taken as the target mesh. In step 502, a mapping mesh corresponding to each target mesh is constructed based on the projection plane coordinates and equivalent water depth value corresponding to each target mesh, resulting in a mapping mesh set. In step 503, the energy of the beam spot of at least one beam and the overall depth-dose information of the beam spot of at least one beam are obtained. In step 504, based on the energy and overall depth-dose information, a preset truncation configuration information for the beam spot of at least one beam is determined; the truncation configuration information includes configuration parameters for multiple different truncation sizes. In step 505, the energy and beam spot projection coordinates corresponding to the beam spot of at least one beam are obtained; the beam spot projection coordinates are the coordinates of the beam spot projection of at least one beam onto the isocenter plane. In step 506, by combining the energy and beam spot projection coordinates corresponding to the beam spot of at least one beam, and the configuration parameters based on each cutoff size, the cutoff regions of different sizes for the beam spot of at least one beam in the mapping grid set are determined, resulting in a reference grid corresponding to each size of cutoff region. In step 507, based on the target regions of different sizes for the beam spot of at least one beam in the spatial grid set, different optimization strategies are applied to optimize the beam spot dose. It should be noted that the specific limitations of the above steps can be found in the specific limitations of a beam spot processing method described above, and will not be repeated here.

[0100] It should be understood that although the steps in the flowcharts of the embodiments described above are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the embodiments described above may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages of other steps.

[0101] Based on the same inventive concept, this application also provides a beam spot processing system for implementing the beam spot processing method described above. The solution provided by this system is similar to the solution described in the above method; therefore, the specific limitations in one or more beam spot processing system embodiments provided below can be found in the limitations of the beam spot processing method described above, and will not be repeated here.

[0102] In one embodiment, such as Figure 6 As shown, a beam spot processing system is provided, comprising:

[0103] The target mesh acquisition module 601 is used to acquire a set of spatial meshes generated based on medical images of a target object, and to acquire the beam direction of at least one beam for the target object, and to take the spatial mesh in the beam direction of the at least one beam as the target mesh.

[0104] The mapping mesh acquisition module 602 is used to construct a mapping mesh corresponding to each target mesh based on the projection plane coordinates and equivalent water depth values ​​corresponding to each target mesh, thereby obtaining a mapping mesh set; the mapping mesh set is used to characterize the material density fluctuation in the beam direction of the at least one beam;

[0105] The truncated region determination module 603 is used to determine the truncated region of the beam spot of the at least one beam in the mapping grid set according to the preset truncated configuration information, and to use the mapping grid in the truncated region as the reference grid.

[0106] The target region acquisition module 604 is used to obtain the target region of the beam spot of the at least one beam in the spatial grid set according to the target grid that has a mapping relationship with the reference grid; the target region is used to participate in the beam spot dose optimization process.

[0107] In one embodiment, the target mesh acquisition module 601 includes:

[0108] The beam direction acquisition submodule is used to acquire a set of spatial grids generated from medical images of a target object, and to acquire the beam direction of at least one beam for the target object;

[0109] A grid-determining submodule is used to determine the spatial grid through which each ray in the beam direction passes, based on the beam direction of the at least one beam.

[0110] The target mesh obtaining submodule is used to obtain a spatial mesh in the beam direction of the at least one beam based on the spatial mesh through which each of the rays passes, as the target mesh.

[0111] In one embodiment, the system further includes:

[0112] The projection plane coordinate acquisition module is used to determine the coordinates of the target grid projected onto the isocenter plane for each target grid, and use them as the projection plane coordinates corresponding to the target grid;

[0113] The equivalent water depth value acquisition module is used to obtain the equivalent water depth value corresponding to the target grid based on the material density of the target grid and the length of the ray path passing through the target grid.

[0114] In one embodiment, the truncation region determination module 603 includes:

[0115] The beam spot projection coordinate acquisition submodule is used to acquire the energy and beam spot projection coordinates corresponding to the beam spot of the at least one beam; the beam spot projection coordinates are the coordinates of the beam spot projection of the at least one beam onto the isocenter plane;

[0116] The truncated region submodule is used to combine the energy and beam spot projection coordinates corresponding to the beam spot of the at least one beam, as well as the truncated configuration information, to obtain the truncated region of the beam spot of the at least one beam in the mapping grid set.

[0117] In one embodiment, the system further includes:

[0118] An energy and overall depth dose information acquisition module is used to acquire the energy of the beam spot of the at least one beam, and to acquire the overall depth dose information of the beam spot of the at least one beam;

[0119] The cutoff configuration information determination module is used to determine preset cutoff configuration information for the beam spot of the at least one beam based on the energy and the overall depth dose information.

[0120] In one embodiment, the truncation configuration information includes configuration parameters for multiple different truncation sizes, and the truncation region determination module 603 includes:

[0121] Different cut-off regions are obtained by sub-modules, which are used to determine the cut-off regions of different sizes of the beam spot of the at least one beam in the mapping grid set based on the configuration parameters of each cut-off size, and obtain the reference grid corresponding to each size of cut-off region;

[0122] The target region acquisition module 604 includes:

[0123] Different target regions are obtained by sub-modules, which are used to obtain the beam spot of the at least one beam of different sizes in the spatial grid set based on the truncated region of each size and the target grid that has a mapping relationship with the reference grid in the truncated region.

[0124] In one embodiment, the system further includes:

[0125] The separate dose optimization module is used to perform beam spot dose optimization processing based on the beam spot of the at least one beam in target regions of different sizes in the spatial grid set, using different optimization strategies respectively.

[0126] Each module in the aforementioned beam spot processing system can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in the processor of a computer device in hardware form or independent of it, or stored in the memory of the computer device in software form, so that the processor can call and execute the corresponding operations of each module.

[0127] In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as follows: Figure 7 As shown, the computer device includes a processor, memory, communication interface, display screen, and input device connected via a system bus. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, mobile cellular networks, NFC (Near Field Communication), or other technologies. When executed by the processor, the computer program implements a beam pattern processing method.

[0128] Those skilled in the art will understand that Figure 7 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0129] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to perform the following steps:

[0130] Obtain a set of spatial grids generated from medical images of a target object, and obtain the beam direction of at least one beam directed toward the target object. Use the spatial grids located in the beam direction of the at least one beam as the target grid.

[0131] Based on the projection plane coordinates and equivalent water depth values ​​corresponding to each target grid, a mapping grid corresponding to each target grid is constructed to obtain a mapping grid set; the mapping grid set is used to characterize the material density fluctuation in the beam direction of the at least one beam;

[0132] Based on the preset truncation configuration information, the truncation region of the beam spot of the at least one beam in the mapping grid set is determined, and the mapping grid in the truncation region is used as the reference grid.

[0133] Based on the target grid that has a mapping relationship with the reference grid, the target region of the beam spot of the at least one beam in the spatial grid set is obtained; the target region is used to participate in the beam spot dose optimization process.

[0134] In one embodiment, the processor, when executing a computer program, also implements the steps of the beam speckle processing method in the other embodiments described above.

[0135] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon, the computer program performing the following steps when executed by a processor:

[0136] Obtain a set of spatial grids generated from medical images of a target object, and obtain the beam direction of at least one beam directed toward the target object. Use the spatial grids located in the beam direction of the at least one beam as the target grid.

[0137] Based on the projection plane coordinates and equivalent water depth values ​​corresponding to each target grid, a mapping grid corresponding to each target grid is constructed to obtain a mapping grid set; the mapping grid set is used to characterize the material density fluctuation in the beam direction of the at least one beam;

[0138] Based on the preset truncation configuration information, the truncation region of the beam spot of the at least one beam in the mapping grid set is determined, and the mapping grid in the truncation region is used as the reference grid.

[0139] Based on the target grid that has a mapping relationship with the reference grid, the target region of the beam spot of the at least one beam in the spatial grid set is obtained; the target region is used to participate in the beam spot dose optimization process.

[0140] In one embodiment, when the computer program is executed by a processor, it also implements the steps of the beam spot processing method in the other embodiments described above.

[0141] In one embodiment, a computer program product is provided, including a computer program that, when executed by a processor, performs the following steps:

[0142] Obtain a set of spatial grids generated from medical images of a target object, and obtain the beam direction of at least one beam directed toward the target object. Use the spatial grids located in the beam direction of the at least one beam as the target grid.

[0143] Based on the projection plane coordinates and equivalent water depth values ​​corresponding to each target grid, a mapping grid corresponding to each target grid is constructed to obtain a mapping grid set; the mapping grid set is used to characterize the material density fluctuation in the beam direction of the at least one beam;

[0144] Based on the preset truncation configuration information, the truncation region of the beam spot of the at least one beam in the mapping grid set is determined, and the mapping grid in the truncation region is used as the reference grid.

[0145] Based on the target grid that has a mapping relationship with the reference grid, the target region of the beam spot of the at least one beam in the spatial grid set is obtained; the target region is used to participate in the beam spot dose optimization process.

[0146] In one embodiment, when the computer program is executed by a processor, it also implements the steps of the beam spot processing method in the other embodiments described above.

[0147] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0148] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0149] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A method for treating beam spots, characterized in that, The method includes: Obtain a set of spatial grids generated from medical images of a target object, and obtain the beam direction of at least one beam directed toward the target object. Use the spatial grids located in the beam direction of the at least one beam as the target grid. Based on the material density of the target grid and the ray path length passing through the target grid, the equivalent water depth value corresponding to the target grid is obtained; based on the projection plane coordinates and equivalent water depth value corresponding to each target grid, a mapping grid corresponding to each target grid is constructed to obtain a mapping grid set; the mapping grid set is used to characterize the material density fluctuation in the beam direction of the at least one beam. Based on the preset truncation configuration information, the truncation region of the beam spot of the at least one beam in the mapping grid set is determined, and the mapping grid in the truncation region is used as the reference grid. Based on the target grid that has a mapping relationship with the reference grid, the target region of the beam spot of the at least one beam in the spatial grid set is obtained; the target region is used to participate in the beam spot dose optimization process.

2. The method according to claim 1, characterized in that, The step of obtaining the beam direction of at least one beam targeting the target object, and using the spatial grid along the beam direction of the at least one beam as the target grid, includes: Based on the beam direction of the at least one beam, determine the spatial grid through which each ray in the beam direction passes; Based on the spatial grid through which each of the rays passes, a spatial grid in the beam direction of the at least one beam is obtained, which serves as the target grid.

3. The method according to claim 1, characterized in that, After the step of taking the spatial grid located in the beam direction of the at least one beam as the target grid, the method further includes: For each target mesh, determine the coordinates of the target mesh projected onto the isocenter plane, and use these coordinates as the projection plane coordinates corresponding to the target mesh.

4. The method according to claim 1, characterized in that, The step of determining the truncated region of the beam spot of the at least one beam in the mapping grid set according to preset truncation configuration information includes: Obtain the energy and beam projection coordinates corresponding to the beam spot of the at least one beam; By combining the energy and beam projection coordinates corresponding to the beam spot of the at least one beam, and the truncation configuration information, the truncation region of the beam spot of the at least one beam in the mapping grid set is obtained.

5. The method according to claim 1, characterized in that, Before the step of determining the truncated region of the beam spot of the at least one beam in the mapping grid set according to preset truncation configuration information, the method further includes: The energy of the beam spot of the at least one beam is obtained, and the overall depth dose information of the beam spot of the at least one beam is obtained; Based on the energy and the overall depth dose information, a preset cutoff configuration for the beam spot of the at least one beam is determined.

6. The method according to claim 1, characterized in that, The truncation configuration information includes configuration parameters for multiple different truncation sizes. The step of determining the truncation region of the beam spot of the at least one beam in the mapping grid set according to the preset truncation configuration information, and using the mapping grid within the truncation region as a reference grid, includes: Based on the configuration parameters for each cutoff size, the cutoff regions of different sizes for the beam spots of the at least one beam in the mapping grid set are determined respectively, and the reference grid corresponding to each size cutoff region is obtained.

7. The method according to claim 6, characterized in that, After the step of obtaining the target region of the beam spot of the at least one beam in the spatial grid set based on the target grid having a mapping relationship with the reference grid, the method further includes: Based on the beam spot of the at least one beam, different optimization strategies are used to optimize the beam spot dose in target regions of different sizes in the spatial grid set.

8. A beam spot treatment system, characterized in that, The system includes: The target mesh acquisition module is used to acquire a set of spatial meshes generated based on medical images of a target object, and to acquire the beam direction of at least one beam directed toward the target object, and to take the spatial meshes located in the beam direction of the at least one beam as the target meshes; The mapping mesh acquisition module is used to obtain the equivalent water depth value corresponding to the target mesh based on the material density of the target mesh and the ray path length passing through the target mesh; and to construct the mapping mesh corresponding to each target mesh based on the projection plane coordinates and equivalent water depth value of each target mesh, thereby obtaining a mapping mesh set; the mapping mesh set is used to characterize the material density fluctuation in the beam direction of the at least one beam; The truncation region determination module is used to determine the truncation region of the beam spot of the at least one beam in the mapping grid set according to the preset truncation configuration information, and to use the mapping grid in the truncation region as the reference grid. The target region acquisition module is used to obtain the target region of the beam spot of the at least one beam in the spatial grid set based on the target grid that has a mapping relationship with the reference grid; the target region is used to participate in the beam spot dose optimization process.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, 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 to 7.