蔣先剛，丘赟立，范德營

(華東交通大學(xué)基礎(chǔ)科學(xué)學(xué)院，江西南昌330013)

0 引言

基于真彩圖像的二維和三維圖像更能顯示出真實(shí)的器官分布細(xì)節(jié)和層次，而目前的計算機(jī)仿真醫(yī)學(xué)模型基本上是灰色基調(diào)的。由于醫(yī)學(xué)上常用的人體成像技術(shù)如MRI的成像結(jié)果是灰度圖像，由這樣的圖像構(gòu)造出來的三維模型無法清晰有效地反映人體的實(shí)際情況，且簡單的灰度圖像彩色化方法的效果不好而復(fù)雜的彩色化方法運(yùn)算效率不高而無法快速獲得能夠有效地反映人體實(shí)際的彩色切片。2001年Reinhard首次在《Color Transfer Between Images》提出Lαβ顏色空間，次年 Welsh在《Transferring Color to Grayscale Images》提出了基于Lαβ顏色空間的顏色遷移算法，2009年由陳倩在《一種改進(jìn)的Welsh灰度圖像彩色化算法》中對Welsh顏色遷移算法進(jìn)化了改進(jìn)優(yōu)化。本文在他們所作的研究的基礎(chǔ)上結(jié)合多種搜索尋優(yōu)算法并作適當(dāng)?shù)母倪M(jìn)以提高遷移顏色算法的效率和質(zhì)量。本文討論多種灰色圖像彩色化技術(shù)，在顏色遷移的匹配點(diǎn)搜索過程中，著重研究和對比多種諸如窮舉法、隨機(jī)法、GA和PSO算法，并重點(diǎn)分析和改進(jìn)PSO算法。最后對彩色化的圖片進(jìn)行三維重構(gòu)，得到的三維重構(gòu)模型還可以反過來對彩色化算法進(jìn)行檢驗。

1 灰色圖像的顏色遷移

灰色圖片的彩色化技術(shù)主要包括偽彩色和顏色遷移方法。偽彩色圖像中的顏色并非反映器官切片組織的分布顏色。顏色遷移是首先根據(jù)一幅灰色圖像和一幅彩色圖像的特征建立映射關(guān)系，然后根據(jù)建立的映射將彩色圖像相應(yīng)的彩色值遷移到灰色圖像中，最后獲得相應(yīng)的彩色化圖像。經(jīng)過合理遷移而得到的彩色化圖像能反映組織構(gòu)造的真實(shí)的顏色效果。本文使用的是Welsh顏色遷移算法，是以灰色MRI切片和真彩色切片的亮度和紋理特征進(jìn)行配準(zhǔn)而實(shí)現(xiàn)顏色遷移技術(shù)。該算法的處理步驟主要包括對灰色圖片進(jìn)行逐個像素塊掃描，而且對彩色圖片也進(jìn)行逐個像素塊掃描，并在彩色圖像中選擇與灰色圖像相應(yīng)像素塊最相似的像素塊的顏色，然后將選擇的顏色遷移至灰色圖片中。圖1是人體腿部灰色圖像彩色化的過程圖。圖1(a)是作為彩色源圖像的人體腿部的真彩色切片;圖1(b)是人體腿部的MRI切片;圖1(c)是圖1(b)經(jīng)過灰度變換拉伸后的圖片，該圖片的亮度和紋理特征更符合人體器官的實(shí)際情況;圖1(d)是利用顏色遷移技術(shù)將圖1(a)的顏色遷移至圖1(c)的彩色化效果圖。

圖1 人體腿部灰度圖片彩色遷移過程

常用的顏色空間中各分量的相關(guān)性較大，在進(jìn)行顏色遷移前，應(yīng)先將圖像從RGB空間轉(zhuǎn)換到一個各分量相對獨(dú)立的空間中。Lαβ顏色空間是Reinhard等人在1998年提出的，其中L是亮度通道，α是黃綠通道;β是黃藍(lán)通道。該空間的三個通道相互垂直，即亮度與顏色分量是相對獨(dú)立的。技術(shù)實(shí)現(xiàn)中首先將兩幅圖像從RGB顏色空間轉(zhuǎn)換到Lαβ顏色空間，然后建立映射關(guān)系并對其執(zhí)行顏色遷移的相關(guān)操作，最后再將遷移所得的圖像從Lαβ顏色空間轉(zhuǎn)換回到RGB顏色空間。RGB顏色空間轉(zhuǎn)換到Lαβ顏色空間的具體步驟如 (1)至 (3)所示，而Lαβ顏色空間轉(zhuǎn)換到RGB顏色空間可以由該過程的逆過程得到。

