基于Retinex理論的彩色眼底圖像增強(qiáng)方法研究 通信工程專業(yè)

第二章主要眼底圖像增強(qiáng)方法彩色眼底圖像增強(qiáng)對(duì)醫(yī)學(xué)診斷具有重要的作用，目前主要的彩色眼底圖像增強(qiáng)方法有：直方圖均衡化、對(duì)比度受限自適應(yīng)直方圖均衡化以及Hessian矩陣增強(qiáng)方法。2.1直方圖均衡化方法一般來說，圖像對(duì)比度的可用較為常見的兩種方法進(jìn)行增強(qiáng)處理，分別為間接對(duì)比度增強(qiáng)方法是直方圖拉伸方法和直方圖均衡化(HistogramEqualization，簡(jiǎn)稱HE)方法。對(duì)于直方圖均衡化而言，圖像灰度改變的是通過累積函數(shù)來實(shí)現(xiàn)的，以此達(dá)到增強(qiáng)對(duì)比度的效果。其基本的操作步驟的核心思路即，對(duì)原始圖像的非均質(zhì)化拉伸處理，使其像素值間距擴(kuò)張，均勻化各灰度范圍的像素量。這種方法也存在一些缺點(diǎn)：(1)增強(qiáng)后圖像的灰度級(jí)會(huì)變少，部分細(xì)節(jié)會(huì)消失；(2)當(dāng)輸入圖像的直方圖有非常密集的部分時(shí)，增強(qiáng)后的圖像的對(duì)比度會(huì)增強(qiáng)過度。通過直方圖均衡化，圖像的亮度可以更好地分布在直方圖上，讓圖像更易于觀察。用這種方法來增強(qiáng)圖像局部的對(duì)比度就不會(huì)使圖像整體的對(duì)比度產(chǎn)生影響，直方圖均衡化通過有均衡亮度密集的區(qū)域來實(shí)現(xiàn)這種功能。直方圖均衡化對(duì)增強(qiáng)背景太亮或者前景太暗的圖像有很好的效果，尤其是增強(qiáng)X光圖像中清晰度較差的骨骼結(jié)構(gòu)以及曝光過度和曝光不足的圖像中的細(xì)節(jié)信息。這種方法具有一個(gè)特殊優(yōu)勢(shì)是它的直觀性和可逆操作性，若均衡化的函數(shù)是已知的，則可以構(gòu)造出初始的直方圖。但該方法的缺點(diǎn)也很明顯，即必須對(duì)所有的數(shù)據(jù)進(jìn)行分析，這就可能會(huì)增加背景的對(duì)比度并且降低有用信息的對(duì)比度。圖像的直方圖可以表現(xiàn)出圖像像素值的分布規(guī)律。由于圖像是由大量像素組建而成，因而可以將像素分布的直方圖進(jìn)行列表統(tǒng)計(jì)來對(duì)其特征進(jìn)行分析研究。直方圖對(duì)圖像特征的提取和確定其相似度上都具有巨大的貢獻(xiàn)，它能通過對(duì)不同區(qū)間的像素值分布特征進(jìn)行整體上的調(diào)整，優(yōu)化其灰度分度，進(jìn)而達(dá)到增強(qiáng)圖像的視覺感。直方圖與圖像清晰度的有如下關(guān)系：(1)亮度不足，即代表其在直方圖中主要位于像素值較小區(qū)間；(2)亮度高，即表示其在直方圖中主要位于像素較大區(qū)間；(3)灰度級(jí)隨對(duì)比度的降低而降低，且中間水平的灰度級(jí)是主要信息的儲(chǔ)存區(qū)；(4)灰度級(jí)隨對(duì)比度的升高而升高，且主要信息呈均勻化分布。直方圖均衡化的基本思想是使輸入圖像的直方圖分布變的均勻，這樣就會(huì)使圖像的灰度級(jí)增加，從而可達(dá)到圖像對(duì)比度整體增強(qiáng)的效果。假設(shè)未處理前的圖像在二維坐標(biāo)系當(dāng)中，令其灰度值在(x,y)處為f，處理后其灰度值為g，則可認(rèn)為當(dāng)f轉(zhuǎn)變?yōu)間時(shí)，就實(shí)現(xiàn)了對(duì)圖像效果的增強(qiáng)。綜上，可將其轉(zhuǎn)換函數(shù)表示為：g=EQ(f)，這個(gè)映射函數(shù)EQ(f)必須滿足兩個(gè)條件(其中L為圖像的灰度級(jí)數(shù))：(1)為使輸入圖像的灰度排列不被打亂，在區(qū)間0≤f≤L-1內(nèi)，EQ(f)必為單調(diào)遞增函數(shù)。(2)為使轉(zhuǎn)換前后灰度不發(fā)生動(dòng)態(tài)變化，當(dāng)0≤f≤L-1時(shí)，必有0≤g≤L-1。HE方法如下：(1)設(shè)f、g分別為輸入圖像和增強(qiáng)處理后的圖像。HE方法第一步如圖2-1所示。計(jì)算輸入圖像f的灰度直方圖，設(shè)為h。h為一個(gè)256維的向量。圖2-1HE方法第一步(2)HE方法第二步如圖2-2所示。求出輸入圖像f的像素總數(shù)：(2-1)式(2-1)(m，n分別為圖像的長(zhǎng)和寬)不同灰度級(jí)范圍內(nèi)的像素量占整個(gè)圖像的比例可表示為：(2-2)圖2-2HE方法第二步(3)HE方法第三步如圖2-3所示。計(jì)算圖像的灰度級(jí)的累計(jì)分布hp：(2-3)圖2-3HE方法第三步(4)HE方法第四步如圖2-4所示。求出增強(qiáng)圖像g的灰度值：(2-4)圖2-4HE方法第四步圖2-5展示了直方圖均衡化的結(jié)果圖，上述方法處理前后的直方圖對(duì)比。從圖中可以發(fā)現(xiàn)，前后的直方圖有很大的變化，增強(qiáng)效果非常明顯，清晰程度與對(duì)比度增強(qiáng)。直方圖均衡化方法是一種十分有效的增強(qiáng)方法，它對(duì)于景物中輪廓背景鮮明，即亮度或暗度明顯的圖像非常有效果。但是這種方法對(duì)所處理的圖像的信息不做區(qū)分，這會(huì)導(dǎo)致增加復(fù)雜背景的信息，并且對(duì)需要增強(qiáng)的部分也有影響，會(huì)降低其對(duì)比度。將直方圖均衡化方法應(yīng)用于圖像去霧問題，雖然可以取得一定的增強(qiáng)效果，但是去霧后的整體視覺效果并不算理想。(a)(b)(c)(d)圖2-5直方圖均衡化前后及直方圖彩色眼底圖像及其直方圖均衡化處理結(jié)果示例如圖2-6所示。(a)(b)圖2-6HE增強(qiáng)(a)輸入圖像；(b)HE方法增強(qiáng)結(jié)果2.2對(duì)比度受限自適應(yīng)直方圖均衡化方法圖像增強(qiáng)技術(shù)可具體為時(shí)域和頻域。時(shí)域增強(qiáng)往往被應(yīng)用于提高圖像的對(duì)比度且改進(jìn)其灰度級(jí)，其原理是在灰度映射轉(zhuǎn)化的基礎(chǔ)上，對(duì)像素進(jìn)行直接性地處理；頻域增強(qiáng)的作用是以強(qiáng)化圖像的低、高頻來達(dá)到改善圖像的平滑性和邊緣為主，其原理是通過傅里葉轉(zhuǎn)換的方式提升興趣區(qū)的頻率分量。圖像灰度分布不均的問題對(duì)于圖像處理技術(shù)而言是一個(gè)影響較大的問題，為了能解決這個(gè)問題，基于時(shí)域增強(qiáng)理論，采用映射函數(shù)當(dāng)中的非線性函數(shù)，通過轉(zhuǎn)化后使其像素分布均勻。直方圖均衡化方法通過不同的灰度分布方案使轉(zhuǎn)變?yōu)榛叶戎稻鶆蚧植嫉膱D像，以眼底圖為例，利用直方圖提取灰度映射曲線，再對(duì)其進(jìn)行轉(zhuǎn)換，從而提升亮度。但該方法仍有瑕疵，即無法增強(qiáng)其對(duì)比度，且在此過程當(dāng)中，一并將噪聲信號(hào)也同時(shí)放大了。以AHE(AdaptiveHistogtamEqualization)為例，該方法將圖像利用網(wǎng)格線切割為數(shù)量眾多的小格子區(qū)域，并對(duì)每一個(gè)格子進(jìn)行均衡化處理。但該方法會(huì)導(dǎo)致圖像失真，且其微噪也會(huì)因此而放大。因此，這里提出了對(duì)比度受限制自適應(yīng)直方圖(ContrastlimitedAdaptiveHistogtamEqualization，CLAHE)，即對(duì)每一個(gè)劃分單元進(jìn)行對(duì)比度的限制處理，這也是該方法與傳統(tǒng)的AHE的不同之處[22]。CLAHE通過限制局部直方圖的高度來限制噪聲放大和局部對(duì)比度增強(qiáng)。該方法將圖像劃分為多個(gè)子區(qū)域；然后對(duì)每個(gè)子區(qū)域的直方圖進(jìn)行分類。再對(duì)每個(gè)子區(qū)域分別進(jìn)行直方圖均衡化，最后通過對(duì)每個(gè)像素進(jìn)行插值運(yùn)算來獲得變換后的灰度值，以此實(shí)現(xiàn)對(duì)比度受限自適應(yīng)直方圖均衡化圖像增強(qiáng)。CLAHE方法原理：CLAHE與AHE不同的地方是增加對(duì)比度限幅，這就可以克服AHE的過度放大噪聲的問題；(1)設(shè)置M×M大小的正方形滑動(dòng)窗口，局部映射函數(shù)可表示為：(2-5)(2)滑動(dòng)窗口局部直方圖的累積分布函數(shù)(cumulativedistributionfunction)：(2-6)從上式可看出，控制器直方圖高度，就能直接對(duì)映射函數(shù)的斜率進(jìn)行限制，以此控制對(duì)比度強(qiáng)度；(3)令限定的最大斜率為Smax，可得到直方圖最大高度為：(2-7)(2-8)(4)對(duì)高度大于Hmax的直方圖應(yīng)截去多余的部分：(2-9)實(shí)際處理中，設(shè)定閾值T對(duì)直方圖進(jìn)行截?cái)鄤澐郑撻撝挡坏扔贖max，將其截掉的部分在整體的灰階范圍內(nèi)分布均勻，達(dá)到總面積不變的情況下提升高度。(5)最后改進(jìn)的直方圖為：(2-10)綜上所述，改變最大的映射函數(shù)斜率Smax及相應(yīng)的最大直方圖高度Hmax，可獲得不同增強(qiáng)效果的圖像。CLAHE方法處理結(jié)果如圖2-7所示。(a)(b)圖2-7CLAHE方法(a)CLAHE處理前的直方圖；(b)CLAHE處理后的直方圖如上所述的處理方法，據(jù)綜合分析對(duì)比后，發(fā)現(xiàn)無論是否有無對(duì)比度的限制處理，都必須要對(duì)任一像素去進(jìn)行計(jì)算其直方圖和變換函數(shù)，這就導(dǎo)致該方法極其耗時(shí)。而插值方法大大降低了上述方法的運(yùn)算速度，并且增強(qiáng)的圖像的質(zhì)量并不會(huì)下降。插值方法如圖2-8所示。首先，將圖像均勻分成多分相同大小的矩形，如下圖的右側(cè)部分所示。然后計(jì)算一個(gè)矩形的直方圖和對(duì)應(yīng)的變換函數(shù)。對(duì)于其他像素而言，可通過其周圍矩形中的變換函數(shù)進(jìn)行插值處理得到。其中，采用雙線性插值區(qū)域?yàn)樗{(lán)色矩形部分，采用線性插值區(qū)域?yàn)闇\綠色矩形部分，而對(duì)于淺紅色的矩形，可直接利用其變換函數(shù)。圖2-8插值方法按照這種過程就可以大大地降低變換函數(shù)所需的計(jì)算次數(shù)，這種方法只增加了部分雙線性插值的計(jì)算。采用CLAHE方法對(duì)彩色眼底圖像進(jìn)行增強(qiáng)，結(jié)果示例如圖2-9所示，彩色眼底圖像的飽和度更好，血管紋路更清晰。(a)(b)圖2-9CLAHE方法增強(qiáng)(a)輸入圖像；(b)對(duì)比度受限制自適應(yīng)直方圖處理2.3Hessian矩陣增強(qiáng)方法Hessian矩陣增強(qiáng)的核心是采用求解高階微分的方式來對(duì)圖像特征進(jìn)行提取。它認(rèn)為，最大模量的特征向量和圖像特征在方向性上具有垂直關(guān)系。此外，它可由具有線正交的較大學(xué)對(duì)峙的二階導(dǎo)和沿線的具有小絕對(duì)值的二階導(dǎo)表示出由高斯函數(shù)構(gòu)造的線性函數(shù)，這即是二維Hessian的幾何意義。Hessian矩陣已經(jīng)被應(yīng)用于多種檢測(cè)和分析，并且還被用于分割和重建醫(yī)學(xué)圖像中的曲線結(jié)構(gòu)。關(guān)于判斷圖像中出現(xiàn)的點(diǎn)是否為角點(diǎn)，即像素密集程度激烈變化的點(diǎn)，可利用Hessian矩陣的特征值來進(jìn)行判斷。Hessian矩陣可表示為：(2-11)X方向上的二階偏微分：(2-12)Y方向上的二階偏微分：(2-13)X，Y方向上的混合偏微分：(2-14)對(duì)于一幅圖像I而言，它只有x，y兩個(gè)方向，所以其Hessian矩陣是一個(gè)二元矩陣，對(duì)應(yīng)的有：(2-15)分別對(duì)圖像卷積運(yùn)算，然后構(gòu)成圖像Hessian矩陣：(2-16)(2-17)2.4本章小結(jié)為彌補(bǔ)直方圖均衡化會(huì)使灰度級(jí)降低，丟失部分圖像細(xì)節(jié)等缺陷，可利用CLAHE方法進(jìn)行處理。該方法極大程度上提高了對(duì)比度，但仍然有著缺點(diǎn)，即同時(shí)放大了噪聲，且無法精準(zhǔn)識(shí)別病灶區(qū)域。圖像中亮、暗及其余細(xì)節(jié)部分的提取，可利用數(shù)學(xué)形態(tài)學(xué)濾波，但該方法處理對(duì)象僅為灰度圖，對(duì)彩色眼底圖像增強(qiáng)有限。然而，想要實(shí)現(xiàn)對(duì)視網(wǎng)膜圖像的全方位增強(qiáng)，利用Hessian矩陣顯示是不適用的。首先，傳統(tǒng)的眼底圖增強(qiáng)僅僅是對(duì)血管信息的強(qiáng)化，鮮有能夠同時(shí)對(duì)諸如視覺神經(jīng)、病灶等全方位的信息進(jìn)行全面強(qiáng)化；此外，傳統(tǒng)的增強(qiáng)方法絕大多數(shù)的增強(qiáng)原理是對(duì)彩色圖像對(duì)灰度圖的轉(zhuǎn)化，極其難以使得處理后的圖像在視覺上的保真性。因此，必須考慮多因素的影響條件，以達(dá)到突出圖像細(xì)節(jié)的基礎(chǔ)上，還能對(duì)其動(dòng)態(tài)范圍及相關(guān)色彩性質(zhì)加以完善。

