1.淮北師范大學(xué) 計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院，安徽 淮北 235000

2.南陽(yáng)師范學(xué)院 計(jì)算機(jī)與信息技術(shù)學(xué)院，河南 南陽(yáng) 473061

3.連云港市第二人民醫(yī)院 醫(yī)學(xué)影像科，江蘇 連云港 222005

1.淮北師范大學(xué) 計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院，安徽 淮北 235000

2.南陽(yáng)師范學(xué)院 計(jì)算機(jī)與信息技術(shù)學(xué)院，河南 南陽(yáng) 473061

3.連云港市第二人民醫(yī)院 醫(yī)學(xué)影像科，江蘇 連云港 222005

1 引言

調(diào)強(qiáng)放射治療是目前對(duì)惡性腫瘤進(jìn)行治療的有效手段之一，這種治療手段通過(guò)可調(diào)的照射角度、可調(diào)的各照射野射線強(qiáng)度，在待治療區(qū)域上形成可調(diào)的放射線劑量分布，保證腫瘤殺滅區(qū)域的劑量達(dá)到治療水平，而非腫瘤區(qū)域的放射劑量控制在發(fā)生并發(fā)癥的水平以下，危險(xiǎn)區(qū)域的放射劑量則控制在最低水平。這樣強(qiáng)度可調(diào)的三維形狀適應(yīng)性很強(qiáng)的放療手段改變了過(guò)去傳統(tǒng)放療并發(fā)癥較高的弊端，已經(jīng)在歐美等發(fā)達(dá)國(guó)家對(duì)頭頸部癌、前列腺癌、肺癌、腦癌等癌癥的治療中廣泛使用，目前國(guó)內(nèi)的一些大中城市的醫(yī)院也在推廣普及這種治療技術(shù)。

調(diào)強(qiáng)放療醫(yī)療器械的生產(chǎn)廠商主要還是國(guó)外的知名公司，如通用、瓦里安等。該類器械的一個(gè)關(guān)鍵技術(shù)就是如何從醫(yī)生設(shè)定計(jì)劃逆向求解得到放療規(guī)劃結(jié)果，即前期通過(guò)醫(yī)生對(duì)待治療區(qū)域器官、照射靶區(qū)的勾勒后，形成劑量要求和劑量約束條件，送入放療逆向求解系統(tǒng)中進(jìn)行優(yōu)化求解，最終得到硬件設(shè)備可以施行的照射方案。其中逆向優(yōu)化過(guò)程的一個(gè)關(guān)鍵環(huán)節(jié)就是得到各照射野的照射分布圖[1-3]。一個(gè)好的照射分布圖應(yīng)該具備以下兩個(gè)方面的要求。（1）待照射靶區(qū)的劑量接收值應(yīng)與醫(yī)生設(shè)定的計(jì)劃劑量值盡可能相等，且保證劑量分布的一致性。（2）待照射靶區(qū)周邊組織的劑量接收值應(yīng)越低越好，從而在靶區(qū)三維輪廓上形成劑量快速下降的分布趨勢(shì)，特別是一些危險(xiǎn)的串行組織需要嚴(yán)格限制劑量上限，而一些并行組織則可以對(duì)劑量約束適當(dāng)放松。

目前已經(jīng)有大量的文獻(xiàn)正關(guān)注于物理劑量模型的建立和相應(yīng)的求解算法[4-5]，考慮到求解算法的高效性要求，目前的逆向計(jì)劃大多轉(zhuǎn)化為線性或二次規(guī)劃模型[6-8]。在這些模型中，劑量逼近程度一般采用實(shí)際劑量與設(shè)定劑量的絕對(duì)差值或絕對(duì)差值的平方進(jìn)行描述；并廣泛采用劑量體積曲線的多點(diǎn)限制模式，即通過(guò)劑量體積圖曲線的調(diào)節(jié)指導(dǎo)放療優(yōu)化，讓正常組織滿足一定的劑量體積比要求的情況下，盡可能地讓待治療區(qū)域達(dá)到放射劑量要求。

非?？上У氖牵瑤в袆┝矿w積限制的數(shù)學(xué)規(guī)劃模型為非凸規(guī)劃問(wèn)題[9-11]，為了求解這樣一類較難的數(shù)學(xué)問(wèn)題，近幾年來(lái)國(guó)內(nèi)外科學(xué)家提出了許多求解方法，取得了一定的成效，但離獲得最優(yōu)解仍有較大距離。

