技術(shù)特征：

1.一種血管造影中介下先心病術(shù)中血管辨識的裝置，包括術(shù)前CT圖像預(yù)處理模塊、術(shù)前CT圖像分割模塊、術(shù)前CT圖像三維重建模塊、術(shù)前CT圖像二維平面化模塊、術(shù)中造影成像模塊、術(shù)中配準(zhǔn)初始化模塊、術(shù)中量子更新模塊、術(shù)中個(gè)體變異模塊、術(shù)中變異后處理模塊、圖像配準(zhǔn)模塊、血管辨識結(jié)果展示模塊，其中：

所述的術(shù)前CT圖像預(yù)處理模塊，用于從術(shù)前心臟造影CT圖像中選擇一幅主動(dòng)脈瓣CT圖像作為初始圖像；標(biāo)記初始圖像中心臟和血管所在的位置，并將其作為圖割方法中無向圖的權(quán)值；

所述的術(shù)前CT圖像分割模塊，是利用圖割方法，對術(shù)前CT圖像預(yù)處理模塊得到的初始圖像進(jìn)行分割，得到分割后的參考模板；用參考模板中心臟和血管所在的位置作為圖割方法的新的權(quán)值，采用圖割方法，對當(dāng)前參考模板分割圖像時(shí)的前一幅圖像進(jìn)行分割，得到并保存更新后的參考模板；判斷當(dāng)前參考模板更新時(shí)所選取的圖像是否為所有CT圖像的最前一幅圖像，若是，則執(zhí)行下一步，否則，繼續(xù)分割圖像；用參考模板中的心臟和血管位置作為圖割方法的新的權(quán)值，采用圖割方法，對當(dāng)前參考模板分割圖像時(shí)的后一幅圖像進(jìn)行分割，得到并保存更新后的參考模板；判斷當(dāng)前參考模板更新時(shí)所選取的圖像是否為所有CT圖像的最后一幅圖像，若是，則將所有的參考模板傳入三維重建模塊，否則，繼續(xù)分割圖像；

所述的術(shù)前CT圖像三維重建模塊，是將所有的參考模板加載到具有三維可視化功能的軟件中；利用軟件的體繪制選項(xiàng)，對所有的參考模板中的心臟和血管，進(jìn)行三維可視化重建；

所述的術(shù)前CT圖像二維平面化模塊，是將先心病患者手術(shù)時(shí)擬平躺于水平手術(shù)臺(tái)且面向醫(yī)生時(shí)的心臟位置作為基準(zhǔn)，以先心病患者心臟的主動(dòng)脈竇的中心為坐標(biāo)原點(diǎn)；以過坐標(biāo)原點(diǎn)且垂直于手術(shù)臺(tái)平面的直線為Z軸，將先心病患者面部所在的方向確定為Z軸正方向；以過原點(diǎn)且垂直于脊柱中心線所在平面的直線為X軸，將指向先心病患者左臂的方向確定為X軸正方向；以過原點(diǎn)且垂直于X軸和Z軸所在平面的直線為Y軸，將指向先心病患者頭部的方向確定為Y軸正方向，完成術(shù)前心臟和血管標(biāo)準(zhǔn)坐標(biāo)系的構(gòu)建；構(gòu)建術(shù)前心臟和血管標(biāo)準(zhǔn)坐標(biāo)系完成后，保持建立標(biāo)準(zhǔn)坐標(biāo)系時(shí)分割的術(shù)前心臟和血管的位置不變，從Z軸正半軸看向負(fù)半軸的方向上，得到一張分割的術(shù)前心臟和血管的二維圖像；以Y軸為旋轉(zhuǎn)軸，將分割的術(shù)前心臟和血管分別向X軸的正方向和反方向偏移五度，從Z軸正半軸看向負(fù)半軸的方向上，分別得到兩張分割的術(shù)前心臟和血管的二維圖像；以X軸為旋轉(zhuǎn)軸，將分割的術(shù)前心臟和血管分別向Y軸的正方向和反方向偏移五度，從Z軸正半軸看向負(fù)半軸的方向上，分別得到兩張分割的術(shù)前心臟和血管的二維圖像；

所述的術(shù)中造影成像模塊，用于獲取先心病患者心臟和血管的術(shù)中造影像；

所述的術(shù)中配準(zhǔn)初始化模塊，是將從術(shù)前CT圖像二維平面化模塊中獲得的五張分割的術(shù)前心臟和血管的二維圖像中任選一張圖像；分別將配準(zhǔn)所需要的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群中的每個(gè)種群初始化為5個(gè)個(gè)體，每個(gè)個(gè)體在其取值范圍內(nèi)取一個(gè)隨機(jī)值；利用個(gè)體組合的互信息計(jì)算方法，計(jì)算所選的分割的術(shù)前心臟和血管的二維圖像與術(shù)中造影像配準(zhǔn)的互信息，得到當(dāng)前的極大互信息和當(dāng)前的極大互信息的個(gè)體組合；將水平位移種群、垂直位移種群的復(fù)制規(guī)模設(shè)定為50，旋轉(zhuǎn)角度種群的復(fù)制規(guī)模設(shè)定為20，根據(jù)每種個(gè)體組合的互信息占所有個(gè)體組合的互信息總和的比例，對個(gè)體組合中的個(gè)體按比例進(jìn)行復(fù)制，生成復(fù)制后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群；

所述的術(shù)中量子更新模塊，是利用量子更新方法，對術(shù)中配準(zhǔn)初始化模塊復(fù)制后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群中的個(gè)體進(jìn)行量子更新，得到更新后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群；

所述的術(shù)中個(gè)體變異模塊，是將每個(gè)種群中個(gè)體發(fā)生變異的概率設(shè)定為20％，對術(shù)中量子更新模塊更新后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群中變異的個(gè)體在其取值范圍內(nèi)取隨機(jī)值，得到變異后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群；

所述的術(shù)中變異后處理模塊，是對術(shù)中個(gè)體變異模塊獲得的變異后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群，利用個(gè)體組合的互信息計(jì)算方法，計(jì)算所選的分割的術(shù)前心臟和血管的二維圖像和術(shù)中造影像配準(zhǔn)的互信息，得到當(dāng)前的極大互信息和當(dāng)前的極大互信息的個(gè)體組合；判斷當(dāng)前的極大互信息是否小于配準(zhǔn)初始化模塊的極大互信息，若是，則獲得極大互信息中的最大值，并執(zhí)行下一步，否則，執(zhí)行水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群的復(fù)制步驟；判斷是否選完五張分割的術(shù)前心臟和血管的二維圖像，若是，則執(zhí)行造影像配準(zhǔn)步驟，否則，執(zhí)行得到圖像的最大互信息的步驟；

所述的圖像配準(zhǔn)模塊，是從五張分割的術(shù)前心臟和血管的二維圖像中找出與術(shù)中造影像的最大互信息；從五張分割的術(shù)前心臟和血管的二維圖像中找出得到最大互信息的分割的術(shù)前心臟和血管的二維圖像；從五張分割的術(shù)前心臟和血管的二維圖像中找出得到最大互信息的個(gè)體組合；用最大互信息的個(gè)體組合中的水平位移、垂直位移、旋轉(zhuǎn)角度，對最大互信息的分割的術(shù)前心臟和血管的二維圖像與術(shù)中造影像進(jìn)行配準(zhǔn)；

所述的血管辨識結(jié)果展示模塊，是將術(shù)前三維參考模板的角度調(diào)整到與配準(zhǔn)所用的分割的術(shù)前心臟和血管的二維圖像的角度一致后，將先心病患者術(shù)前血管的位置輸出到顯示器展示。

2.一種血管造影中介下先心病術(shù)中血管辨識的方法，包括如下步驟：

(1)分割術(shù)前心臟造影CT圖像：

(1a)術(shù)前CT圖像預(yù)處理模塊，從術(shù)前心臟造影CT圖像中選擇一幅主動(dòng)脈瓣CT圖像作為初始圖像；

(1b)術(shù)前CT圖像預(yù)處理模塊標(biāo)記初始圖像中心臟和血管所在的位置，并將其作為圖割方法中無向圖的權(quán)值；

(1c)術(shù)前CT圖像分割模塊利用圖割方法，對初始圖像進(jìn)行分割，得到分割后的參考模板；

(1d)術(shù)前CT圖像分割模塊用參考模板中心臟和血管所在的位置作為圖割方法的新的權(quán)值，采用圖割方法，對當(dāng)前參考模板分割圖像時(shí)的前一幅圖像進(jìn)行分割，得到并保存更新后的參考模板；

(1e)術(shù)前CT圖像分割模塊判斷當(dāng)前參考模板更新時(shí)所選取的圖像是否為所有術(shù)前心臟造影CT圖像的最前一幅圖像，若是，則執(zhí)行步驟(1f)，否則，執(zhí)行步驟(1d)；

