基于區(qū)域生長(zhǎng)算法人工心臟瓣膜有效開口面積測(cè)量

1、基于區(qū)域生長(zhǎng)算法人工心臟瓣膜有效開口面積測(cè)量摘 要：人工心臟瓣膜從研制到臨床應(yīng)用必須通過體外血流動(dòng)力學(xué)性能檢測(cè)和評(píng)價(jià)，有效開口面積是重要指標(biāo)之一。傳統(tǒng)有 效開口面積測(cè)量方法自動(dòng)化程度低，人工干預(yù)度高。為了提高人工心臟瓣膜有效開口面積的測(cè)量效率，本文提出了基于自動(dòng) 種子點(diǎn)設(shè)置的區(qū)域生長(zhǎng)算法有效開口面積計(jì)算方法。通過建立搜索窗格，改進(jìn)的區(qū)域生長(zhǎng)算法對(duì)開口面積的分割有較好的 適應(yīng)性。結(jié)果表明：該方法能夠準(zhǔn)確計(jì)算出人工心臟瓣膜不同開合狀態(tài)下開口面積，為分析開口面積與流量、跨瓣壓差的關(guān) 系提供數(shù)據(jù)支撐。關(guān)鍵詞：人工心臟瓣膜；有效開口面積；自動(dòng)種子點(diǎn)；區(qū)域生長(zhǎng)Measurement of effectiv

2、e orifice area of artificial heart valvebased on region growing algorithm【Abstract】From the development to clinical application of artificial heart valve, it is necessary to measure and evaluate the hemodynamic performance in vitro, and the effective orifice area is one of the important indexes. The

3、 traditional effective orifice area measurement method has low automation and high manual intervention. In order to improve the efficiency of measuring the effective orifice area of artificial heart valve, a region growing algorithm based on automatic seed point setting is proposed. Through the esta

4、blishment of search pane, the improved region growing algorithm has better adaptability to the segmentation of orifice area. The results show that this method can accurately calculate the o rif ice area of the pro sthetic heart valve under dif f erent o pening and clo sing conditions, and provide da

5、ta support for the analysis of the relationship between the orifice area, flow rate, and cross valve pressure difference.【Key words 】Artificial heart valve; Effective orifice area; Auto seed point; Region growingo引言人工心臟瓣膜置換手術(shù)是治療瓣膜性心臟病患 者最有效的醫(yī)療措施。從研制到投入臨床使用，人 工心臟瓣膜需要經(jīng)過體外流體動(dòng)力學(xué)性能檢測(cè)和評(píng) 價(jià)，有效開口面積則是該檢測(cè)和評(píng)價(jià)的重

6、要指標(biāo)之 一。人工心臟瓣膜結(jié)果和功能有其固有的特殊性， 心臟瓣膜開口速度和大小受到博出量、跨瓣壓差的 影響，延用Gorlin公式計(jì)算開口面積,通過理論分析 與試驗(yàn)之后發(fā)現(xiàn)，流出系數(shù)改變較大，實(shí)際值與理論 值有較大誤差。隨著試驗(yàn)設(shè)備的完善，圖像采集處 理系統(tǒng)與人工心臟瓣膜體外脈動(dòng)脈測(cè)試系統(tǒng)相結(jié) 合，在采集跨瓣壓差、流量的同時(shí)，實(shí)時(shí)拍攝人工心 臟瓣膜開閉過程中的圖像，進(jìn)一步通過圖像處理計(jì) 算得到人工心臟瓣膜的有效開口面積，并就有效開 口面積與流量、跨瓣壓差的關(guān)系進(jìn)行分析。目前人工心臟瓣膜有效開口面積的計(jì)算主要依 賴于人工后期處理，步驟繁瑣，人工干預(yù)程度高，效 率低下。因此,探討提高人工心臟瓣膜的有

7、效開口 面積的測(cè)量效率和自動(dòng)化程度的新方法有重要意 義。本文通過圖像處理對(duì)采集到的心臟瓣膜開合運(yùn) 動(dòng)過程進(jìn)行分析，通過自動(dòng)種子點(diǎn)選取的區(qū)域生長(zhǎng) 算法計(jì)算人工心臟瓣膜有效開口面積。1研究?jī)?nèi)容和方法1.1系統(tǒng)組成及實(shí)驗(yàn)環(huán)境人工心臟瓣膜體外脈動(dòng)測(cè)試系統(tǒng)主要包括血液 循環(huán)體外模擬裝置、信號(hào)放大儀、電磁流量?jī)x、CCD 相機(jī)、圖像采集卡等，如圖1所示。圖1測(cè)試系統(tǒng)結(jié)構(gòu)圖Fig. 1 Test system structure diagram血液循環(huán)體外模擬裝置是與上海微創(chuàng)醫(yī)療器械 （集團(tuán)）有限公司合作研發(fā)，該裝置通過模擬人體全 身循環(huán)脈動(dòng)模型，可以產(chǎn)生近似生理?xiàng)l件下的脈動(dòng) 壓力和流量波形，其裝置如圖2所示

8、。圖像采集使 用的是OLYMPUS i-SPEED LT高速攝像機(jī),采集速 率為500幀/s,分辨率為800 x600;圖像處理環(huán)境為 Windows 10操作系統(tǒng)與MATLAB R2016P軟件平 臺(tái)。實(shí)驗(yàn)測(cè)試的人工心臟瓣膜由上海微創(chuàng)醫(yī)療器械 （集團(tuán)）有限公司研發(fā)提供。圖2血液循環(huán)體外模擬裝置Fig. 2 Extracorporeal simulation device of blood circulation 1.2圖像預(yù)處理工心臟瓣膜開閉速度較快，同時(shí)受攝像機(jī)分辨 率及光源照射等因素的影響，攝入的圖像存在幾何 畸變、邊緣模糊、圖像亮度分布不勻等問題，影響有 效開口面積的計(jì)算。通過圖像的預(yù)

9、處理可以消除無 關(guān)信息，突出環(huán)境背景下的特征參數(shù),增強(qiáng)目標(biāo)信息 的可識(shí)別性，有效降低了識(shí)別算法的復(fù)雜度。1.2.1圖像灰度化圖像灰度化處理是將三通道的RGB轉(zhuǎn)化為一 通道的灰度圖，灰度變化范圍縮減為028,大大減 少圖像信息數(shù)據(jù)量，更好的凸顯圖像形態(tài)學(xué)特性。 針對(duì)人工心臟瓣膜拍攝的圖片特性，采用加權(quán)平均 灰度法對(duì)圖像進(jìn)行處理，具有比較好的適應(yīng)性。1.2.2 形態(tài)學(xué)處理形態(tài)學(xué)圖像處理是利用數(shù)學(xué)形態(tài)學(xué)從給定圖像 中提取目標(biāo)有用圖像分量。其基本的運(yùn)算包括:腐 蝕、膨脹、開運(yùn)算、閉運(yùn)算、頂帽操作、擊中或擊不中 變化等。在灰度級(jí)圖像中，腐蝕和膨脹是利用結(jié)構(gòu)元原 點(diǎn)，在每個(gè)像素位置上與灰度級(jí)圖像重合區(qū)域進(jìn)

10、行 最小值、最大值運(yùn)算，最后的運(yùn)算效果取決于 結(jié)構(gòu)元大小內(nèi)容以及運(yùn)算性質(zhì)?；叶燃?jí)圖像的開操作，先利用結(jié)構(gòu)元進(jìn)行腐蝕， 隨后對(duì)所得結(jié)果膨脹操作，灰度級(jí)圖像的閉操作則 與之相反。開操作可以平滑圖像輪廓，斷開較窄的 區(qū)域，消除細(xì)小的孤島和突刺。閉操作操作可以平 滑圖像輪廓，融合狹窄的間斷和細(xì)長(zhǎng)的溝壑,消除小 孔，填補(bǔ)斷裂。原始目標(biāo)灰度級(jí)圖像中包括開口陰影、鈦合金 金屬支架、人工瓣膜、瓣膜安裝橡膠圈以及其他結(jié) 構(gòu)，如圖3所示。其中,鈦合金金屬支架與開口陰影 部分有部分重合且灰度級(jí)相近，直接使用圖像分割 不能達(dá)到預(yù)期的分割效果。圖3原始灰度級(jí)圖像Fig. 3 Original gray level im

