CN103520843B

CN103520843B - A kind of Ziye processing method of reverse intensity modulated radiation therapy and device

Info

Publication number: CN103520843B
Application number: CN201310409693.5A
Authority: CN
Inventors: 姚鹏; 鞠光亮
Original assignee: Neusoft Medical Systems Co Ltd
Current assignee: Shenyang Neusoft Zhirui Radiotherapy Technology Co ltd
Priority date: 2013-09-10
Filing date: 2013-09-10
Publication date: 2016-03-30
Anticipated expiration: 2033-09-10
Also published as: CN103520843A

Abstract

The invention discloses a kind of Ziye processing method and device of reverse intensity modulated radiation therapy, described method comprises: from all matrixes be partitioned in advance, extract the Ziye matrix that machine jumping figure is identical; One or more the first Ziye matrix setting threshold value and the second Ziye matrix being less than setting threshold value of being greater than is chosen from the Ziye matrix of described extraction, and reduce Ziye area corresponding to described first Ziye matrix and increase Ziye area corresponding to described second Ziye matrix, be all more than or equal to setting threshold value with the Ziye area of the second Ziye matrix of the first Ziye matrix and described increase that make described reduction; Wherein, the Ziye area that described first Ziye matrix is corresponding is greater than Ziye area corresponding to described second Ziye matrix, the equal and opposite in direction of the Ziye area of described minimizing and the Ziye area of described increase and position is identical.

Description

A kind of Ziye processing method of reverse intensity modulated radiation therapy and device

Technical field

The present invention relates to field of medical technology, particularly relate to a kind of Ziye processing method and device of reverse intensity modulated radiation therapy.

Background technology

The vitality of tumor is eager to excel than normal organ usually, if high-energy ray can damage tumor, the injury so caused normal organ will be larger, so the principle formulating radiotherapy planning is exactly allow tumor receive dosage as much as possible and to allow the organ of surrounding receive the least possible dosage simultaneously.At present, the technology realizing this function is reverse intensity modulated radiation therapy (IMRT, Intensity-modulatedradiationtherapy), and it utilizes the specific region in computer-controlled accelerator whereabouts malignant tumor or tumor to launch accurate radiation dose.IMRT be exactly by doctor according to the input of the situation of patient oneself think suitable patient dose distribution and calculated by computer, finally draw the optimized algorithm that can realize doctor's target, this requires to optimize dose distribution of drawing and accelerator, and to form the dose distribution that dose distribution will wish with doctor consistent as far as possible.Be implemented as follows:

First, will speed up device irradiation field to be divided into one by one " square " according to identical size, each " square " represents the weight of a ray pencil beam, and each pencil beam can produce dose contribution to all organs, and the dosage of each point (" square " that representative marks off) is:

D_{i} = Σ_{k = 1}^{N_{beamlet}} B_{ik} w_{k} (i = 1,2 . . . N_{point}) - - - (1)

In formula (1), D _ithe dosage of i-th point, B _ikbe a kth pencil beam to the dose contribution of i-th point, w _kthe dose weights of a kth pencil beam, thus the dosage of each point be all pencil beams to the contribution of this dose point and, and w _kbe exactly the amount that IMRT algorithm needs to optimize, w _kthe matrix formed is exactly flux matrix to be split.

Then, the dose constraint that doctor inputs is transformed into object function:

f = \frac{1}{N_{T}} Σ_{i = 1}^{N_{T}} {(D_{i} - D_{0}^{T})}^{2} + \frac{1}{N_{OAR}} Σ_{i = 1}^{N_{OAR}} H (D_{i} - D_{0}^{OAR}) {(D_{i} - D_{0}^{OAR})}^{2} - - - (2)

In formula (2),

H (D_{i} - D_{0}^{OAR}) = \{\begin{matrix} 1 & D_{i} > D_{0}^{OAR} \\ 0 & D_{i} \leq D_{0}^{OAR} \end{matrix},

N _tthe quantity of target area (tumor region) mid point, N _oARthe number of critical organ region (organ of tumor periphery) mid point, D _ithe dosage of i-th point, the target dose of target area, it is the target dose of critical organ.It should be noted that, actual object function should be the object function of all organs and, formula (2) merely illustrates the object function sum of two organs (tumor and a critical organ).　

usually can compare much larger, could realize like this dosage of target area large as far as possible while make the dosage of other critical organ little as far as possible, obviously, if in target area dosage a little more close the value of first functional expression is less, exceedes in critical organ the dosage of point less, the value of second function is less, if whole object function is 0, then illustrates that all points meet dose constraint completely.Therefore, in order to make D _iclose need to use optimized algorithm (as gradient method, simulated annealing etc.) to draw the weight w of the pencil beam of optimization _k, and by w _kthe flux matrix formed carries out splitting to form each Ziye, and segmentation Ziye is the matrix that described flux matrix will be decomposed into some same sizes, such as:

(\begin{matrix} 1 & 3 & 3 & 0 \\ 0 & 2 & 4 & 1 \\ 1 & 1 & 4 & 4 \\ 3 & 3 & 1 & 0 \end{matrix}) = 2 (\begin{matrix} 0 & 1 & 1 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 1 & 1 \\ 1 & 1 & 0 & 0 \end{matrix}) + (\begin{matrix} 0 & 1 & 1 & 0 \\ 0 & 1 & 1 & 1 \\ 1 & 1 & 1 & 1 \\ 0 & 0 & 0 & 0 \end{matrix}) + (\begin{matrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 1 & 1 \\ 0 & 0 & 1 & 1 \\ 1 & 1 & 1 & 0 \end{matrix})

As Fig. 1, show the accelerator launched field that each matrix above is corresponding.Divided matrix number represents the number of Ziye under certain direction." 1 " in each Ziye matrix represents does not have blade to block (representative has roentgenization target area) herein, and " 0 " represents has blade to block (representative does not have roentgenization target area) herein.Need to follow following constraints when splitting Ziye:

1, the region that 1 in each Ziye matrix is formed will be connected (can ensure accuracy of dose) as far as possible;

2, the number of Ziye matrix the least possible (radiation duration can be reduced);

3, all Ziye matrixes factor in front (i.e. machine jumping figure " 2,1,1 ") is that be added and little as far as possible (can reduce radiation duration and wear and tear of machinery);

4, the number of " 1 " many as far as possible (can accuracy of dose be ensured) in each subdivision matrix.　

In order to constraints more than As soon as possible Promising Policy, nearly all partitioning algorithm is all make the Ziye area being partitioned into matrix large as far as possible, but do like this and the number of " 1 " in last remaining Ziye matrix will be caused fewer and feweri, last several Ziye matrixes may only have little several " 1 ", such as following matrix:

(\begin{matrix} 1 & 1 & 1 & 1 \\ 0 & 1 & 1 & 0 \\ 0 & 2 & 1 & 0 \\ 1 & 1 & 1 & 1 \end{matrix}) = (\begin{matrix} 1 & 1 & 1 & 1 \\ 0 & 1 & 1 & 0 \\ 0 & 1 & 1 & 0 \\ 1 & 1 & 1 & 1 \end{matrix}) + (\begin{matrix} 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{matrix})

As Fig. 2, show the accelerator launched field that each matrix above is corresponding.For two Ziyes on the right in Fig. 2, a rear Ziye directly can be neglected when controlling accelerator exposure radiation dose, because this Ziye area usually very little (a few mm in width), forms so little launched field by the leaf position controlling accelerator and almost can not realize; Namely the field area that enable realization is so little, the dosage produced also is completely can be uncared-for; And for the launched field of small size like this, the precision of Rapid Dose Calculation is difficult to ensure, the error between the dosage namely calculated and actual exposure dose is very large.But, if given up by this small size launched field, IMRT accuracy of dose will be reduced.

Summary of the invention

In view of this, the main purpose of the embodiment of the present invention is the Ziye processing method and the device that provide a kind of reverse intensity modulated radiation therapy, to realize the object improving IMRT accuracy of dose.

For achieving the above object, embodiments provide a kind of Ziye area equalization methods of reverse intensity modulated radiation therapy, comprising:

The Ziye matrix that machine jumping figure is identical is extracted from all matrixes be partitioned in advance;

One or more the first Ziye matrix setting threshold value and the second Ziye matrix being less than setting threshold value of being greater than is chosen from the Ziye matrix of described extraction, and reduce Ziye area corresponding to described first Ziye matrix and increase Ziye area corresponding to described second Ziye matrix, be all more than or equal to setting threshold value with the Ziye area of the second Ziye matrix of the first Ziye matrix and described increase that make described reduction;

Wherein, the Ziye area that described first Ziye matrix is corresponding is greater than Ziye area corresponding to described second Ziye matrix, the equal and opposite in direction of the Ziye area of described minimizing and the Ziye area of described increase and position is identical.

Preferably, described choosing from the Ziye matrix of described extraction is one or morely greater than the first Ziye matrix of setting threshold value and is less than the second Ziye matrix of setting threshold value, and reduces Ziye area corresponding to described first Ziye matrix and the Ziye area increasing described second Ziye matrix corresponding comprises:

Determine the Ziye area that in the Ziye matrix of described extraction, each Ziye matrix is corresponding;

Described Ziye area is greater than the Ziye matrix of Ziye matrix as the first kind of setting threshold value, described Ziye area is less than or equal to the Ziye matrix of Ziye matrix as Second Type of setting threshold value;

From the Ziye matrix of the described first kind, choose the maximum Ziye matrix of Ziye area as the first Ziye matrix, and from the Ziye matrix of described Second Type, choose the minimum Ziye matrix of Ziye area as the second Ziye matrix;

Exchange the Ziye area of described first Ziye matrix and described second Ziye matrix line by line, and retain the matrix condition before exchange according to the area variance size of exchanging between front and back the first Ziye matrix and the second Ziye matrix or retain the matrix condition after exchanging;

Judge whether the Ziye area that the second Ziye matrix after exchanging line by line is corresponding is less than described setting threshold value, if so, then judges whether there is the Ziye matrix be not selected in the Ziye matrix of the described first kind;

If the Ziye matrix be not selected described in existing, from the described Ziye matrix be not selected, then choose the maximum Ziye matrix of Ziye area as the first Ziye matrix, continue to perform and exchange the step of the Ziye area of described first Ziye matrix and described second Ziye matrix line by line, until Ziye area corresponding to the second Ziye matrix is not less than described setting threshold value.

Preferably, the Ziye matrix be not selected if there is no and Ziye area corresponding to the second Ziye matrix is still less than described setting threshold value, then all Ziye matrixes of the first kind and the second Ziye matrix are all returned to divided rear time state.

Preferably, the described Ziye area exchanging described first Ziye matrix and described second Ziye matrix line by line, and retain the matrix condition before exchange according to the area variance size of exchanging between front and back the first Ziye matrix and the second Ziye matrix or retain the matrix condition after exchanging, comprising:

The current numerical value of going together of described first Ziye matrix and described second Ziye matrix is exchanged, to exchange Ziye area corresponding to described current colleague's numerical value;

Calculate the first area variance between Ziye area corresponding to the first Ziye matrix before the exchanging Ziye area corresponding with the second Ziye matrix, and calculate the second area variance between Ziye area corresponding to the first Ziye matrix after the exchanging Ziye area corresponding with the second Ziye matrix;

If described second area variance is less than described first area variance, and there is next numerical value of going together of described first Ziye matrix and described second Ziye matrix, then using next colleague's numerical value described as current colleague's numerical value, continue to perform the step that the described current numerical value of going together by described first Ziye matrix and described second Ziye matrix carries out exchanging;

If described second area variance is more than or equal to described first area variance, then described first Ziye matrix and described second Ziye matrix are returned to the state before exchange, after being successfully recovered, when next that there is described first Ziye matrix and described second Ziye matrix goes together numerical value, using next colleague's numerical value described as current colleague's numerical value, continue to perform the step that the described current numerical value of going together by described first Ziye matrix and described second Ziye matrix carries out exchanging.

Preferably, after the current numerical value of going together of described first Ziye matrix and described second Ziye matrix is exchanged, also comprise:

In any one matrix in described first Ziye matrix and described second Ziye matrix, judge that the right leaf position of the accelerator whether left leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping, and/or, judge that the left leaf position of the accelerator whether right leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping;

If, then perform the described step described first Ziye matrix and described second Ziye matrix being returned to the state before exchange, if not, then the step of the first area variance between Ziye area corresponding to the first Ziye matrix before the performing described calculating exchange Ziye area corresponding with the second Ziye matrix is continued;

When next that there is not described first Ziye matrix and described second Ziye matrix goes together numerical value, judge whether to exist the area after the exchange be retained, if, then using the first row numerical value of described first Ziye matrix and described second Ziye matrix as current numerical value of going together, continue to perform the step of being carried out by the current numerical value of going together of described first Ziye matrix and described second Ziye matrix exchanging.

The embodiment of the present invention additionally provides a kind of Ziye area shape method for repairing and mending of reverse intensity modulated radiation therapy, comprising:

The Ziye matrix identical to the machine jumping figure in all matrixes be partitioned in advance carries out area equilibrium;

In Ziye matrix after area equilibrium, choose each Ziye matrix successively, and carry out shape repair operation to choosing Ziye matrix;

Described shape repair operation specifically comprises:

Exchange and choose the part Ziye area of Ziye matrix and the 3rd Ziye matrix correspondence position, to reduce the complex-shaped degree choosing Ziye matrix and the 3rd Ziye matrix, described 3rd Ziye matrix be not with the described Ziye matrix chosen Ziye matrix and participate in exchanging;

A Ziye matrix is chosen as the 3rd Ziye matrix from described having neither part nor lot in the Ziye matrix of exchange, continue to perform described exchange and choose the step of the part Ziye area of Ziye matrix and the 3rd Ziye matrix correspondence position, until described in choose Ziye matrix and other each Ziye matrix all carries out exchange operation.

Preferably, the part Ziye area that Ziye matrix and the 3rd Ziye matrix correspondence position are chosen in described exchange comprises:

Choose Ziye matrix and in the 3rd Ziye matrix, choosing mutually corresponding at least one line number value, and exchange all or part of Ziye area corresponding to every line number value of choosing line by line;

After each exchange terminates, calculate the second shape complexity of the first shape complexity choosing Ziye matrix and the 3rd Ziye matrix before exchanging respectively, and choose the 3rd shape complexity of Ziye matrix and the 4th shape complexity of the 3rd Ziye matrix after calculating exchange respectively;

Judge whether the 3rd shape complexity and the 4th shape complexity sum are less than the first shape complexity and the second shape complexity sum, if, then retain the state of exchanging latter two Ziye matrix, if not, then two the Ziye matrixes participating in exchanging are reverted to the state before exchange.

Preferably, the computing formula of described shape complexity is:

C_{i} = Σ_{j = 2}^{R - 1} (I L_{j} + {IR}_{j}), i = 1, . . ., N;

{IL}_{j} = \{\begin{matrix} (l_{j} - l_{j - 1}) + (l_{j} - l_{j + 1}) & l_{j} > l_{j + 1}, l_{j} > l_{j - 1} \\ 0 & else \end{matrix}

{IR}_{j} = \{\begin{matrix} (r_{j - 1} - r_{j}) + (r_{j + 1} - r_{j}) & r_{j - 1} > r_{j}, r_{j + 1} > r_{j} \\ 0 & else \end{matrix}

Wherein, C _ifor shape complexity, R is the line number of Ziye matrix, l _i-1, l _j, l _i+1be respectively the lobus sinister sheet extended position of the i-th-1 row, the i-th row and the i-th+1 row; r _i-1, r _j, r _i+1be respectively the lobus dexter sheet extended position of the i-th-1 row, the i-th row and the i-th+1 row.

Preferably, the described all or part of Ziye area exchanging every line number value of choosing line by line corresponding comprises:

Exchange line by line choose Ziye matrix and the 3rd Ziye matrix current line numerical value to exchange Ziye area corresponding to current line numerical value;

Or, exchange line by line and choose first choose area and second choose area corresponding with the current line numerical value of the 3rd Ziye matrix of Ziye matrix, wherein, first chooses area chooses the equal and opposite in direction of area with second and position is identical, when described first to choose area be the blade openings area choosing Ziye matrix, described second to choose area be that the blade of the 3rd Ziye matrix is remained silent area, when described first choose area be the blade choosing Ziye matrix remain silent area time, described second chooses the blade openings area that area is the 3rd Ziye matrix.

Preferably, after described each exchange terminates, also comprise:

Choosing in any one matrix in Ziye matrix and the 3rd Ziye matrix, judge that the right leaf position of the accelerator whether left leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping, and/or, judge that the left leaf position of the accelerator whether right leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping;

If, then perform the described step two the Ziye matrixes participating in exchanging being reverted to the state before exchange, if not, then continue to perform the described step calculating the second shape complexity of the first shape complexity choosing Ziye matrix before exchanging and the 3rd Ziye matrix respectively.

The embodiment of the present invention additionally provides a kind of Ziye area balancer of reverse intensity modulated radiation therapy, comprising:

Matrix chooses module, for extracting the identical Ziye matrix of machine jumping figure from all matrixes be partitioned in advance;

Area balance module, for choosing one or more the first Ziye matrix setting threshold value and the second Ziye matrix being less than setting threshold value of being greater than from the Ziye matrix of described extraction, and reduce Ziye area corresponding to described first Ziye matrix and increase Ziye area corresponding to described second Ziye matrix, be all more than or equal to setting threshold value with the Ziye area of the second Ziye matrix of the first Ziye matrix and described increase that make described reduction; Wherein, the Ziye area that described first Ziye matrix is corresponding is greater than Ziye area corresponding to described second Ziye matrix, the equal and opposite in direction of the Ziye area of described minimizing and the Ziye area of described increase and position is identical.

Preferably, described area balance module comprises:

Area determining unit, for determine described extraction Ziye matrix in Ziye area corresponding to each Ziye matrix;

Type determining units, for described Ziye area being greater than the Ziye matrix of Ziye matrix as the first kind of setting threshold value, is less than or equal to the Ziye matrix of Ziye matrix as Second Type of setting threshold value using described Ziye area;

Matrix chooses unit, for choosing the maximum Ziye matrix of Ziye area as the first Ziye matrix from the Ziye matrix of the described first kind, and from the Ziye matrix of described Second Type, chooses the minimum Ziye matrix of Ziye area as the second Ziye matrix;

Exchange unit line by line, for exchanging the Ziye area of described first Ziye matrix and described second Ziye matrix line by line, and retain the matrix condition before exchange according to the area variance size of exchanging between front and back the first Ziye matrix and the second Ziye matrix or retain the matrix condition after exchanging;

First judging unit, for judging whether the Ziye area that the second Ziye matrix after exchanging line by line is corresponding is less than described setting threshold value;

Second judging unit, for judge to obtain at described first judging unit exchange line by line after Ziye area corresponding to the second Ziye matrix be less than described setting threshold value time, judge in the Ziye matrix of the described first kind, whether to there is the Ziye matrix be not selected;

Continue to choose unit, for be not selected described in the existence that judges to obtain at described second judging unit Ziye matrix time, the maximum Ziye matrix of Ziye area is chosen as the first Ziye matrix from the described Ziye matrix be not selected, described in continuing to utilize, exchange unit realizes exchanging line by line the function of the Ziye area of described first Ziye matrix and described second Ziye matrix line by line, until Ziye area corresponding to the second Ziye matrix is not less than described setting threshold value.

Preferably, described area balance module also comprises:

Recovering state unit, for described second judging unit judge obtain do not exist described in be not selected Ziye matrix time, if the Ziye area that the second Ziye matrix is corresponding is still less than described setting threshold value, then all Ziye matrixes of the first kind and the second Ziye matrix are all returned to divided rear time state.

Preferably, described exchange unit line by line, comprising:

Colleague exchanges subelement, is exchanged by the current numerical value of going together of described first Ziye matrix and described second Ziye matrix, to exchange Ziye area corresponding to described current colleague's numerical value;

Variance computation subunit, for calculating the first area variance between Ziye area corresponding to the first Ziye matrix before the exchange Ziye area corresponding with the second Ziye matrix, and calculate the second area variance between Ziye area corresponding to the first Ziye matrix after the exchanging Ziye area corresponding with the second Ziye matrix;

Variance judgment sub-unit, for judging whether described second area variance is less than described first area variance;

First chooses subelement, for when described variance judgment sub-unit judges that the described second area variance obtained is less than described first area variance, if there is next numerical value of going together of described first Ziye matrix and described second Ziye matrix, then using next colleague's numerical value described as current colleague's numerical value, continue to utilize described first colleague to exchange subelement and realize the function that the described current numerical value of going together by described first Ziye matrix and described second Ziye matrix carries out exchanging;

Recovering state subelement, for when described variance judgment sub-unit judges that the described second area variance obtained is more than or equal to described first area variance, returns to the state before exchange by described first Ziye matrix and described second Ziye matrix;

Second chooses subelement, after the state before described first Ziye matrix and described second Ziye matrix return to exchange, if there is next numerical value of going together of described first Ziye matrix and described second Ziye matrix, then using next colleague's numerical value described as current colleague's numerical value, continue to utilize described colleague to exchange subelement and realize the function that the described current numerical value of going together by described first Ziye matrix and described second Ziye matrix carries out exchanging.

Preferably, described exchange unit line by line also comprises:

First overlapping judgment sub-unit, for exchanging after the current numerical value of going together of described first Ziye matrix and described second Ziye matrix exchanges by subelement described colleague, in any one matrix in described first Ziye matrix and described second Ziye matrix, judge that the right leaf position of the accelerator whether left leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping, and/or, the left leaf position of the accelerator whether right leaf position of the accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping,

Described recovering state subelement, also for the result that judges to obtain in described first overlapping judgment sub-unit as being time, described first Ziye matrix and described second Ziye matrix are returned to the step of the state before exchange, continue to call described second and choose subelement;

First continues to perform subelement, when judging that the result obtained is no for described first overlapping judgment sub-unit, continue to utilize described variance computation subunit to realize the described function calculating the first area variance between Ziye area corresponding to the first Ziye matrix before the exchanging Ziye area corresponding with the second Ziye matrix;

Again verify subelement, for when next that there is not described first Ziye matrix and described second Ziye matrix goes together numerical value, judge whether to exist the area after the exchange be retained, if, then using the first row numerical value of described first Ziye matrix and described second Ziye matrix as current numerical value of going together, continue to call described colleague and exchange subelement.

The embodiment of the present invention additionally provides a kind of Ziye area shape repair apparatus of reverse intensity modulated radiation therapy, comprising:

The Ziye area balancer of reverse intensity modulated radiation therapy, carries out area equilibrium for the Ziye matrix identical to the machine jumping figure in all matrixes be partitioned in advance;

Shape repairs module, in the Ziye matrix after area equilibrium, chooses each Ziye matrix successively, and carries out shape repair operation to choosing Ziye matrix;

Described shape repairs module, specifically for exchanging the part Ziye area choosing Ziye matrix and the 3rd Ziye matrix correspondence position, to reduce the complex-shaped degree of Ziye matrix and the 3rd Ziye matrix chosen, described 3rd Ziye matrix be not with the described Ziye matrix chosen Ziye matrix and participate in exchanging; A Ziye matrix is chosen as the 3rd Ziye matrix from described having neither part nor lot in the Ziye matrix of exchange, continue to perform described exchange and choose the step of the part Ziye area of Ziye matrix and the 3rd Ziye matrix correspondence position, until described in choose Ziye matrix and other each Ziye matrix all carries out exchange operation.

Preferably, described shape is repaired module and is comprised area exchange unit, and described area exchange unit chooses the part Ziye area of Ziye matrix and the 3rd Ziye matrix correspondence position for exchanging;

Described area exchange unit comprises:

Exchange subelement line by line, for choosing Ziye matrix and choosing mutually corresponding at least one line number value in the 3rd Ziye matrix, and exchange all or part of Ziye area corresponding to every line number value of choosing line by line;

Complicated dynamic behaviour subelement, for after each exchange terminates, calculate the second shape complexity of the first shape complexity choosing Ziye matrix and the 3rd Ziye matrix before exchanging respectively, and choose the 3rd shape complexity of Ziye matrix and the 4th shape complexity of the 3rd Ziye matrix after calculating exchange respectively;

Complexity judgment sub-unit, for judging whether the 3rd shape complexity and the 4th shape complexity sum are less than the first shape complexity and the second shape complexity sum;

State retains subelement, when being less than the first shape complexity and the second shape complexity sum for the 3rd shape complexity that judges to obtain in described complexity judgment sub-unit and the 4th shape complexity sum, and the state of latter two Ziye matrix of reservation exchange;

State returns subelement, when being not less than the first shape complexity and the second shape complexity sum for the 3rd shape complexity that judges to obtain in described complexity judgment sub-unit and the 4th shape complexity sum, two the Ziye matrixes participating in exchanging are reverted to the state before exchange.

Preferably, the computing formula of described shape complexity is:

C_{i} = Σ_{j = 2}^{R - 1} (I L_{j} + {IR}_{j}), i = 1, . . ., N;

{IL}_{j} = \{\begin{matrix} (l_{j} - l_{j - 1}) + (l_{j} - l_{j + 1}) & l_{j} > l_{j + 1}, l_{j} > l_{j - 1} \\ 0 & else \end{matrix}

{IR}_{j} = \{\begin{matrix} (r_{j - 1} - r_{j}) + (r_{j + 1} - r_{j}) & r_{j - 1} > r_{j}, r_{j + 1} > r_{j} \\ 0 & else \end{matrix}

Preferably, describedly exchange subelement line by line, specifically for exchange line by line choose Ziye matrix and the 3rd Ziye matrix current line numerical value to exchange Ziye area corresponding to current line numerical value; Or, exchange line by line and choose first choose area and second choose area corresponding with the current line numerical value of the 3rd Ziye matrix of Ziye matrix, wherein, first chooses area chooses the equal and opposite in direction of area with second and position is identical, when described first to choose area be the blade openings area choosing Ziye matrix, described second to choose area be that the blade of the 3rd Ziye matrix is remained silent area, when described first choose area be the blade choosing Ziye matrix remain silent area time, described second chooses the blade openings area that area is the 3rd Ziye matrix.

Preferably, described area exchange unit also comprises:

Second overlapping judgment sub-unit, for after described each exchange terminates, choosing in any one matrix in Ziye matrix and the 3rd Ziye matrix, judge that the right leaf position of the accelerator whether left leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping, and/or the left leaf position of the accelerator whether right leaf position of the accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping;

Described state returns subelement, also for the result that judges to obtain in described second overlapping judgment sub-unit as being time, two the Ziye matrixes participating in exchanging are reverted to the state before exchange;

Second continues to perform subelement, when judging that the result obtained is no for described second overlapping judgment sub-unit, continue to utilize described complicated dynamic behaviour subelement to realize the described function calculating the second shape complexity of the first shape complexity choosing Ziye matrix before exchanging and the 3rd Ziye matrix respectively.

The Ziye processing method of the reverse intensity modulated radiation therapy that the embodiment of the present invention provides and device, the second Ziye matrix with less Ziye area is filled up by utilizing one or more first Ziye matrixes with larger Ziye area, make Ziye area corresponding to all Ziye matrixes all meet threshold requirement as far as possible, eliminate the Ziye matrix with less Ziye area, avoid the uncared-for phenomenon of Ziye that area corresponding in Ziye matrix is less, make the actual exposure dose of IMRT closer to plan exposure dose, thus improve IMRT accuracy of dose.In addition, after area equilibrium is carried out to Ziye area, also further obviate the situation of outstanding blade, thus further increase the accuracy of dose of IMRT.

Accompanying drawing explanation

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.

Fig. 1 is a kind of first schematic diagram of accelerator launched field;

Fig. 2 is a kind of second schematic diagram of accelerator launched field;

Fig. 3 is a kind of 3rd schematic diagram of accelerator launched field;

Fig. 4 is the schematic flow sheet of the Ziye area equalization methods of the reverse intensity modulated radiation therapy of the embodiment of the present invention;

Fig. 5 is a kind of idiographic flow schematic diagram of the Ziye area equalization methods of the reverse intensity modulated radiation therapy of the embodiment of the present invention;

Fig. 6 is the schematic diagram of exchange operating process line by line of the embodiment of the present invention first Ziye matrix and the second Ziye matrix;

Fig. 7 is the first schematic diagram of the overlapping situation of embodiment of the present invention blade;

Fig. 8 is the second schematic diagram of the overlapping situation of embodiment of the present invention blade;

Fig. 9 is the original segmentation schematic diagram of the embodiment of the present invention first;

Figure 10 is that first of embodiment of the present invention the first row square exchanges schematic diagram;

Figure 11 is that first of the embodiment of the present invention second row square exchanges schematic diagram;

Figure 12 is that first of embodiment of the present invention the third line square exchanges schematic diagram;

Figure 13 is that first of embodiment of the present invention fourth line square exchanges schematic diagram;

Figure 14 is the original segmentation schematic diagram of the embodiment of the present invention second;

Figure 15 is that second of embodiment of the present invention the first row square exchanges schematic diagram;

Figure 16 is that second of embodiment of the present invention the third line square exchanges schematic diagram;

Figure 17 is that second of embodiment of the present invention fourth line square exchanges schematic diagram;

Figure 18 is that the embodiment of the present invention gives prominence to blade schematic diagram;

Figure 19 is the schematic flow sheet of the Ziye area shape method for repairing and mending of the reverse intensity modulated radiation therapy of the embodiment of the present invention;

Figure 20 is the schematic flow sheet of embodiment of the present invention shape repair operation;

Figure 21 is that the embodiment of the present invention adopts exchange mode to realize the schematic flow sheet of shape repairing;

Figure 22 is the original segmentation schematic diagram of the embodiment of the present invention the 3rd;

Figure 23 is the rear segmentation figure of exchange of embodiment of the present invention partial blade position;

Figure 24 is the structural representation of the Ziye area balancer of the reverse intensity modulated radiation therapy of the embodiment of the present invention;

Figure 25 is the structural representation of the Ziye area shape repair apparatus of the reverse intensity modulated radiation therapy of the embodiment of the present invention.

Detailed description of the invention

For making the object of the embodiment of the present invention, technical scheme and advantage clearly, below in conjunction with the accompanying drawing in the embodiment of the present invention, technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.

In the prior art, after flux matrix being split to be divided into multiple Ziye matrix, there will be the Ziye matrix (namely in Ziye matrix the quantity of " 1 " very little) of corresponding less Ziye area, Given this, the main purpose of the embodiment of the present invention is that the number of " 1 " in each Ziye matrix is average as far as possible, affects IMRT accuracy of dose to avoid occurring " 1 " number Ziye matrix very little.The embodiment of the present invention mainly carries out area equilibrium to improve the object of IMRT accuracy of dose to each Ziye be partitioned into, that is, " 0 " in the Ziye matrix with less Ziye area is filled up in main utilization " 1 " had in the Ziye matrix of larger Ziye area.

Such as: first matrix below in matrix is flux matrix, actual flux matrix is generally the matrix of 20*20, but conveniently the flux matrix of illustrating in the embodiment of the present invention is the matrix of 4*4.Second matrix and the 3rd matrix are the Ziye matrix after utilizing existing Factorization algorithm algorithm to split flux matrix, visible, in 3rd matrix, " 1 " place is the Ziye area that this Ziye matrix is corresponding, because this Ziye area is very little, the dose irradiation at this Ziye area place will be left in the basket, thus cause the reduction of actual dose precision, in order to not lose the 3rd exposure dose that matrix is corresponding, " 1 " in second matrix and " 0 " at same position place in the 3rd matrix are exchanged, to increase the Ziye area of the 3rd matrix.

(\begin{matrix} 1 & 1 & 1 & 1 \\ 0 & 1 & 1 & 0 \\ 0 & 2 & 1 & 0 \\ 1 & 1 & 1 & 1 \end{matrix}) = (\begin{matrix} 1 & 1 & 1 & 1 \\ 0 & 1 & 1 & 0 \\ 0 & 1 & 1 & 0 \\ 1 & 1 & 1 & 1 \end{matrix}) + (\begin{matrix} 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{matrix}) + (\begin{matrix} 1 & 1 & 1 & 1 \\ 0 & 1 & 1 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{matrix}) + (\begin{matrix} 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 1 & 1 & 0 \\ 1 & 1 & 1 & 1 \end{matrix})

After carrying out exchange operation, second matrix and the 3rd matrix distinguish the 4th matrix after corresponding exchange and the 5th matrix, shown in Figure 3, the block chart on the left side is that the 4th roentgenization that matrix is corresponding is wild, the block chart on the right is that the 5th roentgenization that matrix is corresponding is wild, wherein, white represents the launched field part not having blade to block, and black represents the part not having ray to pass through having blade to block.

Embodiment one

Shown in Figure 4, show the schematic flow sheet of the Ziye area equalization methods of the reverse intensity modulated radiation therapy that the embodiment of the present invention provides, specifically comprise:

Step 401: extract the Ziye matrix that machine jumping figure is identical from all matrixes be partitioned in advance.

Wherein, described machine jumping figure be each divided go out the coefficient of Ziye matrix.

Step 402: choose one or more the first Ziye matrix setting threshold value and the second Ziye matrix being less than setting threshold value of being greater than from the Ziye matrix of described extraction, and reduce Ziye area corresponding to described first Ziye matrix and increase Ziye area corresponding to described second Ziye matrix, be all more than or equal to setting threshold value with the Ziye area of the second Ziye matrix of the first Ziye matrix and described increase that make described reduction; Wherein, the Ziye area that described first Ziye matrix is corresponding is greater than Ziye area corresponding to described second Ziye matrix, the equal and opposite in direction of the Ziye area of described minimizing and the Ziye area of described increase and position is identical.

It should be noted that, first the Ziye area determining each Ziye matrix that machine jumping figure is identical is needed, wherein, the Ziye area being less than setting threshold value may be left in the basket, in order to ensure the accuracy of dose of IMRT, the embodiment of the present invention under the prerequisite not changing total radiation dose, by changing Ziye area corresponding to each Ziye matrix, the Ziye area making each Ziye matrix corresponding is more than or equal to setting threshold value, and each Ziye area is not all left in the basket.Namely, reduce Ziye area corresponding to one or more first Ziye matrix and ensure reduce after Ziye area be not less than setting threshold value prerequisite under, utilize the Ziye area reduced to fill up to be less than the Ziye area corresponding to the second Ziye matrix of setting threshold value, thus make Ziye area corresponding to the second Ziye matrix be more than or equal to setting threshold value, to ensure that its field area is not left in the basket.

Embodiment two

Shown in Figure 5, show a kind of idiographic flow schematic diagram of the Ziye area equalization methods of the reverse intensity modulated radiation therapy that the embodiment of the present invention provides, specifically comprise:

Step 501: extract the Ziye matrix that machine jumping figure is identical from all matrixes be partitioned in advance.

Step 502: determine the Ziye area that in the Ziye matrix of described extraction, each Ziye matrix is corresponding.

Step 503: the Ziye matrix of Ziye matrix as the first kind described Ziye area being greater than setting threshold value, is less than or equal to the Ziye matrix of Ziye matrix as Second Type of setting threshold value using described Ziye area.

Step 504: choose the maximum Ziye matrix of Ziye area as the first Ziye matrix from the Ziye matrix of the described first kind, and from the Ziye matrix of described Second Type, choose the minimum Ziye matrix of Ziye area as the second Ziye matrix.

Step 505: the Ziye area exchanging described first Ziye matrix and described second Ziye matrix line by line, and retain the matrix condition before exchange according to the area variance size of exchanging between front and back the first Ziye matrix and the second Ziye matrix or retain the matrix condition after exchanging.

Step 506: judge whether the Ziye area that the second Ziye matrix after exchanging line by line is corresponding is less than described setting threshold value, if so, then perform step 507, if not, perform step 510.

After utilizing the maximum Ziye area of a Ziye matrix (the first Ziye matrix) to fill up the minimum Ziye area of another Ziye matrix (the second Ziye matrix), the Ziye area that second Ziye matrix is corresponding also may not meet area threshold requirement, namely the second Ziye matrix also small size Ziye matrix after filling up, in actual use, the irradiation launched field that second Ziye matrix is corresponding is still left in the basket, thus reduces the accuracy of dose of IMRT.

Step 507: judge whether there is the Ziye matrix be not selected in the Ziye matrix of the described first kind, if so, then performs step 508, if not, then performs step 509.

Step 508: choose the maximum Ziye matrix of Ziye area as the first Ziye matrix from the described Ziye matrix be not selected, continues to perform step 505.

Step 509: all Ziye matrixes of the first kind and the second Ziye matrix are all returned to divided rear time state.

Step 510: process ends.

For the embodiment shown in Fig. 5, the Ziye matrix of the first kind may have one or more, and the Ziye matrix of Second Type also may have one or more.For the Ziye matrix not meeting the Second Type setting threshold value, first the embodiment of the present invention is selected to fill up minimum Ziye area, namely the Ziye matrix of the first kind is chosen successively from big to small according to the size of Ziye area, utilize the Ziye matrix selected to fill up the Ziye matrix of the minimum Second Type of Ziye area successively, after satisfied setting threshold requirement, retain the state of each Ziye matrix after filling up, after filling up, still not meet threshold requirement, then each matrix is returned to the state before not carrying out area equilibrium.For the Ziye matrix of other Second Type, carrying out minimum Ziye area after area fills up, continuing the Ziye matrix selecting the less Second Type of next Ziye area, carrying out area equilibrium according to the method shown in Fig. 5 equally.

The embodiment of the present invention mainly realizes step 505 by the following method, and the first Ziye matrix shown in Figure 6 and the schematic diagram of exchange operating process line by line of the second Ziye matrix, specifically comprise:

Step 601: exchanged by the current numerical value of going together of described first Ziye matrix and described second Ziye matrix, to exchange Ziye area corresponding to described current colleague's numerical value.

Step 602: in any one matrix in described first Ziye matrix and described second Ziye matrix, judge whether to occur the overlapping situation of blade, if so, then performs step 606, if not, then performs step 603.

Wherein, following three kinds of modes are had to judge whether to occur the overlapping situation of blade:

If the leaf position that the current line numerical value in 1 first Ziye matrix or the second Ziye matrix is corresponding only has lobus sinister sheet, then judge that the right leaf position of the accelerator whether left leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping;

If the leaf position that the current line numerical value in 2 first Ziye matrixes or the second Ziye matrix is corresponding only has lobus dexter sheet, then judge that the left leaf position of the accelerator whether right leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping;

If the existing lobus sinister sheet of leaf position that the current line numerical value in 3 first Ziye matrixes or the second Ziye matrix is corresponding has again lobus dexter sheet, then judge that the right leaf position of the accelerator whether left leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping, and judge that the left leaf position of the accelerator whether right leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping.

First schematic diagram of the overlapping situation of blade shown in Figure 7, circled portion in the drawings, there is mutual overlapping situation with lobus dexter sheet below in lobus sinister sheet above; Second schematic diagram of the overlapping situation of blade shown in Figure 8, circled portion in the drawings, there is mutual overlapping situation with lobus sinister sheet below in lobus dexter sheet above.

Step 603: calculate the first area variance between Ziye area corresponding to the first Ziye matrix before the exchanging Ziye area corresponding with the second Ziye matrix, and calculate the second area variance between Ziye area corresponding to the first Ziye matrix after the exchanging Ziye area corresponding with the second Ziye matrix.

Step 604: judge whether described second area variance is less than described first area variance, if so, then performs step 605, if not, then performs step 606.

If after exchanging, the area variance of two Ziye matrixes diminishes, and the object reaching area equilibrium is described.

Step 605: judge whether next numerical value of going together that there is described first Ziye matrix and described second Ziye matrix, if so, then performs step 607, if not, then performs step 608.

Step 606: described first Ziye matrix and described second Ziye matrix are returned to the state before exchange, performs step 605.

Step 607: using next colleague's numerical value described as current colleague's numerical value, continue to perform step 601, if not, then perform step 608.

Step 608: judge whether to exist the area after the exchange be retained, if so, then performs step 609, if not, then performs step 610.

Step 609: using the first row numerical value of described first Ziye matrix and described second Ziye matrix as current numerical value of going together, performs step 601.

Step 610: process ends.

The Ziye area equalization methods of the reverse intensity modulated radiation therapy that the embodiment of the present invention provides, the second Ziye matrix with less Ziye area is filled up by utilizing one or more first Ziye matrixes with larger Ziye area, make Ziye area corresponding to all Ziye matrixes all meet threshold requirement as far as possible, eliminate the Ziye matrix with less Ziye area, avoid the uncared-for phenomenon of less Ziye area corresponding in Ziye matrix, make the actual exposure dose of IMRT closer to plan exposure dose, thus improve IMRT accuracy of dose.

In order to understand figure embodiment two of the present invention more easily, illustrate below:

Example 1: the shown in Figure 9 first original segmentation schematic diagram, two segmentation figure in Fig. 9, after utilizing existing partitioning algorithm to be split by 4*4 flux matrix, the segmentation figure that the Ziye matrix that two the machine jumping figures formed are identical is corresponding respectively, wherein, in 16 squares in each segmentation figure, the numerical value " 0 " of the corresponding Ziye matrix of each black bars, the position and the corresponding blade of black bars is remained silent, the numerical value " 1 " of the corresponding Ziye matrix of each white square, and the corresponding blade openings position of white square.But, the corresponding small size launched field of the right segmentation figure, its Ziye area is less than setting threshold value, for needing to carry out segmentation figure corresponding to the second Ziye matrix that area fills up, because the Ziye area of left side segmentation figure is greater than setting threshold value, herein, the Ziye area (white square) shown in left side segmentation figure is utilized to fill up the Ziye area (white square) shown in the right segmentation figure.

Line by line segmentation figure corresponding for two Ziye matrixes is exchanged in the following manner below:

The first step, first of the first row square shown in Figure 10 exchanges schematic diagram, first, the first row of two segmentation figure in original segmentation schematic diagram is exchanged, then, calculate the Ziye area variance of these two segmentation figure, because the Ziye area variance after exchanging is less than the Ziye area variance before exchange, now retain the state after exchanging.

Second step, first of second row square shown in Figure 11 exchanges schematic diagram, first, the second row in two segmentation figure in original segmentation schematic diagram is exchanged, then, calculate the Ziye area variance of these two segmentation figure, because the Ziye area variance after exchanging is less than the Ziye area variance before exchange, now retain the state after exchanging.

3rd step, first of the third line square shown in Figure 12 exchanges schematic diagram, first, the third line in two segmentation figure in original segmentation schematic diagram is exchanged, then, calculate the Ziye area variance of these two segmentation figure, because the Ziye area variance after exchanging is greater than the Ziye area variance before exchange, now do not retain the state after exchange, and two segmentation figure are reverted to the state shown in Figure 11.

4th step, first of fourth line square shown in Figure 13 exchanges schematic diagram, first, fourth line in two segmentation figure in original segmentation schematic diagram is exchanged, then, calculate the Ziye area variance of these two segmentation figure, because the Ziye area variance after exchanging is greater than the Ziye area variance before exchange, now do not retain the state after exchange, and two segmentation figure are reverted to the state shown in Figure 11.

In whole exchange process in above-mentioned example 1, does not occur the overlapping situation of blade, illustrate below and when belong to the overlapping situation of blade, and in the overlapping situation of appearance, the state after exchange is not retained:

Example 2: the shown in Figure 14 second original segmentation schematic diagram, wherein, the corresponding small size launched field of the right segmentation figure, its Ziye area is less than setting threshold value, for needing to carry out segmentation figure corresponding to Ziye matrix that area fills up, the Ziye area of left side segmentation figure is greater than setting threshold value, utilizes the Ziye area (white square) shown in left side segmentation figure to fill up the Ziye area (white square) shown in the right segmentation figure herein.Wherein, black part is divided into blade to block position, and white portion on-bladed blocks.

In exchange operating process, once there is following situation, current exchange operation is not retained:

Second of the first row square shown in Figure 15 exchanges schematic diagram, the first row in two segmentation figure in original segmentation schematic diagram is exchanged, visible, in right side segmentation figure, if the second row blade is lobus sinister sheet, then there is the overlapping situation of blade in the second row blade and the first row blade, if the second row blade is lobus dexter sheet, then there is the overlapping situation of blade in the second row blade and the third line blade, now do not retain the state after exchange.

Second of the third line square shown in Figure 16 exchanges schematic diagram, the third line in two segmentation figure in original segmentation schematic diagram is exchanged, visible, in left side segmentation figure, there is the overlapping situation of blade in the second row lobus dexter sheet and the third line lobus sinister sheet, does not now retain the state after exchange.

Second of fourth line square shown in Figure 17 exchanges schematic diagram, fourth line in two segmentation figure in original segmentation schematic diagram is exchanged, visible, there is not the overlapping situation of blade, and the area variance of two segmentation figure after exchanging is greater than the area variance before exchange, now retains the state after exchanging.

Embodiment three

If the numerical value of certain a line is mostly " 0 " in Ziye matrix, and the periphery numerical value of this line number value is " 1 ", then there will be the situation of the outstanding blade be circled shown in similar Figure 18, but due to the impact of the factors such as ray scattering, the described outstanding blade be circled can not stop injecting of ray completely, such launched field is also the accuracy of dose that can have influence on IMRT, namely plan is irradiated precision and is irradiated precision be not inconsistent with actual, so when integrating Ziye matrix, be also noted that situation about reducing as far as possible as shown in figure 18.

Based on this, the embodiment of the present invention also proposed a kind of Ziye area shape method for repairing and mending of reverse intensity modulated radiation therapy, is eliminated the situation of above-mentioned outstanding blade by the method as far as possible, to promote the object of IMRT accuracy of dose further.

Shown in Figure 19, show the schematic flow sheet of the Ziye area shape method for repairing and mending of the reverse intensity modulated radiation therapy that the embodiment of the present invention provides, specifically comprise:

Step 1901: the Ziye matrix identical to the machine jumping figure in all matrixes be partitioned in advance carries out area equilibrium.

Wherein, above-mentioned area equalization methods is adopted to carry out area equilibrium to described Ziye matrix.

Step 1902: in the Ziye matrix after area equilibrium, chooses each Ziye matrix successively, and carries out shape repair operation to choosing Ziye matrix.

" carrying out shape repair operation to choosing Ziye matrix " that the embodiment of the present invention specifically adopts mode below to realize in step 1902:

The schematic flow sheet of shape repair operation shown in Figure 20, described shape repair operation specifically comprises:

Step 2001: exchange the part Ziye area choosing Ziye matrix and the 3rd Ziye matrix correspondence position, to reduce the complex-shaped degree of Ziye matrix and the 3rd Ziye matrix chosen, described 3rd Ziye matrix be not with the described Ziye matrix chosen Ziye matrix and participate in exchanging.

Step 2002: judge whether to exist the Ziye matrix having neither part nor lot in exchange, if so, then performs step 2003, if not, then performs step 2004.

Step 2003: choose a Ziye matrix as the 3rd Ziye matrix from described having neither part nor lot in the Ziye matrix of exchange, continues to perform step 2001.

Step 2004: process ends.

The embodiment of the present invention mainly adopts following method to realize above-mentioned steps 2001, and ginseng employing as shown in Figure 21 mode of exchanging realizes the schematic flow sheet that shape is repaired, and specifically comprises:

Step 2101: choosing Ziye matrix and choosing mutually corresponding at least one line number value in the 3rd Ziye matrix, and exchange all or part of Ziye area corresponding to every line number value of choosing line by line.

Wherein, described in can choosing, choose every line number value of Ziye matrix and described 3rd Ziye matrix, or selected part line number value.

For " exchanging all or part of Ziye area that every line number value of choosing is corresponding line by line " in step 2101, below can adopting, wherein a kind of mode realizes, or first adopt first kind of way to adopt the second way again, or first adopt the second way to adopt first kind of way again:

First kind of way: exchange line by line choose Ziye matrix and the 3rd Ziye matrix current line numerical value to exchange Ziye area corresponding to current line numerical value;

The second way: exchange line by line and choose first choose area and second choose area corresponding with the current line numerical value of the 3rd Ziye matrix of Ziye matrix, wherein, first chooses area chooses the equal and opposite in direction of area with second and position is identical, when described first to choose area be the blade openings area choosing Ziye matrix, described second to choose area be that the blade of the 3rd Ziye matrix is remained silent area, when described first choose area be the blade choosing Ziye matrix remain silent area time, described second chooses the blade openings area that area is the 3rd Ziye matrix.

In order to understand the second way more easily, illustrate below:

Ginseng the 3rd original segmentation schematic diagram as shown in Figure 22, suppose that the segmentation figure on the left side is for choosing Ziye matrix, the segmentation figure on the right side is the 3rd Ziye matrix, known by detecting, significantly outstanding blade is had in the segmentation figure on the right side, in order to eliminate this outstanding blade, this part outstanding blade i.e. this part blade can be remained silent the blade openings location swap of correspondence position in area and left side segmentation figure, second and the 3rd square by the third line in two segmentation figure exchange, thus eliminate outstanding blade, be the segmentation figure after the exchange of partial blade position see Figure 23.Step 2102: after each exchange terminates, choosing in any one matrix in Ziye matrix and the 3rd Ziye matrix, judging whether to occur the overlapping situation of blade, if so, then performs step 2106, if not, then performs step 2103.

If 1 leaf position choosing current line numerical value in Ziye matrix or the 3rd Ziye matrix corresponding only has lobus sinister sheet, then judge that the right leaf position of the accelerator whether left leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping;

If 2 leaf positions choosing current line numerical value in Ziye matrix or the 3rd Ziye matrix corresponding only have lobus dexter sheet, then judge that the left leaf position of the accelerator whether right leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping;

If the 3 existing lobus sinister sheets of leaf position choosing current line numerical value in Ziye matrix or the 3rd Ziye matrix corresponding have again lobus dexter sheet, then judge that the right leaf position of the accelerator whether left leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping, and judge that the left leaf position of the accelerator whether right leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping.

Step 2103: the second shape complexity calculating the first shape complexity choosing Ziye matrix and the 3rd Ziye matrix before exchanging respectively, and choose the 3rd shape complexity of Ziye matrix and the 4th shape complexity of the 3rd Ziye matrix after calculating exchange respectively.

Wherein, the first shape complexity, the second shape complexity, the 3rd shape complexity and the 4th shape complexity all adopt following formulae discovery:

C_{i} = Σ_{j = 2}^{R - 1} (I L_{j} + {IR}_{j}), i = 1, . . ., N;

{IL}_{j} = \{\begin{matrix} (l_{j} - l_{j - 1}) + (l_{j} - l_{j + 1}) & l_{j} > l_{j + 1}, l_{j} > l_{j - 1} \\ 0 & else \end{matrix}

{IR}_{j} = \{\begin{matrix} (r_{j - 1} - r_{j}) + (r_{j + 1} - r_{j}) & r_{j - 1} > r_{j}, r_{j + 1} > r_{j} \\ 0 & else \end{matrix}

Wherein, C _ifor shape complexity, R is the line number of Ziye matrix, l _i-1, l _j, l _i+1be respectively the lobus sinister sheet extended position of the i-th-1 row, the i-th row and the i-th+1 row; r _i-1, r _j, r _i+1be respectively the lobus dexter sheet extended position of the i-th-1 row, the i-th row and the i-th+1 row.Wherein, described lobus sinister sheet extended position and lobus dexter sheet extended position all calculate position from the left side, and certainly, described lobus sinister sheet extended position and lobus dexter sheet extended position also can calculate position from the right, and now, above-mentioned formula carries out the amendment of correspondence.

Step 2104: judge whether the 3rd shape complexity and the 4th shape complexity sum are less than the first shape complexity and the second shape complexity sum, if so, then performs step 2105, if not, then performs step 2106.

Step 2105: retain the state of exchanging latter two Ziye matrix.

Step 2106: two the Ziye matrixes participating in exchanging are reverted to the state before exchange.

The Ziye area shape method for repairing and mending of the reverse intensity modulated radiation therapy that the embodiment of the present invention provides, after carrying out area equilibrium to Ziye area, also further obviates the situation of outstanding blade, thus further increases the accuracy of dose of IMRT.It should be noted that, this Ziye area shape method for repairing and mending can be implemented after carrying out area equilibrium, also directly can repair the Ziye matrix be partitioned into, not limit during specific implementation.

Embodiment four

See Figure 24, show the structural representation of the Ziye area balancer of reverse intensity modulated radiation therapy, this device comprises:

Matrix chooses module 251, for extracting the identical Ziye matrix of machine jumping figure from all matrixes be partitioned in advance;

Area balance module 252, for choosing one or more the first Ziye matrix setting threshold value and the second Ziye matrix being less than setting threshold value of being greater than from the Ziye matrix of described extraction, and reduce Ziye area corresponding to described first Ziye matrix and increase Ziye area corresponding to described second Ziye matrix, be all more than or equal to setting threshold value with the Ziye area of the second Ziye matrix of the first Ziye matrix and described increase that make described reduction; Wherein, the Ziye area that described first Ziye matrix is corresponding is greater than Ziye area corresponding to described second Ziye matrix, the equal and opposite in direction of the Ziye area of described minimizing and the Ziye area of described increase and position is identical.

Wherein, described area balance module 252 comprises:

Matrix chooses unit, for choosing the maximum Ziye matrix of Ziye area as the first Ziye matrix from the Ziye matrix of the described first kind, and from the Ziye matrix of described Second Type, chooses the minimum Ziye matrix of Ziye area as the second Ziye matrix.

Exchange unit line by line, for exchanging the Ziye area of described first Ziye matrix and described second Ziye matrix line by line, and retain the matrix condition before exchange according to the area variance size of exchanging front and back the first Ziye matrix and the second Ziye matrix or retain the matrix condition after exchanging;

Wherein, exchange unit line by line described above, comprising:

In addition, described exchange unit line by line also comprises:

First continues to perform subelement, when judging that the result obtained is no for described first overlapping judgment sub-unit, continue to utilize described variance computation subunit to realize the described function calculating the first area variance between Ziye area corresponding to the first Ziye matrix before the exchanging Ziye area corresponding with the second Ziye matrix.

The Ziye area balancer of the reverse intensity modulated radiation therapy that the embodiment of the present invention provides, the second Ziye matrix with less Ziye area is filled up by utilizing one or more first Ziye matrixes with larger Ziye area, make Ziye area corresponding to all Ziye matrixes all meet threshold requirement as far as possible, eliminate the Ziye matrix with less Ziye area, to avoid in Ziye matrix corresponding to the uncared-for phenomenon of the Ziye of small size, make the actual exposure dose of IMRT closer to plan exposure dose, thus improve IMRT accuracy of dose.

Embodiment five

See Figure 25, show the structural representation of the Ziye area shape repair apparatus of reverse intensity modulated radiation therapy, this device comprises:

The Ziye area balancer 261 of reverse intensity modulated radiation therapy, carries out area equilibrium for the Ziye matrix identical to the machine jumping figure in all matrixes be partitioned in advance; Wherein, the Ziye area balancer of this reverse intensity modulated radiation therapy is the device shown in above-mentioned Figure 24.

Shape repairs module 262, in the Ziye matrix after area equilibrium, chooses each Ziye matrix successively, and carries out shape repair operation to choosing Ziye matrix;

Described shape repairs module 262, specifically for exchanging the part Ziye area choosing Ziye matrix and the 3rd Ziye matrix correspondence position, to reduce the complex-shaped degree of Ziye matrix and the 3rd Ziye matrix chosen, described 3rd Ziye matrix be not with the described Ziye matrix chosen Ziye matrix and participate in exchanging; A Ziye matrix is chosen as the 3rd Ziye matrix from described having neither part nor lot in the Ziye matrix of exchange, continue to perform described exchange and choose the step of the part Ziye area of Ziye matrix and the 3rd Ziye matrix correspondence position, until described in choose Ziye matrix and other each Ziye matrix all carries out exchange operation.

Wherein, described shape is repaired module 262 and is comprised area exchange unit, and described area exchange unit chooses the part Ziye area of Ziye matrix and the 3rd Ziye matrix correspondence position for exchanging;

Described area exchange unit comprises:

Exchange subelement line by line, for choosing Ziye matrix and choosing mutually corresponding at least one line number value in the 3rd Ziye matrix, and exchange all or part of Ziye area corresponding to every line number value of choosing line by line; Specifically for exchange line by line choose Ziye matrix and the 3rd Ziye matrix current line numerical value to exchange Ziye area corresponding to current line numerical value; Or, exchange line by line and choose first choose area and second choose area corresponding with the current line numerical value of the 3rd Ziye matrix of Ziye matrix, wherein, first chooses area chooses the equal and opposite in direction of area with second and position is identical, when described first to choose area be the blade openings area choosing Ziye matrix, described second to choose area be that the blade of the 3rd Ziye matrix is remained silent area, when described first choose area be the blade choosing Ziye matrix remain silent area time, described second chooses the blade openings area that area is the 3rd Ziye matrix.

The computing formula of above-mentioned shape complexity is:

C_{i} = Σ_{j = 2}^{R - 1} (I L_{j} + {IR}_{j}), i = 1, . . ., N;

{IL}_{j} = \{\begin{matrix} (l_{j} - l_{j - 1}) + (l_{j} - l_{j + 1}) & l_{j} > l_{j + 1}, l_{j} > l_{j - 1} \\ 0 & else \end{matrix}

{IR}_{j} = \{\begin{matrix} (r_{j - 1} - r_{j}) + (r_{j + 1} - r_{j}) & r_{j - 1} > r_{j}, r_{j + 1} > r_{j} \\ 0 & else \end{matrix}

In addition, described area exchange unit also comprises:

The Ziye area shape method for repairing and mending of the reverse intensity modulated radiation therapy that the embodiment of the present invention provides, after carrying out area equilibrium to Ziye area, also further obviates the situation of outstanding blade, thus further increases the accuracy of dose of IMRT.

As seen through the above description of the embodiments, those skilled in the art can be well understood to the mode that all or part of step in above-described embodiment method can add required general hardware platform by software and realizes.Based on such understanding, technical scheme of the present invention can embody with the form of software product the part that prior art contributes in essence in other words, this computer software product can be stored in storage medium, as ROM/RAM, magnetic disc, CD etc., comprising some instructions in order to make a computer equipment (can be personal computer, server, or the network communication equipments such as such as WMG, etc.) perform the method described in some part of each embodiment of the present invention or embodiment.

It should be noted that, in this description, each embodiment adopts the mode of going forward one by one to describe, and what each embodiment stressed is the difference with other embodiments, between each embodiment identical similar portion mutually see.For device disclosed in embodiment, because it corresponds to the method disclosed in Example, so description is fairly simple, relevant part illustrates see method part.

Also it should be noted that, in this article, the such as relational terms of first and second grades and so on is only used for an entity or operation to separate with another entity or operating space, and not necessarily requires or imply the relation that there is any this reality between these entities or operation or sequentially.And, term " comprises ", " comprising " or its any other variant are intended to contain comprising of nonexcludability, thus make to comprise the process of a series of key element, method, article or equipment and not only comprise those key elements, but also comprise other key elements clearly do not listed, or also comprise by the intrinsic key element of this process, method, article or equipment.When not more restrictions, the key element limited by statement " comprising ... ", and be not precluded within process, method, article or the equipment comprising described key element and also there is other identical element.

To the above-mentioned explanation of the disclosed embodiments, professional and technical personnel in the field are realized or uses the present invention.To be apparent for those skilled in the art to the multiple amendment of these embodiments, General Principle as defined herein can without departing from the spirit or scope of the present invention, realize in other embodiments.Therefore, the present invention can not be restricted to these embodiments shown in this article, but will meet the widest scope consistent with principle disclosed herein and features of novelty.

Claims

1. a Ziye area equalization methods for reverse intensity modulated radiation therapy, is characterized in that, comprising:

From all matrixes be partitioned in advance, extract the Ziye matrix that machine jumping figure is identical, described machine jumping figure be each divided go out the coefficient of Ziye matrix;

From the Ziye matrix of described extraction, choose one or more Ziye area be greater than the second Ziye matrix that setting the first Ziye matrix of threshold value and Ziye area be less than setting threshold value, and reduce Ziye area corresponding to described first Ziye matrix and increase Ziye area corresponding to described second Ziye matrix, be all more than or equal to setting threshold value with the Ziye area of the second Ziye matrix of the first Ziye matrix and described increase that make described reduction;

2. method according to claim 1, it is characterized in that, describedly from the Ziye matrix of described extraction, choose one or more Ziye area be greater than the second Ziye matrix that setting the first Ziye matrix of threshold value and Ziye area be less than setting threshold value, and reduce Ziye area corresponding to described first Ziye matrix and the Ziye area increasing described second Ziye matrix corresponding comprises:

Exchange the Ziye area of described first Ziye matrix and described second Ziye matrix line by line, and after a line Ziye area often exchanging described first Ziye matrix and described second Ziye matrix, retain the matrix condition before exchange according to the area variance size of exchanging between front and back the first Ziye matrix and the second Ziye matrix or retain the matrix condition after exchanging;

Judge whether the Ziye area that the second Ziye matrix after all exchanging line by line is corresponding is less than described setting threshold value, if so, then judges whether there is the Ziye matrix be not selected in the Ziye matrix of the described first kind;

3. method according to claim 2, is characterized in that,

If there is no the Ziye matrix be not selected described in and Ziye area corresponding to the second Ziye matrix is still less than described setting threshold value, then all Ziye matrixes of the first kind and the second Ziye matrix are all returned to divided rear time state.

4. method according to claim 2, it is characterized in that, the described Ziye area exchanging described first Ziye matrix and described second Ziye matrix line by line, and after a line Ziye area often exchanging described first Ziye matrix and described second Ziye matrix, retain the matrix condition before exchange according to the area variance size of exchanging between front and back the first Ziye matrix and the second Ziye matrix or retain the matrix condition after exchanging, comprising:

5. method according to claim 4, is characterized in that, after being exchanged by the current numerical value of going together of described first Ziye matrix and described second Ziye matrix, also comprises:

6. a Ziye area shape method for repairing and mending for reverse intensity modulated radiation therapy, is characterized in that, comprising:

The Ziye matrix identical to the machine jumping figure in all matrixes be partitioned in advance carries out area equilibrium, wherein, adopt the method described in any one of claim 1-5 to carry out area equilibrium to described Ziye matrix, described machine jumping figure be each divided go out the coefficient of Ziye matrix;

Described shape repair operation specifically comprises:

7. method according to claim 6, is characterized in that, the part Ziye area that Ziye matrix and the 3rd Ziye matrix correspondence position are chosen in described exchange comprises:

8. method according to claim 7, is characterized in that, the computing formula of described shape complexity is:

C_{i} = Σ_{j = 2}^{R - 1} ({IL}_{j} + {IR}_{j}), i = 1, ..., N;

{IL}_{j} = \{\begin{matrix} (l_{j} - l_{j - 1}) + (l_{j} - l_{j + 1}) & l_{j} > l_{j + 1}, l_{j} > l_{j - 1} \\ 0 & e l s e \end{matrix}

{IR}_{j} = \{\begin{matrix} (r_{j - 1} - r_{j}) + (r_{j + 1} - r_{j}) & r_{j - 1} > r_{j}, r_{j + 1} > r_{j} \\ 0 & e l s e \end{matrix}

9. method according to claim 7, is characterized in that, the described all or part of Ziye area exchanging every line number value of choosing line by line corresponding comprises:

10. method according to claim 7, is characterized in that, after described each exchange terminates, also comprises:

The Ziye area balancer of 11. 1 kinds of reverse intensity modulated radiation therapies, is characterized in that, comprising:

Matrix chooses module, for extracting the identical Ziye matrix of machine jumping figure from all matrixes be partitioned in advance, described machine jumping figure be each divided go out the coefficient of Ziye matrix;

Area balance module, the second Ziye matrix that setting the first Ziye matrix of threshold value and Ziye area be less than setting threshold value is greater than for choosing one or more Ziye area from the Ziye matrix of described extraction, and reduce Ziye area corresponding to described first Ziye matrix and increase Ziye area corresponding to described second Ziye matrix, be all more than or equal to setting threshold value with the Ziye area of the second Ziye matrix of the first Ziye matrix and described increase that make described reduction; Wherein, the Ziye area that described first Ziye matrix is corresponding is greater than Ziye area corresponding to described second Ziye matrix, the equal and opposite in direction of the Ziye area of described minimizing and the Ziye area of described increase and position is identical.

12. devices according to claim 11, is characterized in that, described area balance module comprises:

Exchange unit line by line, for exchanging the Ziye area of described first Ziye matrix and described second Ziye matrix line by line, and after a line Ziye area often exchanging described first Ziye matrix and described second Ziye matrix, retain the matrix condition before exchange according to the area variance size of exchanging between front and back the first Ziye matrix and the second Ziye matrix or retain the matrix condition after exchanging;

First judging unit, for judging whether the Ziye area that the second Ziye matrix after all exchanging line by line is corresponding is less than described setting threshold value;

13. devices according to claim 12, is characterized in that, described area balance module also comprises:

14. devices according to claim 12, is characterized in that, described exchange unit line by line, comprising:

First chooses subelement, for when described variance judgment sub-unit judges that the described second area variance obtained is less than described first area variance, if there is next numerical value of going together of described first Ziye matrix and described second Ziye matrix, then using next colleague's numerical value described as current colleague's numerical value, continue to utilize described colleague to exchange subelement and realize the function that the described current numerical value of going together by described first Ziye matrix and described second Ziye matrix carries out exchanging;

15. devices according to claim 14, is characterized in that, described exchange unit line by line also comprises:

First overlapping judgment sub-unit, for exchanging after the current numerical value of going together of described first Ziye matrix and described second Ziye matrix exchanges by subelement described colleague, in any one matrix in described first Ziye matrix and described second Ziye matrix, judge that the right leaf position of the accelerator whether left leaf position of accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping, and/or the left leaf position of the accelerator whether right leaf position of the accelerator that current line numerical value is corresponding is corresponding with lastrow numerical value or next line numerical value is overlapping;

The Ziye area shape repair apparatus of 16. 1 kinds of reverse intensity modulated radiation therapies, is characterized in that, comprising:

The Ziye area balancer of the reverse intensity modulated radiation therapy as described in any one of claim 11 to 15, carry out area equilibrium for the Ziye matrix identical to the machine jumping figure in all matrixes be partitioned in advance, described machine jumping figure be each divided go out the coefficient of Ziye matrix;

17. devices according to claim 16, is characterized in that, described shape is repaired module and comprised area exchange unit, and described area exchange unit chooses the part Ziye area of Ziye matrix and the 3rd Ziye matrix correspondence position for exchanging;

Described area exchange unit comprises:

18. devices according to claim 17, is characterized in that, the computing formula of described shape complexity is:

C_{i} = Σ_{j = 2}^{R - 1} ({IL}_{j} + {IR}_{j}), i = 1, ..., N;

{IL}_{j} = \{\begin{matrix} (l_{j} - l_{j - 1}) + (l_{j} - l_{j + 1}) & l_{j} > l_{j + 1}, l_{j} > l_{j - 1} \\ 0 & e l s e \end{matrix}

{IR}_{j} = \{\begin{matrix} (r_{j - 1} - r_{j}) + (r_{j + 1} - r_{j}) & r_{j - 1} > r_{j}, r_{j + 1} > r_{j} \\ 0 & e l s e \end{matrix}

19. devices according to claim 17, is characterized in that,

Describedly exchange subelement line by line, specifically for exchange line by line choose Ziye matrix and the 3rd Ziye matrix current line numerical value to exchange Ziye area corresponding to current line numerical value; Or, exchange line by line and choose first choose area and second choose area corresponding with the current line numerical value of the 3rd Ziye matrix of Ziye matrix, wherein, first chooses area chooses the equal and opposite in direction of area with second and position is identical, when described first to choose area be the blade openings area choosing Ziye matrix, described second to choose area be that the blade of the 3rd Ziye matrix is remained silent area, when described first choose area be the blade choosing Ziye matrix remain silent area time, described second chooses the blade openings area that area is the 3rd Ziye matrix.

20. devices according to claim 17, is characterized in that, described area exchange unit also comprises: