郭悅，武漢大學(xué)人民醫(yī)院，湖北武漢 430060

審校：許峰，北京大學(xué)第三醫(yī)院 醫(yī)學(xué)工程處，北京 100191

0 前言

肺通過交替的通氣和灌注進(jìn)行氣體交換。許多疾病，如肺功能受損，以及對治療反應(yīng)的監(jiān)測都采用專門的成像技術(shù)。因?yàn)榉位盍康臏y定方法可能無法單方面地描繪出功能上的變化[1]。平面顯像和單光子發(fā)射式計(jì)算機(jī)斷層儀（SPECT）以放射性標(biāo)記的示蹤劑為基礎(chǔ)進(jìn)行三維成像，是一個(gè)廣泛被接受的肺灌注和通氣的臨床測量標(biāo)準(zhǔn)[2-3]。核成像方式的主要缺點(diǎn)是空間和時(shí)間分辨率比較差。通氣評估是吸入氙氣[4-5]后，通過計(jì)算機(jī)斷層掃描（CT）灌注成像，并使用雙能碘增強(qiáng)CT[6]提供高的空間分辨率。但這項(xiàng)檢查是在電離輻射下進(jìn)行的，所以限制了兒童和青少年的后續(xù)檢查次數(shù)。

盡管磁共振（MR）成像是一項(xiàng)眾所周知的物理技術(shù)，低質(zhì)子密度組織的可視化存在困難，但它為專門評價(jià)肺生理學(xué)相對或絕對的狀況提供了廣闊的頻譜相對法。基于非質(zhì)子的肺癌磁共振成像，通過患者吸入作為示蹤劑的極化氣體（氦-3，氙-129），得到高質(zhì)量的通氣的圖像[7-8]。然而，臨床上常規(guī)應(yīng)用的磁共振成像氣體和硬件的成本高，人員需要經(jīng)過訓(xùn)練，使其可用性很有限?；谄渌椒ㄟ€有通過的質(zhì)子磁共振成像使用縮短示蹤時(shí)間的順磁性T1造影劑來增強(qiáng)肺部的MR信號：增氧磁共振成像和霧化釓螯合物磁共振成像[9-10]。

使用動態(tài)對比示蹤材料或順磁性對比劑可增強(qiáng)磁共振成像，灌注測量可以運(yùn)用這項(xiàng)技術(shù)完成。這種方法提供了高空間和時(shí)間分辨率，并能估計(jì)各種血流動力學(xué)參數(shù)[11]。然而，注射釓的造影劑會引起某些相關(guān)的罕見病癥，如嚴(yán)重的過敏反應(yīng)或存在潛在的腎系統(tǒng)纖維化的風(fēng)險(xiǎn)。作為灌注測量的替代技術(shù)，動脈自旋可在血液中使用水的質(zhì)子作為一種內(nèi)源性的造影劑。但結(jié)果會導(dǎo)致信噪比低，時(shí)間分辨率不高以及信號的快速衰變[12-13]。

近日，平掃傅立葉分解（FD）技術(shù)的磁共振成像技術(shù)正在研究中。有學(xué)者提議，通過1個(gè)單一的非門MR圖像系列[14-16]獲得區(qū)域肺灌注和通氣的相關(guān)信息。該技術(shù)是利用肺部組織生理過程中某些區(qū)域MR信號強(qiáng)度的變化來成像的。在吸氣時(shí)，肺容積增大，而實(shí)質(zhì)組織密度和信號強(qiáng)度降低；在呼氣相時(shí)則會相反。肺組織MR信號強(qiáng)度也受到心臟循環(huán)的調(diào)制。在收縮期中獲得的圖像的信號強(qiáng)度較低于在舒張期，因?yàn)槭湛s期時(shí)血流速度快，會導(dǎo)致MR信號[17]相移。呼吸的2個(gè)生理流程分別對應(yīng)不同的頻率，并與光譜時(shí)間分辨率相關(guān)。通過相關(guān)的數(shù)據(jù)可生成通氣加權(quán)（VW）和灌注加權(quán)（QW）圖像[16]。 FD磁共振通過使用1個(gè)平衡的穩(wěn)態(tài)獲得無差脈沖序列，從而獲得來自快速自由呼吸的圖像，會有非剛性圖像配準(zhǔn)用于呼吸補(bǔ)償。MR譜只受到單獨(dú)的呼吸和心臟信號的調(diào)制。因此,1次檢測就可以得到VW和QW的數(shù)據(jù)。

