1.一種血管提取方法,其特征在于,包括如下步驟:
獲取由若干層切片圖像構(gòu)成的三維圖像,選取一層切片圖像,基于所述切片圖像中像素點的灰度信息獲取若干個連通域;
基于所述連通域的圓度確定種子點;
以所述種子點為起點,基于閾值a對所述三維圖像進行區(qū)域生長,獲取第一圖像;
以所述種子點為起點,基于第一閾值對所述第一圖像進行區(qū)域生長,獲取血管掩膜;
將所述第一圖像與血管掩膜相減獲取第一減影圖像;
粗分割所述第一減影圖像,獲取初始掩膜;
形態(tài)學(xué)處理所述初始掩膜,獲取骨骼掩膜;
將所述第一圖像與所述骨骼掩膜相減獲取第二減影圖像,所述第二減影圖像包括若干個連通域;
基于血管生長條件連接所述若干個連通域,獲取血管提取結(jié)果。
2.根據(jù)權(quán)利要求1所述的方法,其特征在于,所述選取一層切片圖像為沿Z軸方向上的終止層切片圖像。
3.根據(jù)權(quán)利要求1所述的方法,其特征在于,所述確定種子點包括:
計算所述連通域中所有像素點的邊界距離場值;
選取最大邊界距離場值為該連通域的半徑,根據(jù)圓面積公式計算該半徑對應(yīng)的圓形面積,并計算所述圓形面積與所述連通域?qū)嶋H面積的比值,所述比值為所述連通域的圓度;
選取所述切片圖像中圓度與1差值的絕對值的最小值所在的連通域,確定該連通域中最大邊界距離場值所在的像素點為所述種子點。
4.根據(jù)權(quán)利要求1所述的方法,其特征在于,所述獲取血管掩膜包括:
計算所述第一圖像中各體素點的邊界距離場值;
根據(jù)所述種子點基于所述第一閾值進行三維區(qū)域生長,獲取所述血管掩膜;
其中,所述第一閾值為大于3.0個體素點的寬度。
5.根據(jù)權(quán)利要求4所述的方法,其特征在于,獲取所述血管掩膜還包括膨脹操作。
6.根據(jù)權(quán)利要求1所述的方法,其特征在于,所述粗分割包括:
計算所述第一減影圖像上所有體素點的邊界距離場值;
選取邊界距離場值大于閾值b的體素點為骨骼種子點,根據(jù)所述第一減影圖像中體素點的灰度信息或邊界距離場值進行區(qū)域生長,獲取所述初始掩膜;
其中,所述閾值b為大于3.0個體素點的寬度;所述體素點的灰度信息為灰度值大于500,所述體素點的邊界距離場值為大于1.5個體素點的寬度。
7.根據(jù)權(quán)利要求1所述的方法,其特征在于,形態(tài)學(xué)處理包括:
以閾值c為寬度腐蝕所述初始掩膜,獲取若干個連通域;
選取最大體積的連通域為骨骼區(qū)域;
以閾值d為寬度膨脹所述骨骼區(qū)域,獲取第一骨骼掩膜。
8.根據(jù)權(quán)利要求7所述的方法,其特征在于,所述獲取骨骼掩膜還包括填實操作,包括如下步驟:
根據(jù)所述種子點基于閾值e對所述三維圖像進行區(qū)域生長,獲取第二圖像;
逐層填實所述第二圖像中各連通域;
將填實后的第二圖像與所述第一骨骼掩膜進行圖像疊加,獲取疊加圖像。
9.根據(jù)權(quán)利要求8所述的方法,其特征在于,所述填實操作還包括:
提取所述疊加圖像中最大體積的連通域,對該連通域經(jīng)閉操作獲取所述骨骼掩膜。
10.根據(jù)權(quán)利要求1所述的方法,其特征在于,所述血管生長條件為:
計算所述第二減影圖像中若干個連通域所有體素點的邊界距離場值;
基于最短路徑算法,并且排除所述骨骼掩膜中各個連通域所在區(qū)域提取血管中心線,獲取所述血管提取結(jié)果。