Dynamic intensity modulation control method and device of multi-leaf collimator and storage medium

Abstract

The application relates to the technical field of intelligent control, and discloses a dynamic intensity modulation control method and device of a multi-leaf collimator and a storage medium. The method comprises the following steps: acquiring dynamic intensity modulation basic parameters and an original intensity matrix; determining the number of subgroups based on a maximum dose rate subgroup control model according to the dynamic intensity modulation basic parameters; performing subgroup processing on the original intensity matrix to obtain a current intensity matrix according to the number of subgroups; and performing dynamic intensity modulation control of the multi-leaf collimator according to the current intensity matrix. The number of subgroups is determined based on the maximum dose rate subgroup control model; the original intensity matrix is subjected to subgroup processing by using the number of subgroups, which facilitates the subsequent subfield segmentation in the dynamic intensity modulation control step, and finally the dynamic intensity modulation of the maximum dose rate can be realized. The dynamic intensity modulation control of the multi-leaf collimator is performed according to the intensity matrix subjected to the subgroup processing, and the execution efficiency of the dynamic intensity modulation of the multi-leaf collimator can be effectively improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of intelligent control, and in particular to a dynamic intensity control method and device of a multileaf collimator and a storage medium. BACKGROUND

[0002] The multileaf collimator is a mechanical moving part for generating a conformal radiation field, commonly known as a multileaf grating, a multileaf diaphragm, etc., and is widely used in the medical field. The multileaf collimator is composed of two groups of closely arranged leaves, each leaf is made of tungsten alloy and is in the form of a long strip, and is driven by a small motor.

[0003] According to the movement mode of the multileaf collimator, the multileaf collimator has two types of manual and electric. The manual multileaf collimator is driven by hand to adjust the profile of the radiation field; the electric multileaf collimator is driven by a computer to independently drive each leaf to move, so as to achieve dynamic or static shaping of the radiation field.

[0004] How to accurately control the multileaf collimator to implement dynamic intensity control in the most efficient way, i.e., to implement the dynamic intensity control plan with the maximum dose rate and the maximum movement speed at the same time, has become a technical problem to be solved. SUMMARY

[0005] The present application relates to the technical field of intelligent control, and in particular to a dynamic intensity control method and device of a multileaf collimator and a storage medium.

[0006] To solve the above technical problems, the present application provides a dynamic intensity control method of a multileaf collimator, comprising: acquiring dynamic intensity basic parameters and an original intensity matrix; determining the number of clusters based on a maximum dose rate cluster control model according to the dynamic intensity basic parameters; performing cluster processing on the original intensity matrix to obtain a current intensity matrix according to the number of clusters; and performing dynamic intensity control of the multileaf collimator according to the current intensity matrix.

[0007] To solve the above technical problems, the present application further provides a dynamic intensity control device of a multileaf collimator, comprising: a data acquisition module, a cluster number determination module, a cluster processing module and a dynamic intensity control module; the data acquisition module is used for acquiring dynamic intensity basic parameters and an original intensity matrix; the cluster number determination module is used for determining the number of clusters based on a maximum dose rate cluster control model according to the dynamic intensity basic parameters; the cluster processing module is used for performing cluster processing on the original intensity matrix to obtain a current intensity matrix according to the number of clusters; and the dynamic intensity control module is used for performing dynamic intensity control of the multileaf collimator according to the current intensity matrix.

[0008] To solve the above technical problems, the application further provides a computer readable storage medium, comprising instructions, which, when executed on a computer, cause the computer to perform the dynamic intensity modulation control method of the multileaf collimator.

[0009] To solve the above technical problems, the application further provides a dynamic intensity modulation control device of a multileaf collimator, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the dynamic intensity modulation control method of the multileaf collimator according to the program.

[0010] The application has the following beneficial effects: the application determines the number of clusters according to the acquired dynamic intensity modulation basic parameters and the pre-constructed maximum dose rate cluster control model; the original intensity matrix is processed by clustering according to the number of clusters, which facilitates the subsequent sub-field segmentation in the dynamic intensity modulation control step, and finally the dynamic intensity modulation of the maximum dose rate can be realized. The dynamic intensity modulation control of the multileaf collimator is performed according to the intensity matrix processed by clustering, which can effectively improve the execution efficiency of the dynamic intensity modulation of the multileaf collimator.