這項(xiàng)研究的目的是，在動物實(shí)驗(yàn)中，比較基于FD磁共振成像技術(shù)的平掃肺VW和QW圖像與作為臨床參考標(biāo)準(zhǔn)的SPECT /CT圖像。

知識進(jìn)展

肺通氣加權(quán)和灌注加權(quán)傅里葉分解磁共振結(jié)果與SPECT/CT的結(jié)果顯示一致。

斯皮爾曼秩相關(guān)系數(shù)（ρ）在左右肺的平均信號強(qiáng)度之間，分別為0.79的VW和0.86的QW圖像；肺通氣和肺灌注圖像的值分別為0.91和0.92。

所有動物的萬有引力分布的信號強(qiáng)度在2種設(shè)備中都得到了證實(shí)。

肺功能組織異常FD的磁共振成像結(jié)果和SPECT/CT的結(jié)果在區(qū)域上和質(zhì)量上都相似，并且同CT形態(tài)學(xué)上結(jié)果也一致。

對患者護(hù)理的影響

對大型動物的研究表明，F(xiàn)D磁共振成像是一種非侵入性的替代方法，可以替代SPECT/CT進(jìn)行局部肺的通氣和灌注，而不使用造影劑和電離輻射。

FD磁共振成像設(shè)備是常規(guī)的臨床設(shè)備。

FD磁共振成像檢測除一般禁忌癥外，不再需要患者遵循過多的指令，如延長呼吸或反復(fù)呼吸。

1 材料和方法

1.1 動物準(zhǔn)備

這些實(shí)驗(yàn)的協(xié)議由教育部研究，得到了農(nóng)業(yè)、環(huán)境和農(nóng)業(yè)地區(qū)（德國基爾）的批準(zhǔn)，均符合動物保護(hù)法。實(shí)驗(yàn)中，共選用了7只成熟健康的母豬，平均體重43kg（最小40kg，最大45kg），用磁共振成像和SPECT / CT評估區(qū)域肺灌注和通氣。動物在整個(gè)實(shí)驗(yàn)過程中保持全身麻醉的狀態(tài)，并被固定為仰臥位然后進(jìn)行氣管插管。呼吸完全依靠呼吸機(jī)完成。所有功能成像是在持續(xù)的、可控的通氣機(jī)器中進(jìn)行并連續(xù)監(jiān)測血壓和心率。實(shí)驗(yàn)完成后，對動物處以安樂死。

1.2 成像過程

實(shí)驗(yàn)的第1部分使用1.5T全身MR成像儀（Magnetom Avanto，西門子，Healthcare Sector, Erlangen, 德國），峰值梯度幅值45mT/m，最大的回轉(zhuǎn)速率200T/m/s。把全身12通道和24通道脊柱陣列線圈組合使用。這些方式用于評估FD MR灌注和通氣系統(tǒng)的成像，對胸腔進(jìn)行一系列冠狀圖像的采集與二維橫截面時(shí)間分辨的無差脈沖穩(wěn)態(tài)序列分析。序列參數(shù)見表1。進(jìn)一步的技術(shù)詳情見鮑曼等人[16]的研究，載于附錄E1（在線）。

表1 Parameters for FD MR Imaging