(1)將RGB顏色空間轉(zhuǎn)換到LMS顏色空間

(2)將LMS顏色空間轉(zhuǎn)換到對數(shù)空間以消除空間的Skew不對稱性

(3)將處理過的LMS分量轉(zhuǎn)換到Lαβ顏色空間

2 基于改進(jìn)PSO的顏色遷移

粒子群優(yōu)化算法 (PSO)是由Kennedy和Eberhart等在1995年提出的一種進(jìn)化計算算法，PSO算法起初是為了模擬鳥類群體進(jìn)行覓食的優(yōu)美且無法預(yù)測的運(yùn)動。后來的進(jìn)一步研究證明群體中的信息共享對群體中的個體的進(jìn)化發(fā)揮著積極作用。最后對其引入多維度的搜索功能和與距離相關(guān)的加速度權(quán)重，以及慣性權(quán)重而形成了PSO的標(biāo)準(zhǔn)版本。PSO同GA類似，都是基于迭代的進(jìn)化尋優(yōu)算法。它比遺傳算法規(guī)則更為簡單，沒有遺傳算法的交叉和變異操作，但它保存每個粒子在進(jìn)化過程中經(jīng)過的最優(yōu)解作為局部最優(yōu)解pBest和所有粒子在進(jìn)化過程中經(jīng)過的最優(yōu)解作為全局最優(yōu)解gBest，并且通過讓每個粒子趨于pBest和粒子群總體趨于gBest而獲得最終的最優(yōu)解。PSO首先產(chǎn)生初始種群，然后記錄pBest和gBest并按照式 (4)來修正每個粒子的速度和式 (5)來更新每個粒子的位置

式中:t——迭代時間參數(shù)，n——粒子的維數(shù)，j——粒子序號，xjk(t)和xjk(t+1)——粒子j在t和t+1時刻的位置，vjk(t)和vjk(t+1)——粒子j在t和t+1時刻的速度，pjk——粒子j的局部最佳位置，pgk——粒子群的全局最佳位置，w——慣性權(quán)重，c1——認(rèn)識系數(shù)，c2——社會系數(shù)，r1，r2∈ (0，1］——隨機(jī)數(shù)

式中:t——迭代時間參數(shù)，n——粒子的維數(shù)，j——粒子序號，xjk(t)和xjk(t+1)——粒子j在t和t+1時刻的位置，vjk(t+1)——粒子j在t+1時刻的速度。

該算法通過反復(fù)調(diào)整粒子的速度和位置而實(shí)現(xiàn)群體的進(jìn)化，當(dāng)進(jìn)化代數(shù)達(dá)到設(shè)定值或適應(yīng)度值達(dá)到一定的允限時，算法終止并將適應(yīng)度最高的粒子輸出，然后根據(jù)它的坐標(biāo)X，Y獲取源圖相應(yīng)的彩色值并將該彩色值遷移至灰色圖像中。

2.1 PSO在顏色遷移中的參數(shù)設(shè)置

為了將PSO應(yīng)用到Welsh顏色遷移算法中，本文設(shè)置PSO中的粒子的位置就是彩色圖像中相應(yīng)像素的坐標(biāo)，粒子的最大速度取為彩色圖像的寬和高的最小值的5%。為了使粒子具有多樣性且能從更大范圍搜索最優(yōu)解，將彩色像素劃分為m×n個子塊并在每個子塊中至少選擇一個粒子。在本實(shí)驗中設(shè)置代數(shù)計數(shù)器初值t=1，總代數(shù)T=10，粒子個數(shù)N=200。

慣性權(quán)重w是使粒子保持運(yùn)動的慣性，而認(rèn)知系數(shù)c1和社會系數(shù)c2是將每個粒子推向局部最優(yōu)解pBest和全局最優(yōu)解gBest的統(tǒng)計加速項的權(quán)重，c1和c2是加速常數(shù)。較小的加速常數(shù)使得粒子在到達(dá)全局最優(yōu)的位置之前搜索更多的位置但需要更多的時間，而較大的加速常數(shù)能使粒子更快地到達(dá)全局最優(yōu)的位置但會導(dǎo)致粒子突然地沖向或者越過全局最優(yōu)的位置。如果沒有加速常數(shù)，即c1=c2=0，這時粒子將保持直線運(yùn)動。在這種情況下，每個粒子能搜索的范圍將大大減少，不利于獲得全局最優(yōu)解。如果沒有慣性權(quán)重，即w=0，粒子的速度沒有記憶性而只取決于粒子當(dāng)前的位置和pBest及gBest的相對位置。如果一個粒子已經(jīng)處于全局最優(yōu)的位置上，它將保持靜止。而其它不在全局最優(yōu)的位置上的粒子則飛向pBest和gBest的加權(quán)中心。而認(rèn)知系數(shù)c1和社會系數(shù)c2分別決定粒子趨向于pBest和gBest的能力或程度。基于全面兼顧，實(shí)驗中設(shè)置慣性權(quán)重w=0.9，認(rèn)知系數(shù)c1=2，社會系數(shù)c2=2。

2.2 PSO在顏色遷移中的適應(yīng)度函數(shù)選擇

適應(yīng)度函數(shù)是用來反映粒子在群體中對目標(biāo)的適應(yīng)程度，即灰色圖像和彩色圖像中相應(yīng)像素塊的相似程度。本文使用圖像的亮度和紋理特征來反映像素塊的相似程度。適應(yīng)度函數(shù)中，亮度參數(shù)的具體定義如下

式中:ls和σs——彩色圖片的像素塊的亮度均值和標(biāo)準(zhǔn)差，lt和σt——灰色圖片的像素塊的亮度均值和標(biāo)準(zhǔn)差，λ1和λ2——彩色圖片和灰色圖片的亮度均值偏差和標(biāo)準(zhǔn)差偏差的權(quán)重。而紋理特征采用差分矩特征u)2P(i，j)，P是像素塊的灰度共生矩陣，u是P(i，j)均值，則相應(yīng)的紋理參數(shù)的定義如下

