本發(fā)明涉及圖像處理領(lǐng)域，具體而言，涉及一種肺裂完整度評估方法、裝置和系統(tǒng)。

背景技術(shù)：

肺葉間裂(簡稱肺裂)作為肺葉的邊界把相鄰肺葉隔開，是體現(xiàn)人體肺部結(jié)構(gòu)的重要標(biāo)識。對肺裂的完整度和結(jié)構(gòu)的全面認(rèn)識，在早期肺部疾病檢測，分類，病情發(fā)展以及疾病的治療上，具有很大的臨床實(shí)用價(jià)值。比如，肺裂的不完整性會減少肺部疾病如肺炎在肺葉之間擴(kuò)散的阻礙，也會影響胸腔積液在肺部的分布形態(tài)；肺裂不完整的肺癌病人在手術(shù)移除肺葉后會產(chǎn)生空氣泄漏。肺裂的完整度對于經(jīng)支氣管鏡肺減容手術(shù)適應(yīng)癥的影像學(xué)判定也具有重要意義。采用支氣管鏡下肺減容手術(shù)(安放支氣管腔內(nèi)活瓣)代替?zhèn)鹘y(tǒng)的外科肺減容手術(shù)對于嚴(yán)重肺氣腫病患具有優(yōu)勢，但并不是所有嚴(yán)重肺氣腫病患適用于這種支氣管鏡下肺減容手術(shù)。病患手術(shù)前必須進(jìn)行肺裂完整度分析，肺裂不完整的病人不能考慮活瓣治療。肺部CT影像可以可視化肺裂，因此臨床上采用標(biāo)準(zhǔn)化CT掃描進(jìn)行肺裂完整度的評估，評估方法分為人工打分和計(jì)算機(jī)自動(dòng)評估兩類。

人工打分方法：醫(yī)師手動(dòng)瀏覽肺部CT影像，通過觀察CT影像中的肺裂(顯示為高亮度曲線)，根據(jù)經(jīng)驗(yàn)給出肺裂完整度的百分比。

計(jì)算機(jī)自動(dòng)評估：肺裂自動(dòng)分割后，無論是否為完整的肺裂，都可以得到肺葉的分割，進(jìn)而得到肺葉邊界(視為完整肺裂)，把檢測到的不完整肺裂與肺葉邊界對比，評估肺裂的完整度。因此肺裂完整度的精確度取決于肺裂和肺葉分割結(jié)果的準(zhǔn)確性。有關(guān)肺裂、肺葉檢測分割的算法很多，可概括為基于解剖結(jié)構(gòu)關(guān)系的方法和基于形狀分析的方法兩類。解剖結(jié)構(gòu)關(guān)系的方法利用了肺裂與血管、支氣管之間在解剖上的關(guān)系來確定肺裂位置，進(jìn)而得到肺葉的分割結(jié)果；也可以圖庫為基礎(chǔ)的模版匹配方法來尋找肺裂?；谛螤罘治龅姆椒ㄍㄟ^在二維空間中確定直線或者在三維空間尋找曲面或者平面來實(shí)現(xiàn)肺裂的檢測分割。

(1)人工打分方法的缺點(diǎn)是：

在肺部CT影像上，肺裂顯示為高亮度曲線，其亮度與其他肺部結(jié)構(gòu)如血管、甚至噪音接近，不易區(qū)分。另外，高分辯率CT影像通常包含大量的圖像(>300張)，讓專家通過手動(dòng)方式一張一張圖像去標(biāo)志肺裂所在空間位置及評估肺裂完整度，非常耗時(shí)。

(2)計(jì)算機(jī)自動(dòng)評估的缺點(diǎn)是：

計(jì)算機(jī)自動(dòng)評估方法的目的就是提高肺裂完整度評估的效率。這些自動(dòng)評估方法首先要先分割肺葉，再對已分割的不完整肺裂與肺葉邊界(作為完整肺裂)比較。

肺裂理論上是平滑的曲面，但實(shí)際上由于肺裂檢測分割算法的缺陷，使得分割的肺裂表面不平滑，并且包含大小不等的“洞”。

用于補(bǔ)“洞”的方法有區(qū)域增長法及基于洪水消退模型的各種圖像分割算法，但這些算法會產(chǎn)生“過分割”現(xiàn)象，會把相鄰其他肺葉的部分組織劃分到當(dāng)前肺葉中，導(dǎo)致肺葉分割不準(zhǔn)確。

肺裂分割的不平滑也會導(dǎo)致作為完整肺裂的肺葉邊界不平滑，完整肺裂的面積發(fā)生變化。

肺裂類型的識別不準(zhǔn)確，即無法準(zhǔn)確識別左肺斜裂，右肺斜裂，及右肺水平裂。

這些缺點(diǎn)都會導(dǎo)致肺裂完整度評估結(jié)果偏離真實(shí)值，使得肺裂完整度評估結(jié)果準(zhǔn)確度低。

針對上述的問題，目前尚未提出有效的解決方案。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明實(shí)施例提供了一種肺裂完整度評估方法、裝置和系統(tǒng)，以至少解決現(xiàn)有技術(shù)中肺裂完整度評估準(zhǔn)確度低的技術(shù)問題。

根據(jù)本發(fā)明實(shí)施例的一個(gè)方面，提供了一種肺裂完整度評估方法，包括：獲取肺部影像；檢測所述肺部影像中的肺裂；使用線性濾波器對所述肺裂進(jìn)行平滑處理；從所述肺裂中識別出左肺斜裂、右肺斜裂、右肺水平裂；根據(jù)所述左肺斜裂、所述右肺斜裂、所述右肺水平裂確定出完整肺裂；根據(jù)識別出的所述左肺斜裂、所述右肺斜裂、所述右肺水平裂和所述完整肺裂確定肺裂完整度。

