自動消除DSA圖像中運(yùn)動偽影的方法

1、    自動消除DSA圖像中運(yùn)動偽影的方法        李 雯，聶生東 時間:2008年07月01日     字 體: 大 中 小        關(guān)鍵詞:        摘要：關(guān)鍵詞： 數(shù)字減影血管造影 運(yùn)動偽影 配準(zhǔn)通常，造影劑注入前的DSA血管圖像被稱為蒙片(Mask Im

2、age)，而同一部位在注入造影劑后的X射線圖像被稱為活片(Live Image)。在假設(shè)血管周圍組織、成像條件完全不變的情況下，將蒙片與活片相減應(yīng)得到僅包含血管的圖像。但在臨床應(yīng)用中，這種理論假設(shè)總是不成立的。在成像的同時，病人會有各種運(yùn)動，其中有些運(yùn)動是不可避免的，如呼吸、肌肉運(yùn)動、心臟運(yùn)動以及在造影劑注入時由于病人的自然反應(yīng)而引起的局部運(yùn)動等。這些運(yùn)動會使活片與蒙片之間存在并非由造影劑造成的差異，從而使減影后的圖像出現(xiàn)偽影，妨礙醫(yī)生的正確診斷。本文提出了一種自動消除數(shù)字減影血管造影DSA(Digital Subtraction Angiography)減影圖像中偽影的方法。通過實驗證實，本

3、方法比其他類似的方法1更快、更有效。1 方法原理DSA序列圖像的配準(zhǔn)目的是通過圖像處理和圖像分析手段找到一種幾何變換方法能消除序列圖像中任意二幅之間因為病人的運(yùn)動帶來的偽影。本文選擇對蒙片作相應(yīng)的變形(warping)，使其在圖像相減之前與造影圖像相匹配。DSA圖像的實質(zhì)就是X射線投影圖像。根據(jù)不同物質(zhì)對X射線的吸收系數(shù)的不同，在三維空間把沿著射線方向的不同物體的吸收系數(shù)疊加，形成二維投影圖像。因此，DSA圖像中的運(yùn)動偽影是由三維空間運(yùn)動造成的。2 實現(xiàn)過程2.1 找控制點2，具有很好的邊緣檢測性能，正得到越來越廣泛的應(yīng)用。Canny邊緣檢測法利用高斯函數(shù)的一階微分，能在噪聲抑制和邊緣檢測之間

4、取得較好的平衡。由于沒必要把所有用Canny算子檢測出的點都作為控制點，所以在Canny邊緣點計算出該點的梯度值后，根據(jù)閾值來判斷該點是否作為控制點。為使控制點的分布更合理，利用比例因子對控制點間的距離作規(guī)定：每2個控制點間的最小距離為Dmin=minmax=maxM。通常在DSA圖像中，有用的信息都顯示在一個圓形區(qū)域內(nèi)，區(qū)域之外的像素值為零。這樣的圓形區(qū)域被稱為“曝光區(qū)域”，其定義如下： R=(x-xc)2+(y-yc)2r2，XminxXmax，YminyYmax其中，(xc，yc)為影像增強(qiáng)器的中心點，通過適當(dāng)校正可與圖像的中心點重合。r為圓形區(qū)域的半徑，而Xmin、Xmax、Ymin、

5、Ymax分別是圖像的左、右和上、下邊界。既然所有有用信息都在圓形區(qū)域內(nèi)，因此控制點也應(yīng)該在R之內(nèi)尋找。圖1為DSA蒙片及其控制點圖像，其中圖1(a)為一腦部DSA正面圖像的蒙片，圖1(b)是從圖1(a)中得到的控制點圖像。由圖可見，控制點都位于邊緣上，且保持著一定的均勻度。2.2 選擇相似性測度本文采用模版匹配法求蒙片與活片中像素的相對位移：=(dx，dy)。模版匹配法利用相似性測度來求一幅圖像I0(x，y，t0)中的某一點P(x，y)相對于另一幅圖像I1(x，y，t1)(t1t0)的位移，以P點為中心，以I0中P點周圍W×W鄰域內(nèi)的點為一窗口，在I1中求出這一鄰域移動后的窗口位置，

6、并利用適當(dāng)?shù)南嗨菩詼y度，求出使相似性測度達(dá)到最大的位移，即為要求的最佳位移。本文選用基于圖像直方圖的相似性測度3，當(dāng)兩幅圖達(dá)到最佳匹配時，兩幅圖像之差的直方圖會集中在很窄的一段灰度范圍內(nèi)，且形成明顯的峰值，而當(dāng)兩幅圖像不是很匹配時，直方圖的集中程度會降低；更重要的是，直方圖的這種性質(zhì)不會受血管的影響，DSA造影圖像及其在不同配準(zhǔn)情況下的直方圖比較如圖2所示。其中圖2(a)是DSA造影圖像，其中的白框表示被選取的區(qū)域，左上方的白框不包括血管，而中間的白框包括血管；圖2(b)和圖2(c)分別表示不包括血管時，最佳和非最佳配準(zhǔn)下差異圖像的直方圖；圖2(d)和圖2(e)分別表示包括血管時，最佳和非最佳