Cotrutz等[12]使用基于體素權(quán)重系數(shù)的方法來(lái)幫助求解，在他們的方法中，首先對(duì)線性規(guī)劃問(wèn)題或二次規(guī)劃問(wèn)題進(jìn)行初始求解，并獲得初始劑量分布，然后按該劑量分布，根據(jù)低劑量體素更加容易抑制為更低劑量、高劑量體素更加容易給予更高劑量這樣的啟發(fā)式信息為指導(dǎo)，重新對(duì)體素的權(quán)重進(jìn)行分配，再次進(jìn)行優(yōu)化；重復(fù)上述過(guò)程直到達(dá)到優(yōu)化限定條件。這樣的啟發(fā)式信息有一定的直觀依據(jù)，但是形而上的想法缺乏科學(xué)依據(jù)。本文提出一種新的更加科學(xué)的啟發(fā)式方法對(duì)上述非凸規(guī)劃問(wèn)題進(jìn)行求解，并在四個(gè)病例上進(jìn)行了相關(guān)實(shí)驗(yàn)。實(shí)驗(yàn)結(jié)果證實(shí)，新的啟發(fā)式方法能夠獲得更優(yōu)的逆向計(jì)劃。

2 方法

2.1 強(qiáng)度照射分布圖優(yōu)化問(wèn)題

整個(gè)待治療區(qū)域的人體組織被分為兩個(gè)部分，且用體素進(jìn)行離散化，一部分為待照射靶區(qū)T={T1，T2，…，Tu}，相應(yīng)的體素?cái)?shù)目為…，STu，另一部分為正常組織和危險(xiǎn)器官組織N={N1，N2，…，Nv}，相應(yīng)的體素?cái)?shù)目為…，。照射分布圖用像素進(jìn)行離散化，按照照射野的順序以及分布圖逐行的順序，用列向量X進(jìn)行描述，F(xiàn)T=為靶區(qū)劑量作用矩陣，為正常組織和危險(xiǎn)器官組織的劑量作用矩陣。任意一個(gè)元素 fi，j表示第j個(gè)像素采用單位強(qiáng)度照射時(shí)對(duì)第i個(gè)體素所產(chǎn)生的劑量沉積值。 DT=FT·X和 DN=FN·X代表整個(gè)待照射區(qū)域在強(qiáng)度分布圖X下的三維劑量分布。設(shè)定醫(yī)生指定的劑量為，…，，正常組織以及危險(xiǎn)區(qū)域組織的劑量上限為，…，

給出物理劑量模型所廣泛采用的線性約束二次規(guī)劃模型如下：

2.2 劑量體積曲線限制條件

放射科臨床醫(yī)生常用的一個(gè)概念為劑量體積曲線，這個(gè)曲線能夠從宏觀上描述某個(gè)器官的劑量分布生物學(xué)特征。劑量體積曲線的橫軸為劑量值，縱軸為體積百分比，曲線上的每一個(gè)點(diǎn)表示某個(gè)器官中低于某個(gè)劑量的組織體積占整個(gè)器官組織體積的百分比。因此該曲線為一條遞減曲線。

而相應(yīng)的劑量體積曲線限制即要求曲線在某個(gè)點(diǎn)的上方或某個(gè)點(diǎn)的下方。以正常組織為例，劑量體積曲線限制條件一般限定曲線在一個(gè)點(diǎn)或幾個(gè)點(diǎn)的下方。一般描述為某個(gè)器官組織接受? Gy劑量或更高劑量的體素?cái)?shù)目不超過(guò)該器官組織總體素?cái)?shù)目的η%，其中“Gy”為放射劑量單位。直觀上，即要求曲線在點(diǎn)( ) ?，η的下方。

圖1 劑量體積曲線及約束點(diǎn)

設(shè)正常組織和危險(xiǎn)器官組織 N1，N2，…，Nv有劑量體積限定條件 DVC1，DVC2，…，DVCv，其中為劑量體積限定條件 DVCi中的約束點(diǎn)個(gè)數(shù)，為相應(yīng)的約束點(diǎn)。由此形成以下的帶有劑量體積約束的數(shù)學(xué)規(guī)劃問(wèn)題：

