[發(fā)明專利]一種CT圖像中脊髓的提取方法

有效

• 申請?zhí)枺?/td>: 201310332806.6
• 申請日: 2013-08-01
• 公開（公告）號: CN104346799B
• 公開（公告）日: 2018-02-02
• 發(fā)明（設計）人: 劉靖;李強
• 申請（專利權）人: 上海聯(lián)影醫(yī)療科技有限公司
• 主分類號: G06T7/143
• 分類號: G06T7/143;G06K9/46;A61B6/03
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種CT圖像中脊髓的提取方法，其特征在于，包括如下步驟：

(1)輸入三維CT圖像，通過二維簡化的椎管模型，得到所述三維CT圖像中所述椎管的中軸線位置的三維廣義霍夫變換空間；

(2)根據(jù)所述三維廣義霍夫變換空間，確定所述椎管中軸線在所述三維CT圖像中的位置；

(3)根據(jù)所述椎管中軸線的位置，通過區(qū)域增長方法，得到所述椎管的分割掩模，即得到所述三維CT圖像中的脊髓；

所述椎管的中軸線位置的三維廣義霍夫變換空間通過如下過程得到：

通過邊緣檢測方法，得到第k幀軸切面圖像中的邊緣點位置和所述邊緣點位置處的正規(guī)化梯度向量其中，k為所述三維CT圖像中第k幀二維軸切面圖像，j＝1,...,N^k為所述第k幀軸切面圖像中的第j個邊緣點，N^k為所述得到的第k幀軸切面圖像中的所有邊緣點數(shù)；

根據(jù)所述邊緣點位置和所述正規(guī)化梯度向量得到所述邊緣點映射到廣義霍夫變換空間中的位置和并進行投票，得到所述廣義霍夫空間中的點分布函數(shù)和r₁為所述二維簡化模型中椎孔的半徑，r₂為所述二維簡化模型中椎體的半徑，c為所述廣義霍夫空間中的任一位置；

將所述點分布函數(shù)和進行合并，得到所述第k幀軸切面的廣義霍夫變換空間H^k(c)；

將所有軸切面的廣義霍夫變換空間在軸向進行級聯(lián)，得到H(c,k)；

將所述H(c,k)進行平滑，得到所述椎管的中軸線的三維廣義霍夫變換空間H(c,k)。

2.如權利要求1所述一種CT圖像中脊髓的提取方法，其特征在于，所述點分布函數(shù)的公式為：

H1k(c)=Σr1Σjv1jk(c,r1)r1]]>

其中，為投票函數(shù)，I[·]為指標函數(shù)；所述r₁為所述二維簡化模型中椎孔的半徑。

3.如權利要求1所述一種CT圖像中脊髓的提取方法，其特征在于，所述點分布函數(shù)的公式為：

H2k(c)=Σr2Σj2v2jk(c,r2)r2]]>

其中，為投票函數(shù)，I[·]為指標函數(shù)，為所述正規(guī)化梯度向量沿著人體冠狀軸的分量；所述r₂為所述二維簡化模型中椎體的半徑；n為所述二維簡化模型中椎孔中軸線與軸切面的交點以及椎體中軸線與軸切面的交點之間的相對位移。