這些動物在完成MR檢查之后,被移動到SPECT /CT系統(tǒng)（Symbia T，西門子，醫(yī)療保健部）。實(shí)驗(yàn)之間的時(shí)間間隔大約是1h。 SPECT / CT系統(tǒng)包括雙頭可變角攝像頭和1個(gè)單螺旋CT。成像中將99mTc標(biāo)記的碳（Technegas，Cyclomedica，薩爾茨吉特，德國）在蒸發(fā)器中進(jìn)行蒸發(fā)，并且讓動物在5次呼吸中吸入蒸發(fā)器中的標(biāo)記物，這些均在通氣機(jī)器中完成。大部分應(yīng)用是活性物質(zhì)，未吸附在肺中的锝的生物半衰期為135h[17-18]，完全能滿足灌注顯像的要求。然而，只有一小部分（約12.5%）的活新物質(zhì)被肺部吸附。因此，活動靜脈注射示蹤比1~2個(gè)程度的吸入示蹤[19]更容易實(shí)現(xiàn)。灌注顯像中，靜脈注射99mTc來標(biāo)記大顆粒聚合人血白蛋白。在自由呼吸時(shí)的衰減校正中，用CT進(jìn)行檢查。在吸氣和呼氣時(shí)須屏氣來進(jìn)行形態(tài)學(xué)的CT成像。SPECT / CT成像參數(shù)見表2。

表2 Parameters for SPECT/CT Acquisitions

1.3 圖像的后期處理與FD

所有時(shí)間分辨的平衡穩(wěn)定狀態(tài)自由差數(shù)據(jù)均通過非剛性圖像配準(zhǔn)算法來對呼吸運(yùn)動進(jìn)行校正，該算法采用獨(dú)立的軟件 （fMRLung3.0西門子公司研究，普林斯頓，新澤西州）。呼吸和心臟周期造成傅立葉分析的MR數(shù)據(jù)信號變化。功率譜的生成包含了（圖1）每個(gè)切片的各個(gè)像素。用呼吸峰下的區(qū)域和心臟頻率計(jì)算生成VW和QW圖像。光譜圖像分析使用軟件MATLAB2009b；數(shù)學(xué)工程，內(nèi)蒂克，馬薩諸塞。圖像后期處理可在附錄E1（在線）中找到，進(jìn)一步的細(xì)節(jié)，見鮑曼等人的發(fā)現(xiàn)[16]。

圖1 機(jī)械通氣后的5號豬的肺部組織的信號強(qiáng)度變化幅度的短時(shí)傅里葉變換曲線。

由通氣和血流引起的最大信號強(qiáng)度變化分別在大約0.35Hz和1.3Hz處。

1.4 圖像評估

MR圖像，以及SPECT/ CT的數(shù)據(jù)，被導(dǎo)入到內(nèi)部開發(fā)多模軟件（DIROlab； 德國海德堡癌癥研究中心，F(xiàn)rauenhofer MEVIS，不來梅，德國）中。通過互信息的算法完成剛性圖像配準(zhǔn)形態(tài)的MR和CT數(shù)據(jù)處理。生成的配準(zhǔn)矩陣應(yīng)用于轉(zhuǎn)化VW和量子阱MR圖像并對應(yīng)相應(yīng)的通氣和灌注SPECT的數(shù)據(jù)。數(shù)據(jù)可視化分析和明顯的形態(tài)學(xué)病理檢查由2位達(dá)成共識的醫(yī)師（J.D.and C.H.）共同完成。在多層面的格式化模式中分析VW和QW FD圖像，通氣和灌注SPECT圖像分析采用的是 （MATLAB2009b；MathWorks公司）MATlAB軟件。對每一個(gè)采集的FD磁共振成像和所有相應(yīng)的格式得到的SPECT數(shù)據(jù)進(jìn)行分析，劃定右肺和左肺，包括整個(gè)肺組織，不包括血管的ROI區(qū)域。 配準(zhǔn)數(shù)據(jù)考慮了ROE區(qū)域并保證了評估的一致性。