7、配準(zhǔn)下差異圖像的直方圖。假設(shè)待配準(zhǔn)的兩幅圖像分別為I0(x，y，t0)和I1(x，y，t1)(t1t0)，求I0中某一點P0的位移。在P0周圍選取一個小窗口W0，該窗以P0為中心，大小為W×W個像素，對于每一個位移在I1中找到相應(yīng)的相同大小的窗口W1，將這兩個窗口相減，得到差異圖像Wd(W×W)，針對Wd計算相似性測度，有研究人員提出基于直方圖的判定標(biāo)準(zhǔn)3：其中min，max?奐Z是差異圖像中可能的灰度范圍，h()是差異圖像的歸一化直方圖，并假設(shè)f是嚴(yán)格凸函數(shù)或嚴(yán)格凹函數(shù)，但可微。相關(guān)文獻(xiàn)提出了多種基于直方圖的f函數(shù)形式4，并指出當(dāng)f(x)=x2時計算簡單且精度高。這樣，式

8、(1)應(yīng)寫成：2.3 亞像素精度某些研究者指出，即使是亞像素級的移動也會造成兩幅圖像間的差異很大5，因此在計算像素位移時應(yīng)設(shè)法達(dá)到亞像素精度，而這要用到圖像插值技術(shù)。本方法通過先對圖像插值和重采樣，然后對插值后的圖像求相似性測度來達(dá)到亞像素精度。期望達(dá)到的亞像素精度為0.1 pixel。圖像插值在整個配準(zhǔn)算法中很重要。它的精度直接影響著位移的精度；它的計算速度影響著整個算法的計算時間。如何選擇一種既精確又快速的插值方法成為此配準(zhǔn)算法成敗的關(guān)鍵。本文采用三次多項式插值法對原圖像I0進(jìn)行插值。為得到理想的效果，本文采用4×4的卷積核。2.4 最優(yōu)化搜索在確定了相似性測度為直方圖能量E()

9、后，要求最佳位移使相似性測度達(dá)到最大。本問題屬于無約束優(yōu)化問題，針對該問題，本文采用直接搜索法。最常用的是鮑威爾(Powell)法6和單純形法。這里采用鮑威爾法，也稱為方向加速法。搜索的速度受計算位移時選用的窗W×W大小影響顯著，這是個矛盾的統(tǒng)一體：一方面為了減少搜索時間，希望采用較小的窗；另一方面，較小的窗將使相似性測度產(chǎn)生許多局部極大點，導(dǎo)致搜索結(jié)果不可靠。因此在選擇窗的大小時應(yīng)該慎重考慮，本文選用的窗口大小為51×51個像素。2.5 位移插值求出所有控制點的最佳位移后，進(jìn)而需要求出圖像I0中所有點相對I1的位移，這要利用線性插值來估算出所有點的位移。本文采用Delau

10、nay三角插值法7，它保證在每一個三角形內(nèi)的最小角度盡可能地最大化，且保證了這種三角形分法的惟一性。獲得圖像I0中所有點的位移后，即可對圖像I0進(jìn)行校正。因為位移的單位是0.1像素，則圖像上一點根據(jù)移動后極可能落入整數(shù)值的網(wǎng)格之間，這里需要用插值法(雙線性插值)得到非整數(shù)像素位置上的像素值。3 實驗結(jié)果及討論拍攝后前位的DSA時間序列圖像共39幅，每幅圖像的大小為512×512個像素，選取其中的第1幅(蒙片)和第28幅(活片)作DSA圖像的配準(zhǔn)研究，其中的蒙片如圖1(a)?？刂泣c提取之后(見圖1(b)，根據(jù)模板匹配法計算每個控制點的最佳位移，亦即使差異圖像的直方圖能量達(dá)到最大時的位移。此時，對窗口大小W的選擇很重要。若采用鮑威爾最優(yōu)化方法搜尋最佳位移，則要求E(d)足夠光滑以至于只有一個極值。窗口太小將增加E(d)的局部極小值，從而影響到最優(yōu)解的精度。經(jīng)多次嘗試和比較，選用的窗口大小為51×51個像素。校正結(jié)果如圖3所示，其中圖3(a)是校正前蒙片與活片直接減影后的圖像，而圖3(b)是經(jīng)過校正后再做減影的圖像。從結(jié)果看，使用本方法能去除減影圖像中