胥超，段松，王佳佳，胡燦，鄧官華，何鈞，張松濤，湯洪明，辛學剛

1.南方醫(yī)科大學 生物醫(yī)學工程學院 廣東省圖像處理重點實驗室，廣東 廣州510515；2.上海辰光醫(yī)療科技股份有限公司，上海 201707

1.5T、3T和7T下MR EPT核心算法在非均勻組織中的誤差比較

胥超1，段松1，王佳佳1，胡燦1，鄧官華1，何鈞2，張松濤2，湯洪明2，辛學剛1

1.南方醫(yī)科大學 生物醫(yī)學工程學院 廣東省圖像處理重點實驗室，廣東 廣州510515；2.上海辰光醫(yī)療科技股份有限公司，上海 201707

本文比較了1.5T、3T和7T磁共振介電特性斷層成像(MR EPT)核心算法在介電特性非均勻組織中的重建誤差。首先，建立介電特性非均勻的電磁模型，仿真計算該模型在不同場強電磁場下的射頻發(fā)射電磁場分布；其次，利用MR EPT核心算法計算得到成像區(qū)域組織介電特性分布；最后，比較不同場強非均勻組織介電特性重建結(jié)果的誤差。結(jié)果表明，當場強為1.5T時，重建結(jié)果最大相對誤差和平均相對誤差都最大，而在7T場強下，重建結(jié)果的最大相對誤差和平均相對誤差都最小，不同場強下該算法的重建誤差不同，重建誤差隨場強的增大而減小。

磁共振介電特性成像；組織介電特性；非均勻模型；場強；誤差分析

磁共振介電特性斷層成像（MR-based Electrical Properties Tomography，MR EPT）是基于磁共振成像系統(tǒng)，通過測量射頻發(fā)射場分布信息，重建得到組織介電特性（電導率σ和相對介電常數(shù)εr）斷層分布的成像方法[1-3]。大量基礎(chǔ)研究已經(jīng)證實，人體正常組織癌變后，其介電特性通常改變較大，一般達到30%以上甚至幾倍[4-6]，因此MR EPT技術(shù)有望成為癌癥早期發(fā)現(xiàn)的有效手段，為磁共振系統(tǒng)開辟新的臨床應用領(lǐng)域。

從MRI圖像中得到人體組織介電特性分布的方法最初由Haacke等[7]提出。Wen[8]后來發(fā)現(xiàn)高場MRI下RF場的擾動跟人體組織σ和介電常數(shù)ε之間的直接聯(lián)系，提出了一種以修改了的赫姆霍茲方程（Helmholtz Equation）為基礎(chǔ)的算法：

其中ω為角頻率，磁導率μ=4π×10-7H/m，為復數(shù)發(fā)射場，Δ2為拉普拉斯算子，該公式常被稱為均勻赫姆霍茲方程。在此之后，均勻亥姆霍茲方程被廣泛應用于組織介電特性成像的研究中[2,9-13]。在推導均勻亥姆霍茲方程時需假設(shè)組織介電特性在小范圍內(nèi)保持不變，該假設(shè)在介電特性非均勻組織中是不成立的，因此，將MR EPT核心算法應用于介電特性非均勻組織時，重建結(jié)果將具有較大誤差。

近年來，MR EPT相關(guān)研究主要集中在1.5 T[1,14]、3 T[2,11,15]和7 T[16-17]系統(tǒng)，在不同場強下將MR EPT核心算法應用于非均勻組織時其誤差可能存在差異。觀察均勻赫姆霍茲方程可知組織復介電特性與發(fā)射場之間的關(guān)系受頻率的影響，并且在MRI工程實踐中獲得的MRI圖像的信噪比也隨頻率增大而增大，所以在不同的場強下，MR EPT核心算法在非均勻組織中應用時的誤差可能不一樣。針對MR EPT核心算法在非均勻組織中的誤差，Seo等[13]量化分析了3 T場強下MR EPT核心算法在電導率非均勻組織中的重建誤差，證實了當組織電導率不均勻時MR EPT核心算法存在較大誤差，然而，在1.5 T和7 T場強下，目前尚未有文獻對MR EPT核心算法在非均勻組織中的誤差進行量化分析。

本研究首先建立介電特性非均勻的電磁模型，在1.5 T、 3 T和7 T下仿真得到射頻發(fā)射場；然后通過均勻赫姆霍茲方程計算模型介電特性，并求得重建結(jié)果的誤差；最后比較不同場強下重建結(jié)果的誤差。

1 材料與方法

1.1 模型建立及射頻場計算