11、age采用形態(tài)學(xué)閉操作，對(duì)原始目標(biāo)灰度級(jí)圖像進(jìn) 行處理，結(jié)構(gòu)元選擇半徑為7的平面圓盤形的結(jié)構(gòu) 元素，如圖4所示。經(jīng)過閉運(yùn)算處理之后，亮部細(xì)節(jié) 和背景未受影響，暗部特征明顯削弱，鈦合金金屬支 架能夠成功的從原始目標(biāo)灰度級(jí)圖像消除。圖4形態(tài)學(xué)閉操作圖像Fig. 4 Morphological closed operation image經(jīng)過形態(tài)學(xué)閉操作之后,開口陰影部分仍然存 在一些亮部噪點(diǎn)沒有消除，繼而采用形態(tài)學(xué)開操作, 結(jié)構(gòu)元選擇半徑為7的平面圓盤形的結(jié)構(gòu)元素，如 圖5所示。經(jīng)過灰度級(jí)圖像的閉操作和開操作，在保留開 口陰影部分的同時(shí),鈦合金金屬支架和開口陰影部 分的亮部噪點(diǎn)可以從原始目標(biāo)灰度級(jí)

12、圖像去除。1.3區(qū)域生長(zhǎng)圖像分割1.3.1算法描述圖像分割一般是基于灰度值的二個(gè)基本性質(zhì)： 不連續(xù)性和相似性。不連續(xù)性是基于灰度突變來分 割圖像，比如圖像的邊緣;相似性是根據(jù)預(yù)定的準(zhǔn)則分割圖像相似區(qū)域1_2（區(qū)域生長(zhǎng)算法就是基于灰 度值的相似性，根據(jù)預(yù)先定義的生長(zhǎng)準(zhǔn)則，將像素或 子區(qū)域組分割合成為更大區(qū)域的過程。分割后 的同一區(qū)域中的像素具有相同的灰度級(jí)，兩鄰接區(qū) 域的灰度值不同。圖5形態(tài)學(xué)開操作圖像Fig. 圖5形態(tài)學(xué)開操作圖像Fig. 5 Morphological open operation image區(qū)域生長(zhǎng)算法主要包括區(qū)域生長(zhǎng)點(diǎn)的選擇、區(qū) 域生長(zhǎng)準(zhǔn)則的確定以及終止生長(zhǎng)規(guī)則的確定。生

13、長(zhǎng)種子點(diǎn)可以基于問題的性質(zhì)選擇一組或者 多組。傳統(tǒng)的區(qū)域生長(zhǎng)算法在選擇種子點(diǎn)時(shí)，需要 人為針對(duì)不同的圖像情況設(shè)置不同位置的種子點(diǎn)， 增加測(cè)試人員的工作量，降低了檢測(cè)效率4項(xiàng)。本 文結(jié)合人工心臟瓣膜圖像特性，可以自動(dòng)選擇種子 點(diǎn)進(jìn)行區(qū)域生長(zhǎng)計(jì)算。通過觀察，首先在圖像固定處建立4x4的搜索 窗格,遍歷其中的所有元素，依次獲取其灰度值，當(dāng) 有像素點(diǎn)滿足作為種子點(diǎn)的灰度要求，選擇該點(diǎn)作 為區(qū)域生長(zhǎng)種子點(diǎn)；若該搜索窗格中沒有滿足作 為種子點(diǎn)的像素點(diǎn)，則搜索窗格向上平移四個(gè)像素 位置至處，依照第一次搜索窗格方法進(jìn)行種子點(diǎn) 獲取;其次搜索窗格平移至、處，直至搜索窗 格內(nèi)有像素點(diǎn)滿足作為種子點(diǎn)的灰度要求的像素

