MRI数据预处理流程资料讲解.docx
- 文档编号:23415462
- 上传时间:2023-05-16
- 格式:DOCX
- 页数:15
- 大小:532.69KB
MRI数据预处理流程资料讲解.docx
《MRI数据预处理流程资料讲解.docx》由会员分享,可在线阅读,更多相关《MRI数据预处理流程资料讲解.docx(15页珍藏版)》请在冰豆网上搜索。
MRI数据预处理流程资料讲解
数据处理基本流程
由于MRI是断层扫描,耗费时间较长,患者在进行MRI扫描的时候不可避免的会头部挪动,导致照射出来的图像不能一一映射;不同人的头颅,脑部大小,形状都会有所差异,获得的MRI图像也千差万别,无法对其进行对比。
所以我们就必须用一种算法将所有的MRI图像进行空间转换到一个比较标准的空间
(目前使用较多的是被神经学家广泛认可的Talairach坐标系)将各个解剖结
构一一对应后,再与标准化图谱或者不同个体之间相互比较(目前使用的是
Talairach-Tournoux图谱)
本文使用的是SPM软件和MRIcro软件处理图像数据,将MRI图像进行数据分析。
数据分析的基本流程:
(1)数据预处理:
。
图像格式转换C2slicetiming获取时间校正◎realign头动校正CCoregister不同成像方法间的图像融合C5nomalize不同被试之间的图像标准化(归一化)Csmooth空间平滑《234统称图像的空间变换》
(2)模型构建与参数估计:
C建立统计模型C2将数据应用于统计模型C3进行参数统计得到单个被试的结果,多个被试的组分析
数据预处理
SPM是一款以MATLAB为平台的软件,所以使用SPM前一定要安装MATLAB。
1•图像格式转换
在进行数据预处理第一步要先将图像格式转换成SPM可以识别的
ANALYZE格式。
转换之前先将原始数据放在MATLAB下面的mriimage文件夹
下,将路径设置成D:
\MATLAB\work\mriimage'设置过程如下:
FkEdrtDi?
bugO«»lil!
apWndawHdp
DQ■丄科吁*CU*WI:
feTKfcfVE^FIr■仲
FC449PI・9E#9l確dlTflrt&ITFUKRIM>i:
r
Cm2:
JD)
山妙北-
pr/OV20l.2
smz:
(me.art
33E3HIK-
Dr/OJ/201.2
ST12:
dpBL_nAEAAlhd^_ (vE.]0b ntMitEd- Dr/OJ/201.2 CaHim*nd ]加昭i阴-DF/oifaaia TodltooEdh£4Ci E>i\Fp»J\twlbOK\Pl^ rlZm亍hgl; 求IE: dpk._lfeM^i? 2-2Q) ]? t»lLD-nr/W2DL2 ^-d: ra-ih ««[§3I&3I§3I&3IS3口: d»l«i|th4XIi =i=DaLd-I); k>; B46—IS->7±£FI2: G2T — 打i>: Vv^rk^itiL口毎命wiTfcx—D: VndIR[! ■! [■»&tk\I—prd «--i-TliJS1—' S-S—S2^>? 上牛乐: M-TC #■8IW+Burr' 点击红色方块所指的按钮,在弹出的窗口中选择工作路径,按确定按钮即可。 设置完工作路径后,利用如下方法,将SPM2及其所有子文件夹添加到 MATLAB的搜索途径中(1.点击file按钮,在下拉菜单选择setpath2在弹出的路 径设置窗口点击"AddFolder"浏览并选择目标文件夹,eg: D: \spm2\3点击save按 钮4.点击close按钮,完成添加) 在打开SPM之前,应先确定默认变量的设置是否准确,具体做法如下: 1. 在matlab命令窗口输入“editspm_defaults‘打开spm_defaults.m文件2.查看 defaults.analyze.flip条目,确认defaults.analyze.fip值是否为1,若不是,改成1 %Fileformatspecific defaultanalyze,multivol=0; defaults,analyze.£lip=1;*、<<=Veryimportant.RelatestoL/R 打开SPM: 在matlab命令窗口输入“spm"回车后出现下面窗口,按黄色长方形 覆盖的按钮,方可打开SPM软件(或者直接输入spmfmri即可打开) 左上角的窗口代表按钮窗口, $叶1: fldminiirir^grj FitEdEVitwT4dkD^tikUipWndiSwH白0匚olgirFCl^ar1Rtfiub-Ff^ 5PM2 Wfl-komhgoSPM2 A丄也西比iiSHEuill•Lteitilu^t£ilt$tx^tapcevidufvtESiotuat Thtee-SLiitEeDceflThtMgviElrt Th卍匕曾(Cite卍jgLfr^eu±fid典LiUMA.titbit1! 0且tOE pUH]flCX- =丄vmwchI-iuaruEWL«irc3>ii: ・hbcLtirpdkud;iEmLpj3B pabtLeulHE.Uu-lbau.t5MT"tcpLeun>fwhlcSillatvnawEn'tUEHB ? PK2. ^urdMixl»C#iria? a4nisybeCc^wda時die5PfW*bvim; bri^pi//mnraCU■ian>ucli屈船 OtMilPQt口r3THCMtl曲酣《n11DC;can2xfijund; blips//miuCU.ictn^ucl.&i? >■ifey^n.'lulp L3dtVC-LOptdOut白出丄丄亡屮«XFUXU•血0白】丄■巾口】上窑&血EPXQ匕1丫 (n.L|iThnVflLlEDME'vphEtaantvtX»/Lnq11・3口aca-incajan^JihtitIhuiu丄住力at^eu1m±Lv! iL^jC^U*gtLduLuii>|UC&]UK. 9HISwasj>eime13tdwtiapedtor>牝圧史工归总工日80*T4>cH>oL&I.. 』口仙JLaklniEnjHX|FI11,XatIrrajxirn|rtl]4Indxnl£aLbve(FIL)a PMliuICTfDFJI■■IIT-.Fionrtl? Hintc«nr: CLhur;l^« CroamascSw■-ns■亍辑CRFE-CIJUj日^群孙,Tn#L凹曲"巧piig”R的萼14和iCKk.'DFfl/fflFJ.Dem.^.FciiweJjT«tBiclwls(Depc»Elio3C4T10Clcs-r01JIl-fbl^EUhHI.UML|ybatxenD^U丄・g(NtiLthMdlltEb Uhiwiflici,1iIlilfi^lIt£ir^-DhfioceTiealinpucMe陰ixhHozsltyi.L'E^'iiE-LaiaiLu£5! _«L13t3-=\3,FIl-^lLiUuLv^s.33Lt.nmLLca丄‘UanadH|. SE*H(bauigUuBciLlflEka-sadEElIbsatwnactha>uiJ.£BStinUm 用以对数据进行处理分析; 左下方的窗口为输入窗口, 右边的 Jhrirr.riwrj|r. 窗口为图形窗口,用以显示图像结果 图像格式转换: 1•点选按钮窗口中Toolboxes菜单中的Dicom选项2.在弹出对 话框中选择所有*IMA文件,点done3转换完毕,删除原文件夹下的*IMA数据4. 结果将在该文件夹下生成三个*img/hdr/mat文件5.结束 图像格式转换方法二: 用MRIcro软件进行格式转换 (1)单击convertforeignto analyze (2)出现第二个窗口,将左边的数据填完后,单击箭头所指“选择” Import HeaderViewROIOverlayEtc Ctd+F Qrl+D Op&nforeign... Display/lhideforeignheader ConvertforeigntoAnalyze...[DICOM/Elsrint/GE/Interfile/PickBr/SieirTiensiA/FF] CcnveirtECATtoAnalyze.^. CorweirtSPMivinVHDtoAralyze..u ConvertTIRTLeiratoAnalyze.i- CorweirtFreesurfertoAnalyze ConvertVista/LipsiatnAnalyze..・ ConvertAnalyzetoVrsta/Upsia.u. ConvertPhilipsPARtoAnalyze... CorwertAFNIHEAD/BR1KtoAnalyze... CcnveirtNlfTI-niito.hdr/Jmg Ccnversicrwizard: ccnvertallAle-sinafplder® Selectrootfoldetthatcontuinwall】.・a輕2[andnoothqyfil^s]' ^udiLion3.□ CTi・《L嵐0■暑DSF冥噓mptlsbmri ■Elimt^e MKIim«.geO听力科硏训练 个文件"*i.img/hdr"。 (如被试者有n次试验,则需要建立n个文件夹,分别转换图像格式) 头动校正: 意义: realignmentoffunctionaltime-series. 1点击按钮窗口中的Realign下拉菜单中的Realign按钮(将同一被试者不同采 Choc*'Reg 」 |Render… Jj |FMF? I PPi= 1 ImC-rsfc: 」 Biascor |Ulils... Defautts 」 Quit1 DIspW 样时间点上的3D脑图像对齐 ◎numbersubjects要处理的被试个数eg: 1]0numbersessions,subj1第一个被 试者的试验次数eg: 1](3images,subj1,sess1选择文件夹中所有i.img文件],点donet^whichoption? [coregister&reslice]⑥createwhat? [*AIIimages+mean image] 结果SPM将更新i.mat文件,并文件夹中生成一个头动参数文件(rp_i.txt),还在文件夹 中生成meani.img/hdr/mat文件跟ri.img/hdr/mat文件。 (如果第一个被试者有n次实验,则头 动校正结果为: 在每个文件夹中SPM均更新i.mat文件,并分别生成一个头动参数文件 (rp_i.txt),还在文件夹中生成meani.img/hdr/mat文件,并在图像窗口中显示n个试验的 的头动曲线图,如下) rotation I艸I IUn soiwi5c image 注: 如试验次数只有一次,就没有上述头动曲线图 Coregiste《图像融合(配准)》[头动校正仅对同一被试的同一种成像方法 (或成像模态)有效,对于同一被试的不同成像方法所的图像,由于它们之间没 有足够的可比性,就需要用图像融合的方法来做空间校正】 基于MRI,健康人的脑组织通常可分为三中组织类型: 灰质,白质,脑脊液。 SPM软 件采用聚类算法中的HARTIGAN混合模型的最大似然方法来分隔这三部分。 (这一方法开 始只用在T加权MRI上,后来也用到PET等其他脑功能成像方法中的脑组织分割)分割完 后,同一种组织不同成像方法所得图像之间就有了足够的可比性。 关键的步骤: 点击按钮窗口中的CoregisterCD.numberofsubjects/session[1] 2whichOption[coregisteronly] SPM2sAdmen冷tytor)Coregister 2Targetimage: 选target文件夹下的T1.mnc或者别的模板] 3Sourceimage: 选功能像目录下的meani.img] ②Otherimage: 选功能像目录下的ri..img] 结果: 产生ri.mat和一个meani.mat文件 NormalisedMutualInformationCoregistration XI=-Q.405*X*0.251*Y-0.192*Z*60.826 Y1=D.291'X*D.lzgnr-D.l7StZ-B.558 zi=-D>i0jM.i01B1FinalJointHistogram ..3pm2yEmpl3te3M1.mnc..spm2M: emplates\T1.mnc NarmallsedMutualInformationCoregistration .Atenipl3te$\Trgmsm.mnc a Y1h0.29DwX*0.37fl,,y-C.13rZ-B.3H bE二ueajEEQAo .vtempiatesMr^n-smmnc FinalJointHistogram bE-EralllUA'c号u =-0405*XW249*Y-0.164*ZWO597 Normalize【spatialnormalizsationofimagevolumestoatemplate空间标准化到一个模板上)】个体大脑在形状,大小等方面存在明显差异,所以我们将不同被试者脑图像标准化到同一个标准脑空间(通常是Talairach标准脑空间),使不同被试者脑图像同一像素代表相同的解剖位置) 为了确保最优的标准化,image和template(模板)应该被放在相似的开始位置,这可以通过选择CheckReg检查。 Normalize前需进行的操作⑴在SPM中选择default,然后在对话框里选择spatial normalize.Odefaultforwritingnormalize.©Preservedconcentration©Boundbox中选: [-90: 90-126: 90-72: 108(Template)]©Voxelsize中选: 3*3*3©Interpolationmethod中选: Trilinearinterpolation⑦Nowarp女口下: Normalize具体步骤: 点击按钮窗口中的Normalize按钮 然后①Whichoption? [*DetermineParameters&WriteNormalised(目的是把image标准化到模板上,然后把normalizeparameter标准化参数》保存在mat文件里)②Templateimage[spm2\template目录下的T1.mnc]③Sourceimage[eg功能像目录下的meani.img文件]④Imagestowrite[功能像目录下的ri.img文件]©当 问第二个被试的Sourceimage如果没有,直接点击Done就可以。 结果产生wri开头的hdr和img文件以及spm2.ps,"meani_sn.maf'文件,并在图形窗口中呈现数据的配准结果。 (如图所示) SpatialNormalisation Image: G: ^nrhjneainii.img Linear(affine)component XI=0・0.15-0J1rz-5386 Y1■0J«31*^^0.99a*Y-0.260*2-2.097 21=0.2BD*X*D.085*¥+1」33*Z*0.001 1Enonlinear 7x'Ox7basteAjnctlun^ Smooth空间平滑就是将数据在空间上用一个光滑的函数<通常为高斯函数>去卷积,这个光滑的函数被称为卷积核函数】 Spatialsmoothing的作用: 1,提高信噪比,脑功能成像探测到的神经生理信号来源于数毫米尺度内脑血流等的改变,图像重建中对应着低空间频率部分,噪声对应着高频部分,空间平滑后,噪声得到很大抑制2•使残差项更符合高斯分布假设, 这对用高斯随机场理论来做统计推断很重要3•消除不同被试脑结构间的细微差 别。 在图像归一化后,不同被试之间的脑结构仍有细微差异,空间平滑后这种差 异将被模糊化 具体步骤: 1,在按钮窗口中点smooth按钮2,smooth{FWHMinmm}[8] 《FWHM全称为FullWaveatHalfMaximum,指高斯滤波的波宽,以毫米为单位。 建议取voxel大小的2~3倍,一般设为8》3.selectscans选择文件夹中所 有"wri.img文件],done结果: 产生swri.img/hdr.mat文件,下图为平滑后的图像: MlAdnniniKralor'召中間血 AleEdtVkwInsswtlo^lsWindo*! I4elpCclo-rsClewSPM-Prit rw翳儈・1開 6pm-)0n«rriMliaedi- rlgMrmmlhnwdHmm|up1'fnnllprtch|rad|>rolgrwJ|rtw■: rMrtsra-14 Oriam^l«: )25 cwCobi.suooo ID1OBDDD>aa C.«Q4O.CltDI.Ht £
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MRI 数据 预处理 流程 资料 讲解