[發明專利]基于CT圖像的巖心目標提取方法

在審

• 申請號: 201310706021.0
• 申請日: 2013-12-20
• 公開（公告）號: CN103700102A
• 公開（公告）日: 2014-04-02
• 發明（設計）人: 方黎勇;陳鵬
• 申請（專利權）人: 電子科技大學
• 主分類號: G06T7/00
• 分類號: G06T7/00
• 代理公司: 成都宏順專利代理事務所(普通合伙) 51227
• 代理人: 王偉;周永宏

1.一種基于CT圖像的巖心目標提取方法，其特征在于：包括以下步驟：

S1.讀取巖心的CT圖像，并將所得CT圖像記為I₀：f(x,y)=g，其中，(x,y)為圖像像素坐標，g為圖像像素坐標(x,y)對應的灰度值，g∈[0,1]；

S2.巖心CT圖像的閾值變換：具體包括以下步驟：

S21.設灰度值變換的閾值為0；

S22.遍歷整張巖心CT圖像的像素，如果像素的灰度值為0，則保持其灰度值不變，如果像素的灰度值大于0，則將其灰度值設置為1，相應的表達式為f1(x,y)=0g=01g>0,]]>所得二值圖像的結果記為I₁；

S3.對圖像I₁進行輪廓提取；

S4.對S3所得的每個輪廓進行填充，即將每個輪廓內部的像素的灰度值都設為1，得到圖像為I₂：f₂(x,y)；

S5.對S4所得每個連通區域進行像素數目統計，記錄每個連通區域包含的像素數目，并對統計結果進行比對，獲得具有最大像素數目的連通區域；

S6.將除具有最大像素數目的連通區域外的連通區域的像素的灰度值都設置為零，得到圖像I₃：f₃(x,y)；

S7.對圖像I₃進行輪廓提取，獲得目標輪廓并標記為I₄：f₄(x,y)；

S8.將圖像I₃的每個像素的灰度值與原始圖像I₀中相應位置的像素的灰度值相乘，即得基于CT圖像的巖心目標圖像，記為I₅：f₅(x,y)＝f₀(x,y)×f₃(x,y)；

S9.將目標圖像I₅與輪廓圖像I₄進行加運算，即得基于CT圖像的巖心目標及巖心目標輪廓。

2.根據權利要求1所述的基于CT圖像的巖心目標提取方法，其特征在于：S3所述的輪廓提取具體包括以下步驟：

S31.遍歷I₁的每個像素，對于灰度值為1的像素，如果該像素的八個鄰域點的灰度值都為1，則將該像素的灰度值設置為零，即可提取I₁中所有輪廓；

S32.按照從左到右、從下到上的順依次對S31所得輪廓進行八鄰域點追蹤，獲得單值化的輪廓。

3.根據權利要求2所述的基于CT圖像的巖心目標提取方法，其特征在于：所述八鄰域點追蹤具體如下：找到第一個灰度值為1的像素，假設為最左下方的邊界點，將該點記為第一邊界點，該像素的右、右上、上、左上四個鄰域中至少有一個是邊界點，記為第二邊界點，從第二邊界點開始，按右、右上、上、左上、左、左下、下和右下的順序找相鄰點中的邊界點，記為第三邊界點，如果第三邊界點就是第一邊界點，則表明閉合輪廓追蹤完成；否則，從第三邊界點繼續，直到找到第一邊界點為止。