l 獲取CT斷層圖像數(shù)據(jù)

　　葉輪斷層圖像的獲取是在江蘇大學(xué)附屬江濱醫(yī)院的多層螺旋CT機(jī)上完成的。該CT機(jī)為美國(guó)GE公司產(chǎn)品;與普通CT機(jī)相比，多層螺旋CT機(jī)具有高清晰度、高分辨率、高速掃描、性能優(yōu)越等特點(diǎn)，而且旋轉(zhuǎn)一圈就能采集出8層圖像，實(shí)現(xiàn)了真正的高速透視成像。

　　根據(jù)葉輪的結(jié)構(gòu)特點(diǎn)，為了重建時(shí)能充分利用葉輪的結(jié)構(gòu)特征，同時(shí)也為了減少三維圖像重建的難度與工作量，采集數(shù)據(jù)時(shí)盡可能使葉輪中心線與CT機(jī)水平軸心線重合(圖1)。采集精度采用512×512 x1.25，掃描后共獲得56幅斷層圖像，如圖2所示為第20、25幅圖像。從圖2中可以看出，經(jīng)多層螺旋CT透視掃描后，葉輪的斷層圖像包含了水泵葉輪的所有內(nèi)部結(jié)構(gòu)，如葉片形狀、輪轂、軸孔、鍵槽等的結(jié)構(gòu)信息，為后續(xù)的數(shù)字圖像處理與三維圖像重建提供了可靠的數(shù)據(jù)與信息基礎(chǔ)。

　　圖1 葉輪在CT中的放置圖

　　圖2 葉輪原始CT圖片

　　由于在CT斷層圖像獲取過程中，CT機(jī)硬件系統(tǒng)、外界干擾以及傳輸過程中的噪聲等都會(huì)對(duì)圖像的質(zhì)量產(chǎn)生影響，使圖像出現(xiàn)不同程度的質(zhì)量下降，必須進(jìn)行濾波去噪;另一方面，為了消除因外界干擾以及斷層圖像重建過程中各種客觀因素的影響可能出現(xiàn)的偽跡現(xiàn)象，也很有必要通過灰度變換等措施使斷層圖像得到增強(qiáng)。

　　2.1 濾波去噪

　　一般在空間域內(nèi)采用鄰域平均即局部平均濾波來減少噪聲，而在頻率域內(nèi)，因噪聲頻譜多在高頻段，故一般采用中值濾波和高斯濾波等低通濾波的辦法。

　　中值濾波是一種非線性的信號(hào)處理方法，可以克服線性濾波器如均值濾波等帶來的圖像細(xì)節(jié)模糊，對(duì)濾除脈沖干擾及圖像掃描噪聲很有效。根據(jù)課題特點(diǎn)，采用中值濾波。

　　使用3×3模板分別對(duì)斷層數(shù)據(jù)第20、25斷層進(jìn)行中值濾波，濾波后的圖像如圖3所示。可以看出，噪聲己基本消除，同時(shí)邊界也很好地得到保護(hù)。

　　圖3 中值濾波后CT圖片

　　2.2 灰度變換

　　灰度變換的目的是增強(qiáng)CT斷層圖像的對(duì)比度，擴(kuò)大圖像灰度范圍，從而達(dá)到增強(qiáng)圖像的目的。為此采用局部分段線性變化法來實(shí)現(xiàn)灰度變化，該算法的數(shù)學(xué)模型為：

　　圖4是處理后的第20、25斷層圖像。需要注意的是，在處理時(shí)還可以通過小心調(diào)整折線拐點(diǎn)的位置并控制分段直線的斜率來動(dòng)態(tài)地對(duì)任一灰度區(qū)間進(jìn)行擴(kuò)展或壓縮，直至達(dá)到滿意的效果。

　　圖4 灰度變換后得到增強(qiáng)的CT圖片

　　3 二維圖像處理

