一種醫(yī)學(xué)圖像的多平面重建成像方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及一種醫(yī)學(xué)圖像成像方法，尤其涉及一種醫(yī)學(xué)圖像的多平面重建成像方法。
【背景技術(shù)】
[0002]多平面重建成像技術(shù)(MPR, Mult1-planer Reformat1n)的圖像是根據(jù)醫(yī)學(xué)成像設(shè)備采集到的初始體數(shù)據(jù)(volumetric data)后，經(jīng)過數(shù)據(jù)處理生成的。通用的MPR圖像處理過程為:在視圖空間內(nèi)獲取坐標(biāo)，將該坐標(biāo)轉(zhuǎn)換成體空間，使用任意形式的插補(bǔ)，對體數(shù)據(jù)進(jìn)行再采樣，并生成針對離散視圖空間坐標(biāo)的新的MPR數(shù)據(jù)值。MPR片層(二維圖像序列)通過在固定z值下對多個(gè)(x，y)坐標(biāo)執(zhí)行所述處理而形成。針對z的多個(gè)值重復(fù)執(zhí)行所述處理時(shí)，決定多個(gè)MPR片層，將這些片層投影而能夠形成MPR三維圖像序列。因此，MPR圖像具備與投影面平行地對齊、且配設(shè)在沿著z軸的不同位置的一系列MPR片層。
[0003]上述對于三維目標(biāo)的多平面重建成像過程中，一個(gè)關(guān)鍵之處是需要確定多平面成像的重建中心，基于重建中心在三個(gè)正交方向上生成MPR片層，但是目前并無確定準(zhǔn)確的重建中心的有效方法。如果重建中心不在需要成像的三維目標(biāo)正中心處則對MPR圖像的影響很大，因?yàn)镸PR圖像很有可能只是反映了目標(biāo)的局部信息。
[0004]因此，有必要提供一種能夠快速準(zhǔn)確的確定多平面成像的重建中心的方法，能夠得到更準(zhǔn)確反映成像目標(biāo)三維信息的MPR圖像。

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

