一个处理IVIM,ADC,DKI数据的SPM模块

IVIM,ADC和DKI是弥散磁共振图像的三个弥散模型,并已广泛应用于胸腹部脏器的研究中。然而,针对这三种模型的软件网上几乎找不到。为了解决这个问题,本人写了一个spm模块,专门用于处理这三种模型的数据。之所以写成spm模块,是因为spm里面提供了NIFTI数据的读写和其他操作功能,另外spm模块可以非常方便的进行批处理。

准备

软件
  1. MATLAB: 2017b以上的版本
  2. SPM: spm12版本
数据

使用mricron软件中的dcm2niigui.exe,将图像转换成FSL/SPM8(4D NIfTI nii)。得到一个4D的DWI图像,和记录b值的bval文件,bval文件后缀给成.txt。

安装和使用

DWI处理模块,可以从该地址下载。

此处假定,MATLAB安装完成,并且spm12已经正确安装。

  1. 在MATLAB的命令窗口输入以下命令,spm软件的启动窗口:
spm fmri
spm1.png
  1. 点击Quit旁边的Batch,启动一个新的窗口,在这个窗口中,点击File,点击Add Application,选择提供给大家的spm模块,里面的cfg_dwi_master.m,然后点击左下角的Done.
    spm2.png

    spm3.png

    这时对比spm2.png,可以在窗口的最上面多了一个叫做DWI Analysis的功能键。
    spm4.png
  2. 点击DWI Analysis,选择要用的模块,比如,选择IVIM analysis, 选择IVIM module,模块便被加载的窗口中,可以导入数据进行分析:
    spm5.png
  3. 右边框中,<-X指代的是这行需要输入制定的数据,具体需要导入什么数据,可以在最下面灰色框中,看到提示。
IVIM[1][2][3]
  • Image: 选择4D DWI文件
  • Mask: 选择ROI图像,指定处理的解剖范围,可用ITK-SNAP等软件勾画
  • b value: 选择记录b值得text文件
  • Fit type: 这个模块提供了两种处理方法,一种是seg,一种是bayes
  • bval threshold: 指定b值分析的阈值,这里指代简化模型拟合时需要的范围
  • outName: 给一个输出文件的名字
ADC和DKI

ADC和DKI的输入参数与IVIM相比参数会少一些。具体输入什么数据可以参考IVIM。

  1. 当所有的数据输入完成后,左上角的箭头会变成绿色,这时,点击箭头,则开始数据分析。

结果

IVIM输出的是一个4D的参数文件,第四维度分别是D, S0,f, Dstar[4]
ADC输出一个后缀为_ADC.nii的文件。
DKI每一个参数会输出一个文件,后缀分别是_D.nii_K.nii
下面介绍一下如何提取IVIM的参数值:
将生成的参数图像导入到ITK-SNAP软件中,在下图箭头所指的地方点击,选择2/4,也就是显示S0图像,这个图像通常解剖位置清晰,可以用来勾画ROI。勾画想要计算的区域。勾画完成后,选择Segmentation, 选择Volumes and Statistics,在弹出的窗口中,就会显示ROI内每个参数的数值。

ITKSNAP2.png

ITKSNAP1.png

ITKSNAP3.png

  1. IVIM模型采用的别人的代码,如果使用该模块处理数据务必引用下面的两篇文章:

  2. Gustafsson, Oscar et al. Impact of prior distributions and central tendency measures on Bayesian intravoxel incoherent motion model fitting. Magn Reson Med . 2018;79:1674–1683

  3. Jalnefjord, Oscar et al. Comparison of methods for estimation of the intravoxel incoherent motion (IVIM) diffusion coefficient (D) and perfusion fraction (f). MAGMA . 2018;31:715–723

  4. S0图像是b=0的图像

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 227,663评论 6 531
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 98,125评论 3 414
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 175,506评论 0 373
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 62,614评论 1 307
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 71,402评论 6 404
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 54,934评论 1 321
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 43,021评论 3 440
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 42,168评论 0 287
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 48,690评论 1 333
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 40,596评论 3 354
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 42,784评论 1 369
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 38,288评论 5 357
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 44,027评论 3 347
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 34,404评论 0 25
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 35,662评论 1 280
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 51,398评论 3 390
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 47,743评论 2 370

推荐阅读更多精彩内容