進(jìn)一步地，從所述肺裂中識別出左肺斜裂、右肺斜裂、右肺水平裂包括：將所述肺裂分解為三個(gè)不同類型的肺裂；利用解剖學(xué)結(jié)構(gòu)將所述三個(gè)不同類型的肺裂分別識別為左肺斜裂、右肺斜裂、右肺水平裂。

進(jìn)一步地，將所述肺裂分解為三個(gè)不同類型的肺裂包括：對檢測出的所述肺裂進(jìn)行表面建模，得到肺裂表面，并使用平滑算法平滑所述肺裂表面；計(jì)算所述肺裂表面每一點(diǎn)的曲率張量；移除曲率張量大于等于預(yù)設(shè)數(shù)值的區(qū)域，其中，在移除曲率張量大于等于所述預(yù)設(shè)數(shù)值的區(qū)域之后，肺裂表面分解為多個(gè)面片；將所述多個(gè)面片進(jìn)行聚類，得到三個(gè)不同類型的肺裂。

進(jìn)一步地，將所述多個(gè)面片進(jìn)行聚類，得到三個(gè)不同類型的肺裂包括：確定第一面片，其中，所述第一面片是所述多個(gè)面片中面積最小的面片；沿所述第一面片的法線方向?qū)γ嫫M(jìn)行標(biāo)記；將具有相同標(biāo)記的面片作為一個(gè)類型的肺裂。

進(jìn)一步地，根據(jù)所述左肺斜裂、所述右肺斜裂、所述右肺水平裂確定出完整肺裂包括：計(jì)算單一肺裂的平均法線向量；將肺裂面片的三角形網(wǎng)格頂點(diǎn)投影到與法線向量垂直的平面上；計(jì)算投影后的肺裂面片的包圍盒；將投影平面中具有最小曲率的網(wǎng)格頂點(diǎn)作為目標(biāo)網(wǎng)格頂點(diǎn)；使用所述目標(biāo)網(wǎng)格頂點(diǎn)進(jìn)行表面擬合；使用同一網(wǎng)格內(nèi)所有頂點(diǎn)的法線向量的平均代替頂點(diǎn)各自的法線向量。

進(jìn)一步地，使用線性濾波器對所述肺裂進(jìn)行平滑處理包括：使用各向異性線性濾波器在所述肺裂所在的二維圖像上沿與肺裂曲線相切的方向線性擴(kuò)展。

進(jìn)一步地，使用各向異性線性濾波器在所述肺裂所在的二維圖像上沿與肺裂曲線相切的方向線性擴(kuò)展包括：給定第一像素點(diǎn)，其中所述第一像素點(diǎn)位于檢測出的肺裂上；將目標(biāo)線段繞所述第一像素點(diǎn)旋轉(zhuǎn)多個(gè)預(yù)設(shè)角度，并且，每旋轉(zhuǎn)到一個(gè)預(yù)設(shè)角度，計(jì)算所述目標(biāo)線段上所有像素亮度的平均值，其中，所述目標(biāo)線段的長度為預(yù)設(shè)長度；將像素亮度的平均值最大時(shí)對應(yīng)的預(yù)設(shè)角度作為第一預(yù)設(shè)角度；將目標(biāo)線段繞所述第一像素點(diǎn)旋轉(zhuǎn)所述第一預(yù)設(shè)角度；在所述第一預(yù)設(shè)角度指示的方向上確定第二像素點(diǎn)，其中，所述第二像素點(diǎn)位于檢測出的肺裂上；將所述目標(biāo)線段延伸至所述第二像素點(diǎn)。

進(jìn)一步地，根據(jù)識別出的所述左肺斜裂、所述右肺斜裂、所述右肺水平裂和所述完整肺裂確定肺裂完整度包括：根據(jù)以下公式計(jì)算肺裂完整度：C_f＝A_i/A_c，其中，C_f為肺裂完整度，A_c為所述完整肺裂，A_i為識別出的所述左肺斜裂、所述右肺斜裂、或所述右肺水平裂。

進(jìn)一步地，檢測所述肺部影像中的肺裂包括：將所述肺部影像劃分為多個(gè)細(xì)分體；利用三維平面擬合方法在所述多個(gè)細(xì)分體中尋找肺裂；將尋找到的肺裂作為所述肺部影像中的肺裂。

根據(jù)本發(fā)明實(shí)施例的另一方面，還提供了一種肺裂完整度評估裝置，該裝置包括：獲取單元，用于獲取肺部影像；檢測單元，用于檢測所述肺部影像中的肺裂；平滑處理單元，用于使用線性濾波器對所述肺裂進(jìn)行平滑處理；識別單元，用于從所述肺裂中識別出左肺斜裂、右肺斜裂、右肺水平裂；第一確定單元，用于根據(jù)所述左肺斜裂、所述右肺斜裂、所述右肺水平裂確定出完整肺裂；第二確定單元，用于根據(jù)識別出的所述左肺斜裂、所述右肺斜裂、所述右肺水平裂和所述完整肺裂確定肺裂完整度。

根據(jù)本發(fā)明實(shí)施例的另一方面，還提供了一種肺裂完整度評估系統(tǒng)，該系統(tǒng)包括：CT影像輸入裝置，用于獲取肺部CT影像；肺裂完整度評估裝置，用于根據(jù)所述肺部CT影像計(jì)算肺裂完整度，所述肺裂完整度評估裝置為權(quán)利要求10所述的肺裂完整度評估裝置；結(jié)果輸出裝置，用于輸出計(jì)算出的肺裂完整度。