2 統(tǒng)計(jì)學(xué)分析

通過雙肺ROI區(qū)域的分析可以對個(gè)體平均信號強(qiáng)度及其變化進(jìn)行評估。通過檢查每個(gè)ROI區(qū)域的同質(zhì)信號分布可以估計(jì)變化系數(shù)。用檢測Wilcoxon符號秩的方法來比較左右肺的變化系數(shù)分布。被計(jì)算出的左右肺之間信號強(qiáng)度的斯皮爾曼等級相關(guān)系數(shù)用于所有動物每個(gè)部位的MR圖像和SPECT圖像，為了證明FD MR圖像的個(gè)體魯棒性。

對于視覺插值，采用梯度前后插值法。合并后的最外層部分像素分布的評估采用Mann-Whitney U檢驗(yàn)法。統(tǒng)計(jì)測試執(zhí)行顯著水平= 0.05。P值<0.05指有統(tǒng)計(jì)學(xué)意義。所有統(tǒng)計(jì)分析均使用 OriginPro8；OriginLab軟件，北安普頓，馬薩諸塞。

3 結(jié)果

VW和QW的FDMR成像和通氣及灌注SPECT / CT圖像在所有的動物中有相似的結(jié)果。因?yàn)樵?種動物中，均由多層面重格式化讀取的數(shù)據(jù)來證明區(qū)分通氣和灌注病變。從相對應(yīng)的肺換氣不足和灌注不足的區(qū)域讀取針對所有模型的橫向和背向數(shù)據(jù)，以確保配準(zhǔn)的數(shù)據(jù)一致性。

如在圖2和圖3（圖所示E1的（在線））所示，在健康的肺組織中，肺血管顯示QW圖像的信號強(qiáng)度高，而在VW圖像上卻不可見。在所有SPECT圖像中，肺氣管中無核素的沉積。

MR和SPECT/ CT圖像顯示平均的CV在0.1345~0.1844和0.1151~0.1917。雙肺部所有的FD MR影像和SPECT/ CT圖像的CV估計(jì)，見表3。CV分布在雙肺部在MR成像時(shí)用Wilcoxon符號秩檢驗(yàn)，顯示差異有統(tǒng)計(jì)學(xué)意義，見表4。在左右肺平均信號強(qiáng)度間的斯皮爾曼秩相關(guān)系數(shù)ρ，VW和QW圖像分別為0.79和 0.86。在SPECT / CT圖像中，通氣和灌注圖像的ρ分別為0.91和0.92。用2種設(shè)備獲得的通氣和灌注圖像中，沿前后方向測得的重力對信號強(qiáng)度的影像情況如圖4所示。

表3 Statistical ROI Analysis for FD MR and SPECT Images

表4 CV Intramodal of Left-to-Right Lung Comparison

所有動物最前切片和最后切片的通氣和灌注的VW和QW圖像中，像素分布顯示相當(dāng)?shù)汀＃≒<0.001，采用Mann Whitney曼-魏特萊U檢驗(yàn)）。

圖2顯示了6號豬的VW和QW的FDMR成像和通氣和灌注SPECT/ CT檢查獲得的數(shù)據(jù)，支氣管阻塞的中下葉最可能造成突然的空氣滯留，右肺的CT圖像和VW圖像顯示信號強(qiáng)度相對于周圍比較健康的肺組織的變化有所降低。然而，在相應(yīng)的QW圖像中也顯示了信號強(qiáng)度區(qū)域的血液的流動減少。這些現(xiàn)象在SPECT的采集的數(shù)據(jù)中也可觀測到。

圖2 6號豬的VW和QW的FDMR成像和通氣及灌注SPECT/CT檢查數(shù)據(jù)

從3號豬獲取的圖片顯示出肺不張，見圖3。病變處位于右肺上葉的下段。從SPECT圖上可以得知，受影響的區(qū)域經(jīng)過通氣定影和灌注定影。從FD MR圖像上可以得到相同的功能信息。在VW圖上可以看到肺部密度的變化，同時(shí)，QW圖上顯示出受損肺部區(qū)域沒有血液的流動。

圖3 3號豬的VW和QW的FDMR成像和通氣及灌注SPECT/ CT檢查數(shù)據(jù)

圖4給出了雙肺上5個(gè)相同區(qū)域的FD MR和SPECT/CT圖的信號強(qiáng)度的平均標(biāo)準(zhǔn)偏差。圖中也給出了從病理學(xué)上畸形區(qū)域（3號和6號豬）探測得到的FD MR和SPECT/CT圖中ROI的信號強(qiáng)度。如圖4所示，這些ROI的信號強(qiáng)度比正常組織要低，并在單一標(biāo)準(zhǔn)偏差的范圍之外。

圖4 雙肺上5個(gè)相同區(qū)域的FD MR和SPECT/CT圖的信號強(qiáng)度的平均標(biāo)準(zhǔn)偏差

4 討論

在核醫(yī)學(xué)中，使用SPECT和SPECT/CT進(jìn)行通氣和灌注成像是一項(xiàng)已被認(rèn)可的技術(shù)，通過這項(xiàng)技術(shù)可以發(fā)現(xiàn)許多種肺部疾病所造成的生理改變。在這項(xiàng)研究中，通過與臨床已建立的參考成像手段進(jìn)行對比，我們對一種非增強(qiáng)的肺功能MR成像方法進(jìn)行了評估。在大量動物身上我們做了通氣和灌注的對比。

通過對雙肺的測量，統(tǒng)計(jì)學(xué)分析表明，平均信號強(qiáng)度（VW和QW）與某些量（通氣與灌注）之間存在很好的相關(guān)性。2種模式下的個(gè)別單側(cè)CV表現(xiàn)出相同的均一性；而威爾科克森標(biāo)志排列測試表明雙肺部的QW成像存在很大的差異。從早先的研究中可以得知，心旁區(qū)域的搏動偽影是造成這種差異的可能原因。

這項(xiàng)工作的實(shí)驗(yàn)環(huán)境中，測量肺部嚴(yán)重肺功能改變的FD MR成像的臨床意義能得到證實(shí)。由于測量時(shí)間的限制，F(xiàn)D MR成像被限制為5s，前后方向上為6cm。通氣和灌注中，被廣為接受的重力依賴的差異能被復(fù)制再現(xiàn)并通過FD MR和SPECT成像證明符合統(tǒng)計(jì)學(xué)規(guī)律。對FD MR成像而言，尚未發(fā)現(xiàn)參考信號，目前還不能計(jì)算對比VW和QW的比率。盡管如此，碰巧發(fā)現(xiàn)的2例通氣與灌注缺陷部位的信號強(qiáng)度明顯降低，其值已在正常組織的信號強(qiáng)度標(biāo)準(zhǔn)偏差之外。該結(jié)果也通過SPECT得到了類似的證實(shí)。29%的肺部患病率似乎很高，也有一些報(bào)告給出更高的發(fā)病率和患病率。此外，通常會發(fā)展為急性或慢性感染。通過與鎮(zhèn)靜、侵入式插管法與仰臥位的常規(guī)操作對比，或許可以解釋我們的發(fā)現(xiàn)。

該項(xiàng)研究選用了在良好環(huán)境下生長的豬。2個(gè)偶然的病理性發(fā)現(xiàn)表明通氣-灌注FD MR成像與標(biāo)準(zhǔn)SPECT/CT成像之間具有一致性。但是，仍需要對大量的患者進(jìn)行調(diào)查，來確認(rèn)這種研究方法的魯棒性和重復(fù)性。另外，這些方面，如觀測者的變化也需要考慮。