第三章基于Retinex理論的方法彩色眼底圖像增強(qiáng)在醫(yī)學(xué)領(lǐng)域具有重要的作用，研究發(fā)現(xiàn)可以采用基于Retinex理論的方法來對(duì)彩色眼底圖像進(jìn)行增強(qiáng)，例如單尺度Retinex方法、多尺度Retinex方法、帶色彩恢復(fù)的多尺度Retinex方法以及基于迭代計(jì)算的Retinex方法。3.1單尺度Retinex方法(SSR)單尺度Retinex方法和多尺度Retinex方法會(huì)造成顏色細(xì)節(jié)丟失；帶色彩恢復(fù)的多尺度Retinex方法與MSR相比引入色彩恢復(fù)因子，這雖然起到了抑制光暈的作用，但是會(huì)導(dǎo)致圖像顏色失真。而且這種方法雖然處理彩色眼底圖像效果較好，但是計(jì)算量卻非常大，在處理上需要耗費(fèi)較多時(shí)間，不能達(dá)到醫(yī)學(xué)上的實(shí)時(shí)性需求。SSR和MSR方法雖然效率高，但往往由于圖像本身或者處理時(shí)對(duì)尺度的把控不合理，使得圖像增強(qiáng)后對(duì)比度低。為了克服以上的各種問題，眾多學(xué)者不斷進(jìn)行深入的研究。如今，Retinex方法己經(jīng)越來越豐富，越來越完善。如圖3-1所示，其Retinex理論可認(rèn)為是入射光經(jīng)物體反射后形成反射光。圖3-1Retinex理論示意圖光線經(jīng)物體反射后形成的反射光被人眼所接收，經(jīng)大腦處理后形成圖像，用公式可表示為：(3-1)式(3-1)中，L(x,y)表示入射光，決定像素的動(dòng)態(tài)范圍；R(x,y)表示物體的反射光，表示為圖像的本性，即內(nèi)在屬性；S(x,y)為人眼所能接收到的經(jīng)反射后的圖像。Retinex最核心的思想即消除或降低源于入射光所造成的影響，以此來保留物體本質(zhì)的反射光圖像。其公式為：(3-2)式(3-2)中：i=1，2，3分別代表彩色圖像的三原色R，G，B；*表示卷積運(yùn)算；G(x,y)是高斯環(huán)繞函數(shù)，可表示為：(3-3)式(3-3)中：c表示高斯環(huán)繞函數(shù)的尺度參數(shù)；K為歸一化因子。高斯環(huán)繞函數(shù)需滿足下式的成立：(3-4)Retinex方法的基本流程如圖3-2所示。圖3-2Retinex方法的基本流程SSR方法的具體實(shí)現(xiàn)步驟：(1)讀取輸入圖像數(shù)據(jù)S(x,y)，對(duì)其數(shù)據(jù)進(jìn)行類型變換，將整型變換為double型；(2)確定尺度c的大小，并根據(jù)式(3-4)確定歸一化因子K的取值；(3)根據(jù)式(3-2)來獲得圖像R(x,y)；(4)對(duì)R(x,y)進(jìn)行線性校正處理(因?yàn)镽(x,y)的范圍不在0-255內(nèi))，校正后的圖像即為方法增強(qiáng)圖像。彩色眼底圖像的SSR處理示例圖像如圖3-3所示。(b)圖3-3彩色眼底圖像的SSR處理圖像(a)輸入圖像；(b)SSR處理后圖像3.2多尺度Retinex方法(MSR)MSR是在SSR的基礎(chǔ)上發(fā)展而來的。在SSR圖像增強(qiáng)中，由于單尺度對(duì)圖像增強(qiáng)影響太大，經(jīng)其處理后的動(dòng)態(tài)壓縮性、保真性和細(xì)節(jié)均無法達(dá)到最優(yōu)化。為了消除這種增強(qiáng)方式帶來的負(fù)面影響，可結(jié)合不同尺度下經(jīng)SSR處理所得結(jié)果，即采用MSR方法，這樣就能使圖像在顏色和細(xì)節(jié)信息方面達(dá)到平衡。MSR方法是經(jīng)SSR進(jìn)行加權(quán)處理后的結(jié)果，其可表示為：(3-5)式(3-5)中，為MSR在第i個(gè)顏色通道的輸出；Wn為不同尺度的加權(quán)求和系數(shù)；N為選擇的尺度個(gè)數(shù)。在選擇尺度時(shí)，需要選擇一個(gè)小尺度、一個(gè)中尺度和一個(gè)大尺度，用來保證實(shí)驗(yàn)結(jié)果能有較好的增強(qiáng)效果。而選擇尺度有如下規(guī)律：(1)選擇圖像大小的1%~5%，使實(shí)驗(yàn)結(jié)果能夠得到很好的邊緣細(xì)節(jié)，作為小尺度；(2)選擇圖像大小的10%~15%，使實(shí)驗(yàn)結(jié)果能夠得到邊緣細(xì)節(jié)和色彩，作為中尺度；(3)選擇圖像大小的30%~50%，使實(shí)驗(yàn)結(jié)果能夠得到比較均衡的色彩，作為大尺度。多尺度Retinex的基本步驟如下：(1)初始化處理。將S(x,y)分解為R、G、B，轉(zhuǎn)化其數(shù)據(jù)類型；(2)將輸入的數(shù)據(jù)轉(zhuǎn)換至對(duì)數(shù)域，由于對(duì)數(shù)的特殊性，即無法取0作為其變量，則可在原有數(shù)據(jù)基礎(chǔ)上加上一個(gè)修正值，一般可取為1；(3)選擇不同尺寸的參數(shù)，根據(jù)式3-4來確定c值；再根據(jù)式3-3計(jì)算高斯環(huán)繞函數(shù)，并對(duì)原始對(duì)象進(jìn)行卷積計(jì)算，得到濾波處理后的圖像I(x,y)；(4)進(jìn)行加權(quán)求和計(jì)算，然后經(jīng)過對(duì)數(shù)處理；(5)計(jì)算步驟(4)與原始圖像之間的相對(duì)關(guān)系，可理解為是二者對(duì)數(shù)域當(dāng)中的數(shù)值相減；(6)將圖像從對(duì)數(shù)域轉(zhuǎn)換到實(shí)數(shù)域；(7)通過gain/offset對(duì)圖像進(jìn)一步調(diào)整并作線性拉伸，將R、G、B三色圖像重疊后再進(jìn)行輸出。彩色眼底圖像的MSR處理示例如圖3-4所示。(a)(b)圖3-4彩色眼底圖像的MSR處理圖像(a)輸入圖像；(b)MSR處理后圖像3.3帶色彩恢復(fù)的多尺度Retinex方法(MSRCR)盡管對(duì)于灰度圖像的處理，MSR方法表現(xiàn)出了較好的效果，但對(duì)于彩色圖像而言卻不盡人意，往往會(huì)出現(xiàn)失真、發(fā)灰和泛白等現(xiàn)象。在之前的增強(qiáng)過程中，圖像可能在處理時(shí)增加了噪聲，使圖像的局部細(xì)節(jié)顏色失真，失去了彩色眼底圖像真正的顏色，整體增強(qiáng)效果較差。為了克服MSR對(duì)彩色圖像增強(qiáng)的不足，MSRCR應(yīng)運(yùn)而生。MSRCR方法利用色彩恢復(fù)因子c，調(diào)整輸入的彩色眼底圖像中三個(gè)顏色通道之間的比例關(guān)系，把較暗區(qū)域的信息增強(qiáng)，達(dá)到了消除圖像顏色失真的缺陷。處理后的圖像局部對(duì)比度提高，亮度與真實(shí)場(chǎng)景相似，在人們視覺感知下，圖像顯得更加逼真。但是MSRCR方法處理圖像后，彩色眼底圖像像素值一般會(huì)出現(xiàn)負(fù)值。Jobson等提出了一個(gè)利用色彩恢復(fù)函數(shù)與多尺度MSR的輸出結(jié)果相乘的方法，其可表示為：(3-6)式(3-6)中Ci(x,y)稱為彩色恢復(fù)函數(shù)，可以表示為：(3-7)式(3-7)中，β為增益常數(shù)；α為非線性強(qiáng)度控制因子。彩色眼底圖像的MSRCR處理示例如圖3-5所示。(a)(b)圖3-5彩色眼底圖像的MSRCR處理圖像(a)輸入圖像；(b)MSRCR處理圖像3.4基于迭代計(jì)算的Retinex方法基于迭代計(jì)算的Retinex方法有Frankle-McCannRetinex方法和McCann99Retinex方法。1983年A.Frankle和J.McCann發(fā)明了Frankle-McCannRetinex方法，簡(jiǎn)稱McCannRetinex方法，這種方法在采樣上精度更高，選擇的是一條螺旋形結(jié)構(gòu)的迭代路徑，這種路徑選擇包含了整個(gè)圖像的全局明暗關(guān)系，F(xiàn)rankle-McCannRetinex方法路徑選擇圖如圖3-6所示。由圖3-6可知，坐標(biāo)零點(diǎn)附近的選取點(diǎn)數(shù)較多，這是由于零點(diǎn)附近的點(diǎn)與中心點(diǎn)的相關(guān)性較遠(yuǎn)處相比，更容易獲得精確的光照估計(jì)。圖3-6Frankle-McCannRetinex方法路徑選擇圖McCannRetinex方法首先利用預(yù)先規(guī)定的距離D分離首先處理的像素，然后順時(shí)針旋轉(zhuǎn)比較方向90度，則D為初始的一半，如此反復(fù)直至到單位像素的距離。其迭代步數(shù)是人為規(guī)定的。在對(duì)彩色眼底圖像處理時(shí)的具體步驟如下：將R、G、B三個(gè)通道進(jìn)行對(duì)數(shù)化處理，將數(shù)據(jù)類型轉(zhuǎn)化為double；(2)確定shift變量，進(jìn)行路徑比較；(3)重新對(duì)shift賦值，即shift=1/2shift；(4)重復(fù)(2)、(3)步驟，直到shift<1；(5)設(shè)定迭代次數(shù)n，重復(fù)(2)、(3)和(4)步驟n次；(6)將輸入圖像有對(duì)數(shù)域轉(zhuǎn)化到實(shí)數(shù)域，顯示圖像。Funt等人提出McCann99Retinex方法，這種方法是多分辨率的，它利用金字塔模型來對(duì)圖像進(jìn)行描述，通過從頂部向下逐漸迭代，以此來達(dá)到效果強(qiáng)化，但這種方式對(duì)