14、 點(diǎn);若五個(gè)位置搜索窗格都沒有符合要求像素點(diǎn)，則 說明該時(shí)刻人工心臟瓣膜處于閉合狀態(tài)，開口面積 為0,如圖6和圖7所示。區(qū)域生長(zhǎng)準(zhǔn)則以像素值的相近程度即相似性為 判斷標(biāo)準(zhǔn)，可基于區(qū)域的灰度范圍、基于區(qū)域形狀或 性質(zhì)、基于區(qū)域內(nèi)灰度分布統(tǒng)計(jì)。開口陰影與瓣膜 有明顯的灰度差，因此設(shè)置一定的區(qū)域的灰度范圍 為相似標(biāo)準(zhǔn)來確定有效開口陰影區(qū)域生長(zhǎng)準(zhǔn)則。當(dāng)不再有像素滿足加入?yún)^(qū)域生長(zhǎng)準(zhǔn)則時(shí)，區(qū)域 生長(zhǎng)就會(huì)停止。這就是終止生長(zhǎng)規(guī)則的確定，代表 著區(qū)域生長(zhǎng)結(jié)束。具體算法流程如圖8所示，區(qū)域生長(zhǎng)算法如圖圖6像素坐標(biāo)系表示Fig. 6 Pixel coordinate system representation圖7

15、搜索窗格建立Fig. 7 圖6像素坐標(biāo)系表示Fig. 6 Pixel coordinate system representation圖7搜索窗格建立Fig. 7 Search pane creation圖8區(qū)域生長(zhǎng)算法流程Fig. 8 Region growing algorithm flow圖9區(qū)域生長(zhǎng)算法結(jié)果Fig. 9 Results of region growing algorithm2結(jié)果與分析經(jīng)過圖像處理以及區(qū)域生長(zhǎng)算法分割后的人工 心臟瓣膜開口二值化圖像可以方便地求出心臟瓣膜 開口陰影部分的像素點(diǎn)數(shù)。通過提前對(duì)相機(jī)的標(biāo)定 處理，可以求出實(shí)際情況下每一個(gè)像素點(diǎn)對(duì)應(yīng)的面 積項(xiàng)。人

16、工心臟瓣膜的實(shí)際開口面積計(jì)算公式 (1)為：0EOA =01 pixels.( !)其中，0GO-為人工心臟瓣膜開口面積;01為單 個(gè)像素的對(duì)應(yīng)的實(shí)際面積；pixels為心臟瓣膜開口 陰影部分的像素點(diǎn)數(shù)。在一個(gè)心動(dòng)周期0.85 s中，每0.1 s對(duì)人工心臟 瓣膜開口面積進(jìn)行計(jì)算，見表1。心臟瓣膜有效開口面積的參考標(biāo)準(zhǔn)0.7 cm2。進(jìn)一 步通過計(jì)算出的有效開口面積，可以分析與流量、跨 瓣壓差的對(duì)應(yīng)關(guān)系。H結(jié)束語為了能夠精確、完整地分割出心臟瓣膜開口陰 影部分，本文利用區(qū)域生長(zhǎng)算法，準(zhǔn)確計(jì)算出人工心 臟瓣膜開口陰影部分的像素個(gè)數(shù)，實(shí)現(xiàn)了對(duì)有效開 口面積的計(jì)算;通過建立搜索窗格，實(shí)現(xiàn)了種子生長(zhǎng) 點(diǎn)的自動(dòng)設(shè)置。實(shí)驗(yàn)結(jié)果表明，基于自動(dòng)種子點(diǎn)設(shè) 置的區(qū)域生長(zhǎng)算法對(duì)人工心臟瓣膜有效開口面積測(cè) 量有很好的適應(yīng)性,相較于人工分割方法，效率和自 動(dòng)化程度有著顯著的提高。表1不同時(shí)刻瓣膜開口面積Tab. 1