在本發(fā)明實(shí)施例中，檢測肺部影像中的肺裂，使用線性濾波器對肺裂進(jìn)行平滑處理，使得到的肺裂更加平滑，肺裂類型的識別更準(zhǔn)確，即，從肺裂中識別出左肺斜裂、右肺斜裂、右肺水平裂的準(zhǔn)確度更高，根據(jù)左肺斜裂、右肺斜裂、右肺水平裂確定出的完整肺裂也更準(zhǔn)確，從而，根據(jù)識別出的左肺斜裂、右肺斜裂、右肺水平裂和完整肺裂計(jì)算出的肺裂完整度更準(zhǔn)確，達(dá)到了提高肺裂完整度評估準(zhǔn)確度的技術(shù)效果，進(jìn)而解決了現(xiàn)有技術(shù)中肺裂完整度評估準(zhǔn)確度低的技術(shù)問題。

附圖說明

此處所說明的附圖用來提供對本發(fā)明的進(jìn)一步理解，構(gòu)成本發(fā)明的一部分，本發(fā)明的示意性實(shí)施例及其說明用于解釋本發(fā)明，并不構(gòu)成對本發(fā)明的不當(dāng)限定。在附圖中：

圖1是根據(jù)本發(fā)明實(shí)施例的一種肺裂完整度評估方法的流程圖；

圖2-a是根據(jù)本發(fā)明實(shí)施例的線性形態(tài)濾波器應(yīng)用前的肺裂表面形態(tài)的示意圖；

圖2-b是根據(jù)本發(fā)明實(shí)施例的線性形態(tài)濾波器應(yīng)用后的肺裂表面形態(tài)的示意圖；

圖3是根據(jù)本發(fā)明實(shí)施例的另一種肺裂完整度評估方法的流程圖；

圖4是根據(jù)本發(fā)明實(shí)施例的單一肺裂識別的流程圖；

圖5是根據(jù)本發(fā)明實(shí)施例的肺裂在完成肺裂分類后的肺裂識別的示意圖；

圖6是根據(jù)本發(fā)明實(shí)施例的一種肺裂完整度評估裝置的示意圖；

圖7是根據(jù)本發(fā)明實(shí)施例的一種肺裂完整度評估系統(tǒng)的示意圖。

具體實(shí)施方式

為了使本技術(shù)領(lǐng)域的人員更好地理解本發(fā)明方案，下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分的實(shí)施例，而不是全部的實(shí)施例?；诒景l(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都應(yīng)當(dāng)屬于本發(fā)明保護(hù)的范圍。

需要說明的是，本發(fā)明的說明書和權(quán)利要求書及上述附圖中的術(shù)語“第一”、“第二”等是用于區(qū)別類似的對象，而不必用于描述特定的順序或先后次序。應(yīng)該理解這樣使用的數(shù)據(jù)在適當(dāng)情況下可以互換，以便這里描述的本發(fā)明的實(shí)施例能夠以除了在這里圖示或描述的那些以外的順序?qū)嵤?。此外，術(shù)語“包括”和“具有”以及他們的任何變形，意圖在于覆蓋不排他的包含，例如，包含了一系列步驟或單元的過程、方法、系統(tǒng)、產(chǎn)品或設(shè)備不必限于清楚地列出的那些步驟或單元，而是可包括沒有清楚地列出的或?qū)τ谶@些過程、方法、產(chǎn)品或設(shè)備固有的其它步驟或單元。

根據(jù)本發(fā)明實(shí)施例，提供了一種肺裂完整度評估方法的實(shí)施例，需要說明的是，在附圖的流程圖示出的步驟可以在諸如一組計(jì)算機(jī)可執(zhí)行指令的計(jì)算機(jī)系統(tǒng)中執(zhí)行，并且，雖然在流程圖中示出了邏輯順序，但是在某些情況下，可以以不同于此處的順序執(zhí)行所示出或描述的步驟。

圖1是根據(jù)本發(fā)明實(shí)施例的一種肺裂完整度評估方法的流程圖。如圖1所示，該方法包括如下步驟：

步驟S102，獲取肺部影像。

步驟S104，檢測肺部影像中的肺裂。

步驟S106，使用線性濾波器對肺裂進(jìn)行平滑處理。

步驟S108，從肺裂中識別出左肺斜裂、右肺斜裂、右肺水平裂。

步驟S110，根據(jù)左肺斜裂、右肺斜裂、右肺水平裂確定出完整肺裂。

步驟S112，根據(jù)識別出的左肺斜裂、右肺斜裂、右肺水平裂和完整肺裂確定肺裂完整度。

在本發(fā)明實(shí)施例中，檢測肺部影像中的肺裂，使用線性濾波器對肺裂進(jìn)行平滑處理，使得到的肺裂更加平滑，肺裂類型的識別更準(zhǔn)確，即，從肺裂中識別出左肺斜裂、右肺斜裂、右肺水平裂的準(zhǔn)確度更高，根據(jù)左肺斜裂、右肺斜裂、右肺水平裂確定出的完整肺裂也更準(zhǔn)確，從而，根據(jù)識別出的左肺斜裂、右肺斜裂、右肺水平裂和完整肺裂計(jì)算出的肺裂完整度更準(zhǔn)確，解決了現(xiàn)有技術(shù)中肺裂完整度評估準(zhǔn)確度低的技術(shù)問題，達(dá)到了提高肺裂完整度評估準(zhǔn)確度的技術(shù)效果。

可選地，從肺裂中識別出左肺斜裂、右肺斜裂、右肺水平裂包括：將肺裂分解為三個(gè)不同類型的肺裂；利用解剖學(xué)結(jié)構(gòu)將三個(gè)不同類型的肺裂分別識別為左肺斜裂、右肺斜裂、右肺水平裂。

