A radiotherapy treatment plan optimization method and device

Abstract

The present invention discloses a radiotherapy plan optimization method and device. Among them, the radiotherapy plan optimization method includes: obtaining an initial intensity matrix of a radiotherapy beam; calculating according to the initial intensity matrix to obtain the irradiation dose MU of each sub-beam under this beam; performing a descending order arrangement operation on the irradiation dose MU of each sub-beam, and performing a filtering operation on the rows of the irradiation dose MU before sorting for the top N%, where N is a set value. After the present invention performs a descending order arrangement operation on the irradiation dose MU of each sub-beam, it only performs a filtering operation on the rows of the irradiation dose MU before sorting for the top N%. On the one hand, it realizes filtering operation on the rows of larger irradiation dose MU after sorting, reduces the irradiation dose MU, and improves the treatment efficiency; on the other hand, it only performs a filtering operation on the rows of the irradiation dose MU before sorting for the top N%, avoiding that the excessive proportion of the filtering operation may cause a large difference between the intensity matrix after the filtering operation and the initial intensity matrix, affecting the radiotherapy accuracy.

Description

Technical Field

[0001] The present invention relates to the technical field of radiotherapy equipment, and particularly to a radiotherapy plan optimization method and device. Background Art

[0002] In the radiotherapy industry, in order to solve the problem of low execution efficiency of radiotherapy plans, it solves this problem through a grid method. However, due to the design principle of the intensity matrix, this method also has problems and disadvantages such as excessive merging of intensity values, resulting in inaccurate intensity matrices.

[0003] In view of this, the present invention patent is proposed. Summary of the Invention

[0004] The present invention provides a radiotherapy plan optimization method and device, which solves the problem of low execution efficiency of radiotherapy plans, and at the same time overcomes the problems and disadvantages of the intensity matrix merging grid scheme.

[0005] Specifically, the following technical solutions are adopted:

[0006] A radiotherapy plan optimization method includes:

[0007] Obtain the initial intensity matrix of the radiotherapy beam;

[0008] Calculate according to the initial intensity matrix to obtain the irradiation dose MU of each sub-beam under this beam;

[0009] Perform a descending order operation on the irradiation dose MU of each sub-beam, and perform a filtering operation on the first N% of the irradiation dose MU rows before sorting, where N is a set value.

[0010] As an optional implementation manner of the present invention, in a radiotherapy plan optimization method of the present invention, the calculating according to the initial intensity matrix to obtain the irradiation dose MU of each sub-beam under this beam includes:

[0011] According to the initial intensity matrix, compare the intensity value of the element [i] with the intensity value of its left element [i - 1] row by row;

[0012] If the intensity value of the element [i] is less than the intensity value of the element [i - 1], record the intensity value of the element [i] as 0. If the intensity value of the element [i] is greater than or equal to the intensity value of the element [i - 1], record the intensity value of the element [i] as the difference between the two, to obtain a new intensity matrix;

[0013] Traverse and accumulate the sum of the intensity values of each element in the new intensity matrix to obtain the irradiation dose MU of the corresponding sub-beam for each row.

[0014] As an alternative implementation of the present invention, in a radiotherapy plan optimization method of the present invention, after obtaining the irradiation dose MU of each sub-beam under the beam, traverse the irradiation dose MU of all sub-beams to obtain the maximum irradiation dose MU-max.

[0015] As an alternative implementation of the present invention, in a radiotherapy plan optimization method of the present invention, the operation of sorting the irradiation dose MU of each sub-beam in descending order, and the filtering operation for the irradiation dose MU of the first N% before sorting includes:

[0016] Sort the irradiation dose MU of each sub-beam in descending order, and obtain the rows of the irradiation dose MU of the first N% before sorting as the target rows for the filtering operation. The total number of target rows is rounded.

[0017] Judge whether each target row in the initial intensity matrix meets the filtering operation strategy. If it meets, execute the corresponding filtering operation strategy. If it does not meet, take the average value of the intensities of the three adjacent elements [i-1], [i], and [i+1] in each target row of the initial intensity matrix as the new intensity value of the middle position element [i].

[0018] As an alternative implementation of the present invention, in a radiotherapy plan optimization method of the present invention, the filtering operation strategy includes a first filtering operation strategy. Judging whether each target row in the initial intensity matrix meets the first filtering operation strategy includes:

[0019] Judge whether there is an element [i] in each target row of the initial intensity matrix whose intensity value is equal to the maximum irradiation dose MU-max. If it exists, it meets the first filtering operation strategy, update the element [i] to k1*MU-max, and jump out of this loop, where k1 is a preset value.

[0020] As an alternative implementation of the present invention, in a radiotherapy plan optimization method of the present invention, the filtering operation strategy includes a second filtering operation strategy. If the judgment result of the first filtering operation strategy is that it does not exist, then judge whether each target row in the initial intensity matrix meets the second filtering operation strategy. The second filtering operation strategy includes:

[0021] Judge whether there is a situation where the intensity values of the elements [i-1], [i], and [i+1] in each target row of the initial intensity matrix are less than k2*MU-max. If it exists, it meets the second filtering operation strategy, and the intensity value of the element [i] remains unchanged, and jump out of this loop, where k2 is a preset value.

[0022] As an alternative implementation of the present invention, in a radiotherapy plan optimization method of the present invention, the filtering operation strategy includes a third filtering operation strategy. If the judgment result of the second filtering operation strategy is negative, it is determined whether each target row in the initial intensity matrix satisfies the third filtering operation strategy. The third filtering operation strategy includes:

[0023] Determine whether each target row in the initial intensity matrix satisfies |diff1| > 1 / 2 * MU - max or |diff2| > 1 / 2 * MU - max, where diff1 = Original[i] - Original[i - 1]; diff2 = Original[i + 1] - Original[i], diff1 is the left difference, diff2 is the right difference, Original is the intensity matrix, and Original[i] is the intensity value of element [i] in the initial intensity matrix;

[0024] If the judgment result is positive, the third filtering operation strategy is satisfied, the intensity value of element [i] remains unchanged, and this loop is exited. If the judgment result is negative, it is determined that each target row in the initial intensity matrix does not satisfy the filtering operation strategy.

[0025] As an alternative implementation of the present invention, in a radiotherapy plan optimization method of the present invention, according to the intensity matrix processed by the filtering operation, the vane movement strategy of the multi - leaf collimator is executed:

[0026] For each element [i] in the intensity matrix processed by the filtering operation, control the initial position of the vane to be in the closed state;

[0027] Control the left vane to remain stationary at the initial position, and move the right vane to the shaped sub - field position, and maintain the position states of the left and right vanes for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of element [i];

[0028] Control the right vane to remain stationary at the initial position, and move the left vane to the shaped sub - field position, and maintain the position states of the left and right vanes for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of element [i];

[0029] Then the vane movement strategy for element [i] is completed.

[0030] The present invention also provides a radiotherapy plan optimization device, including:

[0031] An intensity acquisition module, which acquires the initial intensity matrix under the radiotherapy field;

[0032] An irradiation dose calculation module, which calculates according to the initial intensity matrix to obtain the irradiation dose MU of each sub - field under the field;

[0033] The filtering module performs a descending order arrangement operation on the irradiation doses MU of each sub - beam, and performs a filtering operation on the irradiation doses MU of the top N% before sorting, where N is a set value.

[0034] As an optional implementation manner of the present invention, a radiotherapy plan optimization device of the present invention includes a leaf movement module:

[0035] According to the intensity matrix processed by the filtering operation, execute the leaf movement strategy of the multi - leaf collimator:

[0036] For each element [i] in the intensity matrix processed by the filtering operation, control the initial position of the leaf to be in a closed state;

[0037] Control the left leaf to remain stationary at the initial position, and move the right leaf to the shaping sub - beam position, and maintain the position states of the left leaf and the right leaf for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of the element [i];

[0038] Control the right leaf to remain stationary at the initial position, and move the left leaf to the shaping sub - beam position, and maintain the position states of the left leaf and the right leaf for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of the element [i];

[0039] Then the leaf movement strategy for the element [i] is completed.

[0040] Compared with the prior art, the beneficial effects of the present invention:

[0041] In a radiotherapy plan optimization method of the present invention, after performing a descending order arrangement operation on the irradiation doses MU of each sub - beam, only perform a filtering operation on the irradiation doses MU of the top N% before sorting. On the one hand, it realizes filtering the larger irradiation doses MU after sorting, reduces the irradiation dose MU, and improves the treatment efficiency. On the other hand, only performing a filtering operation on the irradiation doses MU of the top N% before sorting avoids the situation that the proportion of the filtering operation is too large, which may cause a large difference between the intensity matrix after filtering and the initial intensity matrix, affecting the radiotherapy accuracy.

[0042] A radiotherapy plan optimization method of the present invention sets the leaf movement strategy of the multi - leaf collimator according to the intensity matrix after filtering, and realizes the movement trajectory of the corresponding leaf. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] Figure 1 Schematic diagram of the execution of the leaf movement strategy in a radiotherapy plan optimization method according to an embodiment of the present invention;

[0044] Figure 2 Block diagram of the modules of a radiotherapy plan optimization device according to an embodiment of the present invention. DETAILED DESCRIPTION OF THE INVENTION

[0045] To make the objectives, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part rather than all of the embodiments of the present invention.

[0046] Therefore, the following detailed description of the embodiments of the present invention is not intended to limit the scope of the claimed invention, but merely represents some embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the scope of protection of the present invention.

[0047] It should be noted that, without conflict, the embodiments in the present invention and the features and technical solutions in the embodiments may be combined with each other.

[0048] It should be noted that like reference numerals and letters denote like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.

[0049] In the description of the present invention, it should be noted that the orientation or positional relationship indicated by the terms "upper", "lower", etc. is based on the orientation or positional relationship shown in the drawings, or the orientation or positional relationship when the product of the present invention is normally placed, or the orientation or positional relationship commonly understood by those skilled in the art. Such terms are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus should not be construed as limiting the present invention. In addition, the terms "first", "second", etc. are only used for descriptive distinction and cannot be construed as indicating or implying relative importance.

[0050] A radiotherapy plan optimization method according to this embodiment includes:

[0051] Obtaining an initial intensity matrix of a radiotherapy beam;

[0052] Calculating according to the initial intensity matrix to obtain the irradiation dose MU of each sub-beam under this beam;

[0053] Performing a descending order arrangement operation on the irradiation dose MU of each sub-beam, and performing a filtering operation on the rows of the irradiation dose MU before sorting with the top N%, where N is a set value.

[0054] In a radiotherapy plan optimization method of this embodiment, after arranging the irradiation doses MU of each sub-beam in descending order, a filtering operation is only performed on the irradiation doses MU of the top N% before sorting. On the one hand, it realizes filtering the larger irradiation doses MU after sorting, reducing the irradiation dose MU and improving the treatment efficiency. On the other hand, only filtering the irradiation doses MU of the top N% before sorting can avoid the large proportion of the filtering operation, which may lead to a large difference between the intensity matrix after filtering and the initial intensity matrix, affecting the radiotherapy accuracy.

[0055] In addition, in this embodiment, a filtering operation is performed on the top N% of the irradiation doses MU after sorting, and N is set according to experience. Optionally, N is taken as 30.

[0056] As an optional implementation manner of this embodiment, in a radiotherapy plan optimization method of this embodiment, the calculation based on the initial intensity matrix to obtain the irradiation doses MU of each sub-beam in this beam includes:

[0057] According to the initial intensity matrix, compare the intensity value of the element [i] with the intensity value of its left element [i - 1] row by row;

[0058] If the intensity value of the element [i] is less than the intensity value of the element [i - 1], record the intensity value of the element [i] as 0. If the intensity value of the element [i] is greater than or equal to the intensity value of the element [i - 1], record the intensity value of the element [i] as the difference between the two, and obtain a new intensity matrix;

[0059] Traverse and accumulate the sum of the intensity values of each element in the new intensity matrix to obtain the irradiation doses MU of the corresponding sub-beams in each row.

