本發(fā)明涉及，特別是一種基于自適應(yīng)卷積核的磁共振相位圖的背景場(chǎng)去除方法。
背景技術(shù)：
：磁化率定義為物質(zhì)放入外磁場(chǎng)后的磁敏感性反應(yīng)，是物質(zhì)的固有屬性。磁共振成像(magneticresonanceimaging,mri)中，每一種物質(zhì)在放入磁場(chǎng)后都可以得到一定程度的磁化，磁化與磁場(chǎng)大小及組織的磁化率成正比。如果施加一個(gè)足夠長(zhǎng)的te，自旋頻率不同的質(zhì)子間將形成明顯的相位差別，從而在相位圖上區(qū)別出磁敏感性不同的組織。定量磁化率成像(quantitativesusceptibilitymapping,qsm)利用梯度回波數(shù)據(jù)的相位信息產(chǎn)生組織的磁場(chǎng)特性圖。qsm通過(guò)求解關(guān)于感應(yīng)磁場(chǎng)分布的病態(tài)逆問(wèn)題，由測(cè)得的相位圖像導(dǎo)出潛在的組織體磁化率的定量圖，在醫(yī)學(xué)磁共振成像領(lǐng)域具有廣闊的應(yīng)用前景。qsm技術(shù)有助于分辨并量化特定生物標(biāo)記，如鐵、鈣、釓和超順磁性氧化鐵納米顆粒(spio)，已經(jīng)廣泛應(yīng)用于腦損傷，多發(fā)性硬化癥、多種神經(jīng)退行性疾病、骨骼礦化以及動(dòng)脈粥樣硬化的研究中，對(duì)此類疾病的臨床診斷具有重要意義。qsm還能夠?qū)铙w組織神經(jīng)束和腦白質(zhì)纖維進(jìn)行非入侵式檢測(cè)，為神經(jīng)纖維成像研究提供新的對(duì)比機(jī)制，這對(duì)于神經(jīng)影像學(xué)定量的連接性研究和生物物理研究非常重要。利用相位信息計(jì)算定量磁化率圖之前，必須通過(guò)預(yù)處理對(duì)相位圖像進(jìn)行校正，包括相位解纏繞和去除背景場(chǎng)。磁化率成像對(duì)于局部磁場(chǎng)不均勻性特別敏感，在某些磁化率差異特別大的區(qū)域，如顱底的含氣鼻竇部位，會(huì)造成局部特別強(qiáng)的相位偽影，使我們無(wú)法有效觀察及利用感興趣區(qū)的相位信息。傳統(tǒng)的高通濾波器方法假設(shè)背景場(chǎng)是緩慢變化的，雖然可以較好的去除背景低頻相位擾動(dòng)，但是可能從大的解剖結(jié)構(gòu)上將一些生理和病理相關(guān)的相位信息去除，尤其在嚴(yán)重磁場(chǎng)不均勻的區(qū)域，并不能完全濾除背景的相位改變，從而殘留大量的背景相位。因此，在應(yīng)用相位圖之前如何對(duì)其進(jìn)行有效的相位校正是定量磁化率成像的首要關(guān)鍵問(wèn)題。受物理原理啟發(fā)，近來(lái)提出的偶極場(chǎng)投影方法(pdf)建立在hilbert空間的投影理論之上，即voi內(nèi)部和外部的組織磁化率源產(chǎn)生的磁場(chǎng)近似正交，該方法去除背景場(chǎng)的有效性優(yōu)于傳統(tǒng)方法，但是在組織的邊界或是組織—空氣的交界面仍存在問(wèn)題。因?yàn)榉椒ǖ幕炯僭O(shè)的有效性在邊界處受到限制，而且難以區(qū)分邊界處的磁場(chǎng)變化來(lái)自內(nèi)部的組織磁化率還是外部的組織源。相位數(shù)據(jù)復(fù)雜諧波偽影抑制算法(sharp)假設(shè)voi外部的背景在voi內(nèi)滿足調(diào)和函數(shù)性質(zhì)，可用歸一化的球均值卷積核求解。由于采用球均值卷積核，sharp方法在voi邊緣區(qū)域的體素因與voi外的無(wú)效信號(hào)做卷積運(yùn)算而導(dǎo)致邊緣模糊，模糊點(diǎn)的數(shù)量正比于所用球體的直徑，需要被剔除。實(shí)際中smv濾波的過(guò)程中邊界常存在殘余相位誤差，這是由于空氣—組織的交界處高磁化率差異導(dǎo)致的。這個(gè)殘余誤差在去卷積濾波的過(guò)程中被進(jìn)一步放大，放大的尺度與球體的直徑成反比。可見(jiàn)，采用大直徑的球均值卷積核會(huì)減小相位誤差，但是不足在于導(dǎo)致剔除較多邊界點(diǎn)。針對(duì)sharp存在的上述缺陷，研究人員提出變卷積核半徑sharp法(v-sharp)，正則化sharp法(resharp)和拉普拉斯方程邊界值法(lbv)等改進(jìn)方法。磁化率映射的不準(zhǔn)確主要是由感興趣區(qū)域外的磁化率引起的背景場(chǎng)分布不均勻?qū)е碌模畿|干和頭部以及組織—組織和空氣—組織的交界面。相位處理過(guò)程中，常對(duì)變化比較激烈的背景場(chǎng)和組織結(jié)構(gòu)交界處的磁化率差異估計(jì)不足，故上述的背景場(chǎng)去除方法均無(wú)法有效的抑制強(qiáng)磁化率變化引入的相位偏差。由于腦組織與鼻竇或顱骨的交界處存在很強(qiáng)的磁化率變化現(xiàn)象，為了降低對(duì)后續(xù)qsm重建準(zhǔn)確性的影響，現(xiàn)有方法使用掩模板將該區(qū)域的腦組織從voi中去除以減少相位殘余。然而組織結(jié)構(gòu)交界處往往含有重要的醫(yī)學(xué)信息，因算法的局限性而導(dǎo)致腦組織完整性的缺失，將使我們無(wú)法獲取完整的大腦診斷信息，這不利于qsm技術(shù)的臨床應(yīng)用。因此，在去除強(qiáng)磁化率變化區(qū)域的腦組織背景場(chǎng)時(shí)，現(xiàn)有方法仍存在不足，有待進(jìn)一步提高。技術(shù)實(shí)現(xiàn)要素：本發(fā)明的主要目的在于克服現(xiàn)有技術(shù)中的上述缺陷，根據(jù)相位的水平集能量值，構(gòu)造尺度可調(diào)的高斯核函數(shù)，提出一種基于自適應(yīng)卷積核的磁共振相位圖的背景場(chǎng)去除方法。本發(fā)明采用如下技術(shù)方案：一種基于自適應(yīng)卷積核的磁共振相位圖的背景場(chǎng)去除方法，其特征在于：利用水平集函數(shù)創(chuàng)建解纏繞相位圖的能量泛函；根據(jù)求解得到的相位水平集能量提取出磁化率區(qū)域變化的顯著度，并逐個(gè)體素的創(chuàng)建自適應(yīng)高斯卷積核；采用自適應(yīng)高斯卷積核去除背景場(chǎng)。優(yōu)選的，所述利用水平集函數(shù)創(chuàng)建解纏繞相位圖的能量泛函，其模型為其中：ψ為水平集函數(shù)，f1(x)和f2(x)為灰度擬合函數(shù)，αi是權(quán)系數(shù)，y是以x為中心的局部圖像域中任意點(diǎn)的坐標(biāo)；i(y)表示點(diǎn)y的灰度值，s為局部圖像域的尺度，由二維高斯核函數(shù)ks定義；為弧長(zhǎng)約束項(xiàng)，為水平集正則項(xiàng)，μ和ν為權(quán)常數(shù)，和div分別為梯度算子和散度算子，d為單位沖擊函數(shù)。優(yōu)選的，所述的采用能量泛函求解得到的相位水平集函數(shù)值，構(gòu)造自適應(yīng)的高斯卷積核，其定義為其中：標(biāo)準(zhǔn)差σ由當(dāng)前體素的水平集函數(shù)值ψ(x,y,z)和該點(diǎn)所處層面的水平集均值ψ0(z)之比調(diào)制，表示為σ(x,y,z)＝rψ0(z)/cψ(x,y,z)；常數(shù)c用于調(diào)節(jié)比值ψ0(z)/ψ(x,y,z)與卷積核標(biāo)準(zhǔn)差σ之間的關(guān)系。優(yōu)選的，根據(jù)水平集函數(shù)值將感興趣區(qū)域劃分為不均勻背景場(chǎng)和均勻背景場(chǎng)兩個(gè)區(qū)域，采用不同的卷積核去除背景場(chǎng)，其模型為其中：為單位沖擊函數(shù)，bloc為局部場(chǎng)，為中間變量，ρsgk為自適應(yīng)高斯卷積核，ρsmv為球均值核。優(yōu)選的，所述磁化率變化顯著的區(qū)域?yàn)椴痪鶆虮尘皥?chǎng)，采用自適應(yīng)高斯卷積核ρsgk去除背景場(chǎng)；而磁化率變化不明顯的區(qū)域?yàn)榫鶆虮尘皥?chǎng)，采用球均值核ρsmv作卷積去除背景場(chǎng)。由上述對(duì)本發(fā)明的描述可知，與現(xiàn)有技術(shù)相比，本發(fā)明具有如下有益效果：1.本發(fā)明采用水平集函數(shù)創(chuàng)建解纏繞相位圖的能量泛函，根據(jù)求解得到的相位水平集能量提取出磁化率顯著變化的區(qū)域，并逐個(gè)體素的創(chuàng)建自適應(yīng)高斯卷積核。卷積核自適應(yīng)可調(diào)的權(quán)重和半徑能夠有效去除背景場(chǎng)，抑制強(qiáng)磁化率變化導(dǎo)致的偽影，防止邊緣信息丟失。2.本發(fā)明的方法根據(jù)各體素的水平集函數(shù)值將感興趣區(qū)域劃分為不均勻背景場(chǎng)和均勻背景場(chǎng)兩個(gè)區(qū)域，采用不同的卷積核去除背景場(chǎng)。用于人腦定量磁化率成像時(shí)，本方法能夠獲得準(zhǔn)確局部場(chǎng)的同時(shí)有效保留鼻竇周圍組織以及腦組織與顱骨交界處的結(jié)構(gòu)完整性，其效果明顯優(yōu)于現(xiàn)有的方法。附圖說(shuō)明圖1為活體人腦磁共振數(shù)據(jù)圖。(a)解纏繞相位圖，(b)水平集函數(shù)值圖，(c)局部場(chǎng)圖，(d)磁化率圖。圖2為數(shù)值仿真數(shù)據(jù)實(shí)驗(yàn)圖。(a)仿真的局部場(chǎng)圖，(b)v-sharp計(jì)算得到的局部場(chǎng)圖，(c)resharp計(jì)算得到的局部場(chǎng)圖，(d)r-sharp計(jì)算得到的局部場(chǎng)圖，(e)仿真的磁化率圖，(f)由v-sharp得到的局部場(chǎng)重建磁化率圖，(g)由resharp得到的局部場(chǎng)重建磁化率圖，(h)由r-sharp得到的局部場(chǎng)重建磁化率圖。圖3為釓溶液體模數(shù)據(jù)實(shí)驗(yàn)圖。(a)含有相位殘余的局部場(chǎng)圖，(b)v-sharp計(jì)算得到的局部場(chǎng)圖，(c)resharp計(jì)算得到的局部場(chǎng)圖，(d)r-sharp計(jì)算得到的局部場(chǎng)圖。圖4為活體人腦去除背景場(chǎng)后的局部場(chǎng)圖，包含橫截面，矢狀面和冠狀面三個(gè)方向。(a)v-sharp局部場(chǎng)圖，(b)resharp局部場(chǎng)圖，(c)r-sharp局部場(chǎng)圖。圖5為活體人腦重建出的磁化率圖，包含橫截面，矢狀面和冠狀面三個(gè)方向。(a)v-sharp對(duì)應(yīng)的磁化率圖，(b)resharp對(duì)應(yīng)的磁化率圖，(c)r-sharp對(duì)應(yīng)的磁化率圖,(d)多方向cosmos磁化率圖。具體實(shí)施方式以下通過(guò)具體實(shí)施方式對(duì)本發(fā)明作進(jìn)一步的描述。一種基于自適應(yīng)卷積核的磁共振相位圖的背景場(chǎng)去除方法，利用水平集函數(shù)創(chuàng)建解纏繞相位圖的能量泛函，其模型為水平集函數(shù)ψ演化的主要驅(qū)動(dòng)能量是相位圖中各點(diǎn)的局部能量項(xiàng)在圖像域的積分，局部能量項(xiàng)表征局部灰度值與相應(yīng)的灰度擬合函數(shù)f1(x)和f2(x)的近似度。αi是權(quán)系數(shù)，y是以x為中心的局部圖像域中任意點(diǎn)的坐標(biāo)，i(y)表示點(diǎn)y的灰度值，局部圖像域的尺度為s，由二維高斯核函數(shù)ks定義。為了提高水平集的演化速度以及穩(wěn)定性，引入弧長(zhǎng)約束項(xiàng)和水平集正則項(xiàng)其中μ和ν為權(quán)常數(shù)，和div分別為梯度算子和散度算子，δ為單位沖擊函數(shù)。圖像的能量泛函具有很強(qiáng)的克服圖像灰度不均勻性的特點(diǎn)，因此將能量泛函應(yīng)用于解纏繞后的人腦相位圖，當(dāng)由局部能量驅(qū)動(dòng)的水平集曲線停留在同質(zhì)區(qū)域邊界時(shí)，水平集函數(shù)可以準(zhǔn)確區(qū)分出不同磁化率值的組織，克服了相位圖像對(duì)比度低的問(wèn)題。而后，根據(jù)求解得到的相位水平集能量提取出磁化率區(qū)域變化的顯著度，并逐個(gè)體素的創(chuàng)建自適應(yīng)高斯卷積核，采用能量泛函求解得到的相位水平集函數(shù)值，構(gòu)造自適應(yīng)的高斯卷積核，定義為通過(guò)水平集函數(shù)值檢測(cè)出磁化率顯著變化的區(qū)域，逐個(gè)體素的創(chuàng)建高斯卷積核，標(biāo)準(zhǔn)差σ由當(dāng)前體素的水平集函數(shù)值ψ(x,y,z)和該點(diǎn)所處層面的水平集均值ψ0(z)之比調(diào)制，表示為σ(x,y,z)＝rψ0(z)/cψ(x,y,z)；其中常數(shù)c用于調(diào)節(jié)比值ψ0(z)/ψ(x,y,z)與卷積核標(biāo)準(zhǔn)差σ之間的關(guān)系。根據(jù)水平集函數(shù)值將感興趣區(qū)域劃分為不均勻背景場(chǎng)和均勻背景場(chǎng)兩個(gè)區(qū)域，采用不同的卷積核去除背景場(chǎng)，其模型為其中δ為單位沖擊函數(shù)，bloc為局部場(chǎng)，為中間變量。卷積核半徑隨著體素逼近感興趣區(qū)域(volumeofinterest,voi)邊緣而減小，卷積核的權(quán)重由水平集函數(shù)值調(diào)制。卷積核自適應(yīng)可調(diào)的權(quán)重和半徑能夠有效抑制感興趣區(qū)域外的組織或外部信號(hào)對(duì)卷積結(jié)果的影響，從而有效去除區(qū)域的背景場(chǎng)，抑制強(qiáng)磁化率變化導(dǎo)致的偽影，防止邊緣信息丟失，保留由組織產(chǎn)生的局部場(chǎng)信息。磁化率變化顯著的區(qū)域背景場(chǎng)分布不均勻，采用自適應(yīng)高斯卷積核ρsgk去除背景場(chǎng)；而磁化率變化不明顯的區(qū)域背景場(chǎng)分布均勻，采用球均值核ρsmv作卷積去除背景場(chǎng)。利用水平集函數(shù)值逐個(gè)體素的創(chuàng)建自適應(yīng)卷積核，當(dāng)前體素的水平集函數(shù)值越高則核的標(biāo)準(zhǔn)差越小，即卷積核的中心點(diǎn)權(quán)重越大，而周圍點(diǎn)權(quán)重衰減的越快速。卷積核半徑隨著體素逼近感興趣區(qū)域邊緣而減小，卷積核的權(quán)重由水平集函數(shù)值調(diào)制。卷積核自適應(yīng)可調(diào)的權(quán)重和半徑能夠有效抑制感興趣區(qū)域外的組織或信號(hào)對(duì)卷積結(jié)果的影響，從而有效去除背景場(chǎng)，抑制強(qiáng)磁化率變化導(dǎo)致的偽影，防止邊緣信息丟失。例如，人腦中鼻竇附近組織的背景場(chǎng)強(qiáng)度通常比局部場(chǎng)高一到二個(gè)數(shù)量級(jí)，通過(guò)水平集函數(shù)值調(diào)制高斯核的中心權(quán)重使之接近于1，可有效地消除背景場(chǎng)成分，減少相位殘余。另外，腦組織與顱骨交界處，采用調(diào)制的高中心權(quán)值卷積核，因其周圍點(diǎn)權(quán)值接近于0而能夠在卷積運(yùn)算時(shí)很好地抑制voi外部的無(wú)效信號(hào)，防止腦組織邊緣區(qū)域被模糊化，保護(hù)voi的完整性。應(yīng)用舉例解纏繞后的活體人腦磁共振相位圖及其對(duì)應(yīng)的水平集函數(shù)值如圖1所示，水平集函數(shù)值圖在鼻竇周圍以及腦組織與顱骨交界處表現(xiàn)出明顯的等級(jí)差異，這有助于更好地設(shè)計(jì)自適應(yīng)卷積核，針對(duì)區(qū)域不同程度的磁化率變化，獲得不同的高斯核權(quán)重分布，從而取得自適應(yīng)的背景場(chǎng)去除效果，有效抑制磁化率偽影。為驗(yàn)證本發(fā)明提出的基于自適應(yīng)卷積核的背景場(chǎng)去除方法r-sharp能夠在強(qiáng)磁化率變化區(qū)域中有效地去除背景場(chǎng)、獲得準(zhǔn)確的局部場(chǎng)，選取v-sharp和resharp作為背景場(chǎng)去除對(duì)比算法，采用數(shù)值仿真數(shù)據(jù)、體模仿真數(shù)據(jù)和活體人腦數(shù)據(jù)分別進(jìn)行試驗(yàn)驗(yàn)證。對(duì)比實(shí)驗(yàn)中三種算法卷積核的半徑均初始化為6；隨著體素靠近voi邊緣，半徑逐漸減小；當(dāng)體素位于voi邊緣上時(shí)，核半徑等于1。為進(jìn)一步驗(yàn)證結(jié)果的準(zhǔn)確性，將得到的局部場(chǎng)圖作磁化率反演計(jì)算其定量磁化率圖，比較估計(jì)出的磁化率值與參考值的差異程度。采用均方根誤差(rmse)，平均結(jié)構(gòu)相似度(mssim)和組織的磁化率均值來(lái)定量評(píng)價(jià)實(shí)驗(yàn)結(jié)果。首先進(jìn)行數(shù)值仿真實(shí)驗(yàn)。創(chuàng)建一個(gè)128ⅹ128ⅹ64的大橢球體用于仿真人腦，在大橢球體中放置四個(gè)小的橢球體分別仿真鼻竇、血管、蒼白球和尾狀核，相應(yīng)的磁化率值為9.4ppm，0.3ppm，0.1ppm和0.05ppm，如圖2e所示。在仿真圖中加入與幅值圖比值滿足snr＝40的高斯白噪聲。r-sharp算法中水平集函數(shù)的參數(shù)經(jīng)過(guò)優(yōu)化設(shè)置為α1＝1,α2＝1000,μ＝0.5,v＝0.65,s＝5，自適應(yīng)卷積核的調(diào)整參數(shù)c＝1.8；v-sharp算法中奇異值分解的截?cái)嚅撝祎svd＝0.05；resharp算法的正則化參數(shù)λ＝0.05。三種方法估算出來(lái)的局部場(chǎng)與磁化率反演結(jié)果分別如圖2(b-d)及圖2(f-h)所示。在局部場(chǎng)圖中，v-sharp和resharp兩種方法在仿真鼻竇周圍存在很強(qiáng)的相位殘余，如箭頭所示處，這些相位殘余相應(yīng)地在重建的磁化率圖上產(chǎn)生嚴(yán)重偽影。相反，r-sharp方法在抑制相位殘余上具有明顯優(yōu)勢(shì)。實(shí)驗(yàn)結(jié)果的定量評(píng)價(jià)結(jié)果如表1和表2所示，可見(jiàn)本發(fā)明的方法各項(xiàng)評(píng)級(jí)指標(biāo)均顯著優(yōu)于對(duì)比方法。表1：數(shù)值仿真數(shù)據(jù)實(shí)驗(yàn)結(jié)果的評(píng)價(jià)指標(biāo)對(duì)比度量v-sharpresharpr-sharprmse0.0110.1740.070mssim0.6110.7810.926表2：數(shù)值仿真數(shù)據(jù)實(shí)驗(yàn)結(jié)果的磁化率值對(duì)比真實(shí)值v-sharpresharpr-sharp0.050.035±0.0250.042±0.0240.041±0.0090.100.139±0.0420.154±0.0430.107±0.0100.300.304±0.0050.308±0.0030.298±0.007其次進(jìn)行釓溶液體模數(shù)據(jù)實(shí)驗(yàn)驗(yàn)證。實(shí)驗(yàn)數(shù)據(jù)在3t飛利浦人體成像儀上采集得到，采用三維梯度回波序列，成像參數(shù)為te/δte/tr＝3/3/35ms，回波個(gè)數(shù)10，成像視野為120×120mm×120mm，層厚1.5mm。采用同樣的優(yōu)化方法，r-sharp算法中水平集函數(shù)的參數(shù)設(shè)置為α1＝α2＝1,μ＝0.5,v＝0.65,s＝5，自適應(yīng)卷積核的調(diào)整參數(shù)c＝1.8；v-sharp算法中奇異值分解的截?cái)嚅撝祎svd＝0.12；resharp算法的正則化參數(shù)λ＝0.004。實(shí)驗(yàn)分別采集含有兩支釓溶液試管的體模與取出釓溶液試管的體模，通過(guò)laplacian相位解纏繞和多回波線性擬合得到相位數(shù)據(jù)。將兩次采集的相位數(shù)據(jù)做差值，得到的局部場(chǎng)圖作為參照用的標(biāo)準(zhǔn)場(chǎng)圖，如圖3a所示，在釓溶液與空氣交界處由于強(qiáng)烈的磁化率變化，存在很明顯的背景場(chǎng)相位殘余，如箭頭所示。采用v-sharp方法無(wú)法完全將該相位殘余消除，如圖3b所示；采用resharp方法，對(duì)殘余相位的抑制作用優(yōu)于v-sharp方法；然而，使用本發(fā)明方法r-sharp則能很好地消除相位殘余，顯著優(yōu)于其他二者。與標(biāo)準(zhǔn)場(chǎng)圖對(duì)照，三種方法v-sharp，resharp和r-sharp獲得的局部場(chǎng)圖rmse值分別為2.674，2.167和2.155。最后進(jìn)行活體人腦數(shù)據(jù)實(shí)驗(yàn)驗(yàn)證。實(shí)驗(yàn)數(shù)據(jù)在7t飛利浦人體成像儀上，使用32通道的頭部線圈采集得到。采用三維梯度回波序列，成像參數(shù)為tr＝45ms，te1＝2ms，δte＝2ms，回波個(gè)數(shù)8，成像視野為220mm×220mm×110mm，層厚1mm，數(shù)據(jù)矩陣為224×224×110。采用同樣的優(yōu)化方法，r-sharp算法中水平集函數(shù)的參數(shù)設(shè)置為α1＝α2＝1,μ＝0.5,v＝0.65,s＝5，自適應(yīng)卷積核的調(diào)整參數(shù)c＝1.8；v-sharp算法中奇異值分解的截?cái)嚅撝祎svd＝0.12；resharp算法的正則化參數(shù)λ＝0.005。圖4列出了v-sharp，resharp和r-sharp三種方法的去背景場(chǎng)結(jié)果，由三個(gè)方向分別示出。由圖可見(jiàn)，在鼻竇周圍的腦組織(如箭頭所示)，v-sharp方法無(wú)法正確地除去背景場(chǎng)的干擾，表現(xiàn)為高亮的光斑，掩蓋了組織的真實(shí)結(jié)構(gòu)信息。resharp方法較v-sharp雖有提升，但仍存在明顯的相位殘余，而r-sharp方法則很好地去除了相位殘余信息，復(fù)原了鼻竇處組織的真實(shí)信息。采用v-sharp，resharp估計(jì)出的局部場(chǎng)圖進(jìn)行磁化率圖像重建，相位殘余將導(dǎo)致陰影狀的磁化率偽影，從而無(wú)法分辨該區(qū)域的腦組織結(jié)構(gòu)細(xì)節(jié)，如圖5a-b。然而，采用r-sharp的局部場(chǎng)圖重建磁化率圖，抑制偽影優(yōu)勢(shì)明顯，能夠辨別鼻竇處組織結(jié)構(gòu)細(xì)節(jié)，重建結(jié)果更接近于由多方向采樣數(shù)據(jù)計(jì)算出的cosmos定量磁化率圖。三種方法v-sharp，resharp和r-sharp獲得的磁化率圖rmse值分別為0.017,0.015和0.011；mssim值分別為r-sharp0.861，v-sharp0.813和resharp0.775。上述僅為本發(fā)明的具體實(shí)施方式，但本發(fā)明的設(shè)計(jì)構(gòu)思并不局限于此，凡利用此構(gòu)思對(duì)本發(fā)明進(jìn)行非實(shí)質(zhì)性的改動(dòng)，均應(yīng)屬于侵犯本發(fā)明保護(hù)范圍的行為。當(dāng)前第1頁(yè)12