對(duì)于臨床醫(yī)生常用的一個(gè)概念被引入到逆向計(jì)劃中卻帶來(lái)了求解過(guò)程的極大困難，因?yàn)樯鲜瞿Ｐ蜑榻M合優(yōu)化模型。

2.3 正規(guī)化空間內(nèi)的空間距離啟發(fā)信息

上述組合優(yōu)化模型的求解為非凸規(guī)劃問(wèn)題的求解，獲得該問(wèn)題的最優(yōu)解是NP難問(wèn)題。目前國(guó)外一些參考文獻(xiàn)的求解方法是采用基于體素權(quán)重系數(shù)的問(wèn)題轉(zhuǎn)換求解思路，即通過(guò)對(duì)各體素設(shè)定一定的權(quán)重，從而將上述問(wèn)題轉(zhuǎn)換為凸規(guī)劃問(wèn)題進(jìn)行求解，當(dāng)然，還有一些思路則回避劑量體積曲線約束，考慮采用其他約束取代這類非凸約束，從而轉(zhuǎn)換為凸規(guī)劃問(wèn)題。對(duì)于后者，認(rèn)為約束條件的轉(zhuǎn)換已經(jīng)基本違背了臨床醫(yī)生的劑量設(shè)定初衷，因此本文不予考慮。對(duì)于問(wèn)題轉(zhuǎn)換思路，本文予以采納，即將上述問(wèn)題轉(zhuǎn)換為凸規(guī)劃從而獲得轉(zhuǎn)換問(wèn)題的最優(yōu)解，但是對(duì)于體素排序的思路，本文提出新的啟發(fā)式想法，即從原來(lái)的劑量高低排序方式轉(zhuǎn)換為二次規(guī)劃正規(guī)化空間內(nèi)的空間距離遠(yuǎn)近排序方式。具體算法流程如圖2所示。

圖2 優(yōu)化算法流程圖

根據(jù)Cholesky分解，Q=RTR，帶入 x=R-1y，上述模型轉(zhuǎn)換為：

對(duì)于線性約束二次規(guī)劃問(wèn)題：

設(shè)(R-1)Tc=e，BR-1=G，上述線性約束二次規(guī)劃問(wèn)題有以下的正規(guī)化形式：

下面通過(guò)一個(gè)簡(jiǎn)單的示例來(lái)對(duì)劑量高低排序方式以及二次規(guī)劃正規(guī)化空間內(nèi)的空間距離遠(yuǎn)近排序方式進(jìn)行比較。

該示例只有一個(gè)規(guī)劃靶區(qū)，并為一個(gè)凹的形狀，另一個(gè)危險(xiǎn)組織區(qū)域?yàn)閳A形，并一定程度地被規(guī)劃靶區(qū)包圍。照射方向?yàn)槲鍌€(gè)均勻分布的方向，靶區(qū)設(shè)定的劑量為單位劑量，危險(xiǎn)組織區(qū)域的最大劑量不能超過(guò)0.5，正常組織劑量不超過(guò)0.8，危險(xiǎn)組織區(qū)域的劑量體積曲線要求接受高于0.3劑量的體積大小不超過(guò)50%（如圖3所示）。

圖3 示例初始劑量分布示意圖

求解了不帶有劑量體積曲線限制條件的規(guī)劃問(wèn)題，并按照劑量高低排序方式以及二次規(guī)劃正規(guī)化空間內(nèi)的空間距離遠(yuǎn)近排序方式對(duì)危險(xiǎn)組織區(qū)域進(jìn)行了排序，排序結(jié)果如圖4所示，顯然，基于二次規(guī)劃正規(guī)化空間內(nèi)的空間距離遠(yuǎn)近排序方式的排序結(jié)果更加自然，遠(yuǎn)離靶區(qū)的體素具有更高的優(yōu)先級(jí)被納入線性約束中去，從而達(dá)到所要求的劑量體積曲線約束。

圖4 危險(xiǎn)組織區(qū)域的正規(guī)化空間內(nèi)空間距離和劑量排序比較示意圖

2.4 實(shí)驗(yàn)設(shè)定

