[發明專利]一種頸動脈狹窄程度計算方法、裝置及存儲介質

有效

• 申請號: 201811472202.0
• 申請日: 2018-12-04
• 公開（公告）號: CN109754388B
• 公開（公告）日: 2021-03-16
• 發明（設計）人: 毛順億;周建華
• 申請（專利權）人: 眾安信息技術服務有限公司
• 主分類號: G06T7/00
• 分類號: G06T7/00;G06T7/62
• 代理公司: 北京市萬慧達律師事務所 11111
• 代理人: 張慧娟

1.一種頸動脈狹窄程度計算方法，其特征在于，所述方法包括步驟：

S1、接收頸部MRA影像，并對所述頸部MRA影像中的第一序列和第二序列進行層間配準，以形成層間對應關系，其中，所述第二序列為TOF序列，所述第一序列為不同于所述TOF序列的序列；

S2、從配準好的所述第一序列的每層影像中分別提取出橫斷位血管，并計算各個所述橫斷位血管的血管面積；

S3、從配準好的所述第二序列的每層影像中分別提取出橫斷位血液目標，并計算各個所述橫斷位血液目標的血液面積；

S4、根據所述層間對應關系，將各個所述橫斷位血管的血管面積分別與對應的各個所述橫斷位血液目標的血液面積代入頸動脈狹窄程度計算公式，計算得到多個頸動脈狹窄程度；

S5、將所述多個頸動脈狹窄程度中的最大值輸出作為所述頸部MRA影像的頸動脈狹窄程度；

所述頸動脈狹窄程度計算公式為：

其中，R為頸動脈狹窄程度，ε為血管壁厚度的調控因子，A₁為在第一序列中的一個橫斷位血管的血管面積，A_TOF為與該橫斷位血管相對應的TOF序列中的一個橫斷位血液目標的血液面積。

2.根據權利要求1所述的方法，其特征在于，所述步驟S1進一步包括：

S11、針對所述第一序列中的每層影像，分別從所述第二序列中查找與其在Z軸坐標上的差值為最小的TOF影像；

S12、將查找到的多個所述TOF影像構成配準好的所述第二序列，并形成所述層間對應關系。

3.根據權利要求1所述的方法，其特征在于，所述步驟S2進一步包括：

S21、將配準好的所述第一序列的每層影像分別輸入到預先訓練好的橫斷位動脈血管分割模型進行預測，得到所述第一序列的每層影像的血管圖像概率圖；

S22、將所述每層影像的血管圖像概率圖分別輸入到指示函數以得到多個輸出結果圖，并對所述多個輸出結果圖分別轉化為二值圖，以獲取多個所述橫斷位血管；

S23、提取出各個所述橫斷位血管的血管外邊緣，并將各個所述血管外邊緣對應勾畫在所述第一序列的每層影像上；

S24、分別對所述第一序列的每層影像上的所述血管外邊緣內的像素點數進行統計計算，得到各個所述橫斷位血管的血管面積。

4.根據權利要求3所述的方法，其特征在于，所述步驟S23進一步還包括：

對多個所述橫斷位血管分別進行像素掩膜以形成多個掩膜圖像；

將多個所述掩膜圖像分別疊加在對應的所述第一序列的每層影像上，以用于調整所述每層影像上的橫斷位血管的血管外邊緣。

5.根據權利要求3所述的方法，其特征在于，所述橫斷位動脈血管分割模型的訓練過程如下：

a、構建橫斷位血管數據集及對應的標注數據集；

b、構建橫斷位動脈血管分割網絡；

c、使用所述橫斷位血管數據集和所述標注數據集對所述橫斷位動脈血管分割網絡進行訓練，生成所述橫斷位動脈血管分割模型。

6.根據權利要求1所述的方法，其特征在于，所述步驟S3進一步包括：

S31、采用閾值法從配準好的所述第二序列的每層影像中分別提取出初步的橫斷位血液目標；

S32、對所述每層影像中的所述初步的橫斷位血液目標中錯誤橫斷位血液目標進行濾除，得到多個所述橫斷位血液目標；

S33、提取出各個所述橫斷位血液目標的血液外邊緣，并將各個所述血液外邊緣對應勾畫在所述第二序列的每層影像上；

S34、對所述第二序列的每層影像上的所述血液外邊緣內的像素點數進行統計計算，得到各個所述橫斷位血液目標的血液面積。

7.根據權利要求1所述的方法，其特征在于，所述步驟S4進一步包括：

S41、基于所述層間對應關系，將各個所述橫斷位血管分別與各個所述橫斷位血液目標進行位置對應；

S42、將各個所述橫斷位血管的血管面積分別與位置對應的各個所述橫斷位血液目標的血液面積代入頸動脈狹窄程度計算公式，計算得到多個頸動脈狹窄程度。