[發明專利]X射線圖像分割方法及終端

有效

• 申請號: 202010040502.2
• 申請日: 2020-01-15
• 公開（公告）號: CN111311605B
• 公開（公告）日: 2023-07-04
• 發明（設計）人: 藍重洲;劉德建;李華勇;葉超;成富平
• 申請（專利權）人: 深圳市安健科技股份有限公司
• 主分類號: G06T7/11
• 分類號: G06T7/11;G06T5/00
• 代理公司: 深圳市博銳專利事務所 44275
• 代理人: 林棟

1.一種X射線圖像分割方法，其特征在于，包括：

獲取X射線空氣圖；

計算所述X射線空氣圖的噪聲層級比例系數；

獲取X射線臨床圖像；

根據所述X射線臨床圖像和噪聲層級比例系數計算得到CNR圖像；

對所述CNR圖像進行二值化處理，得到CNR圖像掩模；

根據所述CNR圖像計算得到邊界抑制掩模；

根據所述X射線臨床圖像分別計算得到空氣區域掩模和束光器邊緣抑制掩模；

對所述CNR圖像掩模、邊界抑制掩模、空氣區域掩模和束光器邊緣抑制掩模作與計算，得到組織區域掩模；

所述計算所述X射線空氣圖的噪聲層級比例系數具體包括：

對所述X射線空氣圖進行AT變換；

采用圖像金字塔分解AT變換后的X射線空氣圖，得到低頻信號層G₀、G₃，以及細節信號層D₀、D₃；

在第一預設閾值范圍內獲取X射線空氣圖的G₃層中的最大灰度值maxV和最小灰度值minV；

根據公式計算得到G₀層的掩模，其中，T表示第一預設比例；

在G₀層的掩模范圍內，根據公式計算得到標準差圖像S0，其中m表示窗口大小；

計算標準差圖像S₀直方圖的最大值，得到X射線空氣圖的噪聲水平nl₀；

根據公式計算得到標準差圖像S₃；

計算標準差圖像S₃直方圖的最大值，得到X射線空氣圖的噪聲水平nl₃；

根據所述X射線空氣圖的噪聲水平nl₀和噪聲水平nl₃計算得到噪聲層級比例系數；

所述根據所述X射線臨床圖像和噪聲層級比例系數計算得到CNR圖像具體包括：

對所述X射線臨床圖像進行AT變換；

采用圖像金字塔分解AT變換后的X射線臨床圖像，得到低頻信號層G'₀、G'₃以及細節信號層D'₀和D'₃；

在第二預設閾值范圍內獲取X射線臨床圖像的G₃'層中的最大灰度值maxV'和最小灰度值minV'；

根據公式計算得到G'₀層的掩模，其中，T'表示第二預設比例；

在G'₀層的掩模范圍內，根據公式計算得到標準差圖像S'₀，其中m'表示窗口大小；

計算標準差圖像S'₀直方圖的最大值，得到X射線臨床圖像的噪聲水平nl'₀；

根據所述噪聲層級比例系數和X射線臨床圖像的噪聲水平nl'₀計算得到X射線臨床圖像第3層的噪聲水平nl'₃；

根據公式計算得到圖像標準差S'₃；

將圖像標準差S'₃除以噪聲水平nl'₃計算得到CNR圖像；

所述對所述CNR圖像進行二值化處理，得到CNR圖像掩模具體為：

在[a,b]的范圍內，根據公式計算得到CNR圖像掩模mask₀，其中，a、b均為常數；

所述根據所述CNR圖像計算得到邊界抑制掩模具體為：

獲取CNR圖像的寬度w和高度h；

根據公式計算得到邊界抑制掩模mask₁，其中，bor＝floor(w*c)，c為常數，bor為邊界抑制掩膜mask₁的寬度，floor為向下取整函數；

根據所述X射線臨床圖像分別計算得到空氣區域掩模和束光器邊緣抑制掩模具體為：

根據公式airThr＝minV'+airRatio·(maxV'-minV')計算得到空氣區域閾值airThr，其中，airRatio為常數；

根據公式計算得到空氣區域掩模mask₂；

根據公式colThr＝colRatio·(maxV'-minV')計算得到束光器區域閾值colThr，其中colRatio為常數；

計算X射線臨床圖像G'₃層中任一點周圍像素的灰度差異；

若所述灰度差異大于colThr，則將所述任一點的灰度值設為0；若所述灰度差異小于colThr，則將所述任一點的灰度值設為1，得到束光器邊緣抑制掩模mask₃。