為了驗(yàn)證所提方法的有效性，使用四個(gè)臨床實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了測(cè)試，劑量設(shè)定值見(jiàn)表1～4。本文中實(shí)驗(yàn)強(qiáng)度矩陣的計(jì)算基于QIB劑量計(jì)算引擎。所有實(shí)驗(yàn)均在Matlab計(jì)算平臺(tái)下進(jìn)行。

表1 頭頸部測(cè)試?yán)膭┝吭O(shè)定值

表2 前列腺測(cè)試?yán)膭┝吭O(shè)定值

表3 肺部測(cè)試?yán)膭┝吭O(shè)定值

3 實(shí)驗(yàn)結(jié)果和分析

表4 口咽部測(cè)試?yán)膭┝吭O(shè)定值

表5 優(yōu)化結(jié)果比較

圖5 頭頸部測(cè)試?yán)膭┝矿w積曲線

在四個(gè)測(cè)試病例上，對(duì)基于劑量高低排序方式的啟發(fā)信息以及基于正規(guī)化空間內(nèi)的空間距離遠(yuǎn)近排序方式的啟發(fā)信息進(jìn)行了實(shí)驗(yàn)對(duì)比，表5中給出了兩種方法所得到的目標(biāo)函數(shù)值。圖5～8中給出了兩種方法得到的靶區(qū)各照射目標(biāo)以及非靶區(qū)各器官的劑量體積曲線比較結(jié)果（標(biāo)有實(shí)心圖例的實(shí)線是本文方法的結(jié)果，標(biāo)有對(duì)應(yīng)的虛心圖例的虛線是基于劑量高低排序方式方法的結(jié)果）。圖9～12中還給出了兩種方法得到的斷層劑量分布圖的比較結(jié)果。實(shí)驗(yàn)結(jié)果顯示，在達(dá)到同樣的劑量體積曲線限制條件的約束情況下，采用提出的啟發(fā)信息，病例一的目標(biāo)函數(shù)值降低了（（26.9－44.9）/44.9=－40%，病例二降低了（54.5－ 84.0）/84.0=－35%，病例四降低了（17.2－23.5）/23.5=－27%），除病例三以外降低幅度均較大。通過(guò)比較靶區(qū)各目標(biāo)的劑量體積曲線，可以看到實(shí)線比虛線下降得更為陡峭，也就是說(shuō)采用本文提出的啟發(fā)信息，靶區(qū)的劑量分布更加均勻，且更接近于所設(shè)定的處方劑量值，通過(guò)斷層劑量分布圖可以明顯地得到這個(gè)結(jié)論。

圖6 前列腺測(cè)試?yán)膭┝矿w積曲線

圖7 肺部測(cè)試?yán)膭┝矿w積曲線

圖8 口咽部測(cè)試?yán)膭┝矿w積曲線

圖9 頭頸部測(cè)試?yán)哪硵鄬觿┝糠植紙D

圖10 前列腺測(cè)試?yán)哪硵鄬觿┝糠植紙D

圖11 肺部測(cè)試?yán)哪硵鄬觿┝糠植紙D

圖12 口咽部測(cè)試?yán)哪硵鄬觿┝糠植紙D

4 結(jié)束語(yǔ)

本文對(duì)調(diào)強(qiáng)放療數(shù)學(xué)規(guī)劃問(wèn)題中的非凸規(guī)劃提出一種新的啟發(fā)式思想，利用這種思想得到的非凸規(guī)劃問(wèn)題求解方法能夠獲得更優(yōu)的優(yōu)化結(jié)果。在四個(gè)病例上的實(shí)驗(yàn)結(jié)果證實(shí)了所提新方法的有效性。此外，這種啟發(fā)思想能夠在其他非凸數(shù)學(xué)規(guī)劃問(wèn)題的求解中同樣得到廣泛應(yīng)用。

[1]Bortfeld T.IMRT：a review and preview[J].Physics in Medicine and Biology，2006，51：363-379.

[2]Hatano K，Araki H，Sakai M，et al.Current status of Intensity-ModulatedRadiation Therapy（IMRT）[J].International Journal of Clinical Oncology，2007，12：408-415.

[3]Yu C X，Amies C J，Svatos M.Planning and delivery of intensity-modulated radiation therapy[J].Medical Physics，2008，25：5233-5241.