獲取肺部影像，檢測肺部影像中的肺裂，使用線性濾波器對肺裂進(jìn)行平滑處理之后，得到的肺裂在空間中彼此相連，融為一體，為了正確區(qū)分肺裂，首先把這個(gè)相連的肺裂分成三個(gè)單獨(dú)的面片，進(jìn)而判斷識別三個(gè)單獨(dú)肺裂的類型。為了提高單一肺裂識別的穩(wěn)健性，結(jié)合肺裂的解剖學(xué)結(jié)構(gòu)及肺裂曲率分析結(jié)果，使用了一個(gè)具有方向不變性的肺裂分類方法，來識別肺裂。

可選地，將肺裂分解為三個(gè)不同類型的肺裂包括：對檢測出的肺裂進(jìn)行表面建模，得到肺裂表面，并使用平滑算法平滑肺裂表面；計(jì)算肺裂表面每一點(diǎn)的曲率張量；移除曲率張量大于等于預(yù)設(shè)數(shù)值的區(qū)域，其中，在移除曲率張量大于等于預(yù)設(shè)數(shù)值的區(qū)域之后，肺裂表面分解為多個(gè)面片；將多個(gè)面片進(jìn)行聚類，得到三個(gè)不同類型的肺裂。

可以用移動(dòng)立方體算法(Marching Cubes Algorithm)對檢測到的肺裂進(jìn)行表面建模，并使用拉普拉斯平滑算法平滑肺裂表面。此步驟可以相對準(zhǔn)確估計(jì)肺裂表面上給定點(diǎn)的曲率張量。

計(jì)算肺裂表面每一點(diǎn)的曲率張量。曲率張量大于等于預(yù)設(shè)數(shù)值的區(qū)域?yàn)楦咔蕝^(qū)域。高曲率區(qū)域是各個(gè)肺裂之間的相交部分，把高曲率區(qū)域移除后，肺裂表面自動(dòng)就變成了多個(gè)面片。將多個(gè)面片進(jìn)行聚類，得到三個(gè)不同類型的肺裂。預(yù)設(shè)數(shù)值可以為0.2mm^-1。

通過計(jì)算肺裂表面每一點(diǎn)的曲率張量，移除高曲率區(qū)域，把每個(gè)肺裂從一個(gè)整體分解出來。

可選地，將多個(gè)面片進(jìn)行聚類，得到三個(gè)不同類型的肺裂包括：確定第一面片，其中，第一面片是多個(gè)面片中面積最小的面片；沿第一面片的法線方向?qū)γ嫫M(jìn)行標(biāo)記；將具有相同標(biāo)記的面片作為一個(gè)類型的肺裂。

找到具有最小面積的面片(第一面片)，沿面片各頂點(diǎn)的法線方向標(biāo)記關(guān)聯(lián)肺裂，標(biāo)記的肺裂被分配一個(gè)索引(如某一顏色)。通過簡單的聚類把曲面片連接起來作為一個(gè)類型的肺裂。左肺中，最大聚類用來代表左肺斜裂。右肺中，第一和第二個(gè)最大聚類分別用來代表右肺斜裂和右肺水平裂，其它所有的聚類作為非肺裂區(qū)被移除。

通過使用基于法線向量的聚類方法，將法線向量方向相近或相同的小曲面片聚為一類。

可選地，根據(jù)左肺斜裂、右肺斜裂、右肺水平裂確定出完整肺裂包括：計(jì)算單一肺裂的平均法線向量；將肺裂面片的三角形網(wǎng)格頂點(diǎn)投影到與法線向量垂直的平面上；計(jì)算投影后的肺裂面片的包圍盒；將投影平面中具有最小曲率的網(wǎng)格頂點(diǎn)作為目標(biāo)網(wǎng)格頂點(diǎn)；使用目標(biāo)網(wǎng)格頂點(diǎn)進(jìn)行表面擬合；使用同一網(wǎng)格內(nèi)所有頂點(diǎn)的法線向量的平均代替頂點(diǎn)各自的法線向量。

為了評判肺裂的完整度，需要構(gòu)建一個(gè)完整的肺裂，使之與檢測到的不完整肺裂比較。首先，利用曲面擬合的方法對檢測到的肺裂進(jìn)行平滑處理，再進(jìn)行肺葉分割，構(gòu)建一個(gè)準(zhǔn)確的完整肺裂。選擇有穩(wěn)定法線向量的約束點(diǎn)，來進(jìn)行肺裂隱性表面擬合，這個(gè)條件對于構(gòu)建完整肺裂是至關(guān)重要的。使用多個(gè)小面片平均的法線向量，而不是以小面片各頂點(diǎn)作為約束點(diǎn)用以表面擬合，這種方法可以減少與以法線向量表示的相鄰約束點(diǎn)相關(guān)的隨機(jī)擾動(dòng)。具體地，(1)計(jì)算單一肺裂的平均法線向量；(2)把肺裂面片的三角形網(wǎng)格頂點(diǎn)投影到與法線向量垂直的平面上，并計(jì)算投影后的肺裂面片的包圍盒；(3)找到投影平面中具有最小曲率的網(wǎng)格頂點(diǎn)，用這些頂點(diǎn)進(jìn)行表面擬合，并用同一網(wǎng)格內(nèi)所有頂點(diǎn)的法線向量的平均代替頂點(diǎn)各自的法線向量。采用上述策略，即可得到平滑的肺葉表面。

如果得到的肺裂不是完整平滑的，那么獲得的肺葉邊界不準(zhǔn)確，在把肺葉邊界作為完整肺裂來分析肺裂完整度的過程中，會影響結(jié)果的精確度。在本發(fā)明實(shí)施例中，為了解決這個(gè)問題，使用各向異性線性濾波器在肺裂所在的二維圖像上沿與肺裂曲線相切的方向線性擴(kuò)展，達(dá)到對肺裂進(jìn)行平滑處理的效果。

