如何利用freesurfer求取被试的LGI指标

老师,我在wsl子系统里装了matlab和freesurfer,并且用reconall跑数据没有问题,按照求取LGI指标的步骤,在MATLAB\R2016b\toolbox\local文件中创造了startup.m文件,编辑如下:

然后在子系统中输入 recon-all -i a_CDQ_a.nii /mnt/d/fmri_data/T1_cnsp/freesurfer_scn/recon_IGL/ sub001 -localGI命令,出现以下问题,

请问如何解决哇,求取LGI指标的步骤是这样的吗?

你recon-all是怎么跑的?计算LGI就是把-all换成-localGI,比如:

recon-all -s sub001 -localGI

recon-all是这样跑的recon-all -i a_CDQ_a.nii /mnt/d/fmri_data/T1_cnsp/freesurfer_scn/recon sub001 -all -qcache,没有问题,但是按照上图中all换成-localGI就出现问题 ,不知道咋回事,添加的步骤和上图一样

没看懂你的代码,没有-s选项?

源代码这样,
image,运行也能跑,得出的结果文件这样的:




里边也都有东西,老师你的意思是如果recon-all能跑的话,只需要把-all换成-localGI就可以了吗,但是还是不行,换成代码如下:
image
跑结果命令窗:

  1. 这个报错的原因是拼写错误,你写的是-loaclGI
  2. 可以这样写(计算LGI的时候,-i选项可以不用):
export SUBJECTS_DIR=...
recon-all -i input.nii -subject sub001 -all -qcache
recon-all -subject sub001 -localGI
1 个赞

老师,跑了一个,已经跑出来了,在suf文件里多了几个文件,这个就是LGI吗,是不是后期需要自己去平滑?
image

是的,这个文件就相当于是lh.thickness,是个体空间的LGI,如果要做vertex-wise的统计分析,还需要转换到fsaverage空间等。至于要不要平滑,我不太清楚,因为LGI本身已经很平滑了。建议参考一下文献。

好的,谢谢老师,最后一个问题,就是我想用freesurfer根据DK图谱提取每个区域对应的皮层指标(厚度这些)平均值,搜索csdn:


但是我自己输入进去无法实现。是这样直接操作吗

  1. 如果是提取皮层厚度等recon-all默认生成的指标,这个代码是没问题的。唯一可能需要注意的是,代码里的ls.txt表示所有被试文件夹的名字,要确保你有提前准备这个文件。另外,如果你自己的代码没有跑通,应该贴出你自己的代码以及报错,这样才能定位原因。
  2. 如果是提取LGI等不是recon-all默认生成的指标,那么还需要一些额外步骤,可以参考我的博客

好的,老师,参考您的博客应该没啥问题 ,谢谢老师,有问题及时向您请教! :smiling_face:

老师,如何利用freesurfer提取皮层下体积呀,我看了一个文章,下边是这样写的:


,我现在想要提取这个体积,该如何用freesurfer去做呢

  1. FreeSurfer的recon-all计算得到的脑区指标存放在每个被试的stats目录下,其中aseg.stats文件里应该包含你提到的这些脑区体积。
  2. 可以使用asegstats2table命令去合并所有被试的结果到一个文件中,你查看一下该命令的帮助信息(asegstats2table --help)应该就会用了。
  3. 另外,如果是和原帖子不同的话题,最好是新开一个帖子来讨论,这样不用点进帖子就能大概知道内容。

好的,老师,下次提问我注意这个问题,不好意思!我已经使用命令得到皮层下的体积,打开后是这样的:
image
不知道数值范围是不是这样正确的?命令:asegstats2table --subjects `cat list.txt` --meas volume -t volume0.txt --skip;这些提取各个大脑部位的体积都是在同一个标准空间下的吗?都是直接利用的recon-all命令得到的aseg.stats文件。

  1. 具体数值范围我并没有关注过,FreeSurfer能生成这些脑区的体积不表示一定是准确的。可能要结合其他文献来参考。
  2. 命令是没问题的,FreeSurfer的这些指标都是在个体空间计算的。

后期需要用这些指标对每个被试构建矩阵的话,是不是可以直接拿来用呀?而且这篇文中说的是体积测量自动化程序经过MNI配准的,


还是说需要再去配准,应该不需要吧,毕竟前边提取的皮层指标都是基于DK图谱的~~~

FreeSurfer的recon-all流程中会有很多配准的过程,但是这些指标的计算是在个体空间进行的。这些指标可以直接用于统计分析。

好的,我明白啦!谢谢老师 :blush: