眼底圖像的動(dòng)靜脈視網(wǎng)膜血管分割方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及計(jì)算機(jī)輔助診斷技術(shù)領(lǐng)域，具體涉及一種眼底圖像的動(dòng)靜脈視網(wǎng)膜血管分割方法。
【背景技術(shù)】
[0002]隨著計(jì)算機(jī)技術(shù)中的人工智能領(lǐng)域的快速發(fā)展，計(jì)算機(jī)輔助診斷技術(shù)也逐漸發(fā)展。計(jì)算機(jī)輔助診斷技術(shù)是指通過(guò)影像學(xué)、醫(yī)學(xué)圖像處理技術(shù)以及其他可能的生理、生化手段，結(jié)合計(jì)算機(jī)的分析計(jì)算，輔助影像科醫(yī)師發(fā)現(xiàn)病灶，提高診斷的準(zhǔn)確率。
[0003]通常醫(yī)學(xué)影像學(xué)中計(jì)算機(jī)輔助診斷分為三步，具體如下:第一步是把病變從正常結(jié)構(gòu)中提取出來(lái)；第二步是圖像特征的量化；第三步是對(duì)數(shù)據(jù)進(jìn)行處理并得出結(jié)論。
[0004]因?yàn)橛?jì)算機(jī)可以全面利用影像信息進(jìn)行精確的定量計(jì)算，去除人的主觀性，避免因個(gè)人知識(shí)和經(jīng)驗(yàn)的差異而引起的“千差萬(wàn)別”的診斷結(jié)果；所以它的結(jié)果是不含糊的，是確定的，它使診斷變得更為準(zhǔn)確、更為科學(xué)。
[0005]隨著現(xiàn)代高科技的發(fā)展，計(jì)算機(jī)輔助診斷將與圖像處理和PACS系統(tǒng)等技術(shù)融合，變得更易于操作、也更趨于準(zhǔn)確，其臨床應(yīng)用范圍將進(jìn)一步擴(kuò)大。
[0006]在醫(yī)學(xué)檢測(cè)中，眼睛是唯一可無(wú)損檢測(cè)同時(shí)信息豐富的器官。研宄指出視網(wǎng)膜血管病變中的血管局限縮窄、彌漫縮窄、動(dòng)靜脈交叉壓迫、血管行走改變、銅絲動(dòng)脈、出血、棉絮斑、硬性滲出以及視網(wǎng)膜神經(jīng)纖維層缺損與腦卒有顯著的相關(guān)性。且對(duì)于腦卒中的預(yù)測(cè)，眼底檢查僅需40元，而MRI檢查則需要上千元，頸動(dòng)脈超聲也需要140元。相比之下眼底檢查的性價(jià)比最高。眼底圖像計(jì)算機(jī)分析的全自動(dòng)化的方法，包括可以提供即時(shí)的視網(wǎng)膜病變分類，而不需要專家意見(jiàn)，建立以眼底血管視神經(jīng)預(yù)測(cè)三高并發(fā)癥的系統(tǒng)具有其確實(shí)的經(jīng)濟(jì)意義。因此，視網(wǎng)膜血管的病變檢測(cè)在對(duì)腦卒的輔助檢測(cè)具有突出作用。其中構(gòu)建一個(gè)動(dòng)靜脈交叉壓迫視網(wǎng)膜血管病變的自動(dòng)檢測(cè)系統(tǒng)更是其中的關(guān)鍵部分。
[0007]對(duì)眼底圖像進(jìn)行血管分割、視盤(pán)定位和血管分類(動(dòng)靜脈分裂)是視網(wǎng)膜血管的病變檢測(cè)的基礎(chǔ)，現(xiàn)有的血管分割方法需要人工添加標(biāo)注信息，自動(dòng)化程度不高。

【發(fā)明內(nèi)容】

[0008]針對(duì)現(xiàn)有技術(shù)的不足，本發(fā)明提供了一種眼底圖像的動(dòng)靜脈視網(wǎng)膜血管分割方法。
[0009]一種眼底圖像的動(dòng)靜脈視網(wǎng)膜血管分割方法，包括如下步驟:
[0010](1-1)按照預(yù)設(shè)的二值化閾值對(duì)經(jīng)過(guò)預(yù)處理的眼底圖像進(jìn)行二值化處理，并提取二值化處理后的眼底圖像中的中心線和邊緣，得到血管樹(shù)；
[0011](1-2)對(duì)所述的血管樹(shù)分叉處做斷開(kāi)處理得到血管段，并對(duì)各個(gè)血管段進(jìn)行線分割得到血管，得到原始血管集；
[0012](1-3)確定誤分割血管，并從原始血管集中去除得到全局血管集。
[0013]本發(fā)明中對(duì)眼底圖像進(jìn)行小波變換以對(duì)眼底圖像進(jìn)行預(yù)處理。
[0014]二值化閾值為二值化處理后為血管的像素點(diǎn)個(gè)數(shù)占整個(gè)眼底圖像的像素點(diǎn)比例，通常取值為4?20%。二值化閾值越大，則越寬松。作為優(yōu)選，所述的二值化閾值為14%。
[0015]所述步驟(1-3)中的誤分割血管包括分為兩類，第一類誤分割血管基于血管兩側(cè)的背景差異確定，第二類基于血管形狀確定。
[0016]本發(fā)明中所述的背景差異值背景顏色的差異，步驟(1-3)中通過(guò)如下步驟確定第一類誤分割血管:
[0017](al)針對(duì)每個(gè)血管，提取該血管兩側(cè)背景的特征向量；
[0018]任意一側(cè)背景的特征向量根據(jù)該側(cè)距離中心線5?10個(gè)像素點(diǎn)以內(nèi)區(qū)域中所有像素點(diǎn)(即與中心線的距離小于5?10個(gè)像素點(diǎn))的RGB三個(gè)通道上的顏色值得到。
[0019]每側(cè)的特征向量為一個(gè)三維向量，分別表示血管兩側(cè)背景的在RGB三個(gè)通道上的顏色值信息。具體實(shí)現(xiàn)時(shí)獲取該側(cè)距離中心線5?10個(gè)像素點(diǎn)的區(qū)域中的所有像素點(diǎn)在R、G、B三個(gè)通道上的顏色值并分別在每個(gè)通道上求平均，進(jìn)而得到該側(cè)的特征向量。
[0020](a2)根據(jù)所述特征向量的采用聚類法將所有血管聚為兩類，得到的小類即為視盤(pán)周?chē)`分割血管。
[0021]針對(duì)每個(gè)血管，計(jì)算兩側(cè)特征向量的歐氏距離，然后對(duì)所有血管對(duì)應(yīng)的歐式距離進(jìn)行聚類，即完成對(duì)血管的聚類。
[0022]由于K均值聚類(即K-means算法)不需要調(diào)整參數(shù)，且運(yùn)行速度較快。作為優(yōu)選，所述步驟(a2)中采用K均值聚類法將所有血管聚為兩類。
[0023]本發(fā)明中所述的血管形狀實(shí)際上指各個(gè)血管的連接關(guān)系，步驟(1-3)中通過(guò)如下步驟確定第二類誤分割血管:
[0024]確定劃分出原始血管集的眼底圖像中的環(huán)狀結(jié)構(gòu)，針對(duì)各個(gè)環(huán)狀結(jié)構(gòu)，若該環(huán)狀結(jié)構(gòu)中長(zhǎng)度最大的血管的長(zhǎng)度小于預(yù)設(shè)的分割長(zhǎng)度閾值，則認(rèn)為該環(huán)狀結(jié)構(gòu)中所有血管均為第二類誤分割血管，進(jìn)一步進(jìn)行如下操作:
[0025]確定該環(huán)狀結(jié)構(gòu)的中心，并計(jì)算該中心到長(zhǎng)度大于或等于分割長(zhǎng)度閾值的血管的最短距離(即該中心到距離其最近的長(zhǎng)度大于或等于分割長(zhǎng)度閾值的血管的距離)，認(rèn)為以該中心為圓心、最短距離為半徑的圓形區(qū)域內(nèi)所有血管為第二類誤分割血管。
[0026]分割長(zhǎng)度閾值的設(shè)定根據(jù)實(shí)際經(jīng)驗(yàn)值設(shè)定，通常為分割長(zhǎng)度閾值為α = χ/60?x/45，X為眼底圖像的橫向大小，即眼底圖像沿橫軸的長(zhǎng)度。
[0027]未作特殊說(shuō)明，本發(fā)明中對(duì)長(zhǎng)度、距離、圖片大小等參數(shù)進(jìn)行衡量時(shí)統(tǒng)一以像素點(diǎn)為單位。
[0028]與現(xiàn)有技術(shù)相比，本發(fā)明在得到原始血管集后進(jìn)一步還利用血管的背景和形狀確定誤分割的血管，可以有效去除因?yàn)榕恼赵斐傻沫h(huán)狀反光、視盤(pán)周?chē)姆茄艿能S階邊緣、斑狀病變以及出血病變等原因造成的誤分割血管。
【附圖說(shuō)明】
[0029]圖1為本實(shí)施例的眼底圖像；
[0030]圖2為本實(shí)施例中對(duì)眼底圖像進(jìn)行血管分割的流程圖；
[0031]圖3為血管分割得到的原始血管集的示意圖；
[0032]圖4為血管分割得到的全局血管集的示意圖。
【具體實(shí)施方式】
[0033]下面將結(jié)合附圖和具體實(shí)施例對(duì)本發(fā)明進(jìn)行詳細(xì)描述。
[0034]本實(shí)施例以圖1所示的眼底圖像為例來(lái)說(shuō)明眼底圖像的動(dòng)靜脈視網(wǎng)膜血管分割方法，該眼底圖像的大小為3000 X 3000ο由拍照造成的環(huán)狀反光、視盤(pán)周?chē)姆茄艿能S階邊緣、斑狀病變以及出血病變等原因，該眼底圖像中存在亮環(huán)。
[0035]對(duì)該眼底圖像進(jìn)行動(dòng)靜脈視網(wǎng)膜血管分割的分類流程如圖2所示，包括如下步驟:
[0036](1-1)對(duì)眼底圖像進(jìn)行小波變換(IUWT小波)，按照預(yù)設(shè)的二值化閾值對(duì)經(jīng)過(guò)小波變換的眼底圖像進(jìn)行二值化處理，并提取二值化處理后的眼底圖像中的中心線和邊緣，得到血管樹(shù)；
[0037](1-2)對(duì)血管樹(shù)分叉處做斷開(kāi)處理得到血管段，并對(duì)各個(gè)血管段進(jìn)行線分割得到血管，得到原始血管集。
[0038]對(duì)血管樹(shù)分叉處做斷開(kāi)處理時(shí):當(dāng)血管樹(shù)中的血管中心線中多根中心線匯集到一點(diǎn)時(shí)，去除中心點(diǎn)(匯集的交叉點(diǎn))，得到單獨(dú)的多根血管中心線。
[0039]對(duì)各個(gè)血管段進(jìn)行線分割時(shí):以每一根中心線作為一個(gè)血管段。血管段為一條曲線，運(yùn)用圖像處理的線分割的傳統(tǒng)方法，將曲線用多根直線逼近。得到的多根直線，每根直線即代表一根血管，所有直線的集合即為原始血