基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法
【專利摘要】本發(fā)明提出了基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，包括以下步驟：(1)圖像預(yù)處理：獲得冠狀動(dòng)脈分割圖像及中心線；(2)剪枝操作；(3)球心及球空間的定義；(4)左、右冠狀動(dòng)脈的劃分規(guī)則；(5)左冠狀動(dòng)脈解剖學(xué)命名算法；(6)右冠狀動(dòng)脈解剖學(xué)命名算法。該發(fā)明利用心臟形似倒置圓錐體的特征，以血管分叉點(diǎn)為球心建立球坐標(biāo)系，根據(jù)冠狀動(dòng)脈各分支節(jié)段的解剖學(xué)走形和相互之間的幾何結(jié)構(gòu)關(guān)系，在球空間中對(duì)血管進(jìn)行定位剖分，達(dá)到自動(dòng)分段與解剖學(xué)標(biāo)記的目的。解決現(xiàn)有的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法計(jì)算時(shí)間較長，算法復(fù)雜度較高，以及無法窮舉冠狀動(dòng)脈分布類型導(dǎo)致分支匹配不準(zhǔn)確等問題。相比于將所提取的血管與先驗(yàn)?zāi)Ｐ推ヅ涞姆椒?，這種方法的優(yōu)點(diǎn)在于節(jié)省時(shí)間、標(biāo)記準(zhǔn)確，且能對(duì)更多的節(jié)段進(jìn)行識(shí)別和標(biāo)記。
【專利說明】
基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及醫(yī)學(xué)圖像處理技術(shù)領(lǐng)域，特別是用于冠狀動(dòng)脈圖像的基于球空間劃分 的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記。
【背景技術(shù)】
[0002] 近年來心血管疾病的發(fā)病率和死亡率都高居疾病榜首位且呈年輕化趨勢(shì)，其中較 為嚴(yán)重的冠狀動(dòng)脈粥樣硬化等會(huì)直接導(dǎo)致心肌梗死和心臟受損。早期冠狀動(dòng)脈血管壁上微 小的病灶并不容易被發(fā)現(xiàn)，這類病灶之后迅速發(fā)展將引起嚴(yán)重后果，甚至危及生命。因此微 小病灶的早期診療，冠狀動(dòng)脈血管各節(jié)段情況的進(jìn)展管理尤為重要。作為無創(chuàng)的影像獲取 方式，冠脈CTA被廣泛地應(yīng)用于冠狀動(dòng)脈疾病的診斷。按照2009年《SCCT guidelines for the interpretation and reporting of coronary computed tomographic angiography [1]》的要求，放射科醫(yī)生和心外科醫(yī)生通常會(huì)報(bào)告解剖學(xué)的相關(guān)檢查信息，例如冠狀動(dòng)脈 的分形和各節(jié)段情況，硬化斑塊的節(jié)段位置，以及每條動(dòng)脈或每個(gè)分支的鈣化、狹窄和閉塞 等。由于冠狀動(dòng)脈的微小病變不容易被人眼所識(shí)別，人工解讀冠脈CTA圖像時(shí)醫(yī)生要借助多 種成像模式才能診斷，耗時(shí)、費(fèi)勁且容易引起誤差。因此需要高效準(zhǔn)確的方法，可以自動(dòng)為 冠狀動(dòng)脈分段并標(biāo)記其解剖學(xué)名稱。
[0003] 利用解剖學(xué)信息對(duì)冠脈樹進(jìn)行標(biāo)記是冠狀動(dòng)脈樹自動(dòng)標(biāo)記可用的方法，已有的研 究側(cè)重于利用建立先驗(yàn)?zāi)Ｐ蛠肀碚鹘馄市畔?，然后求待?biāo)記模型和先驗(yàn)?zāi)Ｐ椭g的相似性 來確定冠脈節(jié)段解剖學(xué)分類。Claire Chalopin等人[2]設(shè)計(jì)了二維X射線血管造影的冠狀 動(dòng)脈樹標(biāo)記方法，他們利用三維拓?fù)淠Ｐ徒⒘巳S仿真人體形態(tài)模型Coronix，在其不同 視角下進(jìn)行分析得到二維拓?fù)淠Ｐ蛿?shù)據(jù)基礎(chǔ)，再將病人造影圖像中提取的血管骨架與具有 最相似血管網(wǎng)形狀的二維拓?fù)淠Ｐ瓦M(jìn)行比較，以分層的方式從主要?jiǎng)用}到分支逐步標(biāo)記左 冠狀動(dòng)脈，該方法的不足之處在于二維圖像中存在冠狀動(dòng)脈重疊的現(xiàn)象，在匹配解讀時(shí)可 能產(chǎn)生誤差;Guanyu Yang等人[3]構(gòu)建了一個(gè)基于統(tǒng)計(jì)學(xué)結(jié)果先驗(yàn)知識(shí)的右優(yōu)勢(shì)型冠狀動(dòng) 脈樹三維模型，將從冠脈CTA數(shù)據(jù)集中提取的冠狀動(dòng)脈樹與三維模型配準(zhǔn)，主要分支利用關(guān) 鍵點(diǎn)的權(quán)重因子進(jìn)行配準(zhǔn)，所有分段利用代價(jià)函數(shù)進(jìn)行配準(zhǔn)，最后根據(jù)臨床標(biāo)準(zhǔn)對(duì)以上配 準(zhǔn)標(biāo)記結(jié)果進(jìn)行調(diào)整，得到最優(yōu)的冠狀動(dòng)脈樹自動(dòng)標(biāo)記結(jié)果;Roman Goldenberg等人[4]倉ij 建了基于冠脈CTA訓(xùn)練數(shù)據(jù)集人工標(biāo)記的概率解剖模型，根據(jù)最大可能性方法和一系列檢 查主要?jiǎng)用}間的預(yù)期空間關(guān)系的測試，對(duì)冠狀動(dòng)脈四條主要分支（左冠脈主干、前降支、旋 支和右冠脈)進(jìn)行標(biāo)記。以上方法都是根據(jù)一定的解剖學(xué)先驗(yàn)知識(shí)模型，利用數(shù)學(xué)方法尋找 待命名冠狀動(dòng)脈樹與模型之間的最近似方案并進(jìn)行匹配，得到分段和解剖學(xué)標(biāo)記。此類方 法計(jì)算時(shí)間較長，算法復(fù)雜度較高，由于冠狀動(dòng)脈的個(gè)體差異較大，無法窮舉所有的冠狀動(dòng) 脈分布類型，對(duì)與模型不是同一分布類型、分支較為復(fù)雜或因病變而發(fā)生明顯改變的冠狀 動(dòng)脈并不適用，所得結(jié)果中還可能存在較多未能找到相應(yīng)匹配的分支。
[0004] [l]Raff G L,Abidov A,Achenbach S,et al.SCCT guidelines for the interpretation and reporting of coronary computed tomographic angiography.
[J]-Journal of Cardiovascular Computed Tomography,2009,3(2):122-136.
[0005] [2]Chalopin C,Finet G,Magnin I.Modeling the 3D coronary tree for labeling purposes.[J].Medical Image Analysis,2002,5(4):301-315.
[0006] [3]Yang G,Broersen A,Petr R,et al.Automatic coronary artery tree labeling in coronary computed tomographic angiography datasets[C]·Computing in Cardiology,2011,38:109-112.
[0007] [4]Roman G,Dov E,Grigory B,et al.Computer-aided simple triage(CAST) for coronary CT angiography(CCTA)[J]. International Journal of Computer Assisted Radiology & Surgery,2012,7(6):819-827.

【發(fā)明內(nèi)容】

