來(lái)指令相關(guān)的硬件來(lái)完成，所述的程序可存儲(chǔ)于一般計(jì)算機(jī)可讀取存儲(chǔ)介質(zhì)中，該程序在執(zhí)行時(shí)，可包括如上述各方法的實(shí)施例的流程。其中，所述的存儲(chǔ)介質(zhì)可為磁碟、光盤(pán)、只讀存儲(chǔ)記憶體(Read-Only Memory，ROM)或隨機(jī)存儲(chǔ)記憶體(Random AccessMemory，RAM)等。
[0245]本領(lǐng)域技術(shù)人員還可以了解到本發(fā)明實(shí)施例列出的各種功能是通過(guò)硬件還是軟件來(lái)實(shí)現(xiàn)取決于特定的應(yīng)用和整個(gè)系統(tǒng)的設(shè)計(jì)要求。本領(lǐng)域技術(shù)人員可以對(duì)于每種特定的應(yīng)用，可以使用各種方法實(shí)現(xiàn)所述的功能，但這種實(shí)現(xiàn)不應(yīng)被理解為超出本發(fā)明實(shí)施例保護(hù)的范圍。
[0246]本發(fā)明中應(yīng)用了具體實(shí)施例對(duì)本發(fā)明的原理及實(shí)施方式進(jìn)行了闡述，以上實(shí)施例的說(shuō)明只是用于幫助理解本發(fā)明的方法及其核心思想；同時(shí)，對(duì)于本領(lǐng)域的一般技術(shù)人員，依據(jù)本發(fā)明的思想，在【具體實(shí)施方式】及應(yīng)用范圍上均會(huì)有改變之處，綜上所述，本說(shuō)明書(shū)內(nèi)容不應(yīng)理解為對(duì)本發(fā)明的限制。
【主權(quán)項(xiàng)】
1.一種基于胸部橫斷面CT圖像的肺部分割提取方法，其特征是，所述的方法包括: 獲取胸部橫斷面的CT圖像； 對(duì)所述的CT圖像進(jìn)行預(yù)處理； 對(duì)預(yù)處理后的CT圖像進(jìn)行閾值分割； 對(duì)閾值分割后的CT圖像進(jìn)行肺區(qū)域提取。2.根據(jù)權(quán)利要求1所述的方法，其特征是，采用中值濾波與小波去噪的聯(lián)合去噪方法對(duì)所述的CT圖像進(jìn)行預(yù)處理。3.根據(jù)權(quán)利要求1所述的方法，其特征是，所述對(duì)預(yù)處理后的CT圖像進(jìn)行閾值分割包括: 從預(yù)處理后的CT圖像中確定出灰度值小于O的像素； 將灰度值小于O的像素置為O，得到第一圖像； 確定所述第一圖像的灰度直方圖； 確定所述灰度直方圖中兩波峰間的波谷，視為分割閾值； 根據(jù)所述的分割閾值對(duì)所述的第一圖像進(jìn)行初步二值分割，得到第二圖像； 對(duì)所述的第二圖像進(jìn)行取反，得到第三圖像； 采用形態(tài)學(xué)開(kāi)運(yùn)算對(duì)所述第三圖像進(jìn)行去除處理，得到二值圖像； 對(duì)所述的二值圖像中的連通區(qū)域按照從左到右、從上到下的策略進(jìn)行區(qū)域標(biāo)號(hào)，得到閾值分割后的CT圖像。4.根據(jù)權(quán)利要求3所述的方法，其特征是，對(duì)閾值分割后的CT圖像進(jìn)行肺區(qū)域提取包括: 從閾值分割后的CT圖像中提取出區(qū)域標(biāo)號(hào)后的連通區(qū)域； 判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最左端像素所在列是否大于20; 當(dāng)判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最右端像素所在列是否小于290列； 當(dāng)判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的面積是否大于2000; 當(dāng)判斷為是時(shí)，所述區(qū)域標(biāo)號(hào)為2的連通區(qū)域?yàn)樽筮叿螀^(qū)域。5.根據(jù)權(quán)利要求4所述的方法，其特征是，對(duì)閾值分割后的CT圖像進(jìn)行肺區(qū)域提取還包括:從區(qū)域標(biāo)號(hào)不為I和2的連通區(qū)域中選取出面積大于2000的連通區(qū)域，即為右邊肺區(qū)域。6.根據(jù)權(quán)利要求3所述的方法，其特征是，對(duì)閾值分割后的CT圖像進(jìn)行肺區(qū)域提取包括: 從閾值分割后的CT圖像中提取出區(qū)域標(biāo)號(hào)后的連通區(qū)域； 判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最左端像素所在列是否大于20; 當(dāng)判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的面積是否大于2000; 當(dāng)判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最右端像素所在列是否小于290列； 當(dāng)判斷為否時(shí)，所述區(qū)域標(biāo)號(hào)為2的連通區(qū)域?yàn)樽笥曳握尺B肺區(qū)域。7.根據(jù)權(quán)利要求3所述的方法，其特征是，對(duì)閾值分割后的CT圖像進(jìn)行肺區(qū)域提取包括: 從閾值分割后的CT圖像中提取出區(qū)域標(biāo)號(hào)后的連通區(qū)域； 判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最左端像素所在列是否大于20; 當(dāng)判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最右端像素所在列是否小于290列； 當(dāng)判斷為否時(shí)，根據(jù)所述的CT圖像確定所述區(qū)域標(biāo)號(hào)為3的連通區(qū)域是否為左邊肺區(qū)域。8.根據(jù)權(quán)利要求3所述的方法，其特征是，對(duì)閾值分割后的CT圖像進(jìn)行肺區(qū)域提取包括: 從閾值分割后的CT圖像中提取出區(qū)域標(biāo)號(hào)后的連通區(qū)域； 判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最右端像素所在列是否小于290列； 當(dāng)判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最左端像素所在列是否大于20; 當(dāng)判斷為否時(shí)，根據(jù)所述的CT圖像確定所述區(qū)域標(biāo)號(hào)為3的連通區(qū)域是否為左邊肺區(qū)域。9.根據(jù)權(quán)利要求7或8所述的方法，其特征是，根據(jù)所述的CT圖像確定所述區(qū)域標(biāo)號(hào)為3的連通區(qū)域是否為左邊肺區(qū)域包括: 判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域最左端像素所在列是否小于2; 當(dāng)判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為3的連通區(qū)域的最右端像素所在列是否小于290; 當(dāng)判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為3的連通區(qū)域的面積是否大于2000; 當(dāng)判斷為是時(shí)，所述區(qū)域標(biāo)號(hào)為3的連通區(qū)域?yàn)樽筮叿螀^(qū)域。10.根據(jù)權(quán)利要求9所述的方法，其特征是，根據(jù)所述的CT圖像確定所述區(qū)域標(biāo)號(hào)為3的連通區(qū)域是否為左邊肺部還包括: 從區(qū)域標(biāo)號(hào)不為1、2、3和4的連通區(qū)域中選取出面積大于2000的連通區(qū)域，即為右邊肺區(qū)域。11.根據(jù)權(quán)利要求3所述的方法，其特征是，對(duì)閾值分割后的CT圖像進(jìn)行肺區(qū)域提取包括: 從閾值分割后的CT圖像中提取出區(qū)域標(biāo)號(hào)后的連通區(qū)域； 從所述連通區(qū)域中確定出區(qū)域集合，所述區(qū)域集合包括了多個(gè)區(qū)域，每個(gè)區(qū)域均滿(mǎn)足最左端像素所在列小于2且面積大于100000; 從所述區(qū)域集合中確定出區(qū)域標(biāo)號(hào)最大的區(qū)域，稱(chēng)為開(kāi)始區(qū)域； 從所述連通區(qū)域中篩選出疑似區(qū)域，所述疑似區(qū)域的區(qū)域標(biāo)號(hào)大于所述開(kāi)始區(qū)域的區(qū)域標(biāo)號(hào)； 從所述疑似區(qū)域中確定出目標(biāo)區(qū)域，所述目標(biāo)區(qū)域的面積大于2000，所述目標(biāo)區(qū)域即為肺區(qū)域。12.—種基于胸部橫斷面CT圖像的肺部分割提取系統(tǒng)，其特征是，所述的系統(tǒng)包括: CT圖像獲取裝置，用于獲取胸部橫斷面的CT圖像； 預(yù)處理裝置，用于對(duì)所述的CT圖像進(jìn)行預(yù)處理； 閾值分割裝置，用于對(duì)預(yù)處理后的CT圖像進(jìn)行閾值分割； 肺區(qū)域提取裝置，用于對(duì)閾值分割后的CT圖像進(jìn)行肺區(qū)域提取。13.根據(jù)權(quán)利要求12所述的系統(tǒng)，其特征是，采用中值濾波與小波去噪的聯(lián)合去噪方法對(duì)所述的CT圖像進(jìn)行預(yù)處理。14.根據(jù)權(quán)利要求12所述的系統(tǒng)，其特征是，所述的閾值分割裝置包括: 確定模塊，用于從預(yù)處理后的CT圖像中確定出灰度值小于O的像素； 第一圖像確定模塊，用于將灰度值小于O的像素置為O，得到第一圖像； 直方圖確定模塊，用于確定所述第一圖像的灰度直方圖； 分割閾值確定模塊，用于確定所述灰度直方圖中兩波峰間的波谷，視為分割閾值； 第二圖像確定模塊，用于根據(jù)所述的分割閾值對(duì)所述的第一圖像進(jìn)行初步二值分割，得到第二圖像； 第三圖像確定模塊，用于對(duì)所述的第二圖像進(jìn)行取反，得到第三圖像； 二值圖像確定模塊，用于采用形態(tài)學(xué)開(kāi)運(yùn)算對(duì)所述第三圖像進(jìn)行去除處理，得到二值圖像； 區(qū)域標(biāo)號(hào)模塊，用于對(duì)所述的二值圖像中的連通區(qū)域按照從左到右、從上到下的策略進(jìn)行區(qū)域標(biāo)號(hào)，得到閾值分割后的CT圖像。15.根據(jù)權(quán)利要求14所述的系統(tǒng)，其特征是，所述的肺區(qū)域提取裝置包括: 連通區(qū)域提取模塊，用于從閾值分割后的CT圖像中提取出區(qū)域標(biāo)號(hào)后的連通區(qū)域； 第一判斷模塊，用于判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最左端像素所在列是否大于20; 第二判斷模塊，用于當(dāng)所述的第一判斷模塊判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最右端像素所在列是否小于290列； 第三判斷模塊，用于當(dāng)所述的第二判斷模塊判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的面積是否大于2000 ； 第一肺部確定模塊，用于當(dāng)所述第三判斷模塊判斷為是時(shí)，確定出所述區(qū)域標(biāo)號(hào)為2的連通區(qū)域?yàn)樽筮叿螀^(qū)域。16.根據(jù)權(quán)利要求15所述的系統(tǒng)，其特征是，所述的閾值分割裝置還包括: 第二肺部確定模塊，用于從區(qū)域標(biāo)號(hào)不為I和2的連通區(qū)域中選取出面積大于2000的連通區(qū)域，即為右邊肺區(qū)域。17.根據(jù)權(quán)利要求14所述的系統(tǒng)，其特征是，所述的閾值分割裝置還包括: 第四判斷模塊，用于當(dāng)所述的第一判斷模塊判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的面積是否大于2000 ； 第五判斷模塊，用于當(dāng)所述第四判斷模塊判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最右端像素所在列是否小于290列； 第三肺部確定模塊，用于當(dāng)所述第四判斷模塊判斷為否時(shí)，確定出所述區(qū)域標(biāo)號(hào)為2的連通區(qū)域?yàn)樽笥曳握尺B肺區(qū)域。18.根據(jù)權(quán)利要求14所述的系統(tǒng)，其特征是，所述的肺區(qū)域提取裝置包括: 第六判斷模塊，用于判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最右端像素所在列是否小于290列； 第七判斷模塊，用于當(dāng)所述第六判斷模塊判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域的最左端像素所在列是否大于20 ； 當(dāng)所述第七判斷模塊判斷為否時(shí)，所述的肺區(qū)域提取裝置還包括:第四肺部確定模塊，用于根據(jù)所述的CT圖像確定所述區(qū)域標(biāo)號(hào)為3的連通區(qū)域是否為左邊肺區(qū)域。19.根據(jù)權(quán)利要求14所述的系統(tǒng)，其特征是，當(dāng)所述第二判斷模塊判斷為否時(shí)，所述的肺區(qū)域提取裝置還包括: 第四肺部確定模塊，用于根據(jù)所述的CT圖像確定所述區(qū)域標(biāo)號(hào)為3的連通區(qū)域是否為左邊肺區(qū)域。20.根據(jù)權(quán)利要求18或19所述的系統(tǒng)，其特征是，所述的第四肺部確定模塊包括: 第一判斷單元，用于判斷區(qū)域標(biāo)號(hào)為2的連通區(qū)域最左端像素所在列是否小于2; 第二判斷單元，用于當(dāng)所述第一判斷單元判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為3的連通區(qū)域的最右端像素所在列是否小于290; 第三判斷單元，用于當(dāng)所述第二判斷單元判斷為是時(shí)，判斷區(qū)域標(biāo)號(hào)為3的連通區(qū)域的面積是否大于2000; 第一肺部確定單元，用于當(dāng)所述第三判斷單元判斷為是時(shí)，確定出所述區(qū)域標(biāo)號(hào)為3的連通區(qū)域?yàn)樽筮叿螀^(qū)域。21.根據(jù)權(quán)利要求20所述的系統(tǒng)，其特征是，所述的第四肺部確定模塊還包括: 第二部確定單元，用于從區(qū)域標(biāo)號(hào)不為1、2、3和4的連通區(qū)域中選取出面積大于2000的連通區(qū)域，即為右邊肺區(qū)域。22.根據(jù)權(quán)利要求14所述的系統(tǒng)，其特征是，所述的肺區(qū)域提取裝置包括: 區(qū)預(yù)計(jì)和確定模塊，用于從所述連通區(qū)域中確定出區(qū)域集合，所述區(qū)域集合包括了多個(gè)區(qū)域，每個(gè)區(qū)域均滿(mǎn)足最左端像素所在列小于2且面積大于100000; 開(kāi)始區(qū)域確定模塊，用于從所述區(qū)域集合中確定出區(qū)域標(biāo)號(hào)最大的區(qū)域，稱(chēng)為開(kāi)始區(qū)域； 疑似區(qū)域確定模塊，用于從所述連通區(qū)域中篩選出疑似區(qū)域，所述疑似區(qū)域的區(qū)域標(biāo)號(hào)大于所述開(kāi)始區(qū)域的區(qū)域標(biāo)號(hào)； 目標(biāo)區(qū)域確定模塊，用于從所述疑似區(qū)域中確定出目標(biāo)區(qū)域，所述目標(biāo)區(qū)域的面積大于2000，所述目標(biāo)區(qū)域即為肺區(qū)域。
【專(zhuān)利摘要】本發(fā)明提供一種基于胸部橫斷面CT圖像的肺部分割提取方法以及系統(tǒng)，所述方法包括：獲取胸部橫斷面的CT圖像；對(duì)所述的CT圖像進(jìn)行預(yù)處理；對(duì)預(yù)處理后的CT圖像進(jìn)行閾值分割；對(duì)閾值分割后的CT圖像進(jìn)行肺區(qū)域提取。本發(fā)明能夠?qū)崿F(xiàn)對(duì)肺區(qū)域的精準(zhǔn)分割，保證肺實(shí)質(zhì)區(qū)域分割的完整性，避免由于肺區(qū)域的邊緣缺失及區(qū)域的缺失而在后續(xù)診斷過(guò)程中造成漏診的問(wèn)題。
【IPC分類(lèi)】G06T7/00
【公開(kāi)號(hào)】CN105574882
【申請(qǐng)?zhí)枴緾N201511023356
【發(fā)明人】劉記奎, 李燁, 蔡云鵬, 尹麗妍
【申請(qǐng)人】中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院
【公開(kāi)日】2016年5月11日
【申請(qǐng)日】2015年12月30日