[0005]本發(fā)明所要解決的技術(shù)問題是提供一種醫(yī)學(xué)圖像的多平面重建成像方法，能夠更準(zhǔn)確的獲得檢查區(qū)域的有效區(qū)域圖像，提高定位及重建成像的準(zhǔn)確性。
[0006]本發(fā)明為解決上述技術(shù)問題而采用的技術(shù)方案是提供一種醫(yī)學(xué)圖像的多平面重建成像方法，包括如下步驟:a)掃描待成像目標(biāo)的檢查區(qū)域，獲取所述檢查區(qū)域的三維圖像序列，所述三維圖像序列由一組二維圖像序列組成山)選取所述二維圖像序列內(nèi)的有效區(qū)域圖像序列，并確定所述有效區(qū)域圖像序列的三維中心；c)以所述有效區(qū)域圖像序列的三維中心作為重建中心，用于進(jìn)行所述檢查區(qū)域的多平面重建成像。
[0007]上述的醫(yī)學(xué)圖像的多平面重建成像方法，其中，所述步驟b)中有效區(qū)域圖像序列的選取過程包括:對所述三維圖像序列的像素值進(jìn)行均值化；計(jì)算所述二維圖像序列的最大像素連通區(qū)域的面積；比較預(yù)設(shè)的面積閾值和所述二維圖像序列的最大像素連通區(qū)域的面積；選取所述二維圖像序列中最大像素連通區(qū)域面積高于所述面積閾值的像素區(qū)域，將所述像素區(qū)域作為所述二維圖像序列的有效區(qū)域圖像序列。
[0008]上述的醫(yī)學(xué)圖像的多平面重建成像方法，其中，所述計(jì)算二維圖像序列的最大像素連通區(qū)域的面積包括:計(jì)算所述二維圖像序列的像素值平均值；比較所述二維圖像序列各像素的像素值和所述像素值平均值；選取所述二維圖像序列中像素值大于所述像素值平均值的像素；僅對選取的所述二維圖像序列的像素計(jì)算二維圖像序列的最大像素連通區(qū)域的面積。
[0009]上述的醫(yī)學(xué)圖像的多平面重建成像方法，其中，所述計(jì)算二維圖像序列的最大像素連通區(qū)域的面積包括:通過直方圖統(tǒng)計(jì)得到所述二維圖像序列的像素值閾值；比較所述二維圖像序列各像素的像素值和所述像素值閾值；選取所述二維圖像序列中像素值大于所述像素值閾值的像素；僅對選取的所述二維圖像序列的像素計(jì)算二維圖像序列的最大像素連通區(qū)域的面積。
[0010]上述的醫(yī)學(xué)圖像的多平面重建成像方法，其中，所述步驟b)中有效區(qū)域圖像序列的選取過程包括:對所述三維圖像序列的像素值進(jìn)行均值化，所述三維圖像序列由一組二維圖像序列構(gòu)成；計(jì)算所述二維圖像序列的最大像素連通區(qū)域的像素點(diǎn)個(gè)數(shù)；比較預(yù)設(shè)的像素個(gè)數(shù)閾值和所述二維圖像序列的最大像素連通區(qū)域的像素點(diǎn)個(gè)數(shù)；選取所述二維圖像序列中最大像素連通區(qū)域像素點(diǎn)個(gè)數(shù)高于所述像素個(gè)數(shù)閾值的像素區(qū)域，將所述像素區(qū)域作為所述二維圖像序列的有效區(qū)域圖像序列。
[0011]上述的醫(yī)學(xué)圖像的多平面重建成像方法，其中，所述步驟b)在選取所述二維圖像序列的有效區(qū)域圖像序列時(shí)，僅以所述二維圖像序列的中間部分的二維圖像序列進(jìn)行選取。
[0012]上述的醫(yī)學(xué)圖像的多平面重建成像方法，其中，所述步驟a)中所述三維圖像序列通過頭部的時(shí)間飛躍法掃描后得到，并僅以所述三維圖像序列的中間40%?60%的二維圖像序列進(jìn)行選取。
[0013]上述的醫(yī)學(xué)圖像的多平面重建成像方法，其中，所述步驟a)中的三維圖像序列為MRI圖像序列、CT圖像序列或X-ray圖像序列。
[0014]本發(fā)明對比現(xiàn)有技術(shù)有如下的有益效果:本發(fā)明提供的醫(yī)學(xué)圖像的多平面重建成像方法，通過對檢查區(qū)域的三維圖像序列進(jìn)行篩選，選取有效區(qū)域圖像序列，舍去噪聲區(qū)域圖像，然后確定有效區(qū)域圖像序列的三維中心并作為多平面重建成像的重建中心，從而提高定位及重建成像的準(zhǔn)確性。此外，本發(fā)明通過先采用像素值均值或像素點(diǎn)個(gè)數(shù)對三維圖像序列進(jìn)行初步篩選，再利用最大像素連通區(qū)域的面積進(jìn)行判斷，從而保證快速可靠地獲取有效區(qū)域圖像序列。
【附圖說明】
[0015]圖1為本發(fā)明醫(yī)學(xué)圖像的多平面重建成像流程示意圖。
【具體實(shí)施方式】
[0016]下面結(jié)合附圖和實(shí)施例對本發(fā)明作進(jìn)一步的描述。
[0017]圖1為本發(fā)明醫(yī)學(xué)圖像的多平面重建成像流程示意圖。
[0018]請參見圖1，本發(fā)明提供的醫(yī)學(xué)圖像的多平面重建成像方法包括如下步驟:
[0019]步驟S1:掃描待成像目標(biāo)的檢查區(qū)域，獲取所述檢查區(qū)域的三維圖像序列；所述三維圖像序列可以為MRI圖像序列、CT圖像序列或X-ray圖像序列。所述三維圖像序列由一組二維圖像序列構(gòu)成，比如通過頭部時(shí)間飛躍法掃描后得到的一組二維圖像序列。這里的“序列”可等同理解為“數(shù)據(jù)”。
[0020]步驟S2:選取所述二維圖像序列的有效區(qū)域圖像序列，并確定所述有效區(qū)域圖像序列的三維中心；
[0021]有效區(qū)域圖像序列可通過:比較三維圖像序列中各二維圖像序列的最大像素連通區(qū)域的面積是否達(dá)到預(yù)設(shè)的面積閾值或最大像素連通區(qū)域的像素點(diǎn)個(gè)數(shù)是否達(dá)到預(yù)設(shè)的像素點(diǎn)閾值來判斷。
[0022]在面積閾值判斷(或像素點(diǎn)閾值判斷)前還可先采用像素值均值的判斷對三維圖像序列的像素進(jìn)行初步選取，以便快速可靠地獲取有效區(qū)域圖像序列。像素值均值可通過對各二維圖像序列統(tǒng)計(jì)直方圖進(jìn)行快速準(zhǔn)確地選取。此外，為了進(jìn)一步加快判定過程，本發(fā)明可僅僅針對三維圖像序列中間部分的二維圖像序列進(jìn)行選取。假如通過頭部時(shí)間飛躍法共掃描形成100個(gè)二維圖像序列，考慮到新的中心點(diǎn)不會(huì)偏離體數(shù)據(jù)中心點(diǎn)太遠(yuǎn)，所以可僅僅對中間部分的二維圖像序列進(jìn)行選取，前后兩頭的二維圖像序列直接作為噪聲區(qū)域舍棄。如僅對其中的40%-60%的二維圖像序列(第40-60個(gè)二維圖像序列)進(jìn)行處理，即從第40個(gè)二維圖像序列開始，才對圖像進(jìn)行均值化或像素點(diǎn)數(shù)初步選取，再判斷各二維圖像序列中最大像素連通區(qū)域的面積是否達(dá)到預(yù)設(shè)的面積閾值或最大像素連通區(qū)域的像素點(diǎn)個(gè)數(shù)是否達(dá)到像素點(diǎn)閾值。
[