[0011] Additional aspects and advantages of the application will be set forth in part in the description which follows, and in part will become apparent to those skilled in the art upon examination of the following and / or can be learned by practice of the application. BRIEF DESCRIPTION OF DRAWINGS

[0012] Figure 1 A dynamic intensity modulation control method flow chart of a multileaf collimator is provided for the embodiments of the application;

[0013] Figure 2 A functional module structure block diagram of a dynamic intensity modulation control device of a multileaf collimator is provided for the embodiments of the application;

[0014] Figure 3 A structure schematic diagram of a dynamic intensity modulation control device of a multileaf collimator is provided for the embodiments of the application. DETAILED DESCRIPTION

[0015] The embodiments of the present application will be described in detail below with specific reference to the drawings. Those skilled in the art can easily understand other advantages and effects of the present application from the content disclosed in the specification. Obviously, the described embodiments are only a part of the embodiments of the present application, not all. The present application can also be implemented or applied by other different specific embodiments, and the details in the specification can be modified or changed based on different views and applications without departing from the spirit of the present application. It should be noted that the following embodiments and features in the embodiments can be combined with each other without conflict. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0016] It is to be understood that the embodiments described hereinbelow within the scope of the appended claims. It will be apparent to one of ordinary skill in the art that aspects described herein can be implemented in a wide variety of forms, and that any specific structure and / or function described herein is merely illustrative. An aspect described herein can be implemented alone or in combination with any other aspect(s). Two or more aspects described herein can be implemented in any combination. For example, an apparatus can be implemented using any number of the aspects described herein. Additionally, an apparatus can be implemented using other structures and / or functionality in addition to or other than one or more of the aspects described herein.

[0017] As Figure 1 shown, the dynamic intensity modulation control of the multi-leaf collimator provided by the embodiments of the present application comprises the following steps:

[0018] S1, obtaining dynamic intensity modulation basic parameters and an original intensity matrix.

[0019] The dynamic intensity modulation basic parameters can include a maximum leaf motion speed, an irradiation dose of each sub-field, a maximum dose rate, a grid size, a number of non-zero sub-fields, and a grouping parameter. The original intensity matrix refers to an intensity matrix to be optimized. The dynamic intensity modulation basic parameters and the original intensity matrix described above can be pre-configured in an optimization configuration file.

[0020] S2, determining a number of groupings based on a maximum dose rate grouping control model according to the dynamic intensity modulation basic parameters.

[0021] S3, performing grouping processing on the original intensity matrix to obtain a current intensity matrix according to the number of groupings.

[0022] S4, performing dynamic intensity modulation control of the multi-leaf collimator according to the current intensity matrix.

[0023] The embodiments of the present application determine the number of groupings according to the obtained dynamic intensity modulation basic parameters and the pre-constructed maximum dose rate grouping control model; perform grouping processing on the original intensity matrix using the number of groupings, which facilitates sub-field segmentation in a subsequent dynamic intensity modulation control step, and finally realizes dynamic intensity modulation of the maximum dose rate. Performing dynamic intensity modulation control of the multi-leaf collimator according to the intensity matrix after the grouping processing can effectively improve the execution efficiency of the dynamic intensity modulation of the multi-leaf collimator.

[0024] Optionally, the number of groupings is determined based on a maximum dose rate grouping control model according to the dynamic intensity modulation basic parameters, and the formula is as follows:

[0025]

[0026] Wherein, f is the number of groups, Ratio mu is the maximum dose rate, which is a fixed value, unit is mu / min, generally takes the value of 600 mu / min; V_left is the maximum speed of leaf movement, which is a fixed value, unit is mm / s, generally takes the value of 25 mm / s; L Step is the grid size, which is a fixed value, unit is mm, generally takes the value of 2.5 mm; Mu is the total dose of each field, which is a fixed value, unit is mu, generally takes the value of 1-200 mu; Segment_num ≠0 is the number of non-zero segments, which is a constant value, generally takes the value of 1-300, when the irradiated tissue or organ is different, the value of the number of non-zero segments is different; Tau is the group parameter, which is a constant value; when the irradiated tissue or organ is different, the value of the group parameter is different, generally takes the value of 0-10.

[0027] The determination process of the above maximum dose rate group control model is as follows:

[0028] 1) Determine the relationship between the total number of segments and the number of groups, the formula is as follows:

[0029] Segment_num=Segment_num ≠0 +f*tau (1);

[0030] Wherein, Segment_num is the total number of segments, Segment_mum ≠0 is the number of non-zero segments, which is a constant value, f is the number of groups, tau is the group parameter, which is a constant value, it can be seen that Segment_num is positively correlated with f.

[0031] 2) Determine the maximum speed of leaf movement V_leaf and the grid size L Step .

[0032] 3) Determine the total dose of each field Mu, and Mu is a fixed value.

[0033] The calculation formula of the maximum dose rate Ratio mu is as follows:

[0034] Ratio_mu=(Mu / Segment_num) / (L Step / V_leaf) (2);

[0035] Substitute formula (1) into (2) to obtain the relationship between the number of groups f and the maximum dose rate Ratio mu, that is:

[0036] Ratio_mu=[Mu / (Segment_num ≠0 +f*tau)] / (L Step / V_leaf) (3)

[0037] The group number calculation formula can be obtained by transforming formula (3) as follows:

[0038]

[0039] The embodiment of the application pre-constructs a maximum dose rate group control model, and after obtaining the optimized configuration file, the group number can be quickly and accurately determined based on the maximum dose rate group control model. The group method can be applied to the optimization of the intensity matrix of different tissues or organs.

[0040] Optionally, the original intensity matrix is grouped according to the group number to obtain a current intensity matrix B x,y , and the calculation formula is as follows:

[0041] B x,y =A x,y / f;

[0042] wherein, A x,y is the original intensity matrix, f is the group number; x = 1, 2, 3…n, n is the total number of rows of the original intensity matrix; y = 1, 2, 3…m, m is the total number of columns of the original intensity matrix;

[0043] In the above embodiment, the original intensity matrix is grouped according to the group number to obtain the current intensity matrix, and a simplified intensity matrix can be obtained, which can effectively improve the accuracy and efficiency of sub-field segmentation.

[0044] Optionally, the dynamic intensity modulation control of the multi-leaf collimator is performed according to the current intensity matrix, comprising:

[0045] S41, determining the maximum sub-field number according to the beam intensity in the grid of the transition intensity matrix.

[0046] S42, determining a plurality of sub-fields corresponding to the current intensity matrix according to the maximum sub-field number and the beam intensity in each grid of the current intensity matrix.

[0047] The motion trajectory of each leaf of the multi-leaf collimator is determined according to the generation order of the plurality of sub-fields.

[0048] In the above embodiment, the maximum sub-field number is determined according to the beam intensity in the grid of the transition intensity matrix, so that the total number of sub-fields that need to be split according to the maximum sub-field number is determined; and then the plurality of sub-fields corresponding to the current intensity matrix are determined according to the maximum sub-field number and the beam intensity in each grid of the current intensity matrix, so as to ensure accurate irradiation to protect the critical organ and accurate irradiation of the target area.

[0049] The application will be described in detail below with reference to a specific example.

[0050] 1) Assume initial intensity matrix A x,y For

[0051]

[0052] 2) Calculate the number of segments according to the segment number calculation formula:

[0053]

[0054] Where, L Step = 2.5mm; Mu = 10mu;

[0055] Ratio_mu = 600mu / min; V_leaf = 25mm / s; Segment_num ≠0 = 1; Tau = 0.9.

[0056] Calculate f = 10, i.e. the number of segments is 10.

[0057] 3) Group the original intensity matrix matrix by the number of segments to obtain the current intensity matrix B x,y .

[0058]

[0059] 4) Calculate the number of subfields corresponding to each row of the current intensity matrix; in the embodiment of the application, the number of subfields determined according to the first row elements of the current intensity matrix is 5, and the number of subfields determined according to the second row elements is also 5. The subfield segmentation calculation process of calculating the number of subfields corresponding to each row of the intensity matrix is a conventional technical means in the art, and will not be described here.

[0060] 5) Select the maximum value from the number of subfields corresponding to each row of the current intensity matrix as the maximum number of subfields. In the embodiment of the application, the maximum number of subfields Sum_segment is 5.

[0061] 6) Determine a plurality of subfields corresponding to the current intensity matrix according to the maximum number of subfields and the beam intensity in each grid of the current intensity matrix.

[0062] Determine a plurality of subfields corresponding to the maximum number of subfields according to the beam intensity in each grid of the current intensity matrix. The subfield segmentation calculation process of determining the subfield according to the intensity matrix is a conventional technical means in the art, and will not be described here. The plurality of subfields are as follows:

[0063] Initial subfield:

[0064] First intermediate subfield:

[0065] Second intermediate subfield:

[0066] Third intermediate sub-field:

[0067] Termination sub-field:

[0068] 7) determining the motion trajectory of each leaf of the multileaf collimator according to the generation order of the plurality of sub-fields.

[0069] As Figure 2 shown, the dynamic intensity control device of the multileaf collimator provided by the embodiment of the present application comprises a data acquisition module 201, a group number determination module 201, a group processing module 203 and a dynamic intensity control module 204.

[0070] The data acquisition module 201 is used to acquire dynamic intensity basic parameters and an original intensity matrix; the group number setting module 202 is used to determine the group number based on the maximum dose rate group control model according to the dynamic intensity basic parameters; the group processing module 203 is used to perform group processing on the original intensity matrix to obtain a current intensity matrix according to the group number; and the dynamic intensity control module 204 is used to perform dynamic intensity control of the multileaf collimator according to the current intensity matrix.

[0071] The embodiment of the present application further provides a computer readable storage medium, which comprises instructions, when the instructions are run on a computer, the computer performs the dynamic intensity control method of the multileaf collimator provided by the above embodiment.

[0072] As Figure 3 shown, the embodiment of the present application further provides a dynamic intensity control device of a multileaf collimator 3000, which comprises a processor 3001, a memory 3003, a computer program stored in the memory 3003 and capable of running on the processor 30001, and the processor 3001 performs the dynamic intensity control method of the multileaf collimator provided by the above embodiment when the program is executed.

[0073] The processor 3001 and the memory 3003 are connected, such as being connected through a bus 3002. Optionally, the electronic device 3000 can further comprise a transceiver 3003, which can be used for data interaction between the electronic device and other electronic devices, such as data sending and / or data receiving. It should be noted that the transceiver 3003 is not limited to one in actual application, and the structure of the electronic device 3000 does not constitute a limitation on the embodiment of the present application.

[0074] The processor 3001 can be a CPU (Central Processing Unit), a general-purpose processor, a DSP (Digital Signal Processor), an ASIC (Application Specific Integrated Circuit), an FPGA (Field Programmable Gate Array) or other programmable logic device, transistor logic device, hardware component, or any combination thereof. It can implement or execute various exemplary logical blocks, modules and circuits described in connection with the present disclosure. The processor 3001 can also be a combination of computing functions, such as a combination of one or more microprocessors, a combination of a DSP and a microprocessor, and the like.

[0075] The bus 3002 can include a path for transmitting information between the above-mentioned components. The bus 3002 can be a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus, or the like. The bus 3002 can be divided into an address bus, a data bus, a control bus, and the like. For convenience of representation, Figure 3 In the figure, only one thick line is used to represent the bus, but it does not mean that there is only one bus or only one type of bus.

[0076] The memory 3003 can be a ROM (Read Only Memory) or other type of static storage device that can store static information and instructions, a RAM (Random Access Memory) or other type of dynamic storage device that can store information and instructions, an EEPROM (Electrically Erasable Programmable Read Only Memory), a CD-ROM (Compact Disc Read Only Memory) or other optical disk storage, a magnetic disk storage medium or other magnetic storage device, or any other medium that can be used to carry or store desired program code in the form of instructions or data structures and that can be accessed by a computer, but is not limited thereto.

[0077] The memory 3003 is configured to store application code (computer program) for implementing the scheme of the present application, and the processor 3001 is configured to control the execution. The processor 3001 is configured to execute the application code stored in the memory 3003 to implement the content shown in the foregoing method embodiments.

[0078] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described devices and units can refer to the corresponding processes in the foregoing method embodiments, which will not be described here.

[0079] In several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic, for example, the division of units is only a logical function division, and actual implementation can have another division manner, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed.

[0080] The units described as separate components can or can not be physically separate, and the components shown as units can or can not be physical units, that is, they can be located in one place, or can be distributed on a plurality of network units. Part or all of the units can be selected according to actual needs to achieve the purpose of the embodiments of the present application.

[0081] In addition, the functional units in each embodiment of the present application can be integrated in one processing unit, or each unit can be physically present alone, or two or more units can be integrated in one unit. The integrated unit can be realized in the form of hardware or in the form of a software functional unit.

[0082] If the integrated unit is realized in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part that contributes to the prior art, or all or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the embodiments of the method of the present application. The foregoing storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

