[發明專利]基于邊界距離場的血管中心路徑提取方法

無效

• 申請號: 200910022500.4
• 申請日: 2009-05-13
• 公開（公告）號: CN101551862A
• 公開（公告）日: 2009-10-07
• 發明（設計）人: 高新波;呂新榮;鄒華;曹艷艷;梁承志;王向華;查林;李晉舟;寧文博
• 申請（專利權）人: 西安電子科技大學
• 主分類號: G06K9/46
• 分類號: G06K9/46;G06T7/00
• 代理公司: 陜西電子工業專利中心
• 代理人: 王品華;黎漢華

1.一種基于邊界距離場的血管中心路徑提取方法，包括如下過程：

A.對輸入的CTA血管造影體數據進行八叉樹分解，得到每一分塊的信息；

B.根據每個分塊的最大值和最小值與雙閾值之間的關系，遍歷存在血管體素的葉子節點提取血管組織；

C.根據提取的血管組織中所有點的坐標，計算提取后的血管組織邊界距離場：

C1)根據提取的血管組織中所有點的坐標，利用原CTA血管造影體數據中血管體素的梯度值倒數和拉普拉斯變換值之和作為血管體素的邊界距離初始值：

f_start(p)＝1/f_G+f_L

其中，p為提取的血管組織中任意一點的坐標，f_G為原CTA血管數據中點p的梯度值，f_L為原CTA血管數據中點p的拉普拉斯變換值；

C2)根據邊界距離初始值用如下公式計算點p的邊界距離值f_final(p)：

f_final(p)＝min(min(f_start(p₁)+10)，min(f_start(p₂)+14)，min(f_start(p₃)+17))其中，f_start(p₁)為點p的最近鄰鄰域點p₁的邊界距離初始值，f_start(p₂)為點p的面對角鄰域點p₂的邊界距離初始值，f_start(p₃)為點p的體對角鄰域點p₃的邊界距離初始值；

C3)利用八叉樹遍歷所有包含血管體素的葉子節點，完成提取后的血管組織邊界距離場計算：

D.在邊界距離場的基礎上建立最大代價生成樹；

E.提取最大代價生成樹的主干作為血管組織的初始中心路徑；

F.利用如下重心法修正提取的血管組織初始中心路徑，使路徑上的點都處于血管管腔橫截面的重心位置，得到一條符合人類觀察習慣的最終血管中心路徑：

_F1)對于初始中心路徑上的每一點，取與該點法向量垂直，并且與該點相距為10個體素的點組成的圓形二維圖像Ω；

F2)利用如下重心法公式求得Ω的重心C：

C=[xc,yc]=(Σxw(x,y)Σw(x,y),Σyw(x,y)Σw(x,y))]]>

其中，w(x，y)＝α(f_2D(x，y)-m)，α＝-1，m＝max(f_2D(x，y))，f_2D(x，y)為Ω中坐標為(x，y)處的像素值；

F3)連接所有重心點作為最終的血管中心路徑。