smri 预处理脚本

基于ANTs and FSL 的预处理脚本

最近在写一个toolbox,所以尝试了各种主流的预处理方法和算法,通过不断的尝试,的出来了下面的这个预处理的pipeline比较好用,这里记录一下推荐给大家,主要的代码是用python实现,但是其实还是调用的ANTs跟fsl,核心代码如下所示

import os
from multiprocessing import Pool as Pool
Standard_image_T1w = '../Template/MNI152_T1_1mm_brain.nii.gz'
Standard_t1w='../Template/MNI152_T1_1mm.nii.gz'
Standard_mask = '../Template/MNI152_T1_1mm_brain_mask.nii.gz'
#Not BIDS format
def T1w_prep(input_dir,output_dir,sub):
    #T1_directory = os.path.join(input_dir,'anat','')
    os.system('mkdir '+output_dir+'tmp/')
    os.system('mkdir '+output_dir+'tmp/T1/')
    os.system('mkdir '+output_dir+'T1')
    os.system('mkdir '+output_dir+'Seg')
    T1_files = os.listdir(input_dir)
    for i in range(len(T1_files)):
        if T1_files[i].find('T1w.nii') != -1 :
            T1_name = T1_files[i]
            T1_file = os.path.join(input_dir,T1_files[i])
    #preprocess T1 image,核心代码如下
    #Crop
        if os.path.exists(output_dir+'tmp/T1/crop_'+T1_name) == 0:
            os.system('robustfov -i '+T1_file+' -r '+output_dir+'tmp/T1/crop_'+T1_name)
        #Skullstrip
        if os.path.exists(output_dir+'tmp/T1/BEBrainExtractionBrain.nii.gz') == 0:
            if T1_name.find('.gz'):
                os.system('antsBrainExtraction.sh -d 3 -a '+output_dir+'tmp/T1/crop_'+T1_name+' -e '+Standard_t1w+' -m '+Standard_mask+' -o '+output_dir+'tmp/T1/BE')
            else:
                os.system('antsBrainExtraction.sh -d 3 -a '+output_dir+'tmp/T1/crop_'+T1_name+'.gz'+' -e '+Standard_t1w+' -m '+Standard_mask+' -o '+output_dir+'tmp/T1/BE')

        #Registration
        if os.path.exists(output_dir+'T1/T1w2std_'+T1_name) == 0:
            os.system('antsRegistrationSyN.sh -d 3 -f '+Standard_image_T1w+' -m '+output_dir+'tmp/T1/BEBrainExtractionBrain.nii.gz  -o ' \
                +output_dir+'tmp/T1/rega2t -n 20')
            os.system('antsApplyTransforms -d 3 -i '+output_dir+'tmp/T1/crop_'+T1_name+' -o '+output_dir+'T1/T1w2std_'+T1_name+' -r '+Standard_t1w+ \
                '-t '+output_dir+'tmp/T1/rega2t1Warp.nii.gz -t '+output_dir+'tmp/T1/rega2t0GenericAffine.mat')

        #Segmentation
        if os.path.exists(output_dir+'Seg/'+sub+'_CSF_mask.nii.gz') == 0 or os.path.exists(output_dir+'Seg/'+sub+'_WM_mask.nii.gz') == 0:
            os.system('fast -S 1 -t 1 -o '+output_dir+'Seg/seg -g -n 3 -b -I 10 '+output_dir+'tmp/T1/rega2tWarped.nii.gz')
            os.system('mv '+output_dir+'Seg/seg_seg_0.nii.gz '+output_dir+'Seg/'+sub+'_CSF_mask.nii.gz')
            os.system('mv '+output_dir+'Seg/seg_seg_2.nii.gz '+output_dir+'Seg/'+sub+'_WM_mask.nii.gz')
            #os.system('mv '+output_dir+'Seg/seg_seg_1.nii.gz '+output_dir+'Seg/'+sub+'_GM_mask.nii.gz') # mask 灰质的一般不需要

唯一的缺点是,在进行配准的时候是基于迭代优化求解的,所以需要花费比较多的时间,但是效果还是比直接用cat12做出来的好很多的。当然,这个用python调用也是支持并行的,并行代码可以参考下面

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