[0008] 針對(duì)現(xiàn)有技術(shù)的不足，本發(fā)明提出了基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖 學(xué)標(biāo)記方法。該發(fā)明利用心臟形似倒置圓錐體的特征，以血管分叉點(diǎn)為球心建立球坐標(biāo)系， 根據(jù)冠狀動(dòng)脈各分支節(jié)段的解剖學(xué)走形和相互之間的幾何結(jié)構(gòu)關(guān)系，在球空間中對(duì)血管進(jìn) 行定位剖分，達(dá)到自動(dòng)分段與解剖學(xué)標(biāo)記的目的。解決現(xiàn)有的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué) 標(biāo)記方法計(jì)算時(shí)間較長，算法復(fù)雜度較高，以及無法窮舉冠狀動(dòng)脈分布類型導(dǎo)致分支匹配 不準(zhǔn)確等問題。相比于將所提取的血管與先驗(yàn)?zāi)Ｐ推ヅ涞姆椒?，這種方法的優(yōu)點(diǎn)在于節(jié)省 時(shí)間、標(biāo)記準(zhǔn)確，且能對(duì)更多的節(jié)段進(jìn)行識(shí)別和標(biāo)記。
[0009] 基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，該方法包括以下步驟： [0010]步驟1、圖像預(yù)處理:獲得冠狀動(dòng)脈分割圖像及中心線；
[0011]步驟2、剪枝操作；
[0012]步驟3、球心及球空間的定義；
[0013]步驟4、左、右冠狀動(dòng)脈的劃分規(guī)則；
[0014] 步驟5、左冠狀動(dòng)脈解剖學(xué)命名算法；
[0015] 步驟6、右冠狀動(dòng)脈解剖學(xué)命名算法。
[0016] 所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，所述步驟1具體執(zhí) 行以下步驟：
[0017] (1)從原始心臟連續(xù)斷層影像中分割出冠狀動(dòng)脈，再對(duì)分割圖像進(jìn)行血管中心線 (也稱為血管骨架)提??；
[0018] (2)將上述血管骨架轉(zhuǎn)換為由鄰接矩陣、節(jié)點(diǎn)和邊來描述的三維樹形結(jié)構(gòu)。其中節(jié) 點(diǎn)和邊的定義如下：
[0019] 葉節(jié)點(diǎn)：只和一個(gè)節(jié)點(diǎn)相連接的節(jié)點(diǎn)；
[0020] 連接節(jié)點(diǎn):和兩個(gè)節(jié)點(diǎn)相連接的節(jié)點(diǎn)；
[0021 ]分叉節(jié)點(diǎn):和兩個(gè)以上節(jié)點(diǎn)相連接的節(jié)點(diǎn)；
[0022] 游離邊:只有一端和節(jié)點(diǎn)相連接的邊；
[0023] 連接邊:兩端都和節(jié)點(diǎn)相連接的邊。
[0024] 所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，所述步驟2具體執(zhí) 行以下步驟：
[0025] 保留血管骨架中含有體素?cái)?shù)目最大的兩個(gè)連通塊，并將血管骨架中長度小于T的 游離邊剪枝置零(T的取值取決于想保留樹枝細(xì)節(jié)的程度，在18節(jié)段命名方式中T可取值為 15)，再將所得血管骨架采用上述方法重新轉(zhuǎn)換為三維樹形結(jié)構(gòu)。
[0026]所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，所述步驟3具體執(zhí) 行以下步驟：
[0027] (1)以確定的某點(diǎn)0(x0，y0，z0)為球心建立球坐標(biāo)系，將原笛卡爾坐標(biāo)空間中所有 點(diǎn)的坐標(biāo)值轉(zhuǎn)換為球坐標(biāo)值，得到球空間。其中將笛卡爾坐標(biāo)轉(zhuǎn)換為球坐標(biāo)的方法示意圖 如圖1，空間中任意一點(diǎn)P (X，y，Z)的球坐標(biāo)(Θ，私r)的計(jì)算方法為：
[0031] 其中Θ為方位角，是從以點(diǎn)〇為原點(diǎn)的X軸正半軸測得的點(diǎn)p的角位移弧度；
[0032] Φ為仰角，是從以點(diǎn)0為原點(diǎn)的X-Y平面測得的點(diǎn)P的角位移弧度；
[0033] r為半徑，是從點(diǎn)0到點(diǎn)P的距離。
[0034] 三個(gè)參數(shù)的取值范圍為:方位角θ = (-π，π]，仰角
半徑r= [0，+ 00 )。
[0035] (2)根據(jù)方位角Θ的取值方位，球空間中方位的定義如下表：
[0037] B卩：以確定的某點(diǎn)0為球心建立球坐標(biāo)系，空間中任意一點(diǎn)P的球坐標(biāo)(θ，φ:,Γ)。
[0038] 當(dāng)θ= (-π,-引時(shí)，視Ρ點(diǎn)位于〇點(diǎn)左前側(cè)，即區(qū)域^內(nèi)；
[0039] 當(dāng)Θ = (-f, 0]時(shí)，視Ρ點(diǎn)位于〇點(diǎn)左后側(cè)，即區(qū)域如內(nèi)；
[0040] 當(dāng)Θ = 時(shí)，視P點(diǎn)位于0點(diǎn)右后側(cè)，即區(qū)域A3內(nèi)；
[0041] 當(dāng)Θ = (|,知時(shí)，視P點(diǎn)位于〇點(diǎn)右前側(cè)，即區(qū)域A4內(nèi)。
[0042] 所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，所述步驟4具體執(zhí) 行以下步驟：
[0043] 根據(jù)左、右冠狀動(dòng)脈竇口兩點(diǎn)的坐標(biāo)，識(shí)別步驟2所得三維樹形結(jié)構(gòu)中左、右冠狀 動(dòng)脈連通塊，左冠狀動(dòng)脈主干與前降支、旋支、中間支(若存在)的分叉節(jié)點(diǎn)Nl和右冠狀動(dòng)脈 起始節(jié)點(diǎn)Nr。
[0044] 所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，所述步驟5具體執(zhí) 行以下步驟：
[0045] (1)計(jì)算三維樹形結(jié)構(gòu)左冠狀動(dòng)脈的游離邊和葉節(jié)點(diǎn)，并計(jì)算節(jié)點(diǎn)Nl到左冠狀動(dòng) 脈所有葉節(jié)點(diǎn)的最短路徑節(jié)點(diǎn)和路徑上的邊。
[0046] 其中最短路徑節(jié)點(diǎn)的計(jì)算，使用Di jkstra算法;路徑上的邊，通過計(jì)算路徑相鄰兩 節(jié)點(diǎn)共有的邊得到。
[0047] (2)計(jì)算所得所有游離邊、路徑上的邊相對(duì)于其連接的起始節(jié)點(diǎn)的球空間中值角 度。
[0048] 其中某邊的球空間中值角度的計(jì)算方法為:在以某邊連接的起始節(jié)點(diǎn)（近心端的 節(jié)點(diǎn)）為球心建立的球坐標(biāo)系下，計(jì)算邊上各點(diǎn)的球坐標(biāo)值并將方位角弧度值轉(zhuǎn)換為角度 值，再將這組各點(diǎn)方位角角度值取中值即為該邊的球空間中值角度。其中取中值的方法為： 對(duì)于項(xiàng)數(shù)為η的數(shù)組xl，x2,. . .，xn，將其進(jìn)行升序排序得到x(l)，x(2)，. . .，x(n);若η為奇 數(shù)，則數(shù)組的中值為
；若11為偶數(shù)，則數(shù)組的中值為
[0049] (3)左冠狀動(dòng)脈的具體分段：
[0050]查找步驟5(2)所得球空間中值角度中節(jié)點(diǎn)Nl所連各邊（包括游離邊和路徑上的 邊)的中值角度ki，i = 1，2，3，…，η，其中η為與節(jié)點(diǎn)Nl相連接的邊的數(shù)目。則各分支：
[0051 ] a)左冠狀動(dòng)脈主干(LM): 0Li>O，此時(shí)區(qū)域范圍有0Lie (A3 U A4)。
[0052] 若LM是游離邊，則所定義游離邊即為LM。若LM不是游離邊，查找其所在路徑Pu?，查 找Plm路徑中的最長路徑Plmi (若Plm只有一條路徑，則Plmi = Plm)和該路徑上的邊PLlmi，并將 Plmi最末端葉節(jié)點(diǎn)所連最長游離邊置于PLlmi末尾，所得Plmi和PLlmi所有點(diǎn)即為LM。
[0053] b)左前降支(LAD):0Li彡〇n0Li = 0min，此時(shí)0LieAi，其中0min為0Li中的最小值。
[0054] 所得LAD球空間中值角度的具體值為cilad。
[0055] 查找LAD所在路徑PLAD，查找PLAD路徑中的最長路徑P LAD1和路徑上的邊PLLAD1，并將 Pladi最末端葉節(jié)點(diǎn)所連最長游離邊置于PLladi末尾。查找Pladi上第一個(gè)分叉節(jié)點(diǎn)Nladi。則各 節(jié)段：
[0056] 前降支近段(pLAD):節(jié)點(diǎn)Nl到節(jié)點(diǎn)Nladi的邊。
[0057] 前降支中段(mLAD):節(jié)點(diǎn)Nladi到節(jié)點(diǎn)NmLAD的邊，其中NmLAD為PLladi除去pLAD后的中 點(diǎn)（不一定為步驟2所得三維樹形結(jié)構(gòu)中的節(jié)點(diǎn)）。計(jì)算該節(jié)段相對(duì)于節(jié)點(diǎn)Nladi的球空間中 值角度drnLAD。
[0058] 前降支遠(yuǎn)段(dLAD):節(jié)點(diǎn)NmLAD至IjPU皿結(jié)束的邊。
[0059]查找步驟5(2)所得球空間中值角度中節(jié)點(diǎn)Nl/^所連的非PLLAD1上的邊的中值角度 0LADli，i = 1，2，3，…，nl，其中nl為與節(jié)點(diǎn)NLAD1相連接的非PLLAD1上的邊的數(shù)目。則分支：
[0060]第一對(duì)角支(Dl ) : 0LADli^iO Π 9LADli>aLAD，此時(shí)區(qū)域范圍有0LADli G Al。
[0061 ]查找PLAD1上第二個(gè)分叉節(jié)點(diǎn)NLAD2,查找步驟5b所得球空間中值角度中節(jié)點(diǎn)N LAD2所 連的非PLLAD1上的邊的中值角度0LAD2i，i = 1，2，3，…，n2，其中n2為與節(jié)點(diǎn)NLAD2相連接的非 PUadi上的邊的數(shù)目。
[0062] 若上述D1未識(shí)別出，則分支：
[0063] D1 : 0LAD2i^iO Π 0LAD2i>CXLAD，此時(shí)區(qū)域范圍有 9LAD2i G Al。
[0064] 若上述D1已識(shí)別出，則分支：
[0065] 第二對(duì)角支(D2) : 9LAD2i^i0 Π 9LAD2i>amLAD，此時(shí)區(qū)域范圍有 9LAD2ie Al。
[0066] 采用上述方法，在PLAD1分叉節(jié)點(diǎn)上依次進(jìn)行識(shí)別直到D1、D2均識(shí)別出。若Pl/^分叉 節(jié)點(diǎn)較少，則記錄能識(shí)別出的D1、D2。
[0067] c)左旋支(LCx): 0Li彡0 n 0Li = Θ-max，此時(shí) 0Li e A2，其中 Θ-max為0Li 中非正值
[0068] 的最大值。所得LCx球空間中值角度的具體值為cilxx。
[0069] 查找LCx所在路徑PLCx，查找PLCx路徑中的最長路徑P LCxl和路徑上的邊PLLCxl，并將 pLCxl最末端葉節(jié)點(diǎn)所連最長游離邊置于PLLCxl末尾。查找pLCxl上第一個(gè)分叉節(jié)點(diǎn)，若該分叉 節(jié)點(diǎn)到節(jié)點(diǎn)Nl的邊的長度不大于閾值T(參考值80)，則該節(jié)點(diǎn)為N LCxl。查找步驟5(2)所得球 空間中值角度中節(jié)點(diǎn)Nlcxi所連各邊的中值角度θιχχιi，i = 1，2，3，…，η，其中η為與節(jié)點(diǎn)Nlcxi 相連接的邊的數(shù)目。則分支：
[0070] 旋支近段(pCx):節(jié)點(diǎn)Nl到節(jié)點(diǎn)Νιχχ?的邊。
[0071]旋支中遠(yuǎn)段（LCx): 0LCxli 彡 aLCx n ELCxii e PLlcxi，此時(shí)區(qū)域范圍有9LCxli e A3。
[0072] 第一鈍緣支（0M1): eLCxli < aLCx n ELCxli e plLCx1，此時(shí)區(qū)域范圍有 0LCxli e a2，其中 ElXxli為9lXxli所在的邊。
[0073] 若Plcx1存在第二個(gè)分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nl的邊的長度不大于閾值Τ(參考值80)，則該節(jié) 點(diǎn)為第二鈍緣支的起始節(jié)點(diǎn)Νμ 2。則分支：
[0074] 第二鈍緣支(0Μ2):節(jié)點(diǎn)Νομ2發(fā)出的非PLLCxl±的邊。
[0075]若Plcx1后續(xù)分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nl的邊的長度大于閾值T (參考值80)，則該節(jié)點(diǎn)為 NLCx2。查找步驟5(2)所得球空間中值角度中節(jié)點(diǎn)NLCx2所連各邊的中值角度0 LCx2i，i = 1，2， 3，…，n2，其中n2為與節(jié)點(diǎn)NLCx2相連接的邊的數(shù)目。則分支：
[0076] 旋支起源后降支(L-PDA):0LCx2i彡0,此時(shí)區(qū)域范圍有0心 2#(心1^2)。
[0077] 旋支起源后側(cè)支(L-PLB):0LCx2i>O,此時(shí)區(qū)域范圍有9 LCx2ie(A3UA4)。
[0078] 若Plcx1上第一個(gè)分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nl的邊的長度大于閾值T(參考值80)，則該節(jié)點(diǎn)為 NLCx2。查找步驟5(2)所得球空間中值角度中節(jié)點(diǎn)NLCx2所連各邊的中值角度0 LCx2i，i = 1，2， 3，…，n2，其中n2為與節(jié)點(diǎn)NLCx2相連接的邊的數(shù)目。則分支：
[0079] L-PDA: ΘLCx2i 彡 0，此時(shí)區(qū)域范圍有 9LCx2i e (A! U A2)。
[0080] L-PLB:0LCx2i>O,此時(shí)區(qū)域范圍有 9LCx2ie(A3UA4)。
[0081 ] pCx與LCx:節(jié)點(diǎn)Nl到節(jié)點(diǎn)Ν?Χχ2的邊。
[0082]以節(jié)點(diǎn)Nl為球心，計(jì)算節(jié)點(diǎn)Nl到節(jié)點(diǎn)Nlcx2的邊上所有點(diǎn)的球坐標(biāo)，得到該邊上各點(diǎn) 按連接順序的方位角數(shù)組θ2ι和半徑數(shù)組r2l，i = l，2,3,…，n，其中η為該邊上點(diǎn)的總數(shù)目。 則分段：
[0083] p CX:節(jié)點(diǎn)Nl到節(jié)點(diǎn)nlCx2m:in的邊，其中WLCX2min為r 2 i X s i η (Θ 2 i)中最小值所在的點(diǎn) (不一定為步驟2所得三維樹形結(jié)構(gòu)中的節(jié)點(diǎn)）。
[0084] LCx:節(jié)點(diǎn) NLCx2min 到節(jié)點(diǎn) NlcX2 的邊。
[0085] d)中間支(RI):0Li = 0Ln，此時(shí) θυεΑυΑ〗）。
[0086] 若RI是游離邊，則所定義游離邊即為RI。若RI不是游離邊，其下所有分支均記入 RI〇
[0087] 所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，所述步驟6具體執(zhí) 行以下步驟：
[0088] (1)計(jì)算三維樹形結(jié)構(gòu)右冠狀動(dòng)脈的游離邊和葉節(jié)點(diǎn)，并計(jì)算節(jié)點(diǎn)Nr到右冠狀動(dòng) 脈所有葉節(jié)點(diǎn)的最短路徑節(jié)點(diǎn)和路徑上的邊。方法與步驟5(1)相同。
[0089] (2)計(jì)算步驟6(1)所得所有游離邊、路徑上的邊相對(duì)于其連接的起始節(jié)點(diǎn)的球空 間中值角度。方法與步驟5(2)相同。
[0090] (3)右冠狀動(dòng)脈的具體分段：
[0091] a)查找步驟6 (1)所得路徑中的最長路徑PRCA和路徑上的邊PLRCA，并將PRCA最末端葉 節(jié)點(diǎn)所連最長游離邊置于PLRCA末尾。以節(jié)點(diǎn)Nr為球心，計(jì)算PLRCA上所有點(diǎn)的球坐標(biāo)，得到 PLrca上各點(diǎn)按連接順序的方位角數(shù)組0Ri和半徑數(shù)組rRi，i = 1，2，3，…，η，其中η為PLRCA上點(diǎn) 的總數(shù)目。則分段：
[0092] 右冠近段(pRCA):節(jié)點(diǎn)Nr到節(jié)點(diǎn)NmRCA的邊。
[0093] 右冠中段(mRCA):節(jié)點(diǎn)NmRCA到節(jié)點(diǎn)NRCAmin的邊。其中N RCAmiI^rRi X c 〇 s (Θ Ri)中最小 值所在的點(diǎn)（不一定為步驟2所得三維樹形結(jié)構(gòu)中的節(jié)點(diǎn))且區(qū)域范圍內(nèi)NRCA min e (A3 u a4)， NrnRCA為節(jié)點(diǎn)Nr到節(jié)點(diǎn)_NRCAmin的邊的中點(diǎn)。
[0094] b)若右冠沒有分支節(jié)點(diǎn)，或右冠主干上存在分支節(jié)點(diǎn)，且分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nr的邊 的長度不大于閾值T2 (參考值200)，則分段：
[0095] 右冠遠(yuǎn)段（dRCA):節(jié)點(diǎn)..到節(jié)點(diǎn)N Rc.amax.的邊，其中1^釀_為rRixcos(0Ri)中 最大值所在的點(diǎn)（不一定為步驟2所得三維樹形結(jié)構(gòu)中的節(jié)點(diǎn)）且區(qū)域范圍內(nèi)NRC~maxe (A^ U Aq).
[0096] 右冠起源后降支（R-PDA):節(jié)點(diǎn)NRCA_到PLRCA結(jié)束的邊。若N RCAma3^PLRCA結(jié)束的 點(diǎn)，則R-PDA不存在。
[0097] c)若右冠主干上存在分支節(jié)點(diǎn)，且分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nr的邊的長度大于閾值T2(參 考值200)，則該節(jié)點(diǎn)為Nrcm。查找步驟6(2)所得球空間中值角度中節(jié)點(diǎn)Nrcm所連的非PLrca上 的邊的中值角度9此虹，丨=1，2,3，'"，11，其中11為與節(jié)點(diǎn)詠11相連接的非?1^ [^上的邊的數(shù)目。 則分段：
[0098] dRCA:節(jié)點(diǎn) NRCAmin 到節(jié)點(diǎn) Nrcai 的邊。
[0099] R-PDA: 9RCAi < -90,此時(shí)區(qū)域范圍有 0RCAi G ^。
[0100] 右冠起源后側(cè)支(R-PLB):-9O〈0RCAi彡0,此時(shí)區(qū)域范圍有0 RCAieA2。
【附圖說明】
[0101 ]圖1為本發(fā)明方法的實(shí)施流程圖；
[0102] 圖2為將笛卡爾坐標(biāo)轉(zhuǎn)換為球坐標(biāo)的方法示意圖；
[0103] 圖3為左冠狀動(dòng)脈主要分支的分段流程圖；
[0104] 圖4為左冠狀動(dòng)脈前降支(LAD)的具體分段流程圖；
[0105]圖5為左冠狀動(dòng)脈旋支(LCx)的具體分段流程圖；
[0106] 圖6為右冠狀動(dòng)脈的分段流程圖。
【具體實(shí)施方式】
[0107] 按照"2014年SCCT指南"中提出的18分段標(biāo)準(zhǔn)對(duì)冠狀動(dòng)脈進(jìn)行解剖學(xué)命名，冠狀動(dòng) 脈18節(jié)段的名稱及其對(duì)應(yīng)的解剖描述如下表：

