阅读说明：本技术 自适应漏射补偿dmlc强度序列化方法、系统、装置及存储介质 (Adaptive leaky-radiation compensation DMLC (digital multiplex liquid crystal) intensity serialization method, system, device and storage medium ) 是由 吴章文 勾成俊 侯氢 于 2020-12-24 设计创作，主要内容包括：本方面实施例公开了一种自适应漏射补偿DMLC强度序列化方法。所述自适应漏射补偿DMLC强度序列化方法包括如下步骤：设置N个不同角度的射野并对其进行优化,得出优化结果；根据所述优化结果获取给定的第i射野的注量率分布其中,x方向为MLC叶片运动方向,y方向与x方向垂直,形成笛卡尔坐标系；对给定的第y个叶片的注量率分布进行离散化,得到所述第i照射野的MLC漏射自适应后的子野分布序列。通过本发明实施例的自适应漏射补偿DMLC强度序列化方法可以得到所述第i照射野的MLC漏射自适应后的子野分布序列,用该序列进行的照射的剂量分布能最大程度地与优化的结果相符合。(The embodiment of the aspect discloses a DMLC (digital multiplex liquid crystal) intensity serialization method for self-adaptive leakage radiation compensation. The adaptive leaky-radiation compensation DMLC intensity serialization method comprises the following steps of: setting N radiation fields with different angles and optimizing the radiation fields to obtain an optimization result; obtaining the fluence rate distribution of the given ith radiation field according to the optimization result The X direction is the MLC leaf motion direction, and the Y direction is vertical to the X direction to form a Cartesian coordinate system; fluence rate distribution for a given y-th blade Discretizing to obtain the subfield distribution sequence of the ith irradiation field after MLC leakage self-adaption. Adaptive leaky shoot through embodiments of the present inventionThe compensation DMLC intensity serialization method can obtain the subfield distribution sequence after MLC leakage adaptation of the ith irradiation field, and the dose distribution of irradiation by using the sequence can be in accordance with the optimized result to the maximum extent.)

1. An adaptive leakage radiation compensation DMLC intensity serialization method is characterized by comprising the following steps of:

setting N radiation fields with different angles and optimizing the radiation fields to obtain an optimization result;

obtaining the fluence rate distribution of the given ith radiation field according to the optimization resultThe X direction is the MLC leaf motion direction, and the Y direction is vertical to the X direction to form a Cartesian coordinate system;

fluence rate distribution for a given y-th bladeDiscretizing to obtain the subfield distribution sequence of the ith irradiation field after MLC leakage self-adaption.

2. The adaptive leaky ray compensation DMLC intensity serialization method as claimed in claim 1, wherein fluence rate distribution for a given y-th lobeThe discretization method comprises the following steps:

S10：will be provided withGrading fluence rate, irradiating hop number MU each time, and reducing values of fluence rate distribution at all positions by MU to obtain new fluence rateRecordingThe position of which is less than or equal to 0 is stored in the leaf segment sequence groupUp toAll values in (1) are 0;

wherein the content of the first and second substances,

where Δ t denotes the time interval of the control points during the MLC motion, MU denotes the exposure for one time interval,is the accelerator dose rate;respectively representing the A-end leaf sequence and the B-end leaf sequence of the MLC;

s20: grouping leaf segmentsSorting according to the sequence from small to large to obtain the difference value of adjacent positions in the array, and when the difference value is larger than d_maxThen insert a group of closed leaf sequence, i.e. in several groups respectivelyAdding a same position value; wherein d is_maxRepresents the maximum distance that an MLC leaf can move within a time interval;

s30: repeating step S20 until all adjacent position differences in the array are less than or equal to d_maxThen, the initial sub-field sequence of the ith field is obtainedJ sub fields are total, and the irradiation hop count of each sub field is MU;

wherein d is_max＝V_max×Δt；

Wherein d is_maxRepresenting the maximum distance that the MLC leaf can move in a time interval, Vmax representing the maximum moving speed of the MLC leaf, and delta t representing the time interval of the control point when the MLC moves;

s40: calculating the sub-fieldResulting fluence rate distributionWherein the content of the first and second substances,

s50: and carrying out summation operation on the fluence rate distribution of the sequence to obtain the fluence rate distribution adaptive to the MLC miss-fire of the 0 th time of the ith field:

s60: computingAnddivergent ofWherein the content of the first and second substances,

s70: feeding back the divergence to the initial fluence rate distribution to obtain the fluence distribution of the ith field after the first MLC miss-fire adaptationWherein the content of the first and second substances,

s80: repeating the steps S10-S70 to obtain a fluence distribution sequence and a fluence rate divergence distribution sequence after the kth MLC leakage injection self-adaption of the ith field:

s90: as the number of k increases, the number of bits increases,will be smaller and smaller; and when the distribution reaches a preset value, obtaining the distribution sequence of the subfield of the ith irradiation field after MLC leakage self-adaption.

3. An adaptive leaky-radiation compensating DMLC intensity serialization system, comprising:

the field optimization module is used for setting N fields with different angles and optimizing the fields to obtain an optimization result;

a fluence rate distribution obtaining module for obtaining a fluence rate distribution of a given ith field according to the optimization resultThe X direction is the MLC leaf motion direction, and the Y direction is vertical to the X direction to form a Cartesian coordinate system;

and

fluence rate distribution discretization module for fluence rate distribution for a given y-th bladeDiscretizing to obtain the subfield distribution sequence of the ith irradiation field after MLC leakage self-adaption.

4. An adaptive leakage compensation DMLC intensity serialization apparatus, wherein the adaptive leakage compensation DMLC intensity serialization apparatus comprises a processor and a memory; the memory is configured to store instructions that, when executed by the processor, cause the apparatus to perform operations corresponding to the adaptive leaky-radiation-compensating DMLC intensity serialization method of claim 1 or 2.

5. A computer-readable storage medium storing computer instructions which, when read by a computer, cause the computer to perform the adaptive leaky ray compensating DMLC intensity serialization method as claimed in claim 1 or 2.

Technical Field

The invention relates to the technical field related to intensity modulated radiotherapy, in particular to a DMLC (digital multiplex liquid crystal) intensity serialization method, a system, a device and a storage medium for adaptive leaky ray compensation.

Background

The fixed field intensity modulated radiotherapy technology is mainly divided into fixed field static intensity modulation and fixed field dynamic intensity modulation. Static intensity modulation different subfields were formed with MLC and step illumination was performed (Stop and Shot). Static emphasis is characterized in that the MLC movement and beam exit do not occur simultaneously. The dynamic MLC intensity modulation is formed according to the superposition of different dynamic irradiation fields, and is different from the static intensity modulation in that rays are in a beam-emitting state simultaneously in the moving process of the leaves, and the dynamic MLC has the advantage of shortening the total treatment time. After obtaining a more ideal fluence rate distribution, since the blade leakage ray is not considered before obtaining the blade movement pattern, it is necessary to perform compensation irradiation on the leakage ray.

Disclosure of Invention

In order to achieve the purpose, the technical scheme provided by the invention is as follows:

an adaptive leaky-radiation compensation DMLC intensity serialization method, comprising the steps of:

setting N radiation fields with different angles and optimizing the radiation fields to obtain an optimization result;

obtaining the fluence rate distribution of the given ith radiation field according to the optimization resultThe X direction is the MLC leaf motion direction, and the Y direction is vertical to the X direction to form a Cartesian coordinate system;

fluence rate distribution for a given y-th bladeDiscretizing to obtain the subfield distribution sequence of the ith irradiation field after MLC leakage self-adaption.

According to a preferred embodiment of the invention, the fluence rate distribution for a given y-th bladeThe discretization method comprises the following steps:

s10: will be provided withGrading fluence rate, irradiating hop number MU each time, and reducing values of fluence rate distribution at all positions by MU to obtain new fluence rateRecordingThe position of which is less than or equal to 0 is stored in the leaf segment sequence groupUp toAll values in (1) are 0;

wherein the content of the first and second substances,

where Δ t denotes the time interval of the control points during the MLC motion, MU denotes the exposure for one time interval,is the accelerator dose rate;respectively representing the A-end leaf sequence and the B-end leaf sequence of the MLC;

s20: grouping leaf segmentsSorting according to the sequence from small to large to obtain the difference value of adjacent positions in the array, and when the difference value is larger than d_maxThen insert a group of closed leaf sequence, i.e. in several groups respectivelyAdding inA same position value; wherein d is_maxRepresents the maximum distance that an MLC leaf can move within a time interval;

s30: repeating step S20 until all adjacent position differences in the array are less than or equal to d_maxThen, the initial sub-field sequence of the ith field is obtainedJ sub fields are total, and the irradiation hop count of each sub field is MU;

wherein d is_max＝V_max×Δt；

Wherein d is_maxRepresenting the maximum distance that the MLC leaf can move in a time interval, Vmax representing the maximum moving speed of the MLC leaf, and delta t representing the time interval of the control point when the MLC moves;

s40: calculating the sub-fieldResulting fluence rate distributionWherein the content of the first and second substances,

s50: and carrying out summation operation on the fluence rate distribution of the sequence to obtain the fluence rate distribution adaptive to the MLC miss-fire of the 0 th time of the ith field:

s60: computingAnddivergent ofWherein the content of the first and second substances,

s70: feeding back the divergence to the initial fluence rate distribution to obtain the fluence distribution of the ith field after the first MLC miss-fire adaptationWherein the content of the first and second substances,

s80: repeating the steps S10-S70 to obtain a fluence distribution sequence and a fluence rate divergence distribution sequence after the kth MLC leakage injection self-adaption of the ith field:

s90: as the number of k increases, the number of bits increases,will be smaller and smaller; and when the distribution reaches a preset value, obtaining the distribution sequence of the subfield of the ith irradiation field after MLC leakage self-adaption.

An adaptive leaky-radiation compensating DMLC intensity serialization system, comprising:

the field optimization module is used for setting N fields with different angles and optimizing the fields to obtain an optimization result;

a fluence rate distribution obtaining module for obtaining a fluence rate distribution of a given ith field according to the optimization resultWherein, the x direction is the MLC leaf motion direction, the y direction and the x directionVertical to form a Cartesian coordinate system;

and

fluence rate distribution discretization module for fluence rate distribution for a given y-th bladeDiscretizing to obtain the subfield distribution sequence of the ith irradiation field after MLC leakage self-adaption.

An adaptive leaky shoot compensating DMLC intensity serializing device comprises a processor and a memory; the memory is configured to store instructions that, when executed by the processor, cause the apparatus to perform operations corresponding to the adaptive leaky emission compensation DMLC intensity serialization method as described in any one of the above.

A computer readable storage medium storing computer instructions which, when read by a computer, cause the computer to perform a method of adaptive leaky-emission compensation DMLC intensity serialization as claimed in any one of the preceding claims.

Compared with the prior art, the adaptive leaky-radiation compensation DMLC intensity serialization method has the following beneficial effects:

according to the adaptive leaky-radiation compensation DMLC intensity serialization method, N radiation fields with different angles are set and optimized, and an optimization result is obtained; and obtaining the fluence rate distribution of the given ith radiation field according to the optimization resultThe X direction is the MLC leaf motion direction, and the Y direction is vertical to the X direction to form a Cartesian coordinate system; last fluence Rate distribution for a given y bladeDiscretizing to obtain the subfield distribution sequence of the ith irradiation field after MLC leakage self-adaption. The dose distribution of the irradiation with this sequence can be maximally matched to the optimized resultAnd (6) mixing.

Additional features of the invention will be set forth in part in the description which follows. Additional features of some aspects of the invention will become apparent to those of ordinary skill in the art upon examination of the following description and accompanying drawings or may be learned by the manufacture or operation of the embodiments. The features of the present disclosure may be realized and attained by practice or use of various methods, instrumentalities and combinations of the specific embodiments described below.

Drawings

The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention without limiting the invention. Like reference symbols in the various drawings indicate like elements. Wherein the content of the first and second substances,

fig. 1 is a schematic diagram of an adaptive leaky-radiation-compensating DMLC intensity sequencing system according to some embodiments of the invention.

Detailed Description

The following description is presented to enable one of ordinary skill in the art to make and use the application and is provided in the context of a particular application and its requirements. Various modifications to the disclosed embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be applied to other embodiments and applications without departing from the spirit and scope of the present application. Thus, the present application is not limited to the embodiments shown, but is to be accorded the widest scope consistent with the claims.

The terminology used herein is for the purpose of describing particular example embodiments only and is not intended to limit the scope of the present application. As used herein, the singular forms "a", "an" and "the" may include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.