[4]Shepard D M，F(xiàn)erris M C，Olivera G H，et al.Optimizing the delivery of radiation therapy to cancer patients[J].Siam Review，1999，41：721-744.

[5]Lim J，F(xiàn)erris M C，Wright S J，et al.An optimization framework for conformal radiation treatment planning[J].INFORMS Journal on Computing，2007，19：366-380.

[6]Morrill S M，Lane R G，Wong J A，et al.Dose-volume considerations with linear programming optimization[J].Medical Physics，1991，18：1201-1210.

[7]Breedveld S，Storchi P R M，Keijzer M，et al.Fast，multiple optimizations of quadratic dose objective functions in IMRT[J]. Physics in Medicine and Biology，2006，51：3569-3579.

[8]Censor Y，Ben-Israel A，Xiao Y，et al.On linear infeasibility arising in intensity-modulated radiation therapy inverse planning[J].LinearAlgebra and ItsApplications，2008，428：1406-1420.

[9]Deasy J O.Multiple local minima in radiotherapy optimization problems with dose-volume constraints[J].Medical Physics，1997，24：1157-1161.

[10]Wu Q，Mohan R.Multiple local minima in IMRT optimization based on dose-volume criteria[J].MedicalPhysics，2002，29：1514-1527.

[11]Wu C，Jeraj R，Mackie T R.The method of intercepts in parameter space for the analysis of local minima caused by dose-volume constraints[J].Physics in Medicine and Biology，2003，48：149-157.

[12]Cotrutz C，Xing L.Using voxel-dependent importance factors for interactive DVH-based dose optimization[J].Physics in Medicine and Biology，2002，47：1659-1669.

基于啟發(fā)式信息的非凸放療規(guī)劃模型的求解方法

張棟冰1，蘭義華2，萬(wàn)金鑫3

ZHANG Dongbing1,LAN Yihua2,WAN Jinxin3

1.School of Computer Science and Technology,Huaibei Normal University,Huaibei,Anhui 235000,China
2.School of Computer and Information Technology,Nanyang Normal University,Nanyang,Henan 473061,China
3.Medical Imaging Department,Lianyungang Second People’s Hospital,Lianyungang,Jiangsu 222005,China

For the inverse planning process of intensity-modulated radiotherapy—the non-convex mathematical programming for the fluence map optimization with dose volume constraints,this paper presents a novel and scientific heuristic information—distance priority value in normalized space.Compared to the traditional dose ordering heuristic information,this new method leads to better solutions.A phantom example and four test cases show the effectiveness of the proposed method.

non-convex fractional programming;heuristic algorithm;intensity-modulated radiation therapy technology;linear constrained quadratic programming

針對(duì)調(diào)強(qiáng)放療逆向優(yōu)化過(guò)程中的關(guān)鍵環(huán)節(jié)——各照射野的強(qiáng)度照射分布圖在帶有劑量體積曲線限制條件下的非凸數(shù)學(xué)規(guī)劃問(wèn)題，提出了一種新穎的更加科學(xué)的啟發(fā)式信息——正規(guī)化空間內(nèi)的空間距離排序值。與傳統(tǒng)的劑量排序啟發(fā)式信息相比較，新方法可以得到更好的解。一個(gè)簡(jiǎn)單示例和四個(gè)測(cè)試病例表明了該方法的有效性。

非凸數(shù)學(xué)規(guī)劃；啟發(fā)式求解；調(diào)強(qiáng)放療；線性約束二次規(guī)劃

A

TP391

10.3778/j.issn.1002-8331.1302-0021

ZHANG Dongbing,LAN Yihua,WAN Jinxin.Method based on heuristic information for solving non-convex model of radiotherapy planning.Computer Engineering and Applications,2013,49（11）：265-270.

國(guó)家自然科學(xué)基金（No.61102117）。

張棟冰（1974—），男，講師，研究領(lǐng)域?yàn)閳D像處理、算法優(yōu)化；蘭義華，男，博士，講師，研究領(lǐng)域：圖像處理、算法優(yōu)化；萬(wàn)金鑫，男，研究員級(jí)高級(jí)工程師，研究方向：醫(yī)學(xué)圖像診斷方法。E-mail：291508366@qq.com

2013-02-04

2013-04-15

1002-8331（2013）11-0265-06