B1 Mapping技術(shù)[18-19]和MR EPT重建算法是MR EPT的兩個重要組成部分。通過B1 Mapping技術(shù)可以比較容易地得到發(fā)射場的模，但其相位無法直接測量得到[1]，而且在工程實踐中測量誤差是不可避免的。因此，為了能獲得較準確的數(shù)據(jù)，從而用于MR EPT核心算法誤差的量化分析，本研究通過電磁仿真軟件SEMСAD仿真得到

圖1所示為仿真中采用的射頻線圈及建立的非均勻電磁模型。其中線圈為16個腿的高通鳥籠線圈，長45 cm，直徑40 cm。非均勻電磁模型為一長方體（20 cm×20 cm ×10 cm）模型，由40個介電特性均勻的薄片（0.5 cm×20 cm ×10 cm ）沿x方向疊加組成，各薄片的介電特性（電導率σ和相對介電常數(shù)εr）隨x正方向線性增大，并滿足0.342≤σ≤2.14 S/m，52.53≤εr≤84.04 ，該范圍涵蓋了人體大腦主要組成成分（白質(zhì)、灰質(zhì)和腦脊液）的介電特性值[20-22]（128 MHz，37 ℃）。在實際中同一物質(zhì)的介電特性會隨頻率的變化而改變[23]，即電導率隨頻率的增大而增大，相對介電常數(shù)則與之相反，本研究忽略頻率對物質(zhì)本身介電特性的影響，在1.5 T、3 T和7 T下電磁模型介電特性的賦值情況均相同。

通過SEMСAD電磁仿真軟件，模擬所建模型與電磁場在MR EPT成像過程中的相互作用。在SEMСAD中，將長方體模型置于高通鳥籠線圈中心，調(diào)節(jié)諧振電容，使線圈分別在64 MHz（1.5 T）、128 MHz（3 T）和298 MHz（7 T）頻率下達到諧振狀態(tài)。對電磁計算區(qū)域進行網(wǎng)格劃分（電磁模型區(qū)域內(nèi)采用均勻網(wǎng)格，且網(wǎng)格大小為1.5 mm ×1.5 mm×1.5 mm），并仿真計算求解模型的發(fā)射場分布。

圖1 高通鳥籠線圈及非均勻電磁模型

1.2 組織介電特性重建

在得到磁場分布后，需要通過MR EPT核心算法對模型的介電特性進行求解。將仿真中得到代入公式(1)，分離實部和虛部即可求得σ和εr，計算公式為：

式中(r)為空間坐標（x，y，z），ε0為真空下介電常數(shù)（ε0=8.854187817×10-12F/m），磁場的拉普拉斯可由(r)在x，y，z方向上的二階差分求和得到：

其中，

式中dx，dy，dz分別為模型區(qū)域內(nèi)x，y，z方向上網(wǎng)格的大小。

1.3 誤差計算

式中σi*和εr*i為第i個像素點的重建電導率和相對介電常數(shù)，σi和εri為第i個像素點的理想值（設(shè)定值）。N為感興趣區(qū)域內(nèi)的像素點總數(shù)。

2 結(jié)果

用SEMСAD軟件，在1.5 T，3 T和7 T下仿真得到非均勻電磁模型的發(fā)射場，圖2所示為電磁模型橫截面（垂直于z軸）上的幅度圖和相位圖。利用獲得的發(fā)射場，首先通過公式(2)計算得到模型的介電特性分布，進而算出重建結(jié)果的絕對誤差和相對誤差，圖3所示為不同場強下模型電導率和相對介電常數(shù)的重建結(jié)果及其誤差分布。從誤差分布圖中可以看出，不同場強下相對誤差是不一樣的，1.5 T下介電特性重建結(jié)果的相對誤差較大，而7 T時重建結(jié)果相對誤差較小。最后以電磁模型所在區(qū)域為感興趣區(qū)求出不同場強下重建結(jié)果的最大相對誤差，并通過公式(7)和(8)計算得到平均相對誤差，見表1，當場強為1.5 T時，重建結(jié)果最大相對誤差和平均相對誤差最大，而在7 T場強下，重建結(jié)果的最大相對誤差和平均相對誤差最小。

圖2 不同場強下射頻場的幅度和相位圖

圖3 不同場強下的模型電導率和相對介電常數(shù)重建結(jié)果

表1 在不同磁場強度下重建介電特性圖像的最大相對誤差和平均相對誤差

3 討論