[0060] In this embodiment, the irradiation doses MU of each sub-beam in this beam are calculated according to the initial intensity matrix, and descending sorting and filtering operations can be performed according to the irradiation doses MU of each sub-beam.

[0061] At the same time, in a radiotherapy plan optimization method of this embodiment, after obtaining the irradiation doses MU of each sub-beam in this beam, traverse the irradiation doses MU of all sub-beams to obtain the maximum irradiation dose MU-max.

[0062] Further, in a radiotherapy plan optimization method of this embodiment, the operation of arranging the irradiation doses MU of each sub-beam in descending order and performing a filtering operation on the irradiation doses MU of the top N% before sorting includes:

[0063] Arrange the irradiation doses MU of each sub-beam in descending order, and obtain the rows of the irradiation doses MU of the top N% before sorting as the target rows for the filtering operation, and round up the total number of target rows;

[0064] Determine whether each target row in the initial intensity matrix meets the filtering operation strategy. If it meets, execute the corresponding filtering operation strategy. If it does not meet, take the average of the intensities of the three adjacent elements [i-1], [i], and [i+1] in each target row of the initial intensity matrix as the new intensity value of the middle position element [i].

[0065] In the radiotherapy plan optimization method of this embodiment, the total number of target rows is rounded. Specifically, for example, if the total number of sorted irradiation dose MU rows is 3 rows and N is 30, then the total number of target rows that need to be filtered is 3 * 30% = 0.9, which is rounded to 1. Then, select the row with the largest sorted irradiation dose MU for optimization.

[0066] Specifically, in the radiotherapy plan optimization method of this embodiment, the filtering operation strategy includes a first filtering operation strategy. Determining whether each target row in the initial intensity matrix meets the first filtering operation strategy includes:

[0067] Determine whether there is an element [i] in each target row of the initial intensity matrix whose intensity value is equal to the maximum irradiation dose MU-max. If it exists, it meets the first filtering operation strategy, update the element [i] to k1 * MU-max, and jump out of this loop, where k1 is a preset value.

[0068] The first filtering operation strategy of this embodiment performs a filtering operation on the element [i] in the target row whose intensity value is equal to the maximum irradiation dose MU-max, updates it to k1 * MU-max, reduces the irradiation dose MU, and improves the treatment efficiency. Specifically, k1 in this embodiment can be taken as 0.95.

[0069] Further, in the radiotherapy plan optimization method of this embodiment, the filtering operation strategy includes a second filtering operation strategy. If the judgment result of the first filtering operation strategy is that it does not exist, then determine whether each target row in the initial intensity matrix meets the second filtering operation strategy. The second filtering operation strategy includes:

[0070] Determine whether there is a situation where the intensity values of the elements [i-1], [i], and [i+1] in each target row of the initial intensity matrix are less than k2 * MU-max. If it exists, it meets the second filtering operation strategy, and the intensity value of the element [i] remains unchanged, and jump out of this loop, where k2 is a preset value.

[0071] The second filtering operation strategy of this embodiment performs a filtering operation on the element [i] in the target row with too small intensity value. Optionally, k2 in this embodiment can be selected as 1 / 10.

[0072] Further, in a radiotherapy plan optimization method according to this embodiment, the filtering operation strategy includes a third filtering operation strategy. If the judgment result of the second filtering operation strategy is negative, it is judged whether each target row in the initial intensity matrix satisfies the third filtering operation strategy. The third filtering operation strategy includes:

[0073] Judge whether each target row in the initial intensity matrix satisfies |diff1| > 1 / 2 * MU - max or |diff2| > 1 / 2 * MU - max, where diff1 = Original[i] - Original[i - 1]; diff2 = Original[i + 1] - Original[i], diff1 is the left difference, diff2 is the right difference, Original is the intensity matrix, and Original[i] is the intensity value of element [i] in the initial intensity matrix;

[0074] If the judgment result is affirmative, it satisfies the third filtering operation strategy, the intensity value of element [i] remains unchanged, and this loop is exited. If the judgment result is negative, it is determined that each target row in the initial intensity matrix does not satisfy the filtering operation strategy.