可選地，使用線性濾波器對肺裂進(jìn)行平滑處理包括：使用各向異性線性濾波器在肺裂所在的二維圖像上沿與肺裂曲線相切的方向線性擴(kuò)展。

可選地，使用各向異性線性濾波器在肺裂所在的二維圖像上沿與肺裂曲線相切的方向線性擴(kuò)展包括：給定第一像素點(diǎn)，其中第一像素點(diǎn)位于檢測出的肺裂上；將目標(biāo)線段繞第一像素點(diǎn)旋轉(zhuǎn)多個(gè)預(yù)設(shè)角度，并且，每旋轉(zhuǎn)到一個(gè)預(yù)設(shè)角度，計(jì)算目標(biāo)線段上所有像素亮度的平均值，其中，目標(biāo)線段的長度為預(yù)設(shè)長度；將像素亮度的平均值最大時(shí)對應(yīng)的預(yù)設(shè)角度作為第一預(yù)設(shè)角度；將目標(biāo)線段繞第一像素點(diǎn)旋轉(zhuǎn)第一預(yù)設(shè)角度；在第一預(yù)設(shè)角度指示的方向上確定第二像素點(diǎn)，其中，第二像素點(diǎn)位于檢測出的肺裂上；將目標(biāo)線段延伸至第二像素點(diǎn)。

第一像素點(diǎn)可以是檢測出的肺裂上任意一個(gè)像素點(diǎn)。

將目標(biāo)線段繞第一像素點(diǎn)旋轉(zhuǎn)多個(gè)預(yù)設(shè)角度，根據(jù)像素亮度從這多個(gè)預(yù)設(shè)角度中篩選出第一預(yù)設(shè)角度，其中，第一預(yù)設(shè)角度滿足的條件是：目標(biāo)線段上所有像素亮度的平均值最大。

將目標(biāo)線段按照第一預(yù)設(shè)角度指示的方向進(jìn)行延長，與檢測出的肺裂交叉于一個(gè)點(diǎn)，這個(gè)點(diǎn)即為第二像素點(diǎn)。

具體過程如下：

在二維CT影像上肺裂表現(xiàn)為高亮曲線，采用一個(gè)線性形態(tài)濾波器來補(bǔ)“洞”，這個(gè)濾波器在二維影像上使用一組短的線段近似表示肺裂。給定檢測到的肺裂上的一個(gè)像素位置(x,y)(第一像素點(diǎn))，線性濾波器通過執(zhí)行如下的圖像亮度變換來計(jì)算在此像素為中心的局部線段的方向。

其中2r為目標(biāo)線段的長度。

此變換的實(shí)施包括四個(gè)步驟：(1)給定線段繞肺裂上的像素點(diǎn)旋轉(zhuǎn)，旋轉(zhuǎn)角度為θ；(2)計(jì)算此線段所有像素亮度的平均值；(3)線段旋轉(zhuǎn)過程中，當(dāng)像素平均亮度達(dá)到最大值時(shí)，線段旋轉(zhuǎn)的角度設(shè)置為θ_max；(4)在該線段方向上找到已檢測的肺裂像素點(diǎn)，并把線段延伸至該點(diǎn)(第二像素點(diǎn))。設(shè)定目標(biāo)線段的長度(即2r)為15mm，旋轉(zhuǎn)角度增量Δθ為10°。

圖2示出了線性形態(tài)濾波器應(yīng)用前后的肺裂表面形態(tài)。圖2-a中，濾波前肺裂不連續(xù)，在二維空間中表現(xiàn)為不連續(xù)的曲線，在三維空間中表現(xiàn)為曲面上的“洞”；圖2-b中，濾波后表現(xiàn)為連續(xù)的曲線(二維)和曲面(三維)。

這個(gè)線性濾波器在原理上類似于形態(tài)學(xué)運(yùn)算或亮度變換，然而，不像傳統(tǒng)的各項(xiàng)同性的“開”/“關(guān)”形態(tài)學(xué)運(yùn)算均勻地在所有方向擴(kuò)展肺裂表面，各向異性線性濾波器在每個(gè)二維圖像上沿與肺裂曲面相切的方向線性擴(kuò)展，使得肺裂平滑效果更好，補(bǔ)“洞”效果更好。

使用一個(gè)各向異性的線性形態(tài)濾波器來平滑肺裂表面及補(bǔ)“洞”，使得肺裂平滑效果更好，補(bǔ)“洞”效果更好。

本發(fā)明實(shí)施例得到的肺裂是完整平滑的，由此獲得肺葉邊界準(zhǔn)確，進(jìn)而把肺葉邊界作為完整肺裂來分析肺裂完整度的精確度明顯提高。

可選地，根據(jù)識別出的左肺斜裂、右肺斜裂、右肺水平裂和完整肺裂確定肺裂完整度包括：根據(jù)以下公式計(jì)算肺裂完整度：C_f＝A_i/A_c，其中，C_f為肺裂完整度，A_c為完整肺裂，A_i為識別出的左肺斜裂、右肺斜裂、或右肺水平裂。

可選地，檢測肺部影像中的肺裂包括：將肺部影像劃分為多個(gè)細(xì)分體；利用三維平面擬合方法在多個(gè)細(xì)分體中尋找肺裂；將尋找到的肺裂作為肺部影像中的肺裂。

為了提高肺裂完整度的評估效率及精確度，本發(fā)明實(shí)施例提供的肺裂完整度評估方法具有如下優(yōu)點(diǎn)：可準(zhǔn)確識別所有肺裂類型，進(jìn)而得到所有肺裂各自的完整度評估結(jié)果，性能穩(wěn)定，肺裂完整度評估準(zhǔn)確，能夠全自動(dòng)的肺裂完整度評估，提高臨床上評估效率。