式中:Ws和Wt——彩色圖片和灰色圖片的像素塊的差分矩特征，則相應(yīng)的綜合適應(yīng)度函數(shù)為

式中:α1和α2——亮度參數(shù)和紋理參數(shù)的權(quán)重，C=100是常量，這樣F的取值范圍規(guī)約為［0，100］。

由于圖像中相鄰的像素的顏色具有一定的相似性，為了實(shí)現(xiàn)快速顏色遷移需在進(jìn)行PSO搜索前先利用點(diǎn)的彩色局部分布特征的方法搜索匹配點(diǎn)。該方法根據(jù)彩色的局部分布特征而首先搜索待彩色化的像素的左邊像素對應(yīng)于彩色圖像中的對應(yīng)像素的右邊的像素，利用式 (6)定義的適應(yīng)度函數(shù)計算該像素的適應(yīng)度值，如果該值小于一個指定的閾值，則將該像素規(guī)定為待彩色化的像素的最優(yōu)匹配結(jié)果。如果該像素不滿足條件，則依次搜索待彩色化的像素的左上角、上面和右上角的像素對應(yīng)于彩色圖像中的像素的右下角、下面和左下角的像素。如果某一像素滿足搜索條件則輸出其為最優(yōu)結(jié)果。否則采用PSO進(jìn)行搜索。另外，如果待彩色化的像素的左邊、左上角、上面或右上角的像素是無效的，即認(rèn)為其是不在圖像有效范圍內(nèi)的點(diǎn)，則直接跳過而不必考慮。

2.3 改進(jìn)PSO在顏色遷移中的應(yīng)用

由于標(biāo)準(zhǔn)的PSO的遷移效率與隨機(jī)法和GA相比并沒有得到明顯的提高，所以本文嘗試在下面4點(diǎn)對PSO進(jìn)行了改進(jìn)。

(1)設(shè)置粒子的取值范圍。在標(biāo)準(zhǔn)的PSO中，新產(chǎn)生的粒子有可能在圖片的有效點(diǎn)之外，則其適應(yīng)度值為0。該粒子沒有對種群的進(jìn)化起到積極的作用，應(yīng)將這類無效的粒子移至粒子有效范圍的邊界以提高算法的收斂速度。

(2)灰色圖像中每一新的下一個待使用PSO進(jìn)行搜索匹配的點(diǎn)的初始群體使用上一個點(diǎn)的最優(yōu)群體。一般情況下，彩色圖像中相鄰的兩個點(diǎn)的顏色是密切相關(guān)的，所以上一個點(diǎn)的最優(yōu)群體有較大的可能性會適應(yīng)下一個點(diǎn)。

(3)在進(jìn)化的后期對粒子進(jìn)行減速處理，以免其跨過最優(yōu)位置而只能獲得次優(yōu)的結(jié)果，對粒子的最大速度的調(diào)整采用式 (9)