(1f)術(shù)前CT圖像分割模塊用參考模板中心臟和血管所在的位置作為圖割方法的新的權(quán)值，采用圖割方法，對當(dāng)前參考模板分割圖像時(shí)的后一幅圖像進(jìn)行分割，得到并保存更新后的參考模板；

(1g)判斷當(dāng)前參考模板更新時(shí)所選取的圖像是否為所有術(shù)前心臟造影CT圖像的最后一幅圖像，若是，則執(zhí)行步驟(2)，否則，執(zhí)行步驟(1f)；

(2)重建術(shù)前三維參考模板：

(2a)術(shù)前CT圖像三維重建模塊將所有參考模板加載到具有三維可視化功能的軟件中；

(2b)術(shù)前CT圖像三維重建模塊利用軟件的體繪制選項(xiàng)，對參考模板中的心臟和血管，進(jìn)行三維可視化重建；

(3)構(gòu)建術(shù)前心臟和血管標(biāo)準(zhǔn)坐標(biāo)系：

(3a)術(shù)前CT圖像二維平面化模塊將先心病患者手術(shù)時(shí)擬平躺于水平手術(shù)臺(tái)且面向醫(yī)生時(shí)的心臟位置作為基準(zhǔn)，以先心病患者心臟的主動(dòng)脈竇的中心為坐標(biāo)原點(diǎn)；

(3b)術(shù)前CT圖像二維平面化模塊，以過坐標(biāo)原點(diǎn)且垂直于手術(shù)臺(tái)平面的直線為Z軸，將先心病患者面部所在的方向確定為Z軸正方向；

(3c)術(shù)前CT圖像二維平面化模塊，以過原點(diǎn)且垂直于脊柱中心線所在平面的直線為X軸，將指向先心病患者左臂的方向確定為X軸正方向；

(3d)術(shù)前CT圖像二維平面化模塊，以過原點(diǎn)且垂直于X軸和Z軸所在平面的直線為Y軸，將指向先心病患者頭部的方向確定為Y軸正方向，完成術(shù)前心臟和血管標(biāo)準(zhǔn)坐標(biāo)系的構(gòu)建；

(4)將分割的術(shù)前心臟和血管平面化：

(4a)保持建立標(biāo)準(zhǔn)坐標(biāo)系時(shí)分割的術(shù)前心臟和血管的位置不變，術(shù)前CT圖像二維平面化模塊，從Z軸正半軸看向負(fù)半軸的方向上，得到一張分割的術(shù)前心臟和血管的二維圖像；

(4b)以Y軸為旋轉(zhuǎn)軸，術(shù)前CT圖像二維平面化模塊將分割的術(shù)前心臟和血管分別向X軸的正方向和反方向偏移五度，從Z軸正半軸看向負(fù)半軸的方向上，分別得到兩張分割的術(shù)前心臟和血管的二維圖像；

(4c)以X軸為旋轉(zhuǎn)軸，術(shù)前CT圖像二維平面化模塊將分割的術(shù)前心臟和血管分別向Y軸的正方向和反方向偏移五度，從Z軸正半軸看向負(fù)半軸的方向上，分別得到兩張分割的術(shù)前心臟和血管的二維圖像；

(5)獲取先心病患者的術(shù)中造影像：

利用術(shù)中造影成像模塊，獲取先心病患者心臟和血管的術(shù)中造影像；

(6)獲得圖像的最大互信息：

(6a)術(shù)中配準(zhǔn)初始化模塊，從步驟(4a)、步驟(4b)、步驟(4c)中獲得的五張分割的術(shù)前心臟和血管的二維圖像中任選一張圖像；

(6b)術(shù)中配準(zhǔn)初始化模塊，分別將配準(zhǔn)所需要的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群中的每個(gè)種群初始化為5個(gè)個(gè)體，每個(gè)個(gè)體在其取值范圍內(nèi)取一個(gè)隨機(jī)值；

(6c)術(shù)中配準(zhǔn)初始化模塊，利用個(gè)體組合的互信息計(jì)算方法，計(jì)算所選的分割的術(shù)前心臟和血管的二維圖像與術(shù)中造影像配準(zhǔn)的互信息，得到當(dāng)前的極大互信息和當(dāng)前的極大互信息的個(gè)體組合；

(6d)術(shù)中配準(zhǔn)初始化模塊，將水平位移種群、垂直位移種群的復(fù)制規(guī)模設(shè)定為50，旋轉(zhuǎn)角度種群的復(fù)制規(guī)模設(shè)定為20，根據(jù)每種個(gè)體組合的互信息占所有個(gè)體組合的互信息總和的比例，對個(gè)體組合中的個(gè)體按比例進(jìn)行復(fù)制，生成復(fù)制后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群；