The features and characteristics of the present application, as well as the methods of operation and functions of the related elements of structure and the combination of parts and economies of manufacture, will become more apparent upon consideration of the following description of the drawings, which form a part hereof. It is to be understood, however, that the drawings are designed solely for the purposes of illustration and description and are not intended as a definition of the limits of the application.

Flow charts are used herein to illustrate operations performed by systems according to embodiments of the present application. It should be understood that the operations in the flow diagrams are not necessarily performed exactly in order. Rather, various steps may be processed in reverse order or simultaneously. Further, one or more other operations may be added to the flowchart. One or more operations may also be deleted from the flowcharts.

One aspect of the embodiments of the present invention discloses a method for DMLC intensity serialization with adaptive leaky ray compensation.

The adaptive leaky-radiation compensation DMLC intensity serialization method mainly comprises the following steps:

setting N radiation fields with different angles and optimizing the radiation fields to obtain an optimization result;

obtaining the fluence rate distribution of the given ith radiation field according to the optimization resultThe X direction is the MLC leaf motion direction, and the Y direction is vertical to the X direction to form a Cartesian coordinate system;

fluence rate distribution for a given y-th bladeDiscretizing to obtain the subfield distribution sequence of the ith irradiation field after MLC leakage self-adaption.

Specifically, N radiation fields with different angles are set, and the optimization result is obtained through optimization algorithm optimization according to the set optimization conditions. For a given ith field, after the optimization is completed, the original fluence rate distribution is obtainedThe x direction is the MLC leaf motion direction, and the y direction is perpendicular to the x direction to form a Cartesian coordinate system. For a given y-th blade, the fluence rate distribution isThe discretization is performed as follows.

The optimization algorithm can adopt the existing portal angle optimization method.

Wherein, for a given y-th blade, the fluence rate distribution isThe discretization method mainly comprises the following steps:

s10: will be provided withGrading fluence rate, irradiating hop number MU each time, and reducing values of fluence rate distribution at all positions by MU to obtain new fluence rateRecordingThe position of which is less than or equal to 0 is stored in the leaf segment sequence groupUp toAll values in (a) are 0.

S20: grouping leaf segmentsSorting according to the sequence from small to large, judging the difference value of adjacent positions in the array, and when the difference value is larger than d_maxThen insert a group of closed leaf sequence, i.e. in several groups respectivelyAdd an identical position value.

S30: the process in S20 is repeated until all the adjacent position differences in the array are less than or equal to d_maxThen, the initial sub-field sequence of the ith field is obtainedThe number of irradiation hops of each subfield is MU.The sequence of leaves at the A end and the sequence of leaves at the B end of the MLC are respectively shown.

d_max＝V_max×Δt

d_maxRepresents the maximum distance that the MLC leaf can move in a time interval, Vmax represents the maximum moving speed of the MLC leaf, delta t represents the time interval of a control point when the MLC moves, MU represents the irradiation amount of one time interval,is the accelerator dose rate.

S40: calculating the sub-fieldResulting fluence rate distribution

S50: and carrying out summation operation on the fluence rate distribution of the sequence to obtain the fluence rate distribution adaptive to the MLC miss-fire of the 0 th time of the ith field:

s60: due to the fact thatLeakage radiation of MLC is considered, so distributionAndthere is a difference, andthe divergence of the two is calculated:

s70: and feeding back the divergence to the initial fluence rate distribution to obtain the fluence distribution after the first MLC miss-fire self-adaptation of the ith field:

s80: continuing the process from S10 to S70 to obtain a fluence distribution sequence and a fluence rate divergence distribution sequence after the kth MLC leakage self-adaptation of the ith field:

s90: as the number of k increases, the number of bits increases,will be smaller and smaller. When it reaches a preset value, the product is obtainedAnd (4) a subfield distribution sequence after MLC leakage emission self-adaption from the ith irradiation field.

The dose distribution of the irradiation with this sequence can be maximally matched to the optimized result.

In another aspect of an embodiment of the invention, an adaptive leaky ray compensation DMLC intensity serialization system is disclosed.

As shown in fig. 1, the adaptive leaky-radiation-compensating DMLC intensity serialization system 10 includes:

the portal optimization module 100 is configured to set and optimize the portals at N different angles to obtain an optimization result;

a fluence rate distribution obtaining module 200 for obtaining a fluence rate distribution of a given ith field according to the optimization resultThe X direction is the MLC leaf motion direction, and the Y direction is vertical to the X direction to form a Cartesian coordinate system;