[0075] As an alternative implementation of this embodiment, refer to Figure 1 As shown, in a radiotherapy plan optimization method according to this embodiment, according to the intensity matrix processed by the filtering operation, the vane motion strategy of the multi - leaf collimator is executed:

[0076] For each element [i] in the intensity matrix processed by the filtering operation, control the initial position of the vane to be in the closed state;

[0077] Control the left vane to remain stationary at the initial position, and move the right vane to the shaped sub - field position. Keep the position states of the left vane and the right vane for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of element [i];

[0078] Control the right vane to remain stationary at the initial position, and move the left vane to the shaped sub - field position. Keep the position states of the left vane and the right vane for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of element [i];

[0079] Then the vane motion strategy for element [i] is completed.

[0080] A radiotherapy plan optimization method according to this embodiment sets the vane motion strategy of the multi - leaf collimator according to the filtered intensity matrix to realize the motion trajectory of the corresponding vane.

[0081] Further, a radiotherapy plan optimization method according to this embodiment includes successively executing the vane movement strategy of the multi-leaf collimator for each element in the intensity matrix after filtering operation until the intensity values of all elements can be achieved.

[0082] A radiotherapy plan optimization method according to this embodiment is specifically exemplified as follows:

[0083] Given an initial intensity matrix A of a radiation field:

[0084]

[0085]

[0086] Step 1: Calculate the required irradiation dose MU for each row. Pad "0" to the left of the leftmost element in each row. The calculated values are as follows:

[0087] The first row = (3 - 0) + 0 + (3 - 1) = 5;

[0088] The second row = (2 - 0) + (3 - 2) = (2 + 1) = 3;

[0089] The third row = (4 - 0) + 0 + (1 - 1) = 4. Take the maximum value of the irradiation dose MU as the maximum irradiation dose MU-max, that is, MU-max = 5.

[0090] Step 2: Filter the top 30% of the maximum rows, that is, filter the maximum row, that is, filter the "1"st row.

[0091] Step 3: Determine whether there is an element [i] in each target row of the initial intensity matrix whose intensity value is equal to the maximum irradiation dose MU-max = 5. Among them, the maximum value in the "1"st row of the initial intensity matrix A is 4, which does not meet the first filtering operation strategy. Proceed to the next step.

[0092] Step 4: Determine whether there is a situation where the intensity values of elements [i - 1], [i], and [i + 1] in each target row of the initial intensity matrix are less than k2 * MU-max. After judging the "1"st row of the initial intensity matrix A, it does not meet the second filtering operation strategy. Proceed to the next step.

[0093] Step 5: Determine whether each target row of the initial intensity matrix meets |diff1| > 1 / 2 * MU-max or |diff2| > 1 / 2 * MU-max. For the "1"st row of the initial intensity matrix A, calculate diff1 = (1 - 3) and diff2 = (3 - 1). Since |diff1| is equal to 1 / 2 * MU-max or |diff2| is equal to 1 / 2 * MU-max, it does not meet the third filtering operation strategy. Proceed to the next step.

[0094] Step 6: For the intensity values 3, 1, 3 of three adjacent elements in the first row of the initial intensity matrix A, take the average value (3 + 1 + 3) / 3 = 7 / 3 as the new intensity value of the middle position element.

[0095] Step 7: So the final intensity matrix is

[0096] 3 7 / 3 3 2 3 3 4 1 1

[0097] Step 8: For each element [i] in the above intensity matrix, control the initial position of the leaf to be in the closed state;

[0098] Control the left leaf to remain stationary at the initial position, move the right leaf to the shaped sub - field position, and maintain the position states of the left and right leaves for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of element [i];

[0099] Control the right leaf to remain stationary at the initial position, move the left leaf to the shaped sub - field position, and maintain the position states of the left and right leaves for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of element [i];

[0100] Then the leaf motion strategy for element [i] is completed.

[0101] Step 9: Sequentially execute the above - mentioned leaf motion strategy of the multi - leaf collimator for each element in the intensity matrix after the filtering operation until the intensity values of all elements can be achieved.

[0102] See Figure 2 As shown, this embodiment also provides a radiotherapy plan optimization device, including:

[0103] An intensity acquisition module, which acquires the initial intensity matrix under the radiotherapy field;