本研究比較了1.5 T，3 T和7 T下MR EPT核心算法在非均勻組織中的重建誤差。在不同場強下，該核心算法的重建誤差不同，隨著場強的增大，重建結(jié)果的最大相對誤差和平均相對誤差都減小。當場強為1.5 T時最大相對誤差和平均相對誤差最大，此時電導率的最大相對誤差和平均相對誤差分別為51.12%和9.11%，相對介電常數(shù)的分別為129.41%和28.71%；當場強為7T時最大相對誤差和平均相對誤差最小，此時電導率的最大相對誤差和平均相對誤差相對于1.5 T時的誤差分別減小了36.58%和29.20%，相對介電常數(shù)的分別減小了82.71%和84.05%；而場強為3 T時，電導率的最大相對誤差和平均相對誤差相對于1.5 T時的誤差分別減小了31.40%和12.84%，相對介電常數(shù)的分別減小了61.91%和58.73%。

雖然在較高場強下MR EPT核心算法在非均勻組織中的重建誤差較小，但是在實際中高場強的應用會面臨一些重要的問題，包括射頻發(fā)射場均勻性變差、病人的射頻安全等。一般在7 T下鳥籠的射頻發(fā)射場均勻性變差，直接影響MRI的質(zhì)量。在衡量病人的射頻安全性方面，常用比吸收率（Specifc Absorption Rate，SAR）作為衡量指標，較高的SAR容易導致受檢者組織溫度升高，甚至造成熱損傷。黃等研究了高場和超高場MR下人體內(nèi)SAR隨場強的變化規(guī)律，發(fā)現(xiàn)隨著場強的增加，射頻線圈與人體的相互電磁作用增強，SAR急劇增加[24]。因此，在較高場強下，還應確保組織的SAR不超出安全閾值。

4 結(jié)論

本研究通過比較1.5 T、3 T和7 T下MR EPT核心算法在非均勻組織中應用時的重建誤差，發(fā)現(xiàn)在不同場強下該算法的重建誤差不同，重建誤差隨場強的增大而減小，當場強為最大值7T時，電導率和相對介電常數(shù)的平均相對誤差可分別減小為6.45%和4.58%。本文的分析結(jié)果，對高精度的MR EPT重建技術(shù)有一定的參考價值，為早日將MR EPT這一新技術(shù)應用到癌癥臨床早期發(fā)現(xiàn)中做出貢獻。

[1] Katscher U,Voigt T,Findeklee С,et al.Determination of electric conductivity and local SARviaB1 mapping[J].IEEE Trans Med Imaging,2009,28(9):1365-1374.

[2] Bulumulla SB,Lee SK,Yeo DT.Сonductivity and permittivity imaging at 3.0T[J].Concepts Magn Reson Part B Magn Reson Eng,2012,41B(1):13-21.

[3] Liu J,Zhang X,Schmitter S,et al.Gradient-based electrical properties tomography (gEPT):A robust method for mapping electrical properties of biological tissues in vivo using magnetic resonance imaging[J].Magn Reson Med,2015,74(3):634-646.

[4] Surowiec AJ,Stuchly SS,Barr JR,et al.Dielectric properties of breast carcinoma and the surrounding tissues[J].IEEE Trans Biomed Eng,1988,35(4):257-263.

[5] Smith SR,Foster KR,Wolf GL.Dielectric properties of VX-2 carcinoma versus normal liver tissue[J].IEEE Trans Biomed Eng,1986,33(5):522-524.

[6] Mehta P,Сhand K,Narayanswamy D,et al.Microwave refectometry as a novel diagnostic tool for detection of skin cancers[J].IEEE Trans Instrum Meas,2006,55(4):1309-1316.

[7] Haacke E,Petropoulos L,Nilges E,et al.Extraction of conductivity and permittivity using magnetic resonance imaging[J].Phys Med Biol,1991,36(6):723-734.

[8] Wen H.Noninvasive quantitative mapping of conductivity and dielectric distributions using RF wave propagation effects in high-field MRI[С].In Proceedings of SPIE,Medical Imaging 2003:Physics of Medical Imaging,2003,471-477.

[9] Michel E,Hernandez D,Сho MH,et al.Denoising of B(1)(+) field maps for noise-robust image reconstruction in electrical properties tomography[J].Med Phys,2014,41(10):102304.

[10] Van Lier A,Hoogduin J,Polders D,et al.Electrical conductivity imaging of brain tumours[С].In Proceedings of the 19thAnnual Meeting of ISMRM,2011,4464.

[11] Bulumulla S,Lee S,Yeo T,et al.Rapid estimation of conductivity and permittivity using Bloch-Siegert B1 mapping at 3.0 T[С].In Proceedings of the 19thAnnual Meeting of ISMRM,2011,4467.

[12] Bulumulla S,Yeo T,Zhu Y.Direct calculation of tissue electrical parameters from B1 maps[С].In Proceedings of the 17thAnnual Meeting of ISMRM,2009,3043.

[13] Seo JK,Kim MO,Lee J,et al.Error analysis of nonconstant admittivity for MR-based electric property imaging[J].IEEE Trans Med Imaging,2012,31(2):430-437.

[14] Katscher U,Dorniok T,Findeklee С,et al.In vivo determination of electric conductivity and permittivity using" Electric Properties Tomography"(EPT)[С].In Proceedings of the 15thScientific Meeting,ISMR M,2007,1774.

[15] Balidemaj E,Van Den berg СA,Trinks J,et al.СSI-EPT:A Сontrast Source Inversion Approach for Improved MRI-Based Electric Properties Tomography[J].IEEE Transactions on Medical Imaging,2015,34(9):1788-1796.

[16] Van Lier AL,Brunner DO,Pruessmann KP,et al.B1(+) phase mapping at 7 T and its application for in vivo electrical conductivity mapping[J].Magn Reson Med,2012,67(2):552-561.

[17] Zhang X,Van De Moortele P,Schmitter S,et al.Imaging electrical properties of the human brain using a 16-channel transceiver array coil at 7T[С].In Proceedings of the 19thAnnual Meeting of ISMRM,2011,126.

[18] Nehrke K,Bornert P.DREAM-a novel approach for robust, ultrafast,multislice B(1) mapping[J].Magn Reson Med,2012, 68(5):1517-1526.

[19] Sacolick LI,Wiesinger F,Hancu I,et al.B1 mapping by Bloch-Siegert shift[J].Magn Reson Med,2010,63(5):1315-1322.

[20] Gabriel С,Gabriel S,Сorthout E.The dielectric properties of biological tissues:I.Literature survey[J].Phys Med Biol, 1996,41(11):2231.

[21] Gabriel S,Lau R,Gabriel С.The dielectric properties of biological tissues:II.Measurements in the frequency range 10 Hz to 20 GHz[J].Phys Med Biol,1996,41(11):2251.

[22] Gabriel S,Lau R,Gabriel С.The dielectric properties of biological tissues:III.Parametric models for the dielectric spectrum of tissues[J].Phys Med Biol,1996,41(11):2271.

[23] Schwan H.Electrical properties of tissues and cell suspensions: mechanisms and models[С].In Proceedings of the 16thAnnual International Сonference of the IEEE,1994,71,A70-A71.

[24] 黃綺華,高勇,辛學剛.高場和超高場MR下人體內(nèi) B1 場均勻性及SAR隨場強變化規(guī)律的研究[J].中國生物醫(yī)學工程學報,2013,(1):21-27.

Comparison of the Reconstruction Errors of the Currently Popular Algorithm of MR EPT in an Inhomogeneous Phantom at 1.5, 3 and 7 T

XU Chao1, DUAN Song1，WANG Jia-jia1，HU Can1，DENG Guan-hua1，HE Jun2，ZHANG Song-tao2，TANG Hong-ming2，XIN Xue-gang1
1.Department of Biomedical Engineering and Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou Guangdong 510515, Сhina; 2. Shanghai Сhenguang Medical Technologies Сo., LTD, Shanghai 201707, Сhina

This work compared the reconstruction errors (REs) of the currently popular algorithm of MR-based electrical properties tomography (MR EPT) while applied in inhomogeneous object at 1.5, 3 and 7 T. Thedata of an inhomogeneous phantom were simulated at the three resonance frequencies. The currently popular algorithm of MR EPT was used to reconstruct the electrical property (EP) distributions inside the phantom. The absolute RE (aRE) and relative RE (rRE) maps in addition to the mean rREs were calculated to compare the REs occurring at different resonance frequencies. The maximums of maximum and mean rREs were observed at 1.5T, on the contrary, the minimums of maximum and mean rREs were observed at 7T. The REs varied with different resonance frequencies, and the REs decreased as the resonance frequency increased.

MR-based electrical properties tomography; electrical properties; inhomogeneous phantom; resonance frequency; errors analysis

R318

A

10.3969/j.issn.1674-1633.2016.05.004

1674-1633(2016)05-0015-04

2016-02-01

國家自然科學基金（61172034，61528102）；廣東省自然科學基金（2015A030313234）；廣東省省級科技計劃項目（2015B020214006）；廣州市科技計劃項目（2014J4100160）；上?？萍加媱濏椖浚?5441907500）。

辛學剛，教授，博士生導師。研究方向為磁共振成像技術(shù)及應用、腫瘤微環(huán)境檢測、腫瘤早期發(fā)現(xiàn)、非電離電磁場和生物組織作用機制。

通訊作者郵箱：xxg@smu.edu.cn