and

fluence Rate distribution discretization Module 300 for fluence Rate distribution for a given y-th bladeDiscretizing to obtain the subfield distribution sequence of the ith irradiation field after MLC leakage self-adaption.

In another aspect of the embodiments of the invention, an adaptive leaky ray compensation DMLC intensity serialization apparatus is disclosed.

The adaptive leakage radiation compensation DMLC intensity serialization device comprises a processor and a memory; the memory is configured to store instructions that, when executed by the processor, cause the apparatus to perform operations corresponding to the adaptive leaky-radiation compensation DMLC intensity serialization method as described above.

Yet another aspect of an embodiment of the present invention discloses a computer-readable storage medium. The storage medium stores computer instructions, and when the computer instructions in the storage medium are read by the computer, the computer executes the adaptive leakage compensation DMLC intensity serialization method as described above.

Having thus described the basic concept, it will be apparent to those skilled in the art that the foregoing detailed disclosure is to be considered merely illustrative and not restrictive of the broad application. Various modifications, improvements and adaptations to the present application may occur to those skilled in the art, although not explicitly described herein. Such modifications, improvements and adaptations are proposed in the present application and thus fall within the spirit and scope of the exemplary embodiments of the present application.

Also, this application uses specific language to describe embodiments of the application. Reference throughout this specification to "one embodiment," "an embodiment," and/or "some embodiments" means that a particular feature, structure, or characteristic described in connection with at least one embodiment of the present application is included in at least one embodiment of the present application. Therefore, it is emphasized and should be appreciated that two or more references to "an embodiment" or "one embodiment" or "an alternative embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment. Furthermore, some features, structures, or characteristics of one or more embodiments of the present application may be combined as appropriate.

Moreover, those skilled in the art will appreciate that aspects of the present application may be illustrated and described in terms of several patentable species or situations, including any new and useful combination of processes, machines, manufacture, or materials, or any new and useful improvement thereon. Accordingly, various aspects of the present application may be embodied entirely in hardware, entirely in software (including firmware, resident software, micro-code, etc.) or in a combination of hardware and software. The above hardware or software may be referred to as "data block," module, "" engine, "" unit, "" component, "or" system. Furthermore, aspects of the present application may be represented as a computer product, including computer readable program code, embodied in one or more computer readable media.

The computer storage medium may comprise a propagated data signal with the computer program code embodied therewith, for example, on baseband or as part of a carrier wave. The propagated signal may take any of a variety of forms, including electromagnetic, optical, etc., or any suitable combination. A computer storage medium may be any computer-readable medium that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code located on a computer storage medium may be propagated over any suitable medium, including radio, cable, fiber optic cable, RF, or the like, or any combination of the preceding.

Computer program code required for the operation of various portions of the present application may be written in any one or more programming languages, including an object oriented programming language such as Java, Scala, Smalltalk, Eiffel, JADE, Emerald, C + +, C #, VB.NET, Python, and the like, a conventional programming language such as C, Visualbasic, Fortran2003, Perl, COBOL2002, PHP, ABAP, a dynamic programming language such as Python, Ruby, and Groovy, or other programming languages, and the like. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user's computer through any network format, such as a Local Area Network (LAN) or a Wide Area Network (WAN), or the connection may be made to an external computer (for example, through the Internet), or in a cloud computing environment, or as a service, such as a software as a service (SaaS).

Additionally, the order in which elements and sequences of the processes described herein are processed, the use of alphanumeric characters, or the use of other designations, is not intended to limit the order of the processes and methods described herein, unless explicitly claimed. While various presently contemplated embodiments have been discussed in the foregoing disclosure by way of example, it should be understood that such detail is solely for that purpose and that the appended claims are not limited to the disclosed embodiments, but, on the contrary, are intended to cover all modifications and equivalent arrangements that are within the spirit and scope of the embodiments herein. For example, although the system components described above may be implemented by hardware devices, they may also be implemented by software-only solutions, such as installing the described system on an existing server or mobile device.

Similarly, it should be noted that in the preceding description of embodiments of the application, various features are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the embodiments. This method of disclosure, however, is not intended to require more features than are expressly recited in the claims. Indeed, the embodiments may be characterized as having less than all of the features of a single embodiment disclosed above.