圖3是根據(jù)本發(fā)明實(shí)施例的另一種肺裂完整度評估方法的流程圖。

如圖3所示，該方法的實(shí)施分為四個(gè)階段：(1)肺裂的檢測；(2)肺裂形態(tài)分析；(3)完整肺裂的構(gòu)建；(4)肺裂完整度的量化評估。

第一階段檢測的肺裂通常是不平滑的，含有很多“洞”，而且肺裂沒有進(jìn)行分類，即沒有分成左肺斜裂，右肺斜裂，及右肺水平裂。在第二個(gè)階段中，對檢測到的肺裂進(jìn)行了形態(tài)分析，包括(1)線性形態(tài)學(xué)分析：使用一個(gè)線性形態(tài)濾波器來平滑肺裂表面及補(bǔ)“洞”；(2)單一肺裂識別：采用一個(gè)漸進(jìn)形狀分解算法把每個(gè)肺裂從融為一個(gè)整體的肺裂中識別分解出來；(3)肺裂類型確定：基于解剖學(xué)信息對肺裂類型進(jìn)行分類。在第三階段，對已識別并分類的肺裂表面進(jìn)行隱性曲面擬合，分割相鄰肺葉，構(gòu)建完整肺裂。第四個(gè)階段，通過計(jì)算實(shí)際檢測到的肺裂和構(gòu)建的完整肺裂的差異，獲得肺裂完整度的量化結(jié)果。

(1)階段I

對于肺裂的檢測，方法很多，如基于解剖結(jié)構(gòu)關(guān)系的方法，利用肺裂附近沒有大的血管的信息，或通過肺裂與血管、支氣管之間在解剖上的關(guān)系來確定肺裂位置；還有基于形狀分析的方法，這類方法通過在二維空間中確定直線或者在三維空間尋找曲面或者平面來實(shí)現(xiàn)。本發(fā)明實(shí)施例采用一種以三維平面擬合為基礎(chǔ)，高效、準(zhǔn)確地識別CT圖像上肺裂的自動(dòng)檢測與分割方法。該方法將CT影像看作是一組三維空間點(diǎn)云，通過將肺部區(qū)域劃分為一些小的細(xì)分體，然后利用三維平面擬合方法在這些小的細(xì)分體中尋找肺裂平面，從而將具有自由曲面特性的肺裂檢測轉(zhuǎn)換為若干平面的檢測。

(2)階段II

線性形態(tài)學(xué)分析

在二維CT影像上肺裂表現(xiàn)為高亮曲線，采用一個(gè)線性形態(tài)濾波器來補(bǔ)“洞”，這個(gè)濾波器在二維影像上使用一組短的線段近似表示肺裂。給定檢測到的肺裂上的一個(gè)像素位置(x,y)，線性濾波器通過執(zhí)行如下的圖像亮度變換來計(jì)算在此像素為中心的局部線段的方向。

其中2r為線段的長度。

此變換的實(shí)施包括四個(gè)步驟：(1)給定線段(目標(biāo)線段)繞肺裂上的像素點(diǎn)(第一像素點(diǎn))旋轉(zhuǎn)，旋轉(zhuǎn)角度為θ；(2)計(jì)算此線段所有像素亮度的平均值；(3)線段旋轉(zhuǎn)過程中，當(dāng)像素平均亮度達(dá)到最大值時(shí)，線段旋轉(zhuǎn)的角度設(shè)置為θ_max；(4)在該線段方向上找到已檢測的肺裂像素點(diǎn)，并把線段延伸至該點(diǎn)(第二像素點(diǎn))。設(shè)定線段的長度(即2r)為15mm，旋轉(zhuǎn)角度增量Δθ為10°。

圖2-a中，濾波前肺裂不連續(xù)，在二維空間中表現(xiàn)為不連續(xù)的曲線，在三維空間中表現(xiàn)為曲面上的“洞”；圖2-b中，濾波后表現(xiàn)為連續(xù)的曲線(二維)和曲面(三維)。

這個(gè)線性濾波器在原理上類似于形態(tài)學(xué)運(yùn)算或亮度變換，然而，不像傳統(tǒng)的各項(xiàng)同性的“開”/“關(guān)”形態(tài)學(xué)運(yùn)算均勻地在所有方向擴(kuò)展肺裂表面，各向異性線性濾波器在每個(gè)二維圖像上沿與肺裂曲面相切的方向線性擴(kuò)展，使得肺裂平滑效果更好，補(bǔ)“洞”效果更好。

單一肺裂識別

在線性形態(tài)學(xué)分析后，進(jìn)行肺裂的識別。肺裂分為右肺斜裂、右肺水平裂、及左肺斜裂。在第一階段檢測出的肺裂在空間中彼此相連，融為一體，為了正確區(qū)分肺裂，首先把這個(gè)相連的肺裂分成三個(gè)單獨(dú)的面片，進(jìn)而判斷識別三個(gè)單獨(dú)肺裂的類型。為了提高單一肺裂識別的穩(wěn)健性，結(jié)合肺裂的解剖學(xué)結(jié)構(gòu)及肺裂曲率分析結(jié)果，使用了一個(gè)具有方向不變性的肺裂分類方法，來識別肺裂，具體操作如圖4所示。

步驟1：用移動(dòng)立方體算法(Marching Cubes Algorithm)對檢測到的肺裂進(jìn)行表面建模，并使用拉普拉斯平滑算法平滑肺裂表面。此步驟可以相對準(zhǔn)確估計(jì)肺裂表面上給定點(diǎn)的曲率張量。

