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

有效

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

1.一種醫學影像感興趣區域的窗位窗寬計算方法，其特征在于：包括：

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

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

根據頻率直方圖有效組兩端所在組、頻率直方圖的寬度及其Hu值的最小值和最大值進行映射計算得到感興趣區域的Hu值范圍，并據此計算得到窗位窗寬。

2.根據權利要求1所述的醫學影像感興趣區域的窗位窗寬計算方法，其特征在于：所述分割頻率直方圖中波峰的個數具體為：

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

3.根據權利要求2所述的醫學影像感興趣區域的窗位窗寬計算方法，其特征在于：所述頻率閾值F_p設置為0.1F_max。

4.根據權利要求1所述的醫學影像感興趣區域的窗位窗寬計算方法，其特征在于：所述根據剔除無效組后的頻率直方圖計算得到感興趣區域的Hu值范圍具體為：

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

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

P_min=Y_min+(N_min-α)*g

P_max=Y_max+(N_max+β)*g

其中，α與β分別為縮小系數和放大系數；

從而計算得到醫學影像數據感興趣區域的Hu值范圍[P_min, P_max]。

5.根據權利要求4所述的醫學影像感興趣區域的窗位窗寬計算方法，其特征在于：所述感興趣區域為肺部，則所述最大波峰個數設置為2。

6.根據權利要求5所述的醫學影像感興趣區域的窗位窗寬計算方法，其特征在于：所述βα。

7.根據權利要求6所述的醫學影像感興趣區域的窗位窗寬計算方法，其特征在于：所述α=0，所述β=1。

8.根據權利要求1所述的醫學影像感興趣區域的窗位窗寬計算方法，其特征在于：所述設定閾值比例設為10%。

9.根據權利要求1所述的醫學影像感興趣區域的窗位窗寬計算方法，其特征在于：還包括映射至顯示屏的步驟：

將感興趣區域的Hu值范圍變換為對應的體素灰度值范圍[pixel_min, pixel_max]，并通過如下公式映射至顯示屏的顯示范圍：

(pixel-pixel_min)/(pixel_max-pixel_min)*256

其中，pixel表示醫學影像數據中某一體素的灰度值。