　　二維圖像處理主要是對(duì)前面處理過的圖像進(jìn)行邊緣檢測(cè)并跟蹤、提取其輪廓，為二維幾何處理提供數(shù)據(jù)準(zhǔn)備。邊緣檢測(cè)的方法很多，而且新算法還在不斷地出現(xiàn)。根據(jù)需要采用基于梯度的檢測(cè)算子來進(jìn)行邊緣檢測(cè)，用八鄰域輪廓跟蹤來實(shí)現(xiàn)輪廓提取。

　　3.1 邊緣檢測(cè)

　　基于梯度的檢測(cè)算子有很多，本文采用Robert算子對(duì)灰度變換后的第20、25斷層圖像進(jìn)行邊緣檢測(cè)。Robert邊緣檢測(cè)算子是根據(jù)任意一對(duì)互相垂直方向上的差分可用來計(jì)算梯度的原理，采用對(duì)角線方向相鄰像素之差。該算子的數(shù)學(xué)模型為：

　　圖5是處理后的圖像，可以看出：葉片的邊緣已比較完整、清晰地被檢測(cè)出來，說明使用Robert邊緣檢測(cè)算子的定位精度高，效果較好。

　　圖5 邊緣檢測(cè)后的CT圖片

　　3.2 輪廓跟蹤與提取

　　邊緣檢測(cè)后得到的邊緣輪廓是離散、獨(dú)立、互不相關(guān)的各邊界輪廓點(diǎn)的無序組合，需要利用輪廓跟蹤算法，沿著圖像的輪廓進(jìn)行搜索，將搜索到的輪廓上點(diǎn)的坐標(biāo)記錄在點(diǎn)列中存儲(chǔ)，供以后使用，其結(jié)果是，每一條輪廓線都表示成一個(gè)點(diǎn)列。

　　本文采用八鄰域輪廓跟蹤方法，在采用輪廓跟蹤算法時(shí)，邊緣灰度圖中背景像素顏色為白色，邊緣上的像素為黑色。

　　結(jié)果如圖6所示。可以看出，葉輪斷面的輪廓是封閉的，同一切片上的多輪廓也互不相交，實(shí)現(xiàn)輪廓跟蹤后的每一條輪廓線是單像素寬的互不相交的曲線。同時(shí)也得到了每個(gè)輪廓的輪廓點(diǎn)序列，為矢量化處理打下了良好的基礎(chǔ)。

　　圖6 提取輪廓后的CT圖片

　　4 二維幾何處理

　　經(jīng)上述處理后，得到的是葉輪斷層上每一個(gè)輪廓的離散點(diǎn)序列，由于存在大量的冗余信息，不能直觀地描述斷層上輪廓的幾何拓?fù)浣Y(jié)構(gòu)，在商用CAD或自編軟件中進(jìn)行編輯處理時(shí)的工作量很大，必須進(jìn)一步用直線、圓弧來擬合這些用離散點(diǎn)序列表達(dá)的輪廓，實(shí)現(xiàn)用幾何的方式來描述這些輪廓，即對(duì)輪廓進(jìn)行矢量化處理，通過坐標(biāo)變換將以像素為單位的矢量數(shù)據(jù)變換為笛卡爾坐標(biāo)系等以國(guó)際長(zhǎng)度單位為單位的矢量數(shù)據(jù)。

　　采用直線段和圓弧對(duì)提取的輪廓進(jìn)行擬合實(shí)現(xiàn)了二維輪廓的矢量化。基本步驟是：對(duì)跟蹤得到的離散輪廓點(diǎn)序列，先提取輪廓上的特征點(diǎn)，利用特征點(diǎn)作為分界點(diǎn)將平面輪廓?jiǎng)澐譃槿舾煞侄危⑴袛喔骶€段的曲線類型，再應(yīng)用分段一合并算法實(shí)現(xiàn)輪廓的直線段擬合，最后應(yīng)用垂直平分線求交法實(shí)現(xiàn)輪廓的直線圓弧擬合。

　　坐標(biāo)變換的變換公式為：