步驟2：計(jì)算肺裂表面每一點(diǎn)的曲率張量，移除高曲率區(qū)域(各個(gè)肺裂之間的相交部分)，把整個(gè)肺裂表面分解為若干個(gè)面片。曲率閾值定為0.2mm^-1。

步驟3：找到具有最小面積的面片，沿面片各頂點(diǎn)的法線方向標(biāo)記關(guān)聯(lián)肺裂，標(biāo)記的肺裂被分配一個(gè)索引(如某一顏色)。通過簡單的聚類把曲面片連接起來作為一個(gè)類型的肺裂。左肺中，最大聚類用來代表左肺斜裂。右肺中，第一和第二個(gè)最大聚類分別用來代表右肺斜裂和右肺水平裂，其它所有的聚類作為非肺裂區(qū)被移除。

步驟4：重復(fù)步驟1至步驟3，直到所有肺裂均已標(biāo)記。至此，檢測到的肺裂都被分解為不同區(qū)域，如圖5所示。

肺裂類型確定

步驟“單一肺裂識別”可以把融合成一個(gè)整體的每個(gè)肺裂區(qū)分開來，還需要知道每個(gè)肺裂屬于哪個(gè)類型。簡單的基于面積大小(理想情況下，斜裂的表面積大于水平裂表面積)來評判肺裂類型是不充分的，比如由于肺裂分割的不完整，有可能導(dǎo)致某個(gè)斜裂表面積小于水平裂。本發(fā)明提出利用解剖學(xué)結(jié)構(gòu)確定肺裂的類型，分類的結(jié)果如圖5所示。

(3)階段III

為了評判肺裂的完整度，需要構(gòu)建一個(gè)完整的肺裂，使之與檢測到的不完整肺裂比較。首先，利用曲面擬合的方法對檢測到的肺裂進(jìn)行平滑處理，再進(jìn)行肺葉分割，構(gòu)建一個(gè)準(zhǔn)確的完整肺裂。選擇有穩(wěn)定法線向量的約束點(diǎn)，來進(jìn)行肺裂隱性表面擬合，這個(gè)條件對于構(gòu)建完整肺裂是至關(guān)重要的。使用多個(gè)小面片平均的法線向量，而不是以小面片各頂點(diǎn)作為約束點(diǎn)用以表面擬合，這種方法可以減少與以法線向量表示的相鄰約束點(diǎn)相關(guān)的隨機(jī)擾動(dòng)。具體地，(1)計(jì)算單一肺裂的平均法線向量；(2)把肺裂面片的三角形網(wǎng)格頂點(diǎn)投影到與法線向量垂直的平面上，并計(jì)算投影后的肺裂面片的包圍盒；(3)找到投影平面中具有最小曲率的網(wǎng)格頂點(diǎn)，用這些頂點(diǎn)進(jìn)行表面擬合，并用同一網(wǎng)格內(nèi)所有頂點(diǎn)的法線向量的平均代替頂點(diǎn)各自的法線向量。采用上述策略，即可得到平滑的肺葉表面。

(4)階段IV

把檢測到的肺裂投影到完整肺裂所在平面上，并與完整肺裂配準(zhǔn)，采用如下公式計(jì)算肺裂完整度。

C_f為肺裂完整度，A_c為完整肺裂，A_i為檢測到的肺裂。

在本發(fā)明實(shí)施例提供的肺裂完整度評估方法中，在第一個(gè)階段利用三維平面擬合方法在這些小的細(xì)分體(肺部區(qū)域)中尋找肺裂平面，從而將具有自由曲面特性的肺裂檢測轉(zhuǎn)換為若干平面的檢測；在第二個(gè)階段中使用一個(gè)各向異性的線性形態(tài)濾波器來平滑肺裂表面及補(bǔ)“洞”，使得肺裂平滑效果更好，補(bǔ)“洞”效果更好；在第二個(gè)階段使用單一肺裂(所有肺裂)識別策略；在第二個(gè)階段通過計(jì)算肺裂表面每一點(diǎn)的曲率張量，移除高曲率區(qū)域，把每個(gè)肺裂從一個(gè)整體分解出來；在第二個(gè)階段利用解剖學(xué)結(jié)構(gòu)確定所有肺裂的類型；在第三個(gè)階段得到平滑肺葉表面；在肺葉分割前先利用曲面擬合的方法對檢測到的肺裂進(jìn)行平滑處理，并選擇有穩(wěn)定法線向量的約束點(diǎn)，來進(jìn)行肺裂隱性表面擬合；在第四個(gè)階段計(jì)算肺裂完整度。由于得到的肺裂是完整平滑的，由此獲得肺葉邊界準(zhǔn)確，進(jìn)而把肺葉邊界作為完整肺裂來分析肺裂完整度的精確度明顯提高。解決了現(xiàn)有技術(shù)中肺裂完整度評估不準(zhǔn)確的問題，解決了現(xiàn)有技術(shù)中肺裂類型識別分類不準(zhǔn)確的問題，可準(zhǔn)確識別所有肺裂類型，進(jìn)而得到所有肺裂各自的完整度評估結(jié)果，性能穩(wěn)定。達(dá)到了全自動(dòng)的肺裂完整度評估，提高了肺裂完整度的評估效率及精確度。

本發(fā)明實(shí)施例提供的肺裂完整度評估方法能夠應(yīng)用于肺減容手術(shù)術(shù)前規(guī)劃平臺，結(jié)果準(zhǔn)確，并且能進(jìn)行自動(dòng)肺裂完整度分析。

本發(fā)明實(shí)施例提供了一種肺裂完整度評估裝置。該肺裂完整度評估裝置可以執(zhí)行上述肺裂完整度評估方法。

如圖6所示，該肺裂完整度評估裝置包括：獲取單元10、檢測單元20、平滑處理單元30、識別單元40、第一確定單元50、第二確定單元60。

獲取單元10，用于獲取肺部影像。

檢測單元20，用于檢測肺部影像中的肺裂。

平滑處理單元30，用于使用線性濾波器對肺裂進(jìn)行平滑處理。

識別單元40，用于從肺裂中識別出左肺斜裂、右肺斜裂、右肺水平裂。

第一確定單元50，用于根據(jù)左肺斜裂、右肺斜裂、右肺水平裂確定出完整肺裂。

第二確定單元60，用于根據(jù)識別出的左肺斜裂、右肺斜裂、右肺水平裂和完整肺裂確定肺裂完整度。

本發(fā)明實(shí)施例提供了一種肺裂完整度評估系統(tǒng)。

如圖7所示，該肺裂完整度評估系統(tǒng)包括：CT影像輸入裝置72、肺裂完整度評估裝置74、結(jié)果輸出裝置76。肺裂完整度評估裝置74可以為圖6所示的肺裂完整度評估裝置。

CT影像輸入裝置72，用于獲取肺部CT影像。

肺裂完整度評估裝置74，用于根據(jù)肺部CT影像計(jì)算肺裂完整度。

結(jié)果輸出裝置76，用于輸出計(jì)算出的肺裂完整度。

CT影像輸入裝置72可以為：CT機(jī)(計(jì)算機(jī)X線斷層攝影機(jī))、PACS、影像存儲系統(tǒng)等，其與本發(fā)明裝置通過DICOM或TCP/IP接口，自動(dòng)或手動(dòng)上傳肺部CT影像，用戶也可通過互聯(lián)網(wǎng)瀏覽器以HTTP方式上傳影像。

肺裂完整度評估裝置74接收到肺部CT影像后，啟動(dòng)肺裂完整度分析，分析結(jié)果包括：(1)以不同顏色標(biāo)記的所有肺裂的三維可視化結(jié)果；(2)與不同顏色的所有肺裂一一對應(yīng)的肺裂完整度量化結(jié)果(％百分比)(3)所有肺裂的形態(tài)特征量化結(jié)果，如亮度、面積、曲率等。

結(jié)果輸出裝置76將分析結(jié)果以圖像、統(tǒng)計(jì)圖表等形式通過DICOM瀏覽器或網(wǎng)頁的形式輸出。

在本發(fā)明的上述實(shí)施例中，對各個(gè)實(shí)施例的描述都各有側(cè)重，某個(gè)實(shí)施例中沒有詳述的部分，可以參見其他實(shí)施例的相關(guān)描述。

在本發(fā)明所提供的幾個(gè)實(shí)施例中，應(yīng)該理解到，所揭露的技術(shù)內(nèi)容，可通過其它的方式實(shí)現(xiàn)。其中，以上所描述的裝置實(shí)施例僅僅是示意性的，例如所述單元的劃分，可以為一種邏輯功能劃分，實(shí)際實(shí)現(xiàn)時(shí)可以有另外的劃分方式，例如多個(gè)單元或組件可以結(jié)合或者可以集成到另一個(gè)系統(tǒng)，或一些特征可以忽略，或不執(zhí)行。另一點(diǎn)，所顯示或討論的相互之間的耦合或直接耦合或通信連接可以是通過一些接口，單元或模塊的間接耦合或通信連接，可以是電性或其它的形式。

所述作為分離部件說明的單元可以是或者也可以不是物理上分開的，作為單元顯示的部件可以是或者也可以不是物理單元，即可以位于一個(gè)地方，或者也可以分布到多個(gè)單元上?？梢愿鶕?jù)實(shí)際的需要選擇其中的部分或者全部單元來實(shí)現(xiàn)本實(shí)施例方案的目的。

另外，在本發(fā)明各個(gè)實(shí)施例中的各功能單元可以集成在一個(gè)處理單元中，也可以是各個(gè)單元單獨(dú)物理存在，也可以兩個(gè)或兩個(gè)以上單元集成在一個(gè)單元中。上述集成的單元既可以采用硬件的形式實(shí)現(xiàn)，也可以采用軟件功能單元的形式實(shí)現(xiàn)。

所述集成的單元如果以軟件功能單元的形式實(shí)現(xiàn)并作為獨(dú)立的產(chǎn)品銷售或使用時(shí)，可以存儲在一個(gè)計(jì)算機(jī)可讀取存儲介質(zhì)中?；谶@樣的理解，本發(fā)明的技術(shù)方案本質(zhì)上或者說對現(xiàn)有技術(shù)做出貢獻(xiàn)的部分或者該技術(shù)方案的全部或部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，該計(jì)算機(jī)軟件產(chǎn)品存儲在一個(gè)存儲介質(zhì)中，包括若干指令用以使得一臺計(jì)算機(jī)設(shè)備(可為個(gè)人計(jì)算機(jī)、服務(wù)器或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本發(fā)明各個(gè)實(shí)施例所述方法的全部或部分步驟。而前述的存儲介質(zhì)包括：U盤、只讀存儲器(ROM，Read-Only Memory)、隨機(jī)存取存儲器(RAM，Random Access Memory)、移動(dòng)硬盤、磁碟或者光盤等各種可以存儲程序代碼的介質(zhì)。

以上所述僅是本發(fā)明的優(yōu)選實(shí)施方式，應(yīng)當(dāng)指出，對于本技術(shù)領(lǐng)域的普通技術(shù)人員來說，在不脫離本發(fā)明原理的前提下，還可以做出若干改進(jìn)和潤飾，這些改進(jìn)和潤飾也應(yīng)視為本發(fā)明的保護(hù)范圍。