使用SPM_batch进行批量跑脚本（matlab.m）
软件：spm8+matlab2023b+win11
数据格式：
F:\ASL\HC\CBF\HC_caishaoqing\CBF.nii
F:\ASL\HC\CBF\HC_caishaoqing\T1.nii
F:\ASL\HC\CBF\HC_wangdonga\CBF.nii
F:\ASL\HC\CBF\HC_wangdonga\T1.nii
clear
spmdir='D:\AnalysisApps\spm8';
datadir='F:\ASL\HC\CBF';
cd(datadir);
allsubs=dir('HC*');for i=1:numel(allsubs)cursub=allsubs(i).name;matlabbatch{1}.spm.spatial.normalise.estwrite.subj.source = {fullfile(datadir,cursub,'CBF.nii')};matlabbatch{1}.spm.spatial.normalise.estwrite.subj.wtsrc = '';matlabbatch{1}.spm.spatial.normalise.estwrite.subj.resample = {fullfile(datadir,cursub,'CBF.nii')};matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.template = {'D:\AnalysisApps\spm8\templates\PET.nii,1'};matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.weight = '';matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.smosrc = 8;matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.smoref = 0;matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.regtype = 'mni';matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.cutoff = 25;matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.nits = 16;matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.reg = 1;matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.preserve = 0;matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.bb = [-78 -112 -5078 76 85];matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.vox = [2 2 2];matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.interp = 1;matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.wrap = [0 0 0];matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.prefix = 'w';spm_jobman('run',matlabbatch)
end
后面的脚本来自spm-view-show.m code
下面这里是show.m产生的原始脚本。可以看出我们只是把绿色部分用循环语句替换了。
matlabbatch{1}.spm.spatial.normalise.estwrite.subj.source = {'F:\ASL\HC\CBF\HC_caishaoqing\CBF.nii,1'};
matlabbatch{1}.spm.spatial.normalise.estwrite.subj.wtsrc = '';
matlabbatch{1}.spm.spatial.normalise.estwrite.subj.resample = {'F:\ASL\HC\CBF\HC_caishaoqing\CBF.nii,1'};
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.template = {'D:\AnalysisApps\spm8\templates\PET.nii,1'};
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.weight = '';
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.smosrc = 8;
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.smoref = 0;
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.regtype = 'mni';
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.cutoff = 25;
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.nits = 16;
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.reg = 1;
matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.preserve = 0;
matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.bb = [-78 -112 -5078 76 85];
matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.vox = [2 2 2];
matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.interp = 1;
matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.wrap = [0 0 0];
matlabbatch{1}.spm.spatial.normalise.estwrite.roptions.prefix = 'w';
相关文章：
使用SPM_batch进行批量跑脚本（matlab.m）
软件：spm8matlab2023bwin11 数据格式： F:\ASL\HC\CBF\HC_caishaoqing\CBF.nii F:\ASL\HC\CBF\HC_caishaoqing\T1.nii F:\ASL\HC\CBF\HC_wangdonga\CBF.nii F:\ASL\HC\CBF\HC_wangdonga\T1.nii clear spmdirD:\AnalysisApps\spm8; datadirF:\ASL\HC\CBF…...
力扣0124——二叉树的最大路径和
二叉树的最大路径和 难度：困难 题目描述 二叉树中的 路径 被定义为一条节点序列，序列中每对相邻节点之间都存在一条边。同一个节点在一条路径序列中 至多出现一次 。该路径 至少包含一个 节点，且不一定经过根节点。 路径和 是路径中各节点…...
c# 字符串帮助类
public class StringHelper { #region 全角半角互相转换 /// <summary> /// 转全角的函数(SBC case) /// </summary> /// <param name"str">任意字符串</param> /// <returns>全…...
LabVIEW双光子荧光显微成像系统开发
双光子显微成像是一种高级荧光显微技术，广泛用于生物学和医学研究，尤其是用于活体组织的深层成像。在双光子成像过程中，振镜（Galvo镜）扮演了非常关键的角色，它负责精确控制激光束在样本上的扫描路径。以下是…...
Prim模板
通过代码探索Prim算法：最小生成树之旅 在计算机科学领域，图算法占据了至关重要的位置，尤其是在设计高效的网络（无论是社交网络、计算机网络还是交通网）时。在这些算法中，寻找最小生成树（MST&am…...
CSS之盒子模型
盒子模型 01-选择器 结构伪类选择器 基本使用 作用：根据元素的结构关系查找元素。 <!DOCTYPE html> <html lang"en"><head><meta charset"UTF-8"><meta http-equiv"X-UA-Compatible" content"IE…...
Linux系统安装（CentOS Vmware）
学习环境安装 VMware安装 VMware下载&安装 访问官网：https://www.vmware.com 在此处可以选择语言 点击China（简体中文） 点击产品，点击Workstation Pro 下滑，点击下载试用版 下滑找到Workstation 17 Pro for Wi…...
STM32 硬件随机数发生器(RNG)
STM32 硬件随机数发生器 文章目录 STM32 硬件随机数发生器前言第1章 随机数发生器简介1.1 RNG主要特性1.2.RNG应用 第2章 RNG原理框图第3章 RNG相关寄存器3.1 RNG 控制寄存器 (RNG_CR)3.2 RNG 状态寄存器 (RNG_SR)3.3 RNG 数据寄存器 (RNG_DR) 第3章 RNG代码部分第4章 STM32F1 …...
Window环境下使用go编译grpc最新教程
网上的grpc教程都或多或少有些老或者有些问题，导致最后执行生成文件时会报很多错。这里给出个人实践出可执行的编译命令与碰到的报错与解决方法。（ps:本文代码按照煎鱼的教程编写：4.2 gRPC Client and Server - 跟煎鱼学 Go (gitbook.io)&…...
STM32——FLASH（1）简单介绍、分类、读写流程及注意事项
文章目录 FLASH的特点Nor flash和nand flashflash的读写flash 的存储单位 flash的读写过程 FLASH的特点 可擦写数据可修改可重写访问速度<ROM Nor flash和nand flash Nor flash 1、与SDRAM相似，用户可以直接运行装载到NORFLASH里面的代码，减少SRAM…...
MySQL的DML语言
DML：Data Manipulation Language（数据操作语言） DML语言用来对数据库中表的数据记录进行增、删、改操作。 一、添加数据命令 注意： 插入数据时，指定的字段顺序需要与值的顺序是一一对应的。 字符串和日期型数据应该包…...
Vivado-IP核
Vivado-IP核 主程序 timescale 1ns / 1ps ////module ip_clk_wiz(input sys_clk,input sys_rst_n,output clk_out1,output clk_out2,output clk_out3,output clk_out4,output locked);clk_wiz_0 instance_name(// Clock out ports.clk_out1(clk_out1), // output clk_out…...
品牌如何营造生活感氛围？媒介盒子分享
「生活感」简而言之是指人们对生活的感受和意义，它往往没有充斥在各种重要的场合和事件中，而是更隐藏在细碎平凡的生活场景中。在营销越来越同质化的当下，品牌应该如何打破常规模式，洞察消费情绪，找到更能打动消费者心…...
Java 学习和实践笔记（2）
今天的学习进度： 注册并下载安装好了Java 8，之后进行以下配置。 1）path 是一个常见的环境变量，它告诉系统除了在当前的目标下妹寻找此程序外，还可以到path指定的目录下找。这句话是什么意思呢？以下举报例…...
Python：批量url链接保存为PDF
我的数据是先把url链接获取到存入excel中，后续对excel做的处理，各位也可以直接在程序中做处理，下面就是针对excel中的链接做批量处理 excel内容格式如下（涉及具体数据做了隐藏） 标题文件链接文件日期网页标题1http://…...
【LeetCode每日一题】525连续数组 303区域和检索（前缀和的基本概念和3个简单案例）
前缀和 // 构造prefix let prefix [0] arr.forEach(num > {prefix.push(prefix.at(-1) num); })如果想要计算某个区间 i 到 j 这个子数组的和时，可以根据 prefix[j1] - prefix[i] 获得。 例题1：303.区域和检索 - 数组不可变 给定一个整数数组 num…...
形态学算法应用之连通分量提取的python实现——图像处理
原理 连通分量提取是图像处理和计算机视觉中的一项基本任务，旨在识别图像中所有连通区域，并将它们作为独立对象处理。在二值图像中，连通分量通常指的是所有连接在一起的前景像素集合。这里的“连接”可以根据四连通或八连通的邻接关系来定义…...
Kafka系列之：Kafka集群同时设置基于时间和日志大小两种方式保存Topic的数据
Kafka系列之：Kafka集群同时设置基于时间和日志大小两种方式保存Topic的数据 一、基于日志大小二、基于时间大小三、参数设置四、设置命令一、基于日志大小 "log.retention.bytes"是Apache Kafka中的一项配置参数，用于指定每个日志段文件的最大大小。当日志段文件的…...
pytest+allure批量执行测试用例
在 Pytest 中，可以使用装饰器 `@pytest.fixture` 来定义用例级别的前置和后置操作。下面是一个示例代码，演示了如何使用 Pytest 的前置和后置操作： ```python import pytest @pytest.fixture(scope="function") def setup_function(): print("Setup fu…...
SpringBoot和SpringMVC
目录 一、springboot项目 （1）创建springboot项目 （2）目录介绍 （3）项目启动 （4）运行一个程序 （5）通过其他方式创建和运行springboot项目 二、SpringMVC…...
React 第五十五节 Router 中 useAsyncError的使用详解
前言 useAsyncError 是 React Router v6.4 引入的一个钩子，用于处理异步操作（如数据加载）中的错误。下面我将详细解释其用途并提供代码示例。 一、useAsyncError 用途 处理异步错误：捕获在 loader 或 action 中发生的异步错误替…...
地震勘探——干扰波识别、井中地震时距曲线特点
目录 干扰波识别反射波地震勘探的干扰波 井中地震时距曲线特点 干扰波识别 有效波：可以用来解决所提出的地质任务的波；干扰波：所有妨碍辨认、追踪有效波的其他波。 地震勘探中，有效波和干扰波是相对的。例如，在反射波…...
C++：std::is_convertible
C++标志库中提供is_convertible，可以测试一种类型是否可以转换为另一只类型： template <class From, class To> struct is_convertible; 使用举例： #include <iostream> #include <string>using namespace std;struct A { }; struct B : A { };int main…...
3.3.1_1 检错编码（奇偶校验码）
从这节课开始，我们会探讨数据链路层的差错控制功能，差错控制功能的主要目标是要发现并且解决一个帧内部的位错误，我们需要使用特殊的编码技术去发现帧内部的位错误，当我们发现位错误之后，通常来说有两种解决方案。第一…...
视频字幕质量评估的大规模细粒度基准
大家读完觉得有帮助记得关注和点赞！！！ 摘要 视频字幕在文本到视频生成任务中起着至关重要的作用，因为它们的质量直接影响所生成视频的语义连贯性和视觉保真度。尽管大型视觉-语言模型（VLMs）在字幕生成方面…...
Linux云原生安全：零信任架构与机密计算
Linux云原生安全：零信任架构与机密计算 构建坚不可摧的云原生防御体系 引言：云原生安全的范式革命 随着云原生技术的普及，安全边界正在从传统的网络边界向工作负载内部转移。Gartner预测，到2025年，零信任架构将成为超…...
令牌桶 滑动窗口-＞限流 分布式信号量-＞限并发的原理 lua脚本分析介绍
文章目录 前言限流限制并发的实际理解限流令牌桶代码实现结果分析令牌桶lua的模拟实现原理总结： 滑动窗口代码实现结果分析lua脚本原理解析 限并发分布式信号量代码实现结果分析lua脚本实现原理 双注解去实现限流 并发结果分析： 实际业务去理解体会统一注…...
拉力测试cuda pytorch 把 4070显卡拉满
import torch import timedef stress_test_gpu(matrix_size16384, duration300):"""对GPU进行压力测试，通过持续的矩阵乘法来最大化GPU利用率参数:matrix_size: 矩阵维度大小，增大可提高计算复杂度duration: 测试持续时间（秒&…...
Rapidio门铃消息FIFO溢出机制
关于RapidIO门铃消息FIFO的溢出机制及其与中断抖动的关系，以下是深入解析： 门铃FIFO溢出的本质 在RapidIO系统中，门铃消息FIFO是硬件控制器内部的缓冲区，用于临时存储接收到的门铃消息（Doorbell Message）。…...
#Uniapp篇：chrome调试unapp适配
chrome调试设备----使用Android模拟机开发调试移动端页面 Chrome://inspect/#devices MuMu模拟器Edge浏览器：Android原生APP嵌入的H5页面元素定位 chrome://inspect/#devices uniapp单位适配 根路径下 postcss.config.js 需要装这些插件 “postcss”: “^8.5.…...