matlab 读写磁共振影像.nii 数据

2023-10-19 12:52

本文主要是介绍matlab 读写磁共振影像.nii 数据,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

问题描述:

  参考LEiDA方法和周昌松老师的large-scale network 的分离整合的划分[1],脑区划分结束后需要把脑区的结果用脑图的形式呈现出来,但是因为我划分脑区的结果都是保存在mat文件里,Brainnet Viewer 则需要使用 .nii 类型的文件。因此,研究了一下如何将数据写入.nii文件。

在这里插入图片描述

加载BrainNet Viewer的样例数据

nifti_data = load_nii('D:\Program Files\ThirdTools\BrainNetViewer_20191031\Data\ExampleFiles\AAL90\aal.nii');

数据结构

nifti_data = load_nii('D:\Program Files\ThirdTools\BrainNetViewer_20191031\Data\ExampleFiles\AAL90\aal.nii');% Extract the voxel values
voxel_values = nifti_data.img(:);% Filter out zero values
non_zero_values = voxel_values(voxel_values ~= 0);% Calculate basic statistics
mean_value = mean(non_zero_values);
median_value = median(non_zero_values);
min_value = min(non_zero_values);
max_value = max(non_zero_values);

通过代码,也可以判断数据存储在nifti_data.img中,1表示的是脑区1,2表示脑区2,3表示脑区3…,以此类推。

如果只是划分脑网络,那么就只需要将脑区x的值x替换为所属的脑区y,就可以实现脑区划分,

如果需要达到参考文献C图的效果,以参考文献中的第4层8个网络为例,我们需要将img中保存的116个脑区的所属的脑区id替换为该脑区对应的特征向量。

实际操作过程中遇到了数据无法保存的问题,原因是样例aal.nii文件的.img 格式为unit8,表示无符号,因此,我们写入的负数、小数值会被强制转换为无符号整数,因此,修改数据前需要先修改数据类型,数据类型的具体描述在save_nii.m 文件中有详细说明。

%%----*** 保存每个不同level的network***---
function write_tonii(level)disp(level);% Load the original NIfTI filenifti_data = load_nii('D:\Program Files\ThirdTools\BrainNetViewer_20191031\Data\ExampleFiles\AAL90\aal.nii');% 设置数据格式,数据格式在save_nii.m文件中有详细说明,此处选则浮点型,代码为16nifti_data.hdr.dime.datatype=16;nifti_data.hdr.dime.bitpix=16;%将模板文件中的unit8改成浮点型,然后写入数据modified_data = double(nifti_data.img);%加载我的数据,脑区与特征向量值的对照表,这个表格的数据是从mat文件中读取后整理成csv格式了,也可以直接整理成mat格式保存dataTable = readtable(['F:\projects\LEiDA-BD\df_sorted_desc_',num2str(level),'.csv'], 'Delimiter', ',');% 写代码,将相应的脑区的数据更改为对应的值然后保存。for i = 1:116indices = find(modified_data ==i);% 根据脑区值获得对应的特征向量值保存到 nii.img 中eigen_value =  dataTable.eigen_value(dataTable.region_id == i);modified_data(indices) =eigen_value;end% Update the data in the NIfTI structurenifti_data.img = double(modified_data);% Save the modified NIfTI file with the same format and informationsave_file = ['aal_level',num2str(level),'_eigen_network.nii'];save_nii(nifti_data, save_file);

然后,使用BrainNet Viewer查看结果,如下图所示。
在这里插入图片描述

参考文献:
[1] Wang, et al. (2021). Segregation, integration, and balance of large-scale resting brain networks configure different cognitive abilities. Proc Natl Acad Sci U S A, 118(23). doi:10.1073/pnas.2022288118

这篇关于matlab 读写磁共振影像.nii 数据的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/239929

相关文章

Python获取中国节假日数据记录入JSON文件

《Python获取中国节假日数据记录入JSON文件》项目系统内置的日历应用为了提升用户体验,特别设置了在调休日期显示“休”的UI图标功能,那么问题是这些调休数据从哪里来呢?我尝试一种更为智能的方法:P... 目录节假日数据获取存入jsON文件节假日数据读取封装完整代码项目系统内置的日历应用为了提升用户体验,

SpringBoot实现数据库读写分离的3种方法小结

《SpringBoot实现数据库读写分离的3种方法小结》为了提高系统的读写性能和可用性,读写分离是一种经典的数据库架构模式,在SpringBoot应用中,有多种方式可以实现数据库读写分离,本文将介绍三... 目录一、数据库读写分离概述二、方案一:基于AbstractRoutingDataSource实现动态

Java利用JSONPath操作JSON数据的技术指南

《Java利用JSONPath操作JSON数据的技术指南》JSONPath是一种强大的工具,用于查询和操作JSON数据,类似于SQL的语法,它为处理复杂的JSON数据结构提供了简单且高效... 目录1、简述2、什么是 jsONPath?3、Java 示例3.1 基本查询3.2 过滤查询3.3 递归搜索3.4

MySQL大表数据的分区与分库分表的实现

《MySQL大表数据的分区与分库分表的实现》数据库的分区和分库分表是两种常用的技术方案,本文主要介绍了MySQL大表数据的分区与分库分表的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有... 目录1. mysql大表数据的分区1.1 什么是分区?1.2 分区的类型1.3 分区的优点1.4 分

Mysql删除几亿条数据表中的部分数据的方法实现

《Mysql删除几亿条数据表中的部分数据的方法实现》在MySQL中删除一个大表中的数据时,需要特别注意操作的性能和对系统的影响,本文主要介绍了Mysql删除几亿条数据表中的部分数据的方法实现,具有一定... 目录1、需求2、方案1. 使用 DELETE 语句分批删除2. 使用 INPLACE ALTER T

Python Dash框架在数据可视化仪表板中的应用与实践记录

《PythonDash框架在数据可视化仪表板中的应用与实践记录》Python的PlotlyDash库提供了一种简便且强大的方式来构建和展示互动式数据仪表板,本篇文章将深入探讨如何使用Dash设计一... 目录python Dash框架在数据可视化仪表板中的应用与实践1. 什么是Plotly Dash?1.1

Redis 中的热点键和数据倾斜示例详解

《Redis中的热点键和数据倾斜示例详解》热点键是指在Redis中被频繁访问的特定键,这些键由于其高访问频率,可能导致Redis服务器的性能问题,尤其是在高并发场景下,本文给大家介绍Redis中的热... 目录Redis 中的热点键和数据倾斜热点键(Hot Key)定义特点应对策略示例数据倾斜(Data S

Python实现将MySQL中所有表的数据都导出为CSV文件并压缩

《Python实现将MySQL中所有表的数据都导出为CSV文件并压缩》这篇文章主要为大家详细介绍了如何使用Python将MySQL数据库中所有表的数据都导出为CSV文件到一个目录,并压缩为zip文件到... python将mysql数据库中所有表的数据都导出为CSV文件到一个目录,并压缩为zip文件到另一个

Redis中高并发读写性能的深度解析与优化

《Redis中高并发读写性能的深度解析与优化》Redis作为一款高性能的内存数据库,广泛应用于缓存、消息队列、实时统计等场景,本文将深入探讨Redis的读写并发能力,感兴趣的小伙伴可以了解下... 目录引言一、Redis 并发能力概述1.1 Redis 的读写性能1.2 影响 Redis 并发能力的因素二、

SpringBoot整合jasypt实现重要数据加密

《SpringBoot整合jasypt实现重要数据加密》Jasypt是一个专注于简化Java加密操作的开源工具,:本文主要介绍详细介绍了如何使用jasypt实现重要数据加密,感兴趣的小伙伴可... 目录jasypt简介 jasypt的优点SpringBoot使用jasypt创建mapper接口配置文件加密