[0104] An irradiation dose calculation module, which calculates according to the initial intensity matrix to obtain the irradiation dose MU of each sub - field under this field;

[0105] A filtering module, which performs a descending order arrangement operation on the irradiation dose MU of each sub - field, and performs a filtering operation on the first N% of the irradiation dose MU rows before sorting, where N is a set value.

[0106] A radiotherapy plan optimization device according to this embodiment. After the filtering module sorts the irradiation doses MU of each sub-beam in descending order, it only performs filtering operations on the irradiation doses MU of the top N% before sorting. On the one hand, it realizes filtering operations on the rows of larger irradiation doses MU after sorting, reduces the irradiation dose MU, and improves the treatment efficiency. On the other hand, only performing filtering operations on the rows of the top N% of the irradiation doses MU before sorting avoids the situation that the excessive proportion of filtering operations may lead to a large difference between the intensity matrix after filtering and the initial intensity matrix, affecting the radiotherapy accuracy.

[0107] In addition, this embodiment performs filtering operations on the rows of the top N% of the irradiation doses MU after sorting. N is set according to experience. Optionally, 0 < N ≤ 50 is set, and N takes 30.

[0108] Furthermore, a radiotherapy plan optimization device according to this embodiment includes a leaf motion module:

[0109] According to the intensity matrix processed by the filtering operation, execute the leaf motion strategy of the multi-leaf collimator:

[0110] For each element [i] in the intensity matrix processed by the filtering operation, control the initial position of the leaf to be in the closed state;

[0111] Control the left leaf to remain in the initial position without moving, and the right leaf to move to the shaped sub-beam position, and maintain the position states of the left leaf and the right leaf for a certain time to obtain the irradiation dose MU corresponding to the intensity value of the element [i];

[0112] Control the right leaf to remain in the initial position without moving, and the left leaf to move to the shaped sub-beam position, and maintain the position states of the left leaf and the right leaf for a certain time to obtain the irradiation dose MU corresponding to the intensity value of the element [i];

[0113] Then the leaf motion strategy for the element [i] is completed.

[0114] A radiotherapy plan optimization device according to this embodiment. The leaf motion module sets the leaf motion strategy of the multi-leaf collimator according to the filtered intensity matrix, and realizes the motion trajectory of the corresponding leaf.

[0115] This embodiment also provides a computer-readable storage medium storing a computer-executable program. When the computer-executable program is executed, it realizes the radiotherapy plan optimization method as described above.

[0116] The computer-readable storage medium described in this embodiment may include a data signal included in a baseband or propagated as part of a carrier wave, which carries the readable program code. Such a propagated data signal may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the foregoing. The computer-readable storage medium may also be any readable medium other than the readable storage medium, and this computer-readable medium may send, propagate, or transmit a program for use by or in conjunction with an instruction execution system, apparatus, or device. The program code contained on the computer-readable storage medium may be transmitted using any appropriate medium, including but not limited to wireless, wired, optical fiber cable, RF, etc., or any suitable combination of the foregoing.

[0117] This embodiment also provides an electronic device, including a processor and a memory, where the memory is used to store a computer-executable program. When the computer program is executed by the processor, the processor executes the method for optimizing a radiotherapy plan.

[0118] The electronic device is presented in the form of a general-purpose computing device. The processor may be one or multiple and work in cooperation. The present invention does not exclude distributed processing, that is, the processors may be dispersed in different physical devices. The electronic device of the present invention is not limited to a single entity, and may also be the sum of multiple physical devices.

[0119] The memory stores a computer-executable program, usually machine-readable code. The computer-readable program may be executed by the processor so that the electronic device can execute the method of the present invention or at least part of the steps in the method.

[0120] The memory includes a volatile memory, such as a random access storage unit (RAM) and / or a cache storage unit, and may also be a non-volatile memory, such as a read-only storage unit (ROM).

[0121] It should be understood that the electronic device of the present invention may also include elements or components not shown in the above examples. For example, some electronic devices also include a display unit such as a display screen, and some electronic devices also include human-computer interaction elements, such as buttons, keyboards, etc. As long as the electronic device can execute the computer-readable program in the memory to implement the method of the present invention or at least part of the steps of the method, it can be considered as the electronic device covered by the present invention.