(6e)術(shù)中量子更新模塊，利用量子更新方法，對復(fù)制后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群中的個(gè)體進(jìn)行量子更新，得到更新后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群；

(6f)將每個(gè)種群中個(gè)體發(fā)生變異的概率設(shè)定為20％，術(shù)中個(gè)體變異模塊對更新后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群中變異的個(gè)體在其取值范圍內(nèi)取隨機(jī)值，得到變異后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群；

(6g)術(shù)中變異后處理模塊，利用個(gè)體組合的互信息計(jì)算方法，對變異后的水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群，計(jì)算所選的分割的術(shù)前心臟和血管的二維圖像和術(shù)中造影像配準(zhǔn)的互信息，得到當(dāng)前的極大互信息和當(dāng)前的極大互信息的個(gè)體組合；

(6h)判斷當(dāng)前的極大互信息是否小于上一次迭代時(shí)的極大互信息，若是，則獲得極大互信息中的最大值，執(zhí)行步驟(7)，否則，執(zhí)行步驟(6d)；

(7)判斷是否選完五張分割的術(shù)前心臟和血管的二維圖像，若是，則執(zhí)行步驟(8)，否則，執(zhí)行步驟(6)；

(8)造影像配準(zhǔn)：

(8a)圖像配準(zhǔn)模塊，從五張分割的術(shù)前心臟和血管的二維圖像中找出與術(shù)中造影像的最大互信息；

(8b)圖像配準(zhǔn)模塊，從五張分割的術(shù)前心臟和血管的二維圖像中找出得到最大互信息的分割的術(shù)前心臟和血管的二維圖像；

(8c)圖像配準(zhǔn)模塊，從五張分割的術(shù)前心臟和血管的二維圖像中找出得到最大互信息的個(gè)體組合；

(8d)圖像配準(zhǔn)模塊，用最大互信息的個(gè)體組合中的水平位移、垂直位移、旋轉(zhuǎn)角度，對最大互信息的分割的術(shù)前心臟和血管的二維圖像與術(shù)中造影像進(jìn)行配準(zhǔn)；

(9)辨識血管的位置：

(9a)將術(shù)前三維參考模板的角度調(diào)整到與配準(zhǔn)所用的分割的術(shù)前心臟和血管的二維圖像的角度一致，獲得先心病患者術(shù)前血管的位置圖像；

(9b)將先心病患者術(shù)前血管的位置圖像，輸出到血管辨識結(jié)果展示模塊中的顯示器，顯示辨識血管的位置。

3.根據(jù)權(quán)利要求2所述的血管造影中介下先心病術(shù)中血管辨識的方法，其特征在于，步驟(1c)中所述圖割方法的具體步驟如下：

第1步，將標(biāo)記心臟和血管所在的位置后的初始圖像，映射為帶權(quán)值的無向圖；

第2步，利用福特-?？松璅ord-Fulkerson方法，從帶權(quán)值的無向圖中得到最大流；

第3步，根據(jù)最大流最小割定理，將帶權(quán)值的無向圖的最大流作為最小割；

第4步，將最小割分開的帶權(quán)值的無向圖映射回參考圖像，得到分割完成的參考模板。

4.根據(jù)權(quán)利要求2所述的血管造影中介下先心病術(shù)中血管辨識的方法，其特征在于：步驟(6b)中所述的水平位移種群和垂直位移種群的隨機(jī)值的取值范圍均為[-10,10]，旋轉(zhuǎn)角度種群的隨機(jī)值的取值范圍為[-5,5]。

5.根據(jù)權(quán)利要求2所述的血管造影中介下先心病術(shù)中血管辨識的方法，其特征在于：步驟(6c)、步驟(6g)中所述個(gè)體組合的互信息計(jì)算方法的具體步驟如下：

第1步，分別從水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群中任意取出一個(gè)個(gè)體，組成個(gè)體組合；

第2步，按照下式，計(jì)算分割的術(shù)前心臟和血管的二維圖像的熵：

$H\left(M\right)=-\underset{m}{\mathrm{\Σ}}f\left(m\right)\lg f\left(m\right)$

其中，H(M)表示分割的術(shù)前心臟和血管的二維圖像的熵，∑表示求和操作，m表示分割的術(shù)前心臟和血管的二維圖像中的像素灰度值，其取值范圍為[0,255]，f(·)表示像素灰度值的均勻分布概率密度函數(shù)，lg表示以10為底的對數(shù)操作；

