1.一種血管中心線的提取方法,其特征在于,包括如下步驟:
獲取醫(yī)學(xué)圖像;
基于所述醫(yī)學(xué)圖像提取血管,以及若干個連通域;
篩選所述若干個連通域,獲取n個血管連通域;
基于各個血管連通域的方向坐標值降序標記所述n個血管連通域;
根據(jù)所述降序標記,選取第一個血管連通域連接血管,并生長該連通域與血管之間的中心線;
根據(jù)所述降序標記順序依次連接并生長第m+1個和第m個連通域之間的中心線,直至生長完成所述n個血管連通域之間的中心線;
其中,m,n為自然數(shù),m小于或等于n-1。
2.根據(jù)權(quán)利要求1所述的提取方法,其特征在于,所述篩選包括如下步驟:
在所述若干個連通域中識別血管連通域;
計算所述血管連通域中各連通域在x、y和z方向上的平均坐標值;
基于各連通域的y平均坐標值升序標記各連通域;
基于所述升序標記,依次判斷該連通域是否為血管連通域。
3.根據(jù)權(quán)利要求2所述的提取方法,其特征在于,所述識別血管連通域包括:
計算所述若干個連通域的截面積,若該連通域的截面積小于閾值a,則該連通域為血管連通域;否則舍棄該連通域。
4.根據(jù)權(quán)利要求3所述的方法,其特征在于,所述閾值a的選取范圍為大于10。
5.根據(jù)權(quán)利要求2所述的提取方法,其特征在于,所述依次判斷該連通域是否為血管連通域包括:
選取y坐標值最小的連通域為基準連通域,根據(jù)所述升序標記,依次比較各連通域的坐標值與所述基準連通域的坐標值:
若yi-y0小于閾值p,且xi-x0大于閾值q,則該連通域為血管連通域;
否則,舍棄該連通域;
其中,x0和y0為所述基準連通域分別在x、y方向上的平均坐標值,xi和yi為第i個連通域在x、y方向上的平均坐標值,p為自然數(shù),q為整數(shù)。
6.根據(jù)權(quán)利要求5所述的提取方法,其特征在于,所述獲取閾值q的取值范圍為大于10,所述閾值q的取值范圍為大于-15。
7.根據(jù)權(quán)利要求2所述的提取方法,其特征在于,所述篩選還包括,選取z坐標值小于閾值k的連通域進行血管連通域的識別。
8.根據(jù)權(quán)利要求1所述的方法,其特征在于,所述第一個血管連通域連接血管,包括:
選取z坐標值大于所述第一個血管連通域所在圖像的醫(yī)學(xué)圖像,根據(jù)源距離場,選取灰度值為(100,500)范圍內(nèi)的像素點作為種子點生長,連接所述血管連通域與血管。
9.根據(jù)權(quán)利要求1所述的提取方法,其特征在于,生長所述中心線,包括:
以該連通域連接至血管的第一個像素點為起始點,基于最短路徑算法生長所述連通域與血管之間的中心線。
10.根據(jù)權(quán)利要求1~9所述的提取方法,其特征在于,基于連通域的x坐標平均值將各個血管連通域劃分為左、右區(qū)域,如x坐標平均值大于圖像x方向尺寸的一半,則該連通域位于右區(qū)域,否則,該連通域位于左區(qū)域;
對左區(qū)域或右區(qū)域內(nèi)的血管連通域,根據(jù)各個血管連通域的z方向坐標值降序標記所述n個血管連通域;
根據(jù)所述降序標記,選取第一個血管連通域連接血管,并生長該連通域與血管之間的中心線;
根據(jù)所述降序標記順序依次連接并生長第m+1個和第m個連通域之間的中心線,直至生長完成所述n個血管連通域之間的中心線;
根據(jù)所述中心線分別獲取獲取左區(qū)域或右區(qū)域內(nèi)的血管;
其中,m,n為自然數(shù),m小于或等于n-1。