[0122] Through the above description of the embodiments, those skilled in the art can easily understand that the present invention can be implemented by hardware capable of executing a specific computer program, such as the system of the present invention, and the electronic processing unit, server, client, mobile phone, control unit, processor, etc. included in the system. The present invention can also be implemented by computer software that executes the method of the present invention, such as control software executed by a microprocessor, an electronic control unit, a client, a server, etc. However, it should be noted that the computer software that executes the method of the present invention is not limited to being executed in one or specific hardware entities, and it can also be implemented in a distributed manner by unspecified specific hardware. For computer software, the software product can be stored in a computer-readable storage medium (which can be a CD-ROM, USB flash drive, mobile hard disk, etc.), or can be distributed and stored on a network, as long as it can enable an electronic device to execute the method according to the present invention.

[0123] The above embodiments are only used to illustrate the present invention and do not limit the technical solutions described in the present invention. Although the present specification has described the present invention in detail with reference to the above respective embodiments, the present invention is not limited to the above specific embodiments. Therefore, any modification or equivalent replacement of the present invention; and all technical solutions and their improvements that do not depart from the spirit and scope of the invention are covered by the scope of the claims of the present invention.

Claims

1. A method for optimizing a radiotherapy treatment plan, characterized in that, Including: Obtain the initial intensity matrix of the radiotherapy beam; Calculate according to the initial intensity matrix to obtain the irradiation dose MU of each sub - beam under this beam; Perform a descending - order sorting operation on the irradiation dose MUs of each sub - beam, and perform a filtering operation on the rows of the irradiation dose MUs in the first N% before sorting, where N is a set value; The calculating according to the initial intensity matrix to obtain the irradiation dose MU of each sub - beam under this beam includes: According to the initial intensity matrix, compare the intensity value of the element [i] with the intensity value of its left - hand element [i - 1] row by row; If the intensity value of the element [i] is less than the intensity value of the element [i - 1], record the intensity value of the element [i] as 0. If the intensity value of the element [i] is greater than or equal to the intensity value of the element [i - 1], record the intensity value of the element [i] as the difference between the two, to obtain a new intensity matrix; Traverse and accumulate the sum of the intensity values of each element in the new intensity matrix to obtain the irradiation dose MU of each corresponding sub - beam in each row; The performing a descending - order sorting operation on the irradiation dose MUs of each sub - beam and performing a filtering operation on the rows of the irradiation dose MUs in the first N% before sorting includes: Perform a descending - order sorting according to the irradiation dose MUs of each sub - beam, obtain the rows of the irradiation dose MUs in the first N% before sorting as the target rows for the filtering operation, and round up the total number of target rows.

2. The radiotherapy treatment plan optimization method according to claim 1, characterized in that, After obtaining the irradiation dose MU of each sub - beam under this beam, traverse the irradiation dose MUs of all sub - beams to obtain the maximum irradiation dose MU - max.

3. A radiotherapy plan optimization method according to claim 2, characterized in that Judge whether each target row in the initial intensity matrix meets the filtering operation strategy. If it meets, execute the corresponding filtering operation strategy. If it does not meet, take the average value of the intensities of the three adjacent elements [i - 1], [i], and [i + 1] in each target row of the initial intensity matrix as the new intensity value of the middle - position element [i].

4. A radiotherapy treatment plan optimization method according to claim 3, characterized in that The filtering operation strategy includes a first filtering operation strategy. Judging whether each target row in the initial intensity matrix meets the first filtering operation strategy includes: Judge whether there is an element [i] in each target row of the initial intensity matrix whose intensity value is equal to the maximum irradiation dose MU - max. If it exists, it meets the first filtering operation strategy, update the element [i] to k1*MU - max, and jump out of this loop, where k1 is a preset value.