在這項(xiàng)動物研究中，所有檢查都是在人為控制的通氣環(huán)境下進(jìn)行的，該技術(shù)不需要穩(wěn)定受控的呼吸過程。數(shù)據(jù)采集過程中，呼吸和心率的大幅改變會影響到信噪比，但檢查仍可以進(jìn)行。基于MR的通氣和灌注測量不產(chǎn)生電離輻射。實(shí)驗(yàn)中動物處于大劑量的核輻射中，因而需要另外的輻射保護(hù)。在接下來的肺部通氣和灌注實(shí)驗(yàn)中，基于FD的MR成像無此限制。另外，在后面的實(shí)驗(yàn)中，消除造影劑可以減小其他藥理學(xué)上的副作用或相互影響。這表明，F(xiàn)D MR成像在妊娠婦女、兒童及腎衰竭患者的檢測中具有很大的吸引力。至于采集時(shí)間，VW和QW FD MR成像可以縮短成像時(shí)間，從而增加患者的舒適感。不過，這會減少測量的固定時(shí)間。當(dāng)前的技術(shù)已經(jīng)可以在1min之內(nèi)獲取1個(gè)截面厚度為12mm的二維時(shí)間分辨平面圖。研究中，檢測5個(gè)截面用時(shí)5min，同時(shí)獲得了通氣量和灌注量。檢測人肺的測量時(shí)間可能會增加到大約15min,正向覆蓋范圍為18cm。在本文的實(shí)驗(yàn)中，單獨(dú)獲取通氣和灌注的SPECT數(shù)據(jù)持續(xù)了85min，對SPECT來說是可以縮短成像時(shí)間的，但相應(yīng)的計(jì)數(shù)率也會相應(yīng)的下降并會限制空間分辨率。

大血管的直接可視FD灌注MR成像可能是另一個(gè)值得研究的方面。這可應(yīng)用于探測急性肺栓塞患者的大血管內(nèi)的血栓。在灌注SPECT研究中，核素的分布取決于毛細(xì)血管堵塞的比例。核素不會存在于大血管中，因此，這些區(qū)域中存在無效區(qū)域。MR成像的優(yōu)勢在于在檢測中能進(jìn)一步的結(jié)合形態(tài)學(xué)成像，這將使得單一的臨床MR成像機(jī)和一個(gè)擁有高質(zhì)量CT完整的SPECT研究一樣全面。

更多的優(yōu)勢體現(xiàn)在邏輯方面。閃爍掃描術(shù)和SPECT需運(yùn)用99mTc,它來自于鉬99。2008年，所有的核素資源由于各種各樣的原因不可用，導(dǎo)致了世界范圍內(nèi)的同位素短缺。肺部功能評估的前提條件通常是侵入式的，如肺部手術(shù)的先決條件。因此，很多過程都不得不推遲。另一方面，F(xiàn)D MR成像能被廣泛的應(yīng)用而不會遇到任何供應(yīng)鏈的瓶頸。

通氣-灌注FD MR成像的主要限制在于通氣和灌注可視化的間接方法。通氣成像主要是通過持續(xù)測量呼吸引起的胸壁和肺質(zhì)子密度的改變。肺泡和支氣管空間并不直接成像。這項(xiàng)研究中引用的臨床參照具有類似的限制：超細(xì)分布核素標(biāo)記的碳原子被肺泡和支氣管表面吸收。因此，SPECT只能提供呼吸系統(tǒng)的靜態(tài)成像。FD MR成像中，肺灌注成像是通過依賴流量的信號移相。與中心大血管對比或合并，增加肺外部的磁導(dǎo)率是否會影響測量方法的敏感度仍需要進(jìn)一步的評估。在大顆粒白蛋白的穩(wěn)態(tài)分配中，本文又一次運(yùn)用臨床參照模式進(jìn)行了1次靜態(tài)的測量。穩(wěn)態(tài)SPECT成像中核素分布不受胸部或心律的影響，QW成像的信噪比會受心動周期中心律不齊的影響，如果胸部受到機(jī)械的不穩(wěn)定的傷害，VW成像將是不可靠的，見圖5。