[0083] The above merely describes preferred embodiments of the present application, and is not used to limit the present application, and any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A dynamic intensity modulation control method for a multi-leaf collimator, characterized by, The method comprises the following steps: acquiring dynamic intensity modulation basic parameters and an original intensity matrix; determining the number of clusters based on a maximum dose rate cluster control model according to the dynamic intensity modulation basic parameters; the number of clusters is determined based on the maximum dose rate cluster control model according to the dynamic intensity modulation basic parameters, and the formula is as follows: ； wherein f is the number of subgroups, is the maximum dose rate, which is a fixed value, and the unit is mu / min; is the maximum leaf movement speed, which is a fixed value, and the unit is mm / s; L Step is the grid size, which is a fixed value, and the unit is mm; Mu is the total dose of each field, which is a fixed value, and the unit is mu; is the number of non-zero subfields, which is a constant value, is the sub-group parameter, which is a constant value; grouping processing is performed on the original intensity matrix according to the number of clusters to obtain a current intensity matrix; dynamic intensity modulation control of a multileaf collimator is performed according to the current intensity matrix; the dynamic intensity modulation control of the multileaf collimator according to the current intensity matrix comprises the following steps: determining the maximum number of sub-fields according to the beam intensity in a grid of the current intensity matrix; determining a plurality of sub-fields corresponding to the current intensity matrix according to the maximum number of sub-fields and the beam intensity in each grid of the current intensity matrix; determining the motion trajectory of each leaf of the multileaf collimator according to the generation order of the plurality of sub-fields.

2. The method of claim 1, wherein, The current intensity matrix B is obtained by grouping the original intensity matrix according to the group number x,y The formula is as follows: B x,y =A x,y / f; wherein, A x,y is the original intensity matrix, f is the number of groups; x = 1, 2, 3…n, n is the total number of rows of the original intensity matrix; y = 1, 2, 3…m, m is the total number of columns of the original intensity matrix.

3. The method of claim 1, wherein, the maximum number of sub-fields is determined according to the beam intensity in a grid of the current intensity matrix, and the method comprises the following steps: calculating the number of sub-fields corresponding to each row of the current intensity matrix; selecting the maximum value from the number of sub-fields corresponding to each row of the current intensity matrix as the maximum number of sub-fields.

4. A dynamic intensity modulation control device for a multi-leaf collimator, characterized by, The method comprises the following steps: a data acquisition module is configured to acquire dynamic intensity modulation basic parameters and an original intensity matrix; a cluster number setting module is configured to determine the number of clusters based on a maximum dose rate cluster control model according to the dynamic intensity modulation basic parameters; the cluster number setting module determines the number of clusters according to the following formula: ； wherein f is the number of subgroups, is the maximum dose rate, which is a fixed value, and the unit is mu / min; is the maximum leaf movement speed, which is a fixed value, and the unit is mm / s; L Step is the grid size, which is a fixed value, and the unit is mm; Mu is the total dose of each field, which is a fixed value, and the unit is mu; is the number of non-zero subfields, which is a constant value, is the sub-group parameter, which is a constant value; a grouping processing module is configured to perform grouping processing on the original intensity matrix according to the number of clusters to obtain a current intensity matrix; a dynamic intensity modulation control module is configured to perform dynamic intensity modulation control of a multileaf collimator according to the current intensity matrix; the dynamic intensity modulation control of the multileaf collimator according to the current intensity matrix comprises the following steps: determining the maximum number of sub-fields according to the beam intensity in a grid of the current intensity matrix; determining a plurality of sub-fields corresponding to the current intensity matrix according to the maximum number of sub-fields and the beam intensity in each grid of the current intensity matrix; determining the motion trajectory of each leaf of the multileaf collimator according to the generation order of the plurality of sub-fields.

5. The apparatus of claim 4, wherein, The grouping processing module obtains the current intensity matrix B according to the following formula x,y The formula is as follows: B x,y =A x,y / f; wherein, A x,y is the original intensity matrix, f is the number of groups; x = 1, 2, 3…n, n is the total number of rows of the original intensity matrix; y = 1, 2, 3…m, m is the total number of columns of the original intensity matrix.

6. A computer-readable storage medium comprising instructions, characterized in that, When the instructions run on the computer, the computer performs the dynamic intensity modulation control method of the multileaf collimator according to any one of claims 1 to 3.

7. A dynamic intensity modulation control device for a multi-leaf collimator, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein, The processor performs the dynamic intensity modulation control method of the multileaf collimator according to any one of claims 1 to 3 when executing the program.

Images