5. A radiotherapy plan optimization method according to claim 4, characterized in that, The filtering operation strategy includes a second filtering operation strategy. If the judgment result of the first filtering operation strategy is non - existence, judge whether each target row in the initial intensity matrix meets the second filtering operation strategy. The second filtering operation strategy includes: Judge whether there is a situation where the intensity values of the elements [i - 1], [i], and [i + 1] in each target row of the initial intensity matrix are less than k2*MU - max. If it exists, it meets the second filtering operation strategy, and the intensity value of the element [i] remains unchanged, and jump out of this loop, where k2 is a preset value.

6. A radiotherapy treatment plan optimization method according to claim 5, characterized in that, The described filtering operation strategy includes a third filtering operation strategy. If the judgment result of the second filtering operation strategy is non - existent, it is judged whether each target row in the initial intensity matrix meets the third filtering operation strategy. The third filtering operation strategy includes: Judging whether each target row in the initial intensity matrix meets |diff1| > 1 / 2 * MU - max or |diff2| > 1 / 2 * MU - max, where diff1 = Original[i] - Original[i - 1]; diff2 = Original[i + 1] - Original[i], diff1 is the left difference, diff2 is the right difference, Original is the intensity matrix, and Original[i] is the intensity value of element [i] in the initial intensity matrix; If the judgment result is yes, it meets the third filtering operation strategy, the intensity value of element [i] remains unchanged, and this loop is exited. If the judgment result is no, it is determined that each target row in the initial intensity matrix does not meet the filtering operation strategy.

7. A radiotherapy treatment plan optimization method according to any one of claims 1-6, characterized in that, According to the intensity matrix processed by the filtering operation, execute the vane motion strategy of the multi - leaf collimator: For each element [i] in the intensity matrix processed by the filtering operation, control the initial position of the vane to be in the closed state; Control the left vane to remain at the initial position without movement, and the right vane moves to the shaped sub - field position, and maintain the position states of the left and right vanes for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of element [i]; Control the right vane to remain at the initial position without movement, and the left vane moves to the shaped sub - field position, and maintain the position states of the left and right vanes for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of element [i]; Then the vane motion strategy for element [i] is completed.

8. A radiotherapy treatment plan optimization device, characterized in that, Including: An intensity acquisition module that acquires the initial intensity matrix under the radiotherapy field; An irradiation dose calculation module that calculates according to the initial intensity matrix to obtain the irradiation dose MU of each sub - field under this field; A filtering module that performs a descending order arrangement operation on the irradiation dose MU of each sub - field, and performs a filtering operation on the rows of the irradiation dose MU before sorting for the first N%, where N is a set value; The calculation according to the initial intensity matrix to obtain the irradiation dose MU of each sub - field under this field includes: According to the initial intensity matrix, compare the intensity values of element [i] and its left - hand element [i - 1] row by row; If the intensity value of element [i] is less than the intensity value of element [i - 1], record the intensity value of element [i] as 0. If the intensity value of element [i] is greater than or equal to the intensity value of element [i - 1], record the intensity value of element [i] as the difference between the two, to obtain a new intensity matrix; Traverse and accumulate the sum of the intensity values of each element in the new intensity matrix to obtain the irradiation dose MU corresponding to each row of the sub - field; The operation of performing a descending order arrangement on the irradiation dose MU of each sub - field and performing a filtering operation on the rows of the irradiation dose MU before sorting for the first N% includes: Perform a descending order sorting according to the irradiation dose MU of each sub - field, obtain the rows of the irradiation dose MU before sorting for the first N% as the target rows for the filtering operation, and round up the total number of target rows.

9. The radiotherapy treatment plan optimization device according to claim 8, characterized in that, Including a blade movement module: Execute the blade movement strategy of the multi-leaf collimator according to the intensity matrix processed by the filtering operation: For each element [i] in the intensity matrix processed by the filtering operation, control the initial position of the blade to be in the closed state; Control the left blade to remain stationary at the initial position, and the right blade moves to the shaped sub-field position. Keep the position states of the left and right blades for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of the element [i]; Control the right blade to remain stationary at the initial position, and the left blade moves to the shaped sub-field position. Keep the position states of the left and right blades for a certain period of time to obtain the irradiation dose MU corresponding to the intensity value of the element [i]; Then the blade movement strategy for the element [i] is completed.

Images