利用freesurfer6进行海马分割的环境配置和步骤,以及获取海马体积

利用freesurfer6进行海马分割的环境配置和步骤

    • Matlab Runtime 安装
    • 1. 运行recon-all:
    • 2. 利用 recon-all -s subj -hippocampal-subfields-T1 进行海马分割
    • 3. 结束后需要在/$SUBJECTS_DIR/subject/的文件夹/mri路径下输入下面的代码查看分割情况
    • 4. 在文件SUBJECTS_DIR路径下输入quantifyHippocampalSubfields.sh 收集所有被试的海马体积
    • 参考文献

在这里插入图片描述

Matlab Runtime 安装

一些 freesurfer 模块,例如 fsfast 和子场/核分割工具,需要 matlab 运行时包 (MCR)。 MCR 允许用户运行分布式 matlab 编译的程序,而无需支付 matlab 许可证费用。 freesurfer 命令 fs_install_mcr 可用于安装此软件包。 MCR版本要求如下。

在 freesurfer 7.0 以上版本时,请通过运行以下命令安装 MCR 8.4 (R2014b):

fs_install_mcr R2014b

在 freesurfer 6.0 版本时,请通过运行以下命令安装 MCR 8.0 (R2012b):

fs_install_mcr R2012b

注意:如果 fs_install_mcr 脚本在您的 freesurfer 发行版中不可用,可以通过运行以下命令来下载:

cd $FREESURFER_HOME/bin && curl https://raw.githubusercontent.com/freesurfer/freesurfer/dev/scripts/fs_install_mcr -o fs_install_mcr && chmod +x fs_install_mcr

注意如果安装不成功,请切换到 root 用户下。

1. 运行recon-all:

#!/usr/bin/env bash
export SUBJECTS_DIR=/media/data/HC_dataset/ProOutput/freeT1w
T1wpath=/media/data/HC_dataset/HCBIDS
anatpath=/media/data/HC_dataset/ProOutput/anatfor subj in `cat group1.txt`
docp $T1wpath/$subj/anat/*.nii $anatpath
done#并行运行
ls $anatpath/*.nii | parallel --jobs 60 recon-all -s {.} -i {} -all -qcache

2. 利用 recon-all -s subj -hippocampal-subfields-T1 进行海马分割

#!/usr/bin/env bash
export SUBJECTS_DIR=/media/data/HC_dataset/ProOutput/freeT1w
T1wpath=/media/data/HC_dataset/HCBIDS
anatpath=/media/data/HC_dataset/ProOutput/anatfor subj in `cat group1.txt`
dorecon-all -s $subj -hippocampal-subfields-T1
done

3. 结束后需要在/$SUBJECTS_DIR/subject/的文件夹/mri路径下输入下面的代码查看分割情况

freeview -v nu.mgz -v lh.hippoSfLabels-T1.v10.mgz:colormap=lut -v rh.hippoSfLabels-T1.v10.mgz:colormap=lut

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

4. 在文件SUBJECTS_DIR路径下输入quantifyHippocampalSubfields.sh 收集所有被试的海马体积

#!/usr/bin/env bash
export SUBJECTS_DIR=/media/data/HC_dataset/ProOutput/freeT1w
T1wpath=/media/data/HC_dataset/HCBIDS
anatpath=/media/data/HC_dataset/ProOutput/anatquantifyHippocampalSubfields.sh T1 hippocampal_volume.txt

在这里插入图片描述
在这里插入图片描述
Abbreviations: CA, cornu ammonis; DG, dentate gyrus; GC-DG, granule cell layer of dentate
gyrus; HATA, hippocampal-amygdaloid transition region.

参考文献

A computational atlas of the hippocampal formation using ex vivo, ultra-high resolution MRI: Application to adaptive segmentation of in vivo MRI. Iglesias, J.E., Augustinack, J.C., Nguyen, K., Player, C.M., Player, A., Wright, M., Roy, N., Frosch, M.P., McKee, A.C., Wald, L.L., Fischl, B., and Van Leemput, K. Neuroimage, 115, July 2015, 117-137.

Genon, S., Bernhardt, B. C., La Joie, R., Amunts, K., & Eickhoff, S. B. (2021). The many dimensions of human hippocampal organization and (dys) function. Trends in neurosciences, 44(12), 977-989.

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/96113.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

910数据结构(2019年真题)

算法设计题 问题1 有一种排序算法叫做计数排序。这种排序算法对一个待排序的表(采用顺序存储)进行排序,并将排序结果存放到另一个新的表中。必须注意的是,表中所有待排序的关键字互不相同,计数排序算法针对表中的每个元素,扫描待排序的表一趟,统计表中有多少个元素的关…

C++转换函数

什么是转换函数? C转换函数是一种特殊的成员函数,用于将一个类的对象转换为另一个类型。它是通过在类中定义特定的函数来实现的。 转换函数的用途: 类型转换:转换函数可以将一个类的对象从一种类型转换为另一种类型。这样可以方便地在不同…

使用maven框架搭建一个IDEA插件项目

以下是使用 Maven 框架搭建 IDEA 插件项目的步骤: 打开 IDEA,点击 File -> New -> Project,选择 Maven。 在弹出的 New Project 窗口中,选择 Maven,然后选择 Create from archetype,找到 Maven 插件…

Java 时间的加减处理

时间的加减处理 Date date new Date(操作时间(类型Date)-(60000*60*1));600001分钟 60000*60*1 1小时

A*算法和Dijkstra

A*算法 https://www.redblobgames.com/pathfinding/a-star/introduction.html这是个宝藏网页,https://www.redblobgames.com/pathfinding/a-star/introduction.html,里边的图可以一步一步演示! A*算法 个人理解FGH,F是总距离&a…

win10睡眠快捷方式

新建快捷方式 如下图 内容如下 rundll32.exe powrprof.dll,SetSuspendState 0,1,0 下一步 点击完成即可。 特此记录 anlog 2023年10月6日

一年一度的国庆节又结束了

这里写自定义目录标题 欢迎使用Markdown编辑器新的改变功能快捷键合理的创建标题,有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants 创建一个自定义列表如何创建一个…

C++中的static和extern关键字

1 声明和定义 声明就是告诉编译器有这个东西的存在,而定义则是这个东西的实现。 对于变量来说,声明就是告诉编译器存在这个名称的变量,定义则是给这个变量分配内存并赋值: // 变量声明,声明时不能赋值,如…

万字总结网络原理

目录 一、网络基础 1.1认识IP地址 1.2子网掩码 1.3认识MAC地址 1.4一跳一跳的网络数据传输 1.5总结IP地址和MAC地址 二、网络设备及相关技术 2.1集线器:转发所有端口 2.2交换机:MAC地址转换表+转发对应端口 2.3主机:网络分层从上到下封装 2.4主机&路由器:ARP…

(三)正点原子STM32MP135移植——optee移植

一、概述 OP-TEE 是一个开源工程,完整的实现了一个可信执行环境。 主要包括 Secure world OS(optee_os)、normal world client(optee_client)、test suite(optee_test/xtest)以及 Linux 驱动部分…

软件项目和安全项目案例(承接软件和安全项目合作)

公司有专业的软件开发团队和安全研究团队,具备完善的安全测试、安全培训、安全开发、安全服务等安全解决方案,可以助力政企研发专业、高效、安全、稳定的软件产品,欢迎项目咨询、商务合作! 一、软件开发项目咨询 1.承接车载等终…

Iphone文件传到电脑用什么软件,看这里

在数字化时代,文件传输已经成为我们日常生活中不可或缺的一部分。然而,苹果用户在将手机文件传输到电脑时,往往会面临一些困扰。曾经的“文件传输助手”并不能完全满足用户的需求。于是,很多人开始寻找更便捷的解决方案。在本文中…

【CMU15-445 Part-15】Query Planning Optimization II

Part15-Query Planning & Optimization II Selection Statistics 维护每张表中的基本主要信息也就是tuple数量 N R N_R NR​以及每个属性中不同值的数量 V ( A , R ) V(A,R) V(A,R), N R N_R NR​关系R中的元组数量,单独维护,不能用pag…

Python 无废话-基础知识函数详解

函数定义 函数是一段可重复使用的代码块,用于实现特定的功能。 类似的前面已学过print(),len(),input(),str(),list(),tuple()等内置函数 自定义函数 内置函数是python类库的已定义好的函数,如果某些功能不能满足时,需要我们自定义函数&am…

《TCP/IP网络编程》代码实现

文章目录 1. 项目说明1.1 项目特点2. 文件说明2.1 脚本文件2.1.1 `TCP_IP.sln`2.1.2 `xmake.lua`2.1.2.1 编译说明2.1.2.2 运行说明2.1.3 章节说明项目代码已经开源在github上! 微信公众号文章同步发表! 1. 项目说明 根据《TCP/IP网络编程》书籍学习,对其中的代码进行整理,…

alsa pcm接口之在unix环境的传输方法

在unix环境,数据片段响应被接受通过standard I/O call或事件等待路径(poll或select功能),为完成列表,异步通知响应该被列举出来.ALSA实现那些方法被描述在ALSA transfers部分. 标准I/O传输(Standadrd I/O transfers) 这个标准I/O传输常常使用read和write C语言函数集,对于那些函…

网页版”高德地图“如何设置默认城市?

问题: 每次打开网页版高德地图时默认定位的都是“北京”,想设置起始点为目前本人所在城市,烦恼的是高德地图默认的初始位置是北京。 解决: 目前网页版高德地图暂不支持设置起始点,打开默认都是北京,只能将…

华为OD机试真题-补种未成活胡杨(Java/C++/Go/Python)

华为OD机试真题-补种未成活胡杨(Java/C++/Go/Python) 题目描述 近些年来,我国防沙治沙取得显著成果。某沙漠新种植N棵胡杨(编号1-N),排成一排。 一个月后,有M棵胡杨未能成活。现可补种胡杨K棵,请问如何补种(只能补种,不能新种),可以得到最多的连续胡杨树? 输入…

假期题目整合

1. 下载解压题目查看即可 典型的猪圈密码只需要照着输入字符解开即可得到答案 2. 冷门类型的密码题型,需要特意去找相应的解题思路,直接百度搜索天干地支解密即可 3. 一眼能出思路他已经给了篱笆墙的提示提示你是栅栏密码对应解密即可 4. 最简单的社会主…

SLAM面试笔记(8) — 计算机视觉面试题

目录 问题1:目标检测的算法分类 问题2:卷积神经网络的组成 问题3:输入层的作用 问题4:卷积层作用 问题5:卷积核类型 问题6:11卷积核作用 问题7:卷积核是否越大越好 问题8:棋…