(4)為了防止pBest和gBest對粒子的速度方向的影響產(chǎn)生滯后性，隨著迭代次數(shù)的增加而按式 (10)線性減少慣性權(quán)重w

式中:t和T——當(dāng)前代數(shù)和最大進(jìn)化代數(shù)，w(t)——群體進(jìn)化到t代時的慣性權(quán)重，取wmax=0.9和wmin=0.4——設(shè)定的最大值和最小值。

3 實(shí)驗結(jié)果及分析

本文所做實(shí)驗的開發(fā)平臺為Delphi 2010，相應(yīng)的硬件環(huán)境為Intel酷睿i5 480M，2.66GHz的CPU和6GB的內(nèi)存。實(shí)驗所需的彩色源圖像采用中國可視人真彩切片。這些切片真實(shí)客觀地反映了人體各器官組織在各個切面上的分布和顏色特征。實(shí)驗圖片的分辨率均為256×256。窮舉法是遍歷彩色源圖像的每一個像素來搜索灰色圖像中待彩色化的像素的最佳匹配點(diǎn)，而隨機(jī)法則是取彩色源圖像的總像素量的15%的像素進(jìn)行搜索。為了更好地說明改進(jìn)PSO算法的平均遷移效率，本文在依據(jù)點(diǎn)鄰域分布和彩色局部分布特征并采用多種搜索算法如窮舉法、隨機(jī)法、GA、PSO、改進(jìn)GA和改進(jìn)PSO進(jìn)行了對比分析實(shí)驗，而且每組實(shí)驗進(jìn)行20次，實(shí)驗對比結(jié)果如表1所示。

表1 基于不同特征和搜索方法的顏色遷移的效率比較

根據(jù)表1可以得出如下結(jié)論:窮舉法雖然能得到理想最優(yōu)的結(jié)果，但是耗時最長，不宜應(yīng)用于實(shí)際的圖片顏色遷移;跟窮舉法相比，隨機(jī)法的效率有了明顯的提高;與隨機(jī)法相比，GA和PSO也有了進(jìn)一步的提高，但是提高的程度并不明顯;改進(jìn)GA和改進(jìn)PSO的效率是最高的，而改進(jìn)PSO的遷移效率甚至高達(dá)窮舉法的20多倍。圖2是基于GA、PSO及其改進(jìn)算法的進(jìn)化代數(shù)與如式 (8)所示適應(yīng)度的逼近曲線。圖2中表明改進(jìn)PSO的群體的最優(yōu)適應(yīng)度值較快速達(dá)到90而趨于適應(yīng)度值為100的全局最優(yōu)解。

圖2 基于不同搜索算法的進(jìn)化效率

為了突顯基于灰度切片和彩色化切片的三維重構(gòu)圖像的效果的優(yōu)劣，在灰度切片實(shí)現(xiàn)三維重構(gòu)的直接體繪制的過程中使用亮度傳遞函數(shù)，而在彩色化切片實(shí)現(xiàn)三維重構(gòu)的過程中使用與顏色向量和顏色梯度相關(guān)的傳遞函數(shù)如式 (11)

式中:T∈［0，1］——閾值參數(shù)，ω ∈［0，1］——權(quán)衡顏色向量和顏色梯度的參數(shù)，γ∈(1，∞)——指數(shù)權(quán)重參數(shù)，c和g——像素的顏色向量和顏色梯度，而傳遞函數(shù)o反映的是數(shù)據(jù)場中像素點(diǎn)的顏色向量和顏色梯度等因素對三維重構(gòu)模型的影響程度。

圖3(a)是基于灰度切片的三維重構(gòu)模型，圖3(b)是基于強(qiáng)調(diào)白黃紅的彩色化切片的三維彩色仿真模型，圖3(c)是基于顏色遷移后的強(qiáng)調(diào)肌肉和血管顏色的三維彩色仿真模型。通過對圖3的對比分析可知，基于彩色化切片的三維重構(gòu)模型的效果明顯比基于灰度切片的三維重構(gòu)模型的效果更好，而且基于彩色化切片的三維重構(gòu)模型在直接體繪制的過程中選取不同的顏色向量和顏色梯度的權(quán)重能多層次地反映出人體器官組織的三維細(xì)節(jié)和客觀構(gòu)造。

圖3 不同顏色權(quán)重的器官重構(gòu)效果

4 結(jié)束語

本文在傳統(tǒng)的圖像彩色化算法中引入仿生智能算法，并對粒子群優(yōu)化PSO進(jìn)行了改進(jìn)性實(shí)驗，其中通過限定粒子的取值范圍、保留最終進(jìn)化群體作為下一次進(jìn)行PSO搜索的初始群體、隨著進(jìn)化代數(shù)的增加而降低最大速度和慣性權(quán)重這些措施和改進(jìn)對PSO的效率的提高有顯著促進(jìn)。所得的基于改進(jìn)的PSO的Welsh彩色化算法的效率非常高，比一般的窮舉法、隨機(jī)法、GA和PSO的效率要高出很多倍。今后的進(jìn)一步研究將會放在GA和PSO的組合求解尋優(yōu)上。在這些彩色仿真切片的基礎(chǔ)上，通過基于彩色向量和梯度的傳遞函數(shù)可以得到三維仿真彩色模型，而這樣的三維仿真彩色模型的確能有效地反映人體器官組織真正的顏色和多層次空間分布。

［1］LIU Yong，SHAO Changyong.Improvement of Welsh's grey image colorization arithmetics［J］.Modern Electronics Technique，2009(24)141-143(in Chinese).［劉勇，邵常勇.對Welsh灰度圖像彩色化算法的改進(jìn) ［J］.現(xiàn)代電子技術(shù)，2009(24)141-143.］

［2］CHEN Qing.Improvement of Welsh's grey image colorization arithmetic［J］.Journal of Wuhan University of Technology，2009，31(22)151-154(in Chinese).［陳倩.一種改進(jìn)的Welsh灰度圖像彩色化算法 ［J］.武漢理工大學(xué)學(xué)報，2009，31(22)151-154.］

［3］LIU Xiaopei，WAN Liang，LIN S，et al.Intrinsic colorization ［J］.ACM Trans on Graphics，2008，27(5):1521-1529.

［4］XU Zhen，SUN Shaoyuan，DAI Zhonghua，et al.Texture library based color＆ texture transferring for LLL images［J］.Infrared Technology，2011，33(1):49-55(in Chinese).［許真，孫韶媛，代中華，等.基于紋理特征庫的微光圖像色彩紋理傳遞［J］.紅外技術(shù)，2011，33(1):49-55.］

［5］ZHENG Zhi，HUANGJian.Analysis of the application mode of color transfer algorithm in images［J］.Journal of Fujian Agriculture and Forestry University，2010，39(2):206-210(in Chinese).［鄭之，黃鍵.圖像色彩傳遞算法的應(yīng)用模式分析 ［J］.福建農(nóng)林大學(xué)學(xué)報，2010，39(2):206-210.］

［6］FENG Jun，XUE Yuncan，JIANG Jinlong.Study on a new modified particle swarm optimization algorithm［J］.Journal of Hohai University Changzhou，2006，20(1):10-13(in Chinese).［馮駿，薛云燦，江金龍.一種新的改進(jìn)粒子群算法研究［J］.河海大學(xué)常州分校學(xué)報，2006，20(1):10-13.］

［7］LIN Jinping，YU Jinxiang.Parameter estimation of multiple linear regression models based on t he improved particle swarm optimization algorithm［J］.Computer Engineering＆ Science，2010，32(4):101-105(in Chinese).［劉錦萍，郁金祥.基于改進(jìn)的粒子群算法的多元線性回歸模型參數(shù)估計 ［J］.計算機(jī)工程與科學(xué)，2010，32(4):101-105.］

［8］GUO Huiyong，WANG Lei，LI Zhengliang.Two-stage damage detection method based on improved particle swarm optimization algorithm ［J］.Journal of Southwest Jiaotong University，2011，46(6):926-932(in Chinese).［郭惠勇，王磊，李正良.基于改進(jìn)PSO算法的兩階段損傷識別方法［J］.西南交通大學(xué)學(xué)報，2011，46(6):926-932.］

［9］LIU Shicheng，ZHANG Jianhua，LIU Zongqi.Application of parallel adaptive particle swarm optimization algorithm in reactive power optimization of power system［J］.Power System Technology，2012，36(1):108-112(in Chinese).［劉世成，張建華，劉宗岐.并行自適應(yīng)粒子群算法在電力系統(tǒng)無功系統(tǒng)［J］.電網(wǎng)技術(shù)，2012，36(1):108-112.］

［10］Suresh Chandra Satapathy，Gunanidhi Pradhan，Sabyasachi Pattnaik，et al.Performance comparisons of PSO based clustering［J］.InterJRI Computer Science and Networking，2009，1(1):18-23.