第3步，按照下式，計(jì)算術(shù)中造影像的熵：

$H\left(N\right)=-\underset{n}{\mathrm{\Σ}}f\left(n\right)\lg f\left(n\right)$

其中，H(N)表示術(shù)中造影像的熵，∑表示求和操作，n表示術(shù)中造影像中的像素灰度值，其取值范圍為[0,255]，f(·)表示像素灰度值的均勻分布概率密度函數(shù)，lg表示以10為底的對數(shù)操作；

第4步，按照下式，計(jì)算分割的術(shù)前心臟和血管的二維圖像和術(shù)中造影像的聯(lián)合熵：

$K=-\underset{a}{\mathrm{\Σ}}\underset{b}{\mathrm{\Σ}}f(a,b)\lg f(a,b)$

其中，K表示分割的術(shù)前心臟和血管的二維圖像和術(shù)中造影像的聯(lián)合熵，∑表示求和操作，a和b分別表示像素灰度值，其范圍均是[0,255]，f(·)表示像素灰度值的均勻分布概率密度函數(shù)，lg表示以10為底的對數(shù)操作；

第5步，按照下式，計(jì)算分割的術(shù)前心臟和血管的二維圖像和術(shù)中造影像的互信息：

$MI=\frac{H\left(M\right)+H\left(N\right)}{K}$

其中，MI表示分割的術(shù)前心臟和血管的二維圖像和術(shù)中造影像的互信息，H(M)表示分割的術(shù)前心臟和血管的二維圖像的熵，H(N)表示術(shù)中造影像的熵，K表示分割的術(shù)前心臟和血管的二維圖像和術(shù)中造影像的聯(lián)合熵；

第6步，判斷是否獲得所有個(gè)體組合的互信息，若是，則執(zhí)行第7步，否則，執(zhí)行第1步；

第7步，對所有個(gè)體組合的互信息進(jìn)行排序，找出當(dāng)前的極大互信息和得到當(dāng)前的極大互信息的個(gè)體組合。

6.根據(jù)權(quán)利要求2所述的血管造影中介下先心病術(shù)中血管辨識的方法，其特征在于：步驟(6e)中所述量子更新方法的具體步驟如下：

第1步，按照下式，構(gòu)造量子旋轉(zhuǎn)矩陣：

$G=\left(\begin{array}{cc}cos(l\×p)& sin(l\×p)\\ 0& 0\end{array}\right)$

其中，G表示量子旋轉(zhuǎn)矩陣，cos(·)表示求余弦操作，sin(·)表示求正弦操作，l表示步長，取值為0.9，r表示[0,1]中產(chǎn)生的隨機(jī)數(shù)，p表示根據(jù)隨機(jī)數(shù)r得到的方向數(shù)，如果隨機(jī)數(shù)r小于等于1/3，則取方向數(shù)p為-1，如果隨機(jī)數(shù)r大于1/3且小于等于2/3，則取方向數(shù)p為0，如果隨機(jī)數(shù)r大于2/3，則取方向數(shù)p為1；

第2步，按照下式，構(gòu)造個(gè)體矩陣：

$A=\left(\begin{array}{cc}x& a_1\\ a_2& a_3\end{array}\right)$

其中，A表示個(gè)體矩陣，x表示水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群中的任意一個(gè)個(gè)體，a₁,a₂,a₃表示取值不同的三個(gè)量子位，其數(shù)值按照下式產(chǎn)生：

$a_i=\frac{rand}{randlen}$

其中，a_i表示取值不同的三個(gè)量子位中的第i個(gè)，i＝1,2,3,rand表示水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群中的任意一個(gè)個(gè)體的取值范圍內(nèi)的一個(gè)隨機(jī)值，randlen表示水平位移種群、垂直位移種群、旋轉(zhuǎn)角度種群中的任意一個(gè)個(gè)體的取值范圍的區(qū)間長度；

第3步，將量子旋轉(zhuǎn)矩陣G乘以個(gè)體矩陣A，得到結(jié)果矩陣；

第4步，將結(jié)果矩陣中的第一行第一列的元素，作為量子更新操作后的個(gè)體。

7.根據(jù)權(quán)利要求2所述的血管造影中介下先心病術(shù)中血管辨識的方法，其特征在于：步驟(6f)中所述的水平位移種群、垂直位移種群個(gè)體的變異取值范圍為[-10,10]，旋轉(zhuǎn)角度種群個(gè)體的變異取值范圍為[-5,5]。