[發明專利]一種醫學影像感興趣區域的窗位窗寬計算方法

有效

• 申請號: 202210526046.1
• 申請日: 2022-05-16
• 公開（公告）號: CN114663653B
• 公開（公告）日: 2022-08-12
• 發明（設計）人: 劉孝波
• 申請（專利權）人: 南京佗道醫療科技有限公司
• 主分類號: G06V10/25
• 分類號: G06V10/25;G06T5/40;G06T7/11
• 代理公司: 南京瑞弘專利商標事務所(普通合伙) 32249
• 代理人: 梁天彥

本發明公開了一種醫學影像感興趣區域的窗位窗寬計算方法，包括：將醫學影像數據的體素灰度值轉換為Hu值；構建Hu值的頻率直方圖，并根據設定閾值比例剔除頻率直方圖中的無效組；根據剔除無效組后的頻率直方圖計算得到感興趣區域的Hu值范圍，并據此計算得到窗位窗寬。本發明針對不同影像設備采集的醫學影像或經過其他處理后的醫學影像，可以分割計算得到病灶部分的窗位窗寬，并映射至顯示屏上顯示出來。

技術領域

本發明涉及醫學影像技術領域，尤其涉及一種醫學影像感興趣區域的窗位窗寬計算方法。

背景技術

醫學影像大多數來自CT、DR、US等影像設備，位深常見12~16bit，若要供醫生診斷，需針對不同部位調整合適的窗位窗寬，才能將特定色域范圍灰度在顯示屏（多數8bit）顯示。

在患者病灶尤其是肺結節的診斷、治療等方面，醫生關注重點在于能夠清晰、準確顯示病灶部分的紋理及其周圍臟器等，來自不同影像設備采集的醫學影像或經過其他軟件處理后的醫學影像不能很好地在顯示屏上顯示病灶部分的特征。

發明內容

發明目的：本發明針對上述不足，提出一種醫學影像感興趣區域的自適應調窗方法，針對不同影像設備采集的醫學影像或經過其他軟件處理后的醫學影像，可以分割計算得到病灶部分的窗位窗寬，并映射至顯示屏上顯示出來。

技術方案：

一種醫學影像感興趣區域的窗位窗寬計算方法，包括：

將醫學影像數據的體素灰度值轉換為Hu值；

構建Hu值的頻率直方圖，并根據設定閾值比例剔除頻率直方圖中的無效組；

根據剔除無效組后的頻率直方圖計算得到感興趣區域的Hu值范圍，并據此計算得到窗位窗寬。

所述根據設定閾值比例剔除頻率直方圖中的無效組具體為：

分割頻率直方圖中波峰的個數；

若其不小于根據感興趣區域設置的最大波峰個數，則根據設定閾值比例剔除頻率直方圖兩端的干擾波峰，以剔除干擾波峰后得到的有效組中第一個波谷所在組和最后一個波谷所在組作為有效組兩端所在組；

若其小于根據感興趣區域設置的最大波峰個數，則根據設定閾值比例剔除頻率直方圖兩端的干擾組，以剔除干擾組后得到的有效組中開始組和結尾組作為有效組兩端所在組。

所述分割頻率直方圖中波峰的個數具體為：

根據頻率直方圖中的最大頻率F_max設置頻率直方圖的頻率閾值F_p，定義截止頻率F_C，從0開始不斷按迭代頻率步長d_p抬升F_C，根據識別得到的頻率直方圖中波谷的數量分割得到頻率直方圖中孤立的波峰，直到得到的頻率直方圖中波峰的個數n_peak大于設置的最大波峰個數N_peak或者截止頻率F_C超過頻率閾值F_p。

所述頻率閾值F_p設置為0.1F_max。

所述根據剔除無效組后的頻率直方圖計算得到感興趣區域的Hu值范圍具體為：

定義頻率直方圖上的Hu值的最小值Y_min和最大值Y_max，頻率直方圖的組數為N_hist，則直方圖的寬度為g=(Y_max-Y_min)/N_hist；

定義頻率直方圖的有效組兩端所在組為N_min和N_max，則N_min和N_max對應的Hu值為：