專利名稱：：一種最大密度投影圖像數(shù)據(jù)后處理方法
技術(shù)領(lǐng)域：
：本發(fā)明涉及圖像處理
技術(shù)領(lǐng)域：
，具體涉及一種最大密度投影圖像數(shù)據(jù)后處理方法。
背景技術(shù)：
：在圖像處理技術(shù)中，三維圖像的建立可以通iW二維圖像序列的投影完成，最大密度投影(MIP，Maximumintensityprojection)作為一種三維圖像的重建技術(shù)，一直是MR(Magneticresonanceimaging,磁共振)、(computedtomography,計(jì)算機(jī)斷層成像術(shù))凄t據(jù)體得到血管、骨骼和軟組織等三維結(jié)構(gòu)的有效圖〗象重建方法。為了利用計(jì)算機(jī)快速進(jìn)行圖像處理，滿足圖像重建的時(shí)間要求，人們研究、發(fā)展了很多最大密度投影MIP算法，同時(shí)由于在生成MR、CT的原始二維圖像序列時(shí)，是通過數(shù)據(jù)采樣得到二維圖像數(shù)據(jù)的，這些圖像數(shù)據(jù)不可能完全滿足奈奎斯特采樣定律，因此經(jīng)過最大密度投影后不可能完全重現(xiàn)三維物體的原貌，這就需要對二維圖像數(shù)據(jù)在最大密度投影前進(jìn)行插值處理(稱為前插值)，或者對得到的最大密度投影三維圖像數(shù)據(jù)進(jìn)行投影后插值處理(稱為后插值)，以彌補(bǔ)原始圖像數(shù)據(jù)的不足。現(xiàn)有技術(shù)的前插值和后插值方法很多，前插值方法一般計(jì)算量是很大的，加上三維重建的數(shù)據(jù)量本身是"海量，，的，常常直接造成效率的很低，重建速度讓人難以接受。早期的CT、MR的機(jī)器，因掃描間隔很大，要進(jìn)行三維重建，保持采樣一致性5一般要進(jìn)行很大計(jì)算量的前插值處理。但現(xiàn)在多排CT,強(qiáng)場MR所獲取的圖像序列的掃描間隔可以很小，只要最大密度投影重建圖像方法合理，就可以直接重建了。后插值方法雖然效率很高，但現(xiàn)有技術(shù)的的后插值方法僅僅針對最大密度投影三維圖像數(shù)據(jù)進(jìn)行插值處理，而沒有考慮最大密度投影前二維圖像序列各圖像切片之間的空間關(guān)系，造成最終得到三維圖像的組織紋理細(xì)節(jié)模糊，有的甚至變形。
發(fā)明內(nèi)容本發(fā)明要解決的技術(shù)問題是提供一種最大密度投影圖像數(shù)據(jù)后處理方法，克服現(xiàn)有技術(shù)的后插值方法得到三維圖像的組織紋理細(xì)節(jié)模糊，甚至變形的缺陷。本發(fā)明為解決上述技術(shù)問題所采用的技術(shù)方案為一種最大密度投影圖像數(shù)據(jù)后處理方法，包括步驟Al、在對圖像序列進(jìn)行最大密度投影前，選取各幅圖像的同一位置的頂點(diǎn)，沿所述圖^f象序列遞增或遞減的順序計(jì)算并記錄相鄰兩幅圖像的所述頂點(diǎn)在投影平面上的坐標(biāo)偏移量；A2、在對所述圖像序列進(jìn)行最大密度投影后，記錄投影數(shù)據(jù)中每個象素點(diǎn)的坐標(biāo)、灰度值和所在圖像的標(biāo)識及其關(guān)系，對所述象素點(diǎn)進(jìn)行后處理，包括步驟Bl、掃描所述投影數(shù)據(jù)中的象素點(diǎn)，根據(jù)所述象素點(diǎn)的坐標(biāo)和所述坐標(biāo)偏移量計(jì)算所述象素點(diǎn)所在圖像的相鄰圖像在所述投影數(shù)據(jù)中的4殳影匹配象素點(diǎn)的坐標(biāo)，進(jìn)而找到所述投影匹配象素點(diǎn)的灰度值；B2、計(jì)算所述象素點(diǎn)和所述投影匹配象素點(diǎn)之間的距離，若所述距離大于l，則執(zhí)行步驟B3，否則針對下一所述象素點(diǎn)，返回步驟B1;B3、計(jì)算所述象素點(diǎn)和所述"t殳影匹配象素點(diǎn)之間的象素點(diǎn)的坐標(biāo)，進(jìn)而找到其灰度值；B4、根據(jù)公式<formula>formulaseeoriginaldocumentpage7</formula>計(jì)算所述象素點(diǎn)和所述投影匹配象素點(diǎn)之間的象素點(diǎn)的插值灰度值，其中&是所述象素點(diǎn)的灰度值，G&,是所述投影匹配象素點(diǎn)的灰度值，n是所述象素點(diǎn)和所述投影匹配象素點(diǎn)之間的象素點(diǎn)的個數(shù)，i(l《w")是計(jì)算時(shí)當(dāng)前象素點(diǎn)的序數(shù)；B5、若計(jì)算得到的所述插值灰度值大于所述象素點(diǎn)和所述投影匹配象素點(diǎn)之間的象素點(diǎn)的原有灰度4直，則用所述插值灰度值替換所述原有灰度值。所述的最大密度投影圖像數(shù)據(jù)后處理方法，其中所述步驟A1包括步驟All、在統(tǒng)一坐標(biāo)系中建立所述圖像序列的單位法向量；Al2、在統(tǒng)一坐標(biāo)系中建立所述圖像序列中各幅圖像的平面方程；A13、在統(tǒng)一坐標(biāo)系中建立投影面的平面方程；A14、計(jì)算所述圖像序列中相鄰圖像在投影面上的偏移量。所述的最大密度投影圖像數(shù)據(jù)后處理方法，其中所述圖像序列設(shè)為醫(yī)學(xué)領(lǐng)域的數(shù)字圖像和通訊DICOM標(biāo)準(zhǔn)的圖像序列。所述的最大密度投影圖像數(shù)據(jù)后處理方法，其中根據(jù)醫(yī)學(xué)領(lǐng)域的數(shù)字圖像和通訊DICOM標(biāo)準(zhǔn)的圖像數(shù)據(jù)中的標(biāo)簽TAG值確定建立所述圖像序列的單位法向量所需的參數(shù)。所述的最大密度投影圖像數(shù)據(jù)后處理方法，其中根據(jù)醫(yī)學(xué)領(lǐng)域的數(shù)字圖像和通訊DICOM標(biāo)準(zhǔn)的圖像數(shù)據(jù)中的標(biāo)簽TAG值確定所述圖像序列中各幅圖像的平面方程所需的參數(shù)。所述的最大密度投影圖像數(shù)據(jù)后處理方法，其中根據(jù)醫(yī)學(xué)領(lǐng)域的數(shù)字圖像和通訊DICOM標(biāo)準(zhǔn)的圖像數(shù)據(jù)中的標(biāo)簽TAG值確定建立投影面的平面方程所需的參數(shù)。所述的最大密度投影圖像數(shù)據(jù)后處理方法，其中從所述圖像序列中任一選定圖像的圖像數(shù)據(jù)中取出標(biāo)簽TAG值確定相關(guān)的參數(shù)取出TAG(0018，5100)值，確定統(tǒng)一坐標(biāo)系的方向；取出TAG(0020，0032)值，確定所述選定圖4象首個1象素在所述統(tǒng)一坐標(biāo)系中的三維坐標(biāo)值；耳又出TAG(0020，0037)值，確定所述選定圖^f象的單位行向量和單位列向量，計(jì)算所述單位行向量和單位列向量的叉積得到所述圖像序列的單位法向量；取出圖^f象TAG(0028，0030)值，確定所述選定圖l象的l象素的實(shí)際物理寬度和實(shí)際物理高度；:f又出TAG(0028,OOIO)值，確定所述選定圖^f象的高度；取出TAG(0028，OOll)值，確定所述選定圖4象的寬度。所述的最大密度投影圖像數(shù)據(jù)后處理方法，其中根據(jù)所述首個像素的三維坐標(biāo)值和所述單位法向量確定所述選定圖Y象在所述統(tǒng)一坐標(biāo)系中的平面方程。本發(fā)明的有益效果為本發(fā)明最大密度投影圖像數(shù)據(jù)后處理方法，針對圖像序列中相鄰圖像的投影數(shù)據(jù)進(jìn)行灰度值插值處理，并且并不是機(jī)械地用計(jì)算得到的插值灰度值替換原有灰度值，而是先將計(jì)算得到的插值灰度值和原有灰度值進(jìn)行比較，然后再決定是否替換，這樣就克服了現(xiàn)有技術(shù)的前插值和后插值的缺陷，使經(jīng)過后處理的三維投影圖像更加逼真，組織紋理細(xì)節(jié)清晰，不會產(chǎn)生變形。本發(fā)明包括如下附圖圖1為本發(fā)明圖像序列中相鄰圖像在投影面上的投影象素點(diǎn)之間的偏移量示意圖；圖2為本發(fā)明經(jīng)過最大密度投影后得到的才殳影數(shù)據(jù)存儲格式示意圖；圖3為本發(fā)明投影數(shù)據(jù)中象素點(diǎn)與其投影匹配象素點(diǎn)的位置關(guān)系示意圖；圖4為現(xiàn)有技術(shù)棋盤格插值后的三維人體組織圖像；圖5為本發(fā)明對最大密度投影后得到的投影數(shù)據(jù)后處理后得到的三維人體組織圖像。具體實(shí)施方式下面根據(jù)附圖和實(shí)施例對本發(fā)明作進(jìn)一步詳細(xì)說明1.記錄相鄰切片層在投影面上的偏移量如圖1所示，假設(shè)圖像序列中共有iV,幅圖像，記為Q,込，L，Ov,各圖像序列的頂點(diǎn)對應(yīng)到二維平面上的點(diǎn)分別記錄下來為A,^，L，pw。在二維平面上，從這些點(diǎn)可以得到P2相對于《的偏移量(Ax,，Ay,)，g相對于尸2的(AX2，A^),…"直至Pw相對于i^的偏移量(A^-i，A;v,)。2.對MIP圖像進(jìn)4亍處理現(xiàn)有技術(shù)中對圖像序列進(jìn)行最大密度投影的方法很多，本發(fā)明由于是對經(jīng)過最大密度投影后得到的數(shù)據(jù)進(jìn)行后處理，并不局限于特定的最大密度投影算法，因此最大密度投影的方法本申請中不再贅述。2.1最大密度投影得到的圖像數(shù)據(jù)假設(shè)經(jīng)過最大密度投影得到的圖像數(shù)據(jù)的存儲格式及大小如圖二所示，共有&行，Sc列.其中每個數(shù)據(jù)點(diǎn)包含的信息為當(dāng)前數(shù)據(jù)在圖像中的坐標(biāo)(x,y)，該坐標(biāo)下的圖像灰度值gray,以及其來自圖像序列中哪一個切片層k。2.2圖像數(shù)據(jù)進(jìn)行逐行掃描得到匹配點(diǎn)對圖像數(shù)據(jù)進(jìn)行逐行的掃描，假設(shè)掃描到一點(diǎn)&(下標(biāo)sr為source的縮寫)，該點(diǎn)的數(shù)據(jù)信息為圖像中坐標(biāo)為(&,;0,灰度值為(^，以及來增加的圖像序列切片層^.利用第一步驟的得到的偏移量得到要對點(diǎn)&進(jìn)行處理需要的另一個坐標(biāo)點(diǎn)尸血(下標(biāo)des為destination的縮寫),具體如下已知尸w來自于第i^層切片，從第一步得到的圖像序列各切片層相對偏移量可以知道K，,層相對與《,-l層切片的偏移量記為(Ax-AyJ,則坐標(biāo)點(diǎn)i^為知道尸血的坐標(biāo)位置為(^,4)，從而知道該點(diǎn)的灰度值為Gf注意如果(^，^)的值超出圖像的范圍即^<0，4<0，或2.3根據(jù)掃描的點(diǎn)的間距進(jìn)行判斷是否處理從步驟2.2中知道了兩點(diǎn)的坐標(biāo)信息以及灰度值，根據(jù)兩點(diǎn)的距離/>式可以知道它們間的距離d:x血j=xw—7c/es—凡/"AjV^(4-1)10如果dd,即點(diǎn)&與i^相鄰一個像素或重合則不需要對像素?cái)?shù)據(jù)進(jìn)行處理，進(jìn)入下一點(diǎn)的掃描.如果c^2，在不影響數(shù)字圖像運(yùn)算的情況下按照四舍五入取d的整數(shù)值n:這說明&與尸血間的n個點(diǎn)，不妨記為，尸血+2，...，下面對n個像素點(diǎn)進(jìn)行處理.2.4求取掃描點(diǎn)與匹配點(diǎn)間的^^素點(diǎn)的坐標(biāo)經(jīng)過步驟2.3可知，如果兩點(diǎn)距離6/Sl，不需要進(jìn)行處理，直接跳到步驟b，進(jìn)行下一個掃描點(diǎn)判斷處理.如果d22情況下，則需要對^與尸血間的"個點(diǎn)進(jìn)行處理.假設(shè)對該"個點(diǎn)中的第/個點(diǎn)進(jìn)行處理(B"")，但在進(jìn)行處理需要知道第z'個點(diǎn)的坐標(biāo)為(&,wL+,)，從而才能知道該坐標(biāo)點(diǎn)下的像素值為。如圖3所示，在確定了&與戶血坐標(biāo)后，其在才殳影數(shù)據(jù)中的位置有兩種情況1)C與P血在水平方向上這種情況比較簡單，兩點(diǎn)在同一水平方向上，其間點(diǎn)的縱坐標(biāo)值不變的，只有橫坐標(biāo)改變.從而可知^與&s間的第!'個點(diǎn)(1"^)的坐標(biāo)(^+,，:^+,)為(4-5)2)尸5,與尸血在斜線方向上這種情況下，需要對直線尸血C進(jìn)行離散,得到其中的n個點(diǎn)的坐"=(int)(d+0.5)(4-4)標(biāo)。已知點(diǎn)/L(u血)與點(diǎn)di)，則直線的斜率A為(4-6)直線方程為^=n(4-7)&—義血jL—X/"當(dāng)|*>1時(shí)，以縱方向增減為方向，結(jié)合公式(4-7)，可以得到尸血+,的坐標(biāo)(x^，力^)為—一J卞~m當(dāng)|*|<1時(shí)，以橫方向增減為方向，結(jié)合公式(4-7)，可以得到尸血+,的坐標(biāo)&一>^+》為^1、(4-9)至此，掃描點(diǎn)與匹配點(diǎn)間的像素點(diǎn)的坐標(biāo)尸血+1，P血+2，...，Ae自求出出來了。2.5掃描點(diǎn)與匹配點(diǎn)間的像素點(diǎn)的處理從步驟2.4知道了要處理點(diǎn)尸血+,的坐標(biāo)(^+,，:^+,)了，從而也知道了該點(diǎn)下的原始的灰度值G^+,。在e步驟里來確定是否對點(diǎn)進(jìn)行處理以及怎樣處理。已知點(diǎn)&的灰度值&,點(diǎn)尸血的灰度值G^，對兩點(diǎn)進(jìn)行線性插值得到《m+,的灰度植G'G'des+i=&x"+x丄(4_io)比較G'h,與G&,的大小，如果G」+,>G"，則用G:+,賦值為處的灰度值，否則，尸血+,處的灰度值不進(jìn)行任務(wù)改變，仍為Gh,。在經(jīng)過上述步驟的處理后，完成了對MIP圖像的遍歷，得到最后的圖像數(shù)據(jù)。這樣就完成了MIP圖像處理的全部過程。下面給出本發(fā)明方法在醫(yī)學(xué)領(lǐng)域應(yīng)用的實(shí)施例圖像格式約定為通用的遵循DICOM標(biāo)準(zhǔn)的圖像，也就是說，每一幅圖像都是一幅標(biāo)準(zhǔn)的DICOM格式的圖j象。DIC0M圖像文件內(nèi)容由兩個部分組成包存參數(shù)信息的文件頭(Header)和圖點(diǎn)數(shù)據(jù)(PixelData)。DIC0M文件頭(DIC0MFileMetaInformation)包含了標(biāo)識數(shù)據(jù)集合的相關(guān)信息。每個DICOM文件都必須包括該文件頭。文件頭的最開始是文件前言，它由128字節(jié)個00H組成，接下來是DICOM前綴，它是一個長度為4字節(jié)的字符串"DICM",可以根據(jù)該值來判斷一個文件是不是DICOM文件。文件頭中還包括其它一些非常有用的信息，如文件的傳輸格式、生成該文件的應(yīng)用程序等。圖點(diǎn)數(shù)據(jù)則描述圖像的各個點(diǎn)的亮度值。DIC0M包含4個內(nèi)容層次1.Patient(病人)；2.Study(檢驗(yàn))；3.Series(系列)；4.Image(圖像)。盡管前面幾層的內(nèi)容在很多圖像里是相同的，但它們在每個圖像文件里都有。每一層叫一個信息實(shí)體(InformationEntity);每一層又細(xì)分成模塊(Module);每個模塊(Module)里面的最小單元叫^故一個屬性(Attribute)或數(shù)據(jù)元素(Element)。在DICOM文件中，每一個數(shù)據(jù)元素的位置都存放在固定的位置，因此只要知道該文件在內(nèi)存中存;^丈的首地址，就可以^4居存;J文位置的偏移量找到對應(yīng)的數(shù)據(jù)元素。而那些DIC0M中TAG值也正是一個個的數(shù)據(jù)元素。TAG(0018,5100):確定了統(tǒng)一坐標(biāo)系的方向，即所有圖像序列的各個參數(shù)都是在參考該統(tǒng)一坐標(biāo)系計(jì)算給出的(根據(jù)DICOM3.Q的標(biāo)準(zhǔn))。TAG(0020，0032):確定了該張影像首個像素("左上方")的坐標(biāo)X，Y，Z值。它和TAG(0020，0037)的值可以確定整幅圖像的所有點(diǎn)的空間坐標(biāo)。TAG(0020,0037):確定圖Y象的單位行向量和單位列向量，而兩者的外積(叉積)就是單位法向量，因此從TAG(0020,0037)可以獲取3個單位向量。TAG(0028，0030):確定圖像的每個像素代表的實(shí)際物理寬度和代表的實(shí)際物理高度。TAG(0018，0050):確定圖像切片間的物理間距。TAG(0028，OOIO):確定圖像的高度。TAG(0028，0011):確定圖像的寬度。TAG(7EF0,0010):圖像序列數(shù)據(jù)，根據(jù)DICOM標(biāo)準(zhǔn)的規(guī)定把圖像序列解析出來.由TAG(0028，0010)，TAG(0028,0011)，TAG(0028，0030)及TAG(0018，0050)就可以實(shí)現(xiàn)由TAG(0020,0037)確定的3個單位法向量上物理坐標(biāo)(單位毫米)與圖像坐標(biāo)(單位:像素)之間的轉(zhuǎn)換,從而物理坐標(biāo)與圖像對象坐標(biāo)系就可以實(shí)現(xiàn)了變換.TAG(7EF0，0010)里是進(jìn)行投影的圖像切片數(shù)據(jù)，根據(jù)上面的坐標(biāo)法矢量和所述的首個像素的三維坐標(biāo)值確定圖像在統(tǒng)一坐標(biāo)系中的平面方程.最后根據(jù)視角方向矢量可以得到過坐標(biāo)原點(diǎn)的投影平面方程，從而完成了MIP所需的統(tǒng)一坐標(biāo)系下的表示。表l對建立平面方程求取圖像序列偏移量的參數(shù)進(jìn)行了規(guī)定。表l<table>tableseeoriginaldocumentpage14</column></row><table><table>tableseeoriginaldocumentpage15</column></row><table>除了參數(shù)"圖像序列單位法向量"夕卜，其它參數(shù)都已經(jīng)獲取完成。計(jì)算圖〗象序列的單位法向量由于)f=^xr=(x,，j^,zr)x(xc，>>c,zc)，計(jì)算結(jié)果如式(4-11)所示:w=/乂A:AK;A少cZC=-m,zrxc-z入，wc-xc>g(4-12)針對表2-l中所需的參數(shù)，可以得到如式(4-13)的結(jié)果:(4-13)確定圖像序列中每幅圖像的平面方程由于對于圖像序列中的圖像而言，每一幅圖像的首個像素的坐標(biāo)位置都已經(jīng)知道為iUxmJm,^)(第z'幅圖像的首個像素坐標(biāo))。因此序列中所有圖像所在的平面方程就可以計(jì)算出來了，計(jì)算的結(jié)果如式(4-14)所示ox++=oxm+121+cz121(414)M確定投影面的平面方程當(dāng)在進(jìn)行投影，視角方向r"^，，^)是已知的，投影面垂直與視線，即投影面的法矢量就是纊，從而可以得到投影面的平面方程如式(4-15):+e^+ezZ+<i=0(4-15)利用平行投影性質(zhì)，當(dāng)采用投影面過統(tǒng)一坐標(biāo)下的原點(diǎn)不影響投影結(jié)果，則投影平面在空間坐標(biāo)系的平面方程化簡為式(4-16):^+^;+^=0(4-I6)計(jì)算圖l象序列在投影面上的偏移量圖像序列的第/幅圖像的第l個像素的空間坐標(biāo)(B"iV)A,",,即圖像序列的左上頂點(diǎn)，在投影面上的投影位置根據(jù)投影定理，并結(jié)合式(4-14)與(3-5)可以得到其在投影面上的映射點(diǎn)X'm=A/i-^(叫n+m+)Ki=A'i+m化izi)(4—17)其中圖像序列的第/幅圖像的第l個像素的空間坐標(biāo)(i《kw)^",P^,^)是滿足公式(4-14)的。圖像序列的左上頂點(diǎn)在投影面上的二維坐標(biāo)知道后，根據(jù)式(4-17)則第i幅圖像的與前一幅的偏移量(Ax,—i,A;^)為式(4-18):[AxM=(xw-u)+60w-凡—u)+c(z,"-z1M1)][AyM-凡》、[a(xin+6(凡！-凡—)+c(zin-z1M1)](4—18)至此，相鄰切片層在投影面上的偏移量都獲得了.由于該結(jié)構(gòu)包含有圖^f象序列號(SliceNumber),包含有兩個方向上的偏移量(Offset)，所以簡稱該結(jié)構(gòu)為SOffsetT(字母"T"代表Table)。這里以最常用的。++語言表示,該結(jié)構(gòu)的定義如下所示structSOffsetT{intSliceNumjshortOffset—x;shortOffset_y;};使用C++STL中的向量vector,向內(nèi)存申請一個SOffsetT類型的vector,申請的變量如下std::vector<SOffsetT>FSOffsetTable;變量FSOffsetTable中存放進(jìn)行各相鄰序列偏移量的數(shù)據(jù)。注意:圖像序列的第一個切片層的向前無相鄰層，故偏移量為(0,0),即當(dāng)SOffsetT.SliceNum=1時(shí)，SOffsetT.Offset—x=0，SOffsetT.Offset_y=0;最大密度投影得到的圖像數(shù)據(jù)對各圖像序列進(jìn)行投影，形成的MIP圖像的原始數(shù)據(jù)，由于該結(jié)構(gòu)包含有顯示在屏幕上的點(diǎn)的二維坐標(biāo)值(Position),包含投影線穿過體素的最大值(MaximumGrayPixel),還包含有該〗象素來自哪層切片(Slice),所以簡稱該結(jié)構(gòu)為MIP_PMGST(字母"T"代表Table)。該結(jié)構(gòu)的定義如下所示structMIP—PMGSTintCor一x;〃像素橫坐標(biāo)intCor//像素縱坐標(biāo)shortGray;〃像素值intSlice_num;〃像素來自的切片層};最后經(jīng)過最大密度投影后，在投影面上形成MIP圖像原始數(shù)據(jù)的Sc(圖像寬hSr(圖像高)的圖像矩陣。用0++的定義數(shù)據(jù)類型為MIP-PMGST結(jié)構(gòu)的內(nèi)存變量為"MIP一PMGSTFMIPRawData[Sr*Sc];"該數(shù)組存放了最大密度投影后每一個像素點(diǎn)的信息，比如當(dāng)對正幅圖像數(shù)據(jù)進(jìn)行遍歷時(shí)，遍歷到第n個數(shù)據(jù)點(diǎn)，則該點(diǎn)的的兩個坐標(biāo)為Pn(xn，yn):jx=FMIPRawData[n].Cor—x(4—19)l凡=FMIPRawData[n].Cor_y灰度<直為FMIPRawData[n].Gray,i亥灰度<直來自與FMIPRawData[n].幼ce—做w層，/人而為后面的處理^艮方^更地訪問到該點(diǎn)的信息。求取掃描點(diǎn)的處理匹配點(diǎn)的像素點(diǎn)坐標(biāo)已知pn點(diǎn)來自與FMIPRawData[n].幼ce—咖m層，通過查找FSOffsetTable變量中的FOffsetTable[FMIPRawData[n].57/ce—做w]知道了偏移的坐標(biāo)量為jAx'=FOffsetTable[FMIPRawData[n].57/ce_w,].Q^ef一x(4—20)'=FOffsetTable[FMIPRawData[n].S/zce—做w].Q^"—y從而求出與Pn的處理匹配點(diǎn)的像素點(diǎn)P'nOc'n，力)為J"x'=FMIPRawData[n].Cor—x+FOffsetTable[FMIPRawData[n].5//"—做/w].Q^j"—x丄y'=FMIPRawData[n].Cor_y+FOffsdTable[FMIPRawData[n].S"ce_wmm].Q^"—y(4-21)求取掃描點(diǎn)與匹配點(diǎn)間的像素點(diǎn)的坐標(biāo)已知點(diǎn)PJ、,yJ與點(diǎn)i^(x'n，力)，有公式(4-6)直線的斜率A:為&=(4-22)、-V直線方程為H"冗'(4一23)凡—>V當(dāng)間》l時(shí)，以縱方向增減為方向，結(jié)合公式(4-23),可以得到尸州的坐標(biāo)(x"+,，尺J為凡"=凡'+!'—H'—，，(4-24)比如尸"'(v，凡')，另一點(diǎn)為尸"(V+l,凡'+2)，則W=2，所以以y方向的增減為方向，其中的兩個像素點(diǎn)坐標(biāo)為(u+1),(x+2，y+2)為當(dāng)|*|<1時(shí)，以橫方向增減為方向，結(jié)合公式(4-24),可以得到《'的坐標(biāo)(U^)為=氣+fx-x',、，(4-25)K+,="+'—",(凡—凡')+少"'比如iy(V,O，另一點(diǎn)為尸"(V+2,凡'十1),則衿卜%，所以以x方向的增減為方向，其中的兩個像素點(diǎn)坐標(biāo)為(、+1，，(&+2，凡+2)為i,&+1=&+2(4_26)尸+,像素點(diǎn)的處理在計(jì)算出《+,(^+,，凡+,)后，根據(jù)式(4-10)計(jì)算出G州'G'n,^x"+G"'x丄(4—27)如果G肝,'〉FMIPRawData[n+i].Gra;;，則把尸州.，)的像素值替換為C即FMIPRawData[n+i].Gra;;=',否則不變，繼續(xù)掃4笛。在對整幅MIP圖像掃描完成后，則進(jìn)行顯示即可，至此，整個MIP圖像處理過程完成結(jié)束。圖4為現(xiàn)有技術(shù)棋盤格插值后的三維人體組織圖像；圖5為本發(fā)明對最大密度投影后得到的投影數(shù)據(jù)后處理后得到的三維人體組織圖像。通過圖4和圖5的對比可以看出，4吏用本發(fā)明方法進(jìn)行后處理得到的三維人體組織紋理細(xì)節(jié)更加清晰。本領(lǐng)域技術(shù)人員不脫離本發(fā)明的實(shí)質(zhì)和精神，可以有多種變形方案實(shí)現(xiàn)本發(fā)明，以上所述僅為本發(fā)明較佳可行的實(shí)施例而已，并非因此局限本發(fā)明的權(quán)利范圍，凡運(yùn)用本發(fā)明說明書及附圖內(nèi)容所作的等效結(jié)構(gòu)變化，均包含于本發(fā)明的4又利范圍之內(nèi)。權(quán)利要求1、一種最大密度投影圖像數(shù)據(jù)后處理方法，其特征在于，包括步驟A1、在對圖像序列進(jìn)行最大密度投影前，選取各幅圖像的同一位置的頂點(diǎn)，沿所述圖像序列遞增或遞減的順序計(jì)算并記錄相鄰兩幅圖像的所述頂點(diǎn)在投影平面上的坐標(biāo)偏移量；A2、在對所述圖像序列進(jìn)行最大密度投影后，記錄投影數(shù)據(jù)中每個象素點(diǎn)的坐標(biāo)、灰度值和所在圖像的標(biāo)識及其關(guān)系，對所述象素點(diǎn)進(jìn)行后處理，包括步驟B1、掃描所述投影數(shù)據(jù)中的象素點(diǎn)，根據(jù)所述象素點(diǎn)的坐標(biāo)和所述坐標(biāo)偏移量計(jì)算所述象素點(diǎn)所在圖像的相鄰圖像在所述投影數(shù)據(jù)中的投影匹配象素點(diǎn)的坐標(biāo)，進(jìn)而找到所述投影匹配象素點(diǎn)的灰度值；B2、計(jì)算所述象素點(diǎn)和所述投影匹配象素點(diǎn)之間的距離，若所述距離大于1，則執(zhí)行步驟B3，否則針對下一所述象素點(diǎn)，返回步驟B1；B3、計(jì)算所述象素點(diǎn)和所述投影匹配象素點(diǎn)之間的象素點(diǎn)的坐標(biāo)，進(jìn)而找到其灰度值；B4、根據(jù)公式2、根據(jù)權(quán)利要求1所述的最大密度投影圖像數(shù)據(jù)后處理方法，其特征在于所述步驟A1包括步驟All、在統(tǒng)一坐標(biāo)系中建立所述圖^f象序列的單位法向量；A12、在統(tǒng)一坐標(biāo)系中建立所述圖像序列中各幅圖像的平面方程；A13、在統(tǒng)一坐標(biāo)系中建立才殳影面的平面方程；A14、計(jì)算所述圖像序列中相鄰圖像在投影面上的偏移量。3、根據(jù)權(quán)利要求2所述的最大密度投影圖像數(shù)據(jù)后處理方法，其特征在于所述圖像序列設(shè)為醫(yī)學(xué)領(lǐng)域的數(shù)字圖像和通訊DICOM標(biāo)準(zhǔn)的圖像序列。4、根據(jù)權(quán)利要求3所述的最大密度投影圖像數(shù)據(jù)后處理方法，其特征在于根據(jù)醫(yī)學(xué)領(lǐng)域的數(shù)字圖像和通訊DICOM標(biāo)準(zhǔn)的圖像數(shù)據(jù)中的標(biāo)簽TAG值確定建立所述圖像序列的單位法向量所需的參數(shù)。5、根據(jù)權(quán)利要求4所述的最大密度投影圖像數(shù)據(jù)后處理方法，其特征在于根據(jù)醫(yī)學(xué)領(lǐng)域的數(shù)字圖像和通訊DICOM標(biāo)準(zhǔn)的圖像數(shù)據(jù)中的標(biāo)簽TAG值確定所述圖像序列中各幅圖像的平面方程所需的參數(shù)。6、根據(jù)權(quán)利要求5所述的最大密度投影圖像數(shù)據(jù)后處理方法，其特征在于根據(jù)醫(yī)學(xué)領(lǐng)域的數(shù)字圖像和通訊DICOM標(biāo)準(zhǔn)的圖像數(shù)據(jù)中的標(biāo)簽TAG值確定建立投影面的平面方程所需的參數(shù)。7、根據(jù)權(quán)利要求4至6任一所述的最大密度投影圖像數(shù)據(jù)后處理方法，其特征在于從所述圖像序列中任一選定圖像的圖像數(shù)據(jù)中取出標(biāo)簽TAG值確定相關(guān)的參數(shù)取出TAG(0018,5100)值，確定統(tǒng)一坐標(biāo)系的方向；取出TAG(0020,0032)值，確定所述選定圖像首個^f象素在所述統(tǒng)一坐標(biāo)系中的三維坐標(biāo)值；取出TAG(0020,0037)值，確定所述選定圖像的單位行向量和單位列向量，計(jì)算所述單位行向量和單位列向量的叉積得到所述圖像序列的單位法向量；取出圖像TAG(0028，0030)值，確定所述選定圖^f象的像素的實(shí)際物理寬度和實(shí)際物理高度；取出TAG(0028，OOIO)值，確定所述選定圖像的高度；取出TAG(0028，OOll)值，確定所述選定圖像的寬度。8、根據(jù)權(quán)利要求7所述的最大密度投影圖像數(shù)據(jù)后處理方法，其特征在于才艮據(jù)所述首個像素的三維坐標(biāo)值和所述單位法向量確定所述選定圖<象在所述統(tǒng)一坐標(biāo)系中的平面方程。全文摘要本發(fā)明公開了一種最大密度投影圖像數(shù)據(jù)后處理方法，包括步驟A1.在對圖像序列進(jìn)行最大密度投影前，選取各幅圖像的同一位置的頂點(diǎn)，沿所述圖像序列遞增或遞減的順序計(jì)算并記錄相鄰兩幅圖像的所述頂點(diǎn)在投影平面上的坐標(biāo)偏移量；A2.在對所述圖像序列進(jìn)行最大密度投影后，記錄投影數(shù)據(jù)中每個象素點(diǎn)的坐標(biāo)、灰度值和所在圖像的標(biāo)識及其關(guān)系，對所述象素點(diǎn)進(jìn)行后處理。文檔編號G06T11/00GK101295406SQ20081006654公開日2008年10月29日申請日期2008年4月7日優(yōu)先權(quán)日2008年4月7日發(fā)明者呂廣志,陳閩峰申請人:深圳市藍(lán)韻實(shí)業(yè)有限公司