[0110] 基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，該方法實(shí)施流程圖如圖 1，具體實(shí)施流程如下：
[0111] 步驟1、圖像預(yù)處理:獲得冠狀動(dòng)脈分割圖像及中心線；
[0112] 步驟2、剪枝操作；
[0113] 步驟3、球心及球空間的定義；
[0114] 步驟4、左、右冠狀動(dòng)脈的劃分規(guī)則；
[0115] 步驟5、左冠狀動(dòng)脈解剖學(xué)命名算法；
[0116] 步驟6、右冠狀動(dòng)脈解剖學(xué)命名算法。
[0117] 所述步驟1獲得冠狀動(dòng)脈分割圖像及中心線，具體方法如下：
[0118] 首先從原始心臟連續(xù)斷層影像中分割出冠狀動(dòng)脈，并采用形態(tài)學(xué)操作和離散傅里 葉變換等方法對(duì)血管邊緣進(jìn)行平滑處理。接著通過一組3X3X3的匹配模板對(duì)提取出的冠 狀動(dòng)脈血管樹進(jìn)行細(xì)化和單體素化，通過分析體素點(diǎn)的鄰接關(guān)系對(duì)分叉點(diǎn)進(jìn)行標(biāo)記，提取 出符合血管樹實(shí)際情況的冠狀動(dòng)脈中心線（也稱血管骨架）。然后在此基礎(chǔ)上遍歷血管骨 架，同時(shí)構(gòu)造多叉樹，得到血管骨架拓?fù)浣Y(jié)構(gòu)的圖表示，即將血管骨架轉(zhuǎn)換為由鄰接矩陣、 節(jié)點(diǎn)和邊來描述的三維樹形結(jié)構(gòu)。
[0119] 其中節(jié)點(diǎn)和邊的定義如下：
[0120] 葉節(jié)點(diǎn)：只和一個(gè)節(jié)點(diǎn)相連接的節(jié)點(diǎn)；
[0121 ]連接節(jié)點(diǎn):和兩個(gè)節(jié)點(diǎn)相連接的節(jié)點(diǎn)；
[0122] 分叉節(jié)點(diǎn)：和兩個(gè)以上節(jié)點(diǎn)相連接的節(jié)點(diǎn)；
[0123] 游離邊:只有一端和節(jié)點(diǎn)相連接的邊；
[0124] 連接邊:兩端都和節(jié)點(diǎn)相連接的邊。
[0125] 所述步驟2剪枝操作，具體方法步驟如下：
[0126] 步驟2a、判斷步驟1所得血管骨架的連通情況。若血管骨架連通塊的數(shù)目大于2,則 保留含有體素?cái)?shù)目最大的兩個(gè)連通塊，其他連通塊全部置零。所保留的兩個(gè)連通塊即是左、 右冠狀動(dòng)脈血管骨架。
[0127] 步驟2b、計(jì)算三維樹形結(jié)構(gòu)的游離邊，并將步驟2a所得血管骨架中長度小于T的游 離邊置零(T的取值取決于想保留樹枝細(xì)節(jié)的程度，在18節(jié)段命名方式中T可取值為15)。此 處應(yīng)注意游離邊置零后需將該游離邊所連接的節(jié)點(diǎn)重新置一，以避免血管骨架不連續(xù)。
[0128] 步驟2c、與步驟1采用相同方法，將上述處理后的血管骨架重新轉(zhuǎn)換為三維樹形結(jié) 構(gòu)。
[0129] 所述步驟3球心及球空間的定義，具體方法步驟如下：
[0130]步驟3a、以確定的某點(diǎn)0(x0，y0，zO)為球心建立球坐標(biāo)系，將原笛卡爾坐標(biāo)空間中 所有點(diǎn)的坐標(biāo)值轉(zhuǎn)換為球坐標(biāo)值，得到球空間。其中將笛卡爾坐標(biāo)轉(zhuǎn)換為球坐標(biāo)的方法示 意圖如圖2，空間中任意一點(diǎn)P (X，y，z)的球坐標(biāo)(θ，φr)的計(jì)算方法為：
[0134] 其中Θ為方位角，是從以點(diǎn)〇為原點(diǎn)的X軸正半軸測得的點(diǎn)P的角位移弧度；
[0135] Φ為仰角，是從以點(diǎn)0為原點(diǎn)的X-Y平面測得的點(diǎn)P的角位移弧度；
[0136] r為半徑，是從點(diǎn)0到點(diǎn)P的距離。
[0137] 三個(gè)參數(shù)的取值范圍為:方位角θ = (-π，π]，仰角
（半徑r= [0，+ 00)。
[0138] 步驟3b、根據(jù)方位角Θ的取值方位，球空間中方位的定義如下表：