圖5 豬雙肺ROI測量中信號強(qiáng)度的平均差和標(biāo)準(zhǔn)差（誤差線）

總之，我們給出了一個(gè)不受電離輻射和造影劑影響的方法，它能用于肺通氣和灌注的評估，即使是在呼吸困難不順應(yīng)的患者身上也同樣適用。基于這些最初的結(jié)果，F(xiàn)D MR成像需要更進(jìn)一步的實(shí)驗(yàn)和臨床評估來證明該方法的魯棒性，診斷的精確度以及觀測者間的重復(fù)性。

[1]Tetzlaff R,Schwarz T,Kauczor HU,et al.Lung function measurement of single lungs by lung area segmentation on 2D dynamic MRI[J].Acad Radiol,2010,17(4):496-503.

[2]Patz EF,Coleman RE.Nuclear medicine techniques.In:Mason RJ,Murray J,Broaddus VC,Nadel J,eds. Textbook of respiratory medicine[M].3rd ed.Philadelphia,Pa:Saunders Elsevier,2005:594-599.

[3]Roach PJ,Bailey DL,Harris BE.Enhancing lung scintigraphy with single-photon emission computed tomography[J].Semin Nucl Med,2008,38(6):441-449.

[4]Gur D,Drayer BP,Borovetz HS,et al.Dynamic computed tomography of the lung: regional ventilation measurements[J].J Comput Assist Tomogr,1979,3(6):749-753.

[5]Chae EJ,Seo JB,Goo HW,et al.Xenon ventilation CT with a dual-energy technique of dual-source CT: initial experience[J].Radiology,2008,248(2):615-624.

[6]Thieme SF,Hoegl S,Nikolaou K,et al.Pulmonary ventilation and perfusion imaging with dual-energy CT[J].Eur Radiol,2010,20(12):2882 -2889.

[7]Cleveland ZI,Cofer GP,Metz G,et al.Hyperpolarized Xe MR imaging of alveolar gas uptake in humans[J].PLoS ONE,2010:5(8):e12192 .

[8]Fain S,Schiebler ML,McCormack DG,et al.Imaging of lung function using hyperpolarized helium-3 magnetic resonance imaging: review of current and emerging translational methods and applications[J].J Magn Reson Imaging,2010,32(6):1398-1408.

[9]Edelman RR,Hatabu H,Tadamura E,et al.Noninvasive assessment of regional ventilation in the human lung using oxygen-enhanced magnetic resonance imaging[J].Nat Med,1996,2(11):1236-1239.

[10]Berthezène Y,Vexler V,Clément O,et al. Contrastenhanced MR imaging of the lung: assessments of ventilation and perfusion[J].Radiology,1992,183(3):667-672.

[11]Hatabu H,Tadamura E,Levin DL,et al.Quantitative assessment of pulmonary perfusion with dynamic contrast-enhanced MRI[J].Magn Reson Med,1999,42(6):1033-1038.

[12]Grobner T.Gadolinium: a specifi c trigger for the development of nephrogenic fi brosing dermopathy and nephrogenic systemic fi -brosis? [J].Nephrol Dial Transplant,2006,21(4):1104-1108.

[13]Henderson AC,Prisk GK,Levin DL,et al. Characterizing pulmonary blood fl ow distribution measured using arterial spin labeling[J].NMR Biomed,2009,22(10):1025-1035 .

[14]Deimling M,Jellus V,Geiger B,et al .Time resolved lung ventilation imaging by Fourier decomposition[abstr].In:Proceedings of the Sixteenth Meeting of the International Society for Magnetic Resonance in Medicine[M].Berkeley,Calif: International Societyfor Magnetic Resonance in Medicine,2008:2639.

