[發明專利]一種基于廣義S變換和Teager屬性的肺部CT圖像結節篩選方法

有效

• 申請號: 201910065176.8
• 申請日: 2019-01-23
• 公開（公告）號: CN109829902B
• 公開（公告）日: 2022-04-12
• 發明（設計）人: 孫翎馬;彭真明;蒲恬;蒲紅;趙學功;郭璐;王卓然;袁國慧;唐雨瀟;范文瀾;陳江華
• 申請（專利權）人: 電子科技大學
• 主分類號: G06T7/00
• 分類號: G06T7/00;G06T7/11;G06T7/13;G06T7/136;G06T7/168
• 代理公司: 成都弘毅天承知識產權代理有限公司 51230
• 代理人: 楊保剛

1.一種基于廣義S變換和Teager屬性的肺部CT圖像結節篩選方法，其特征在于，包括以下步驟：

步驟1：對輸入的肺部CT圖像進行水平方向和豎直方向上的廣義S變換，得到水平方向時頻譜和豎直方向時頻譜；

步驟2：對所述水平方向時頻譜和豎直方向時頻譜進行Teager主能量的提取，得到水平方向Teager主能量圖和豎直方向Teager主能量圖；

步驟3：將所述水平方向Teager主能量圖和豎直方向Teager主能量圖進行閾值分割，得到疑似結節；

所述步驟1中對輸入的肺部CT圖像進行預處理，所述預處理包括確定肺部CT圖像行變量和列變量的取值范圍和去除肺部以外的像素，具體如下包括如下：

步驟a：確定輸入的肺部CT圖像中行變量和列變量的取值范圍即非零像素，并計算非零像素所在的最小行mir、最大行mar、最小列mic和最大列mac，計算公式如下：

{(x，y)|I(x，y)＞0}

mir＝min(x)，mar＝max(x)，mic＝min(y)，mac＝max(y)

其中，I(x，y)表示肺部CT圖像在點(x，y)處的灰度值，x是水平方向位置變量，y是豎直方向位置變量，min表示最小值，max表示最大值，mir表示最小行數值，mar表示最大行數值，mic表示最小列數值，mac表示最大列數值；

步驟b：初始化水平方向行變量i、豎直方向列變量j和廣義S變換的窗口調節參數p；

步驟c：提取CT圖像水平方向第i行一維信號和第j列一維信號，去除零像素得到水平方向第i行一維信號sigh和豎直方向第j列一維信號sigv。

2.如權利要求1所述的一種基于廣義S變換和Teager屬性的肺部CT圖像結節篩選方法，其特征在于，所述步驟b具體為：

初始化廣義S變換的窗口調節參數p＝0.1；

計算水平方向時頻譜時，初始化水平方向行變量i＝mir，i∈mir～mar；

計算豎直方向時頻譜時，初始化豎直方向列變量j＝mic，j∈mic～mac。

3.如權利要求1所述的一種基于廣義S變換和Teager屬性的肺部CT圖像結節篩選方法，其特征在于，所述步驟c中，計算第i行一維信號sigh和豎直方向第j列一維信號sigv的公式如下：

sigh＝{I(x，：)|I(x，：)＞0，I(x，：)∈I(i，：)}

sigv＝{I(：,y)|I(：,y)O,I(：,y)∈I(：,j)}

其中，I(i，：)表示肺部CT圖像水平方向的第i行一維信號，I(x，：)表示第i行一維信號在處的灰度值；I(：,j)表示肺部CT圖像豎直方向的第j行一維信號，I(：，y)表示第j列一維信號在y處的灰度值；sigh表示水平方向第i行一維信號，sigv表示豎直方向第j列一維信號。