[0141 ] B卩：以確定的某點(diǎn)0為球心建立球坐標(biāo)系，空間中任意一點(diǎn)P的球坐標(biāo)(Θ，r)。
[0142]當(dāng)
時(shí)，視P點(diǎn)位于0點(diǎn)左前側(cè)，即區(qū)域^內(nèi)；
[0143] 當(dāng)
時(shí)，視P點(diǎn)位于〇點(diǎn)左后側(cè)，即區(qū)域如內(nèi)；
[0144] 當(dāng)
|時(shí)，視P點(diǎn)位于0點(diǎn)右后側(cè)，即區(qū)域A3內(nèi)；
[0145] 當(dāng)
I時(shí)，視P點(diǎn)位于0點(diǎn)右前側(cè)，即區(qū)域A4內(nèi)。
[0146] 所述步驟4左、右冠狀動(dòng)脈的劃分規(guī)則，具體方法步驟如下：
[0147] 步驟4a、采用步驟3所述方法，以血管骨架上左、右冠狀動(dòng)脈竇口兩點(diǎn)的中點(diǎn)為球 心建立球坐標(biāo)系，將步驟2處理后的三維樹形結(jié)構(gòu)中各節(jié)點(diǎn)的笛卡爾坐標(biāo)值轉(zhuǎn)換為球坐標(biāo) 值。
[0148] 步驟4b、分別尋找三維樹形結(jié)構(gòu)兩個(gè)連通塊中球坐標(biāo)值半徑最小的節(jié)點(diǎn)。
[0149] 步驟4c、比較步驟4b所得兩個(gè)節(jié)點(diǎn)的方位角。方位角為負(fù)的節(jié)點(diǎn)，即該節(jié)點(diǎn)位于此 球空間區(qū)域六:或如，亦即該節(jié)點(diǎn)位于球心左側(cè)，其所在的連通塊為左冠狀動(dòng)脈;該節(jié)點(diǎn)即為 左冠狀動(dòng)脈主干與前降支、旋支、中間支(若存在）的分叉節(jié)點(diǎn)Nl。方位角為正的節(jié)點(diǎn)，即該 節(jié)點(diǎn)位于此球空間區(qū)域A 3或A4,亦即該節(jié)點(diǎn)位于球心右側(cè)，其所在的連通塊為右冠狀動(dòng)脈； 該節(jié)點(diǎn)即為右冠狀動(dòng)脈起始節(jié)點(diǎn)Nr。
[0150] 所述步驟5左冠狀動(dòng)脈解剖學(xué)命名算法，具體方法步驟如下：
[0151 ] 步驟5a、計(jì)算三維樹形結(jié)構(gòu)左冠狀動(dòng)脈的游離邊和葉節(jié)點(diǎn)，并計(jì)算節(jié)點(diǎn)Nl到左冠 狀動(dòng)脈所有葉節(jié)點(diǎn)的最短路徑節(jié)點(diǎn)和路徑上的邊。
[0152]其中最短路徑節(jié)點(diǎn)的計(jì)算，使用Di jkstra算法;路徑上的邊，通過計(jì)算路徑相鄰兩 節(jié)點(diǎn)共有的邊得到。路徑上的邊一定為連接邊，且所有路徑上的邊與所有連接邊等同。 [0153]步驟5b、計(jì)算步驟5a所得所有游離邊、路徑上的邊相對(duì)于其連接的起始節(jié)點(diǎn)的球 空間中值角度。
[0154]其中某邊的球空間中值角度的計(jì)算方法為:在以某邊連接的起始節(jié)點(diǎn)（近心端的 節(jié)點(diǎn)）為球心建立的球坐標(biāo)系下，計(jì)算邊上各點(diǎn)的球坐標(biāo)值并將方位角弧度值轉(zhuǎn)換為角度 值，再將這組各點(diǎn)方位角角度值取中值即為該邊的球空間中值角度。即，某邊的球空間中值 角度實(shí)質(zhì)上是該邊的大致方位角。其中取中值的方法為:對(duì)于項(xiàng)數(shù)為η的數(shù)組xl，x2, ...， xn，將其進(jìn)行升序排序得到奴1)，以2)，...，以1〇;若11為奇數(shù)，則數(shù)組的中值為1
若η為偶數(shù)，則數(shù)組的中值為
[0155] 步驟5c、左冠狀動(dòng)脈主要分支的分段。流程圖如圖3。
[0156] 由解剖學(xué)先驗(yàn)知識(shí)可知，橫斷面上相對(duì)于節(jié)點(diǎn)Nl，左主干(LM)自右側(cè)走形、左前降 支(LAD)向左前側(cè)走形、左旋支(LCx)向左后側(cè)走形、中間支(RI)走形于LAD與LCx之間。
[0157] 查找步驟5b所得球空間中值角度中節(jié)點(diǎn)Nl所連各邊(包括游離邊和路徑上的邊） 的中值角度ki，i = 1，2，3，…，n，其中η為與節(jié)點(diǎn)見相連接的邊的數(shù)目。則各分支：
[0158] LM:0Li>O,此時(shí)區(qū)域范圍有 0Lie(A3UA4)。
[0159] 1^^〇門0^ = 0_，此時(shí)0^￡厶1，其中0_為叱中的最小值。所得1^球空間中 值角度的具體值為aiAD。
[0160] LCx:0Li<〇n0Li = 0-max，此時(shí)0LieA2,其中Θ-max為0 Li中非正值的最大值。所得LCX球 空間中值角度的具體值為αιχχ。
[0161 ] RI: 0Li = 0Ln，此時(shí)0Li e (A! U Α2)。
[0162] 步驟5d、LM和RI的處理。
[0163] 由于LM、RI下沒有其他定義分段，若步驟5c中定義的LM、RI是游離邊，則所定義游 離邊即為LM、RI。
[0164] 若步驟5c中定義的LM、RI不是游離邊，則步驟5c中定義的LM，查找其所在路徑Pu?， 查找Plm路徑中的最長路徑Plmi (若Plm只有一條路徑，則Plmi = Plm)和該路徑上的邊PLlmi，并將 Plmi最末端葉節(jié)點(diǎn)所連最長游離邊置于PLlmi末尾，所得Plmi和PLlmi所有點(diǎn)即為LM;步驟5c中 定義的RI，其下所有分支均記入RI。
[0165] 步驟5e、LAD的具體分段。流程圖如圖4。
[0166] (1)由"2014年SCCT指南"18分段標(biāo)準(zhǔn)可知，前降支近段(pLAD)為LM末端至第一對(duì) 角支或第一間隔支分叉節(jié)點(diǎn)Nu D1，此分叉節(jié)點(diǎn)至心尖部(LAD末端)按長度二等分，依次分為 前降支中段(mLAD)和前降支遠(yuǎn)段(dLAD)。
[0167] 查找步驟5c中定義的LAD所在路徑PLAD，查找PLAD路徑中的最長路徑P LADdP路徑上 的邊PLladi，并將Pladi最末端葉節(jié)點(diǎn)所連最長游離邊置于PLladi末尾。查找Pladi上第一個(gè)分叉 節(jié)點(diǎn)NlADl。則各節(jié)段：
[0168] pLAD:節(jié)點(diǎn)Nl到節(jié)點(diǎn)Nladi的邊。
[0169] mLAD:節(jié)點(diǎn)Nladi到節(jié)點(diǎn)NmLAD的邊，其中NmLAD為PLladi除去pLAD后的中點(diǎn)（不一定為步 驟2所得三維樹形結(jié)構(gòu)中的節(jié)點(diǎn)）。計(jì)算該節(jié)段相對(duì)于節(jié)點(diǎn)Nladi的球空間中值角度amLAD。
[0170] dLAD:節(jié)點(diǎn) NmLAD 到PLladi 結(jié)束的邊。
[0171] (2)由解剖學(xué)先驗(yàn)知識(shí)可知，橫斷面上第一對(duì)角支(D1)、第二對(duì)角支(D2)均走形于 LAD左側(cè)。
[0172] 查找步驟5b所得球空間中值角度中節(jié)點(diǎn)Nl/^所連的非PLLAD1上的邊的中值角度 0LADli，i = 1，2，3，…，nl，其中nl為與節(jié)點(diǎn)NLAD1相連接的非PLLAD1上的邊的數(shù)目。則分支：
[0173] D1 : 0LADli^iO Π 0LADli>aLAD，此時(shí)區(qū)域范圍有 0LADli G Al。
[0174] (3)查找Pl皿上第二個(gè)分叉節(jié)點(diǎn)NLAD2,查找步驟5b所得球空間中值角度中節(jié)點(diǎn)N LAD2 所連的非PUAD1上的邊的中值角度0LAD2i，i = 1，2，3，…，n2，其中n2為與節(jié)點(diǎn)NLAD2相連接的非 PLladi上的邊的數(shù)目。
[0175] 若上述D1未識(shí)別出，則分支：
[0176] D1 : 0LAD2i^iO Π 0LAD2i>aLAD，此時(shí)區(qū)域范圍有 0LAD2i G Al。
[0177] 若上述D1已識(shí)別出，則分支：
[0178] D2 : 0LAD2i^iO Π 0LAD2i>amLAD，此時(shí)區(qū)域范圍有 0LAD2ie Al。
[0179] (4)采用上述(3)中方法，在PLAD1*叉節(jié)點(diǎn)上依次進(jìn)行識(shí)別直到D1、D2均識(shí)別出。若 Pladi分叉節(jié)點(diǎn)較少，則記錄能識(shí)別出的D1、D2。
[0180] 步驟5f、LCx的具體分段。流程圖如圖5。
[0181] (1)由"2014年SCCT指南"18分段標(biāo)準(zhǔn)可知，旋支近段(pCx)為LM末端至第一鈍緣支 (0M1)分叉節(jié)點(diǎn)N LCxl。由解剖學(xué)先驗(yàn)知識(shí)可知，橫斷面上相對(duì)于分叉節(jié)點(diǎn)NLCxl，0M1走形于LCx 左側(cè)、旋支中遠(yuǎn)段(LCx)向右后側(cè)走形;橫斷面上相對(duì)于LCx與旋支起源后降支(L-PDA)的分 叉節(jié)點(diǎn)N LCx2，L-PDA向左側(cè)走形、旋支起源后側(cè)支(L-PLB)向右側(cè)走形。
[0182] 查找步驟5c中定義的LCx所在路徑PLCx，查找PLCx路徑中的最長路徑P LCxdP路徑上 的邊ΡΙχω，并將P^xi最末端葉節(jié)點(diǎn)所連最長游離邊置于PU&d末尾。
[0183] (2)查找PLCxl上第一個(gè)分叉節(jié)點(diǎn)，若該分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nl的邊的長度不大于閾值T (參考值80)，則該節(jié)點(diǎn)為N LCxl。查找步驟5b所得球空間中值角度中節(jié)點(diǎn)NLCxl*連各邊的中值 角度9ix xii，i = l，2,3,…，n，其中η為與節(jié)點(diǎn)Νιχχ?相連接的邊的數(shù)目。貝分支：
[0184] pCx:節(jié)點(diǎn)Nl到節(jié)點(diǎn)Νιχχ?的邊。
[0185] LCx: 9[Xxii 多 αιχχ Π Eixxii G PLlcxI，此時(shí)區(qū)域范圍有 SiXxii G A3 0
[0186] 0Μ1: θιχχ? < ?lcx η Eixxii ? PLLCx1， 此時(shí)區(qū)域范圍有9LCxli ε Α2,其中ELCxli為9LCxli所 在的邊。
[0187] (3)若PLCxl存在第二個(gè)分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nl的邊的長度不大于閾值T(參考值80)，則 該節(jié)點(diǎn)為第二鈍緣支(0Μ2)的起始節(jié)點(diǎn)Νμ 2。則分支：
[0188] 0Μ2:節(jié)點(diǎn)Νομ2發(fā)出的非PLlcxi上的邊。
[0189]若PLCxl后續(xù)分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nl的邊的長度大于閾值T(參考值80)，則該節(jié)點(diǎn)為 NLCx2。查找步驟5b所得球空間中值角度中節(jié)點(diǎn)NLCx2所連各邊的中值角度0LCx2i，i = 1，2， 3，…，n2，其中n2為與節(jié)點(diǎn)NLCx2相連接的邊的數(shù)目。則分支：
[0190] L-PDA: 0LCx2i彡0，此時(shí)區(qū)域范圍有9LCx2i e (A! U A2)。
[0191 ] L-PLB: 0LCx2i>O，此時(shí)區(qū)域范圍有0LCx2i e (A3 U A4)。
[0192] (4)若PLCxl上第一個(gè)分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nl的邊的長度大于閾值T(參考值80)，則該節(jié) 點(diǎn)為N LCx2。查找步驟5b所得球空間中值角度中節(jié)點(diǎn)NLCx2所連各邊的中值角度0LCx2i，i = 1，2， 3，…，n2，其中n2為與節(jié)點(diǎn)NLCx2相連接的邊的數(shù)目。則分支：
[0193] L-PDA: 0LCx2i彡0，此時(shí)區(qū)域范圍有9LCx2i e (A! U A2)。
[0194] L-PLB: 0LCx2i>O，此時(shí)區(qū)域范圍有9LCx2i e (A3 U A4)。
[0195] pCx與LCx:節(jié)點(diǎn)Nl到節(jié)點(diǎn)Nlcx2的邊。
[0196]以節(jié)點(diǎn)Nl為球心，計(jì)算節(jié)點(diǎn)Nl到節(jié)點(diǎn)NlxX2的邊上所有點(diǎn)的球坐標(biāo)，得到該邊上各點(diǎn) 按連接順序的方位角數(shù)組θ2ι和半徑數(shù)組r2l，i = l，2,3,…，n，其中η為該邊上點(diǎn)的總數(shù)目。 則分段：
[0197] pCx:節(jié)點(diǎn)Nl到節(jié)點(diǎn)NLCX2min.的邊，其中N Lcx2mto..為r2i X s iη (Θ2i)中最小值所在的點(diǎn) (不一定為步驟2所得三維樹形結(jié)構(gòu)中的節(jié)點(diǎn)）。
[0198] LCX:節(jié)點(diǎn) NLCx2min 到節(jié)點(diǎn) NlcX2 的邊。
[0199] 所述步驟6右冠狀動(dòng)脈解剖學(xué)命名算法，具體方法步驟如下：
[0200] 步驟6a、計(jì)算三維樹形結(jié)構(gòu)右冠狀動(dòng)脈的游離邊和葉節(jié)點(diǎn)，并計(jì)算節(jié)點(diǎn)Nr到右冠 狀動(dòng)脈所有葉節(jié)點(diǎn)的最短路徑節(jié)點(diǎn)和路徑上的邊。方法與步驟5a相同。
[0201 ]步驟6b、計(jì)算步驟6a所得所有游離邊、路徑上的邊相對(duì)于其連接的起始節(jié)點(diǎn)的球 空間中值角度。方法與步驟5b相同。
[0202]步驟6c、右冠狀動(dòng)脈的具體分段。流程圖如圖6。
[0203] (1)由"2014年SCCT指南" 18分段標(biāo)準(zhǔn)可知，右冠近段(pRCA)為右冠起始至心臟銳 緣處的一半，右冠中段(mRCA)為pRCA末端至心臟銳緣，右冠遠(yuǎn)段(dRCA)為mRCA末端至右冠 起源后降支(R-PDA)起始，即dRCA與R-TOA、右冠起源后側(cè)支(R-PLB)的分叉節(jié)點(diǎn)NRCA1。由解 剖學(xué)先驗(yàn)知識(shí)可知，橫斷面上相對(duì)于分叉節(jié)點(diǎn)N RCA1，R-PDA向左前側(cè)走形、R-PLB向左后側(cè)走 形。
[0204]查找步驟6a所得路徑中的最長路徑PRCA和路徑上的邊PLRCA，并將P RCA最末端葉節(jié)點(diǎn) 所連最長游離邊置于PLRCA末尾。以節(jié)點(diǎn)Nr為球心，計(jì)算PLRCA上所有點(diǎn)的球坐標(biāo)，得到PLrca上 各點(diǎn)按連接順序的方位角數(shù)組9 Rl和半徑數(shù)組rRl，i = l，2,3,…，n，其中η為PLRCA上點(diǎn)的總數(shù) 目。則分段：
[0205] pRCA:節(jié)點(diǎn)Nr到節(jié)點(diǎn)NmRCA的邊。
[0206] mRCA:節(jié)點(diǎn)NmRCA到節(jié)點(diǎn)NRCAmin的邊。其中NRCAmin為 rRi X cos(0Ri)中最小值所在的點(diǎn) (不一定為步驟2所得三維樹形結(jié)構(gòu)中的節(jié)點(diǎn))且區(qū)域范圍內(nèi)N RCA_ e (A3 u A4)，NmRCA為節(jié)點(diǎn) Nr到節(jié)點(diǎn)NRCAmin的邊的中點(diǎn)。
[0207] (2)若右冠沒有分支節(jié)點(diǎn)，或右冠主干上存在分支節(jié)點(diǎn)，且分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nr的邊 的長度不大于閾值T2 (參考值200)，則分段：
[0208] dRCA:節(jié)點(diǎn) NRCAmiDL到節(jié)點(diǎn)N.R_CAm. ax的邊，其中 NRCAmaj^rRi X cos(0Ri)中最大值所在 的點(diǎn)(不一定為步驟2所得三維樹形結(jié)構(gòu)中的節(jié)點(diǎn))且區(qū)域范圍內(nèi)Nrca_ ε (A2 U As〕:。
[0209] R-PDA:節(jié)點(diǎn)NRCAmax到PLRCA結(jié)束的邊。若N RCAma3^PLRCA結(jié)束的點(diǎn)，則R-PDA不存在。
[0210] (3)若右冠主干上存在分支節(jié)點(diǎn)，且分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nr的邊的長度大于閾值T2(參 考值200)，則該節(jié)點(diǎn)為Nrcm。查找步驟6b所得球空間中值角度中節(jié)點(diǎn)Nrcm所連的非PLrca上的 邊的中值角度Grcm，i = 1，2，3，…，η，其中η為與節(jié)點(diǎn)NRCA1相連接的非PLrca上的邊的數(shù)目。則 分段：
[0211 ] dRCA:節(jié)點(diǎn) NRCAmm|^節(jié)點(diǎn)Nrcai 的邊。
[0212 ] R-PDA: 0RCAi < -90，此時(shí)區(qū)域范圍有 0RCAi e A!。
[0213] R-PLB:-9O〈0RCAi<〇,此時(shí)區(qū)域范圍有 0RCAieA2。
[0214] 本發(fā)明的有益效果：
[0215] 1.本發(fā)明對(duì)冠狀動(dòng)脈樹中各個(gè)分支節(jié)段的空間幾何位置及其毗鄰關(guān)系在球空間 中直接進(jìn)行定義和劃分，有利于準(zhǔn)確的對(duì)各節(jié)段進(jìn)行解剖學(xué)命名標(biāo)記;解決現(xiàn)有的冠狀動(dòng) 脈自動(dòng)分段與解剖學(xué)標(biāo)記方法計(jì)算時(shí)間較長，算法復(fù)雜度較高，以及無法窮舉冠狀動(dòng)脈分 布類型導(dǎo)致分支匹配不準(zhǔn)確等問題；
[0216] 2.易于擴(kuò)展:人體冠脈的個(gè)體差異比較大，本發(fā)明直接對(duì)分布空間定義不需要窮 舉各種類型，相比于將所提取的血管與先驗(yàn)?zāi)Ｐ推ヅ涞姆椒?，這種方法的優(yōu)點(diǎn)在于節(jié)省時(shí) 間、標(biāo)記準(zhǔn)確，且能對(duì)更多的節(jié)段進(jìn)行識(shí)別和標(biāo)記；
[0217] 3.便于規(guī)范化影像診斷報(bào)告和回歸性分析;本發(fā)明的方法是對(duì)空間冠脈樹進(jìn)行解 剖學(xué)標(biāo)記，其標(biāo)記信息是帶有空間三維坐標(biāo)信息的，同時(shí)可與冠脈量化分析信息中其他數(shù) 值(如，管腔半徑，斑塊體積，鈣化和非鈣化的比例等）的空間坐標(biāo)結(jié)合起來，集成在數(shù)據(jù)結(jié) 構(gòu)中的，在醫(yī)生出診斷報(bào)告時(shí)，可以直接提取CCTA指南規(guī)定參數(shù)出診斷報(bào)告，不僅僅規(guī)范了 報(bào)告，同時(shí)便于局部病灶診斷數(shù)據(jù)的回顧性分析；
[0218] 4.有利于動(dòng)態(tài)量化評(píng)估:可根據(jù)解剖標(biāo)記對(duì)PCI手術(shù)后支架植入部分進(jìn)行聚焦，結(jié) 合動(dòng)態(tài)影像技術(shù)，可評(píng)估術(shù)后相應(yīng)節(jié)段的血流暢通情況及其在心臟周期性運(yùn)動(dòng)中的適應(yīng)， 可對(duì)術(shù)后再狹窄等情況進(jìn)行定期的無創(chuàng)監(jiān)測。應(yīng)當(dāng)理解的是，本方法在具體應(yīng)用中并不局 限于"2014年SCCT指南"中提出的18分段標(biāo)準(zhǔn)，凡涉及三維冠狀動(dòng)脈模型解剖學(xué)自動(dòng)分段和 解剖學(xué)命名的應(yīng)用，本發(fā)明均適用。本方法在具體應(yīng)用中并不局限于CTA這種圖像數(shù)據(jù)類 型，凡涉及連續(xù)斷層的圖像該方法均適用。對(duì)本領(lǐng)域普通技術(shù)人員來說，可以根據(jù)上述說明 加以改進(jìn)或變換，而所有這些改進(jìn)和變換都應(yīng)屬于本發(fā)明所附權(quán)利要求的保護(hù)范圍。
【主權(quán)項(xiàng)】
1. 基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，其特征在于，包括W下步 驟： 步驟1、圖像預(yù)處理:獲得冠狀動(dòng)脈分割圖像及中屯、線； 步驟2、剪枝操作； 步驟3、球屯、及球空間的定義； 步驟4、左、右冠狀動(dòng)脈的劃分規(guī)則； 步驟5、左冠狀動(dòng)脈解剖學(xué)命名算法； 步驟6、右冠狀動(dòng)脈解剖學(xué)命名算法。2. 根據(jù)權(quán)利要求1所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，其特 征在于，所述步驟1具體執(zhí)行W下步驟： 從原始屯、臟連續(xù)斷層影像中分割出冠狀動(dòng)脈，再對(duì)分割圖像進(jìn)行血管中屯、線（也稱為 血管骨架)提取;將上述血管骨架轉(zhuǎn)換為由鄰接矩陣、節(jié)點(diǎn)和邊來描述的Ξ維樹形結(jié)構(gòu)。3. 根據(jù)權(quán)利要求1所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，其特 征在于，所述步驟2具體執(zhí)行W下步驟： 保留血管骨架中含有體素?cái)?shù)目最大的兩個(gè)連通塊，并將血管骨架中長度小于T的游離 邊剪枝置零，再將所得血管骨架采用上述方法重新轉(zhuǎn)換為Ξ維樹形結(jié)構(gòu)。4. 根據(jù)權(quán)利要求1所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，其特 征在于，所述步驟3具體執(zhí)行W下步驟： (1似確定的某點(diǎn)〇(x〇，y〇，zO)為球屯、建立球坐標(biāo)系，將原笛卡爾坐標(biāo)空間中所有點(diǎn)的 坐標(biāo)值轉(zhuǎn)換為球坐標(biāo)值，得到球空間；其中將笛卡爾坐標(biāo)轉(zhuǎn)換為球坐標(biāo)的方法示意圖如圖 l，空間中任意一點(diǎn)P(χ，y，z)的球坐標(biāo)(θ，φ,r)的計(jì)算方法為：其中Θ為方位角，是從W點(diǎn)0為原點(diǎn)的X軸正半軸測得的點(diǎn)P的角位移弧度； 為仰角，是從從點(diǎn)0為原點(diǎn)的X-Y平面測得的點(diǎn)P的角位移弧度； r為半徑，是從點(diǎn)0到點(diǎn)P的距離； Ξ個(gè)參數(shù)的取值范圍為:方位角0 = (-JT,ji]，仰角φ = ,半徑r=[0，+°°); (2)根據(jù)方位角Θ的取值方位，球空間中方位的定義如下表：良確定的某點(diǎn)ο為球屯、建立球坐標(biāo)系，空間中任意一點(diǎn)P的球坐標(biāo)(目I奪,r); 當(dāng)白二(-π, - f|時(shí)，視P點(diǎn)位于0點(diǎn)左前側(cè)，即區(qū)域Ai內(nèi)； 當(dāng)β = 別時(shí)，視P點(diǎn)位于0點(diǎn)左后側(cè)，即區(qū)域A2內(nèi)； 當(dāng)曰二批fl時(shí)，視Ρ點(diǎn)位于0點(diǎn)右后側(cè)，即區(qū)域A3內(nèi)； 當(dāng)曰=β，π]時(shí)，視P點(diǎn)位于0點(diǎn)右前側(cè)，即區(qū)域A4內(nèi)。5. 根據(jù)權(quán)利要求1所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，其特 征在于，所述步驟4具體執(zhí)行W下步驟： 根據(jù)左、右冠狀動(dòng)脈竇口兩點(diǎn)的坐標(biāo)，識(shí)別步驟2所得Ξ維樹形結(jié)構(gòu)中左、右冠狀動(dòng)脈 連通塊，左冠狀動(dòng)脈主干與前降支、旋支、中間支(若存在)的分叉節(jié)點(diǎn)化和右冠狀動(dòng)脈起始 節(jié)點(diǎn)Nr。6. 根據(jù)權(quán)利要求1所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，其特 征在于，所述步驟5具體執(zhí)行W下步驟： (1) 計(jì)算Ξ維樹形結(jié)構(gòu)左冠狀動(dòng)脈的游離邊和葉節(jié)點(diǎn)，并計(jì)算節(jié)點(diǎn)化到左冠狀動(dòng)脈所有 葉節(jié)點(diǎn)的最短路徑節(jié)點(diǎn)和路徑上的邊； 其中最短路徑節(jié)點(diǎn)的計(jì)算，使用Dijkstra算法;路徑上的邊，通過計(jì)算路徑相鄰兩節(jié)點(diǎn) 共有的邊得到； (2) 計(jì)算所得所有游離邊、路徑上的邊相對(duì)于其連接的起始節(jié)點(diǎn)的球空間中值角度； 其中某邊的球空間中值角度的計(jì)算方法為:在W某邊連接的起始節(jié)點(diǎn)（近屯、端的節(jié)點(diǎn)） 為球屯、建立的球坐標(biāo)系下，計(jì)算邊上各點(diǎn)的球坐標(biāo)值并將方位角弧度值轉(zhuǎn)換為角度值，再 將運(yùn)組各點(diǎn)方位角角度值取中值即為該邊的球空間中值角度;其中取中值的方法為:對(duì)于 項(xiàng)數(shù)為η的數(shù)組xl，x2,. . .，xn，將其進(jìn)行升序排序得到x(l)，x(2)，. . .，x(n);若η為奇數(shù)，貝1J 數(shù)組的中值為m = 若η為偶數(shù)，則數(shù)組的中值女(3) 左冠狀動(dòng)脈的具體分段： 查找步驟5(2)所得球空間中值角度中節(jié)點(diǎn)化所連各邊(包括游離邊和路徑上的邊)的中 值角度目Li，i = 1，2，3，…，η，其中η為與節(jié)點(diǎn)化相連接的邊的數(shù)目；則各分支： a)左冠狀動(dòng)脈主干(LM):目Li〉0,此時(shí)區(qū)域范圍有目Lie (A3UA4); 若LM是游離邊，則所定義游離邊即為LM;若LM不是游離邊，查找其所在路徑Plm，查找Plm 路徑中的最長路徑Plmi(若Plm只有一條路徑，貝化lmi = Plm)和該路徑上的邊化LM1，并將Plmi最 末端葉節(jié)點(diǎn)所連最長游離邊置于PLlmi末尾，所得Plmi和PLlmi所有點(diǎn)即為LM; b)左前降支化AD):目Li《0 η 0Li=目min,此時(shí)目Lie Al，其中目min為目Li中的最小值;所得LAD 球空間中值角度的具體值為QLAD ; 查找LAD所在路徑Plad，查找Plad路徑中的最長路徑Pl皿和路徑上的邊化LAD1，并將Pl廬最 末端葉節(jié)點(diǎn)所連最長游離邊置于PLladi末尾;查找Pladi上第一個(gè)分叉節(jié)點(diǎn)Nladi ;則各節(jié)段： 前降支近段(pLAD):節(jié)點(diǎn)化到節(jié)點(diǎn)Nladi的邊； 前降支中段(mLAD):節(jié)點(diǎn)化AD1到節(jié)點(diǎn)NmLAD的邊，其中NmLAD為化LAD1除去pLAD后的中點(diǎn)（不 一定為步驟2所得Ξ維樹形結(jié)構(gòu)中的節(jié)點(diǎn)）；計(jì)算該節(jié)段相對(duì)于節(jié)點(diǎn)Nladi的球空間中值角度 日 mLAD; 前降支遠(yuǎn)段(化AD):節(jié)點(diǎn)NmLAD至帆XAD儲(chǔ)束的邊； 查找步驟5(2)所得球空間中值角度中節(jié)點(diǎn)Nladi所連的非PU皿上的邊的中值角度0LADii， i = 1，2，3，…，nl，其中nl為與節(jié)點(diǎn)Nl\di相連接的非PU皿上的邊的數(shù)目；則分支： 第一對(duì)角支化1):目LADli《0 Π 0LADli〉aLAD，此時(shí)區(qū)域范圍有目LADliGAl; 查找Pladi上第二個(gè)分叉節(jié)點(diǎn)化AD2，查找步驟化所得球空間中值角度中節(jié)點(diǎn)Nlad2所連的 非PLladi上的邊的中值角度目LAD2i，i = 1，2，3，…，n2，其中n2為與節(jié)點(diǎn)NLAD2相連接的非PLladi上 的邊的數(shù)目； 若上述D1未識(shí)別出，則分支： D1 :目LAD2i《0 η 0LAD2i〉aLAD，此時(shí)區(qū)域范圍有0LAD2ieAl; 若上述D1已識(shí)別出，則分支： 第二對(duì)角支化2 ) : 0LAD2i《O η 0LAD2i〉amLAD，此時(shí)區(qū)域范圍有目LAD2i e Al; 采用上述方法，在Pladi分叉節(jié)點(diǎn)上依次進(jìn)行識(shí)別直到D1、D2均識(shí)別出；若Pladi分叉節(jié)點(diǎn) 較少，則記錄能識(shí)別出的D1、D2; C)左旋支化Cx):目Li《〇n0Li =目-max,此時(shí)目LiGA2,其中目-max為目Li中非正值的最大值;所 得LCX球空間中值角度的具體值為QlCx ; 查找LCx所在路徑Plcx，查找Plcx路徑中的最長路徑PlcxI和路徑上的邊化LCxl，并將PlcxI最 末端葉節(jié)點(diǎn)所連最長游離邊置于化LGxl末尾;查找PlGxI上第一個(gè)分叉節(jié)點(diǎn)，若該分叉節(jié)點(diǎn)到 節(jié)點(diǎn)化的邊的長度不大于闊值Τ(參考值80)，則該節(jié)點(diǎn)為化Cxi;查找步驟5(2)所得球空間中 值角度中節(jié)點(diǎn)化cxi所連各邊的中值角度目LCxii，i = 1，2，3，…，η，其中η為與節(jié)點(diǎn)化cxi相連接 的邊的數(shù)目；則分支： 旋支近段(pCx ):節(jié)點(diǎn)化到節(jié)點(diǎn)化Cxl的邊； 旋支中遠(yuǎn)段化Cx ):目LCxl冷OLCx η ELCxli e PLlCxI，此時(shí)區(qū)域范圍有目LCxli e A3 ; 第一純緣支(OMl ):日化位〇 ElCxU g PLlCx.1，此時(shí)區(qū)域范圍有目LCxli E A2，其中垃Cxli為 目LCxli所在的邊； 若Plcxi存在第二個(gè)分叉節(jié)點(diǎn)到節(jié)點(diǎn)化的邊的長度不大于闊值Τ(參考值80)，則該節(jié)點(diǎn)為 第二純緣支的起始節(jié)點(diǎn)Ν0Μ2 ;則分支： 第二純緣支(0Μ2):節(jié)點(diǎn)Ν〇μ2發(fā)出的非PUcxi上的邊； 若Plcxi后續(xù)分叉節(jié)點(diǎn)到節(jié)點(diǎn)化的邊的長度大于闊值Τ(參考值80)，則該節(jié)點(diǎn)為Nlcx2;查 找步驟5(2)所得球空間中值角度中節(jié)點(diǎn)化Cx2所連各邊的中值角度0LCx2i，i = l，2,3,…，n2， 其中n2為與節(jié)點(diǎn)化cx細(xì)連接的邊的數(shù)目；則分支： 旋支起源后降支化-PDA):目lXx2i《0,此時(shí)區(qū)域范圍有目LCx2ie (A1UA2); 旋支起源后側(cè)支化-PLB):目LCx2i〉0,此時(shí)區(qū)域范圍有目LCx2ie(A3UA4); 若Plcxi上第一個(gè)分叉節(jié)點(diǎn)到節(jié)點(diǎn)化的邊的長度大于闊值Τ(參考值80)，則該節(jié)點(diǎn)為 化Cx2 ;查找步驟5(2)所得球空間中值角度中節(jié)點(diǎn)化Cx2所連各邊的中值角度白隱，i = 1，2， 3，…，n2，其中n2為與節(jié)點(diǎn)化相連接的邊的數(shù)目；則分支： kPDA :目LCx2i《0，此時(shí)區(qū)域范圍有目LCx2i e ( Al U A2 ); 心口18:目1。曲〉0，此時(shí)區(qū)域范圍有扣械把^3。八4); pCx與LCx:節(jié)點(diǎn)化到節(jié)點(diǎn)化Cx2的邊； W節(jié)點(diǎn)化為球屯、，計(jì)算節(jié)點(diǎn)化到節(jié)點(diǎn)化Gx2的邊上所有點(diǎn)的球坐標(biāo)，得到該邊上各點(diǎn)按連 接順序的方位角數(shù)組目21和半徑數(shù)組。1，1 = 1，2,3^''，11，其中11為該邊上點(diǎn)的總數(shù)目；則分 段： pCx :節(jié)點(diǎn)化到節(jié)點(diǎn)N說X.2.郵η的邊，其中NLCx2mm為r2i X S in (目2i )中最小值所在的點(diǎn)壞一定 為步驟2所得Ξ維樹形結(jié)構(gòu)中的節(jié)點(diǎn)）； LCx :節(jié)點(diǎn)NLCx2>"i。到節(jié)點(diǎn)化Cx2的邊； d)中間支(RI):目Li=0Ln，此時(shí)目Lie(AlUA2); 若RI是游離邊，則所定義游離邊即為RI;若RI不是游離邊，其下所有分支均記入RI。7.根據(jù)權(quán)利要求1所述基于球空間劃分的冠狀動(dòng)脈自動(dòng)分段與解剖學(xué)標(biāo)記方法，其特 征在于，所述步驟6具體執(zhí)行W下步驟： (1) 計(jì)算Ξ維樹形結(jié)構(gòu)右冠狀動(dòng)脈的游離邊和葉節(jié)點(diǎn)，并計(jì)算節(jié)點(diǎn)Nr到右冠狀動(dòng)脈所有 葉節(jié)點(diǎn)的最短路徑節(jié)點(diǎn)和路徑上的邊，方法與步驟5(1)相同； (2) 計(jì)算步驟6(1)所得所有游離邊、路徑上的邊相對(duì)于其連接的起始節(jié)點(diǎn)的球空間中 值角度;方法與步驟5(2)相同； (3) 右冠狀動(dòng)脈的具體分段： a) 查找步驟6(1)所得路徑中的最長路徑扣CA和路徑上的邊化RCA，并將扣CA最末端葉節(jié)點(diǎn) 所連最長游離邊置于PLrca末尾；W節(jié)點(diǎn)Nr為球屯、，計(jì)算化RCA上所有點(diǎn)的球坐標(biāo)，得到PLrca上 各點(diǎn)按連接順序的方位角數(shù)組目Ri和半徑數(shù)組rRi ? = 1，2,3，···，η，其中η為化RCA上點(diǎn)的總數(shù) 目；則分段： 右冠近段(pRCA):節(jié)點(diǎn)Nr到節(jié)點(diǎn)NmRGA的邊； 右冠中段(mRCA ):節(jié)點(diǎn)NmRCA到節(jié)點(diǎn)N.RCAm.i.。的邊;其中Nrow。為rRi X COS (目Ri )中最小值所 在的點(diǎn)（不一定為步驟2所得Ξ維樹形結(jié)構(gòu)中的節(jié)點(diǎn))且區(qū)域范圍內(nèi)NRCAw?E(A3UA4)，NmRCA 為節(jié)點(diǎn)Nr到節(jié)點(diǎn)Nrcawii的邊的中點(diǎn)； b) 若右冠沒有分支節(jié)點(diǎn)，或右冠主干上存在分支節(jié)點(diǎn)，且分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nr的邊的長度 不大于闊值Τ2(參考值200)，則分段： 右冠遠(yuǎn)段(dRCA):節(jié)點(diǎn)NRCAmi,劇節(jié)點(diǎn)NRCAmax的邊，其中NRCAmax為rRi X COS (目Ri )中最大值 所在的點(diǎn)(不一定為步驟2所得立維樹形結(jié)構(gòu)中的節(jié)點(diǎn))且區(qū)域范圍內(nèi)化CAmaxe(A2UA3:.); 右冠起源后降支(R-PDA ):節(jié)點(diǎn)NrcA胃到化RGA結(jié)束的邊;若NRCAmax為PLrGA結(jié)束的點(diǎn)，貝IJ R-PDA不存在； C)若右冠主干上存在分支節(jié)點(diǎn)，且分叉節(jié)點(diǎn)到節(jié)點(diǎn)Nr的邊的長度大于闊值Τ2(參考值 200)，貝喊節(jié)點(diǎn)為Nrcai ;查找步驟6 (2)所得球空間中值角度中節(jié)點(diǎn)Nrcai所連的非PLrca上的邊 的中值角度目RCAi，i = 1，2，3，…，η，其中η為與節(jié)點(diǎn)Nrca袖連接的非化RCA上的邊的數(shù)目；則分 段： dRCA :節(jié)點(diǎn)NRCAmin到節(jié)點(diǎn)Nrcai的邊； R-PDA:目RCAi《-90，此時(shí)區(qū)域范圍有目RCAi e Ai; 右冠起源后側(cè)支(R-化B):-90<目RCAi《0,此時(shí)區(qū)域范圍有目RCAiGAs。
【文檔編號(hào)】G06T7/00GK106097298SQ201610362810
【公開日】2016年11月9日
【申請(qǐng)日】2016年5月27日
【發(fā)明人】李穎, 王如意, 張紹祥, 譚立文
【申請(qǐng)人】中國人民解放軍第三軍醫(yī)大學(xué)