[15]Bauman G,Puderbach M,Deimling M,et al .Non-contrast enhanced perfusion and ventilation MRI of the lung by means of Fourier decomposition[abstr].In: Radiological Society of North America scientifi c assembly and annualmeeting program[M].Oak Brook,Ill: RadiologicalSociety of North America,2008:788.

[16]Bauman G,Puderbach M,Deimling M,et al.Non-contrastenhanced perfusion and ventilation assessment of the human lung by means of fourier decomposition in proton MRI[J].Magn Reson Med,2009,62(3):656-664.

[17]Suga K,Ogasawara N,Okada M,et al.Lung perfusion impairments in pulmonary embolic and airway obstruction with noncontrast MR imaging[J]. J Appl Physiol,2002,92(6):2439-2451.

[18]Kawakami K,Iwamura A,Goto E,et al.Kinetics and clinical application of 99mTctechnegas[J].Kaku Igaku,1990,27(7):725-733.

[19]O'Brodovich HM,Coates G.Quantitative ventilationperfusion lung scans in infants andchildren: utility of a submicronic radiolabeledaerosol to assess ventilation[J].J Pediatr,1984,105(3):377-383.

[20]Kellman P,Chefd'hotel C,Lorenz CH,et al. Fully automatic,retrospective enhancement of real-time acquired cardiac cine MR images using imagebased navigators and respiratory motioncorrected averaging[J].Magn Reson Med,2008,59(4):771-778 .

[21]Chefd' hotel C,Hermosillo G,Faugeras O.A variational approach to multimodal image matching. Proceedings of the IEEE Workshop on Variational and Level Set Methods in Computer Vision (VLSM'2001),ICCV Workshop,Piscataway[M].NJ:Institute of Electricaland Electronics Engineers,2001:21-28.

[22]Chefd'hotel C,Hermosillo G,Faugeras O.Flows of diffeomorphisms for multimodal image registration. Proceedings of the IEEE International Symposium on Biomedical Imagi ng(ISBI'2002),Piscataway[M].NJ: Institute of Electrical and Electronics Engineers,2002:753-756.

[23]Bauman G,Eichinger M,Deimling M,et al.Non-contrastenhanced lung perfusion MRI in comparison to contrastenhanced MRI perfusion in young cystic fi brosis patients[D].Presented at the annual meeting of the European Congress of Radiology,Vienna,Austria,2010.

[24]Bates DV.Measurement of regional ventilation and blood fl ow distribution.In:Fenn WO,Rahn H,eds.Handbook of physiology.Section 3.Respiration [M].Vol 2:Washington,DC:American Physiological Society,1973:1425-1436.

[25]Jirawattanapong P,Stockhofe-Zurwieden N,van Leengoed L,et al. Pleuritis in slaughterslaughter pigs: relations between lung lesions and bacteriology in 10 herds with high pleuritis[J].Res Vet Sci,2010,88(1):11-15 .

[26]Kirschvink N,Reinhold P.Use of alternative animals as asthma models[J].Curr Drug Targets,2008,9(6):470-484.

[27]Bauman G,Puderbach M,Deimling M,et al.Non-contrastenhanced MRI of lung perfusion and ventilation by Fourier decomposition:volunteer study [abstr]. In: Proceed ings of the Seventeenth Meeting of the International Society for Magnetic Resonance in Medicine[J].Berkeley,Calif: International Society for Magnetic Resonance in Medicine,10.

[28]Puderbach M,Hintze C,Ley S,et al. MR imaging of the chest: a practical approach at 1.5T[J].Eur J Radiol,2007,64(3):345-355.

[29]Gheeraert M.Press release 25 August. Association of Imaging Producers & Equipment Suppliers(AIPES).http://www.aipeseeig.org/.Published August25,2008.Accessed October20,2010 .

[30]Gheeraert M.Press information.Associationof Imaging Producers& Equipment Suppliers(AIPES).http://www.aipes-eeig.org/.Published August 28,2008.Accessed October20,2010.