　　式中(Xi，yi)為圖像像素位置，(a，b)為笛卡爾坐標(biāo)的原點(diǎn)相對(duì)于搜索起點(diǎn)的坐標(biāo)，(Xi，Yi)為圖像像素位置在設(shè)定的笛卡爾坐標(biāo)系下的坐標(biāo)，坐標(biāo)的單位均為像素;S為系統(tǒng)標(biāo)定參數(shù)，其含義為圖像中一個(gè)像素長(zhǎng)度代表的實(shí)際長(zhǎng)度為Smm。

　　5 三維模型重建

　　三維重建主要有3種途徑：基于輪廓的表面重建、基于體素的表面重建以及基于體素的直接體視法。前兩種重建方法都屬于表面重建，其中基于輪廓表面重建法是從一組平面輪廓重構(gòu)通過這些輪廓的表面;基于體素的表面重建法首先要確定物體表面在每個(gè)小體素內(nèi)的小面片，然后將這些小面片連接起來構(gòu)成物體的表面。這兩種算法可采用現(xiàn)有的圖形學(xué)顯示方法。直接體視法直接對(duì)體數(shù)據(jù)進(jìn)行顯示，而不是構(gòu)造物體表面的幾何描述，可顯示體數(shù)據(jù)中所包含的所有物體，但由于不包含物體的幾何拓?fù)浣Y(jié)構(gòu)關(guān)系，因而不能很好地表現(xiàn)空間層次。本文先用第一種重建算法實(shí)現(xiàn)了水泵葉輪的三維重建，然后再通過一系列算法將該模型轉(zhuǎn)換成實(shí)體模型。

　　5.1 表面三維重建

　　首先在前面已提取并矢量化的輪廓線上尋找特征點(diǎn)，試圖通過這些點(diǎn)的連線來近似地表達(dá)輪廓線

　　圖7 輪廓點(diǎn)所圍面積計(jì)算

　　由于特征點(diǎn)數(shù)目遠(yuǎn)小于矢量化的輪廓線上的點(diǎn)數(shù)，故構(gòu)造三角面片的工作量大為減小。本文采用了一種基于Delaunay三角剖分思想的表面重建的算法，算法考慮了組成最佳重建表面的三角片的形態(tài)特點(diǎn)，按照剖分三角片集中的最小內(nèi)角在整體上盡可能大的優(yōu)化準(zhǔn)則進(jìn)行剖分，算法簡(jiǎn)單，重建表面質(zhì)量較好。

　　5.2 表面模型轉(zhuǎn)換為實(shí)體模型

　　由于實(shí)體模型的任意曲匾可用多邊形網(wǎng)格來逼近，而邊界模型(B-Rep模型)雖然是用實(shí)體的表面來表示實(shí)體的形狀，但在具體處理時(shí)是將實(shí)體的表面形狀表示為多面體，即用一組三維多邊形所圍成的形體來表示實(shí)體。所以可以在前面已經(jīng)構(gòu)造出的由三角面片形成的表面模型的基礎(chǔ)上用邊界模型來構(gòu)造和表達(dá)重建的三維實(shí)體模型。重建后的葉輪三維模型如圖8所示。可看出，該模型在基本形狀、結(jié)構(gòu)葉輪三維模型尺寸與葉輪原型基本一樣，但是在光順程度上仍有欠缺，需要在今后對(duì)算法與重建過程繼續(xù)加以改進(jìn)。

　　圖8 重建后的水泵葉輪三維模型

　　6 結(jié)束語

　　利用多層螺旋CT對(duì)水泵葉輪進(jìn)行透視掃描，采集得到了含有葉輪設(shè)計(jì)及內(nèi)部構(gòu)造信息的CT斷層圖像;在VC++6.0開發(fā)平臺(tái)上，對(duì)這些斷層圖像進(jìn)行一系列數(shù)字圖像處理，實(shí)現(xiàn)了水泵葉輪的三維表面模型重建以及實(shí)體模型的轉(zhuǎn)換，為水泵葉輪的逆向工程提供了一種新的思路。