当前位置：首页 > news >正文

使用SPM_batch进行批量跑脚本（matlab.m）

news 2025/12/16 16:51:54

软件：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…...

编程日记 2024/2/10 2:13:22

力扣0124——二叉树的最大路径和

二叉树的最大路径和难度：困难题目描述二叉树中的路径被定义为一条节点序列，序列中每对相邻节点之间都存在一条边。同一个节点在一条路径序列中至多出现一次。该路径至少包含一个节点，且不一定经过根节点。路径和是路径中各节点…...

编程日记 2024/2/10 2:07:16

c# 字符串帮助类

public class StringHelper { #region 全角半角互相转换 /// <summary> /// 转全角的函数(SBC case) /// </summary> /// <param name"str">任意字符串</param> /// <returns>全…...

编程日记 2024/2/10 2:06:15

LabVIEW双光子荧光显微成像系统开发

双光子显微成像是一种高级荧光显微技术，广泛用于生物学和医学研究，尤其是用于活体组织的深层成像。在双光子成像过程中，振镜（Galvo镜）扮演了非常关键的角色，它负责精确控制激光束在样本上的扫描路径。以下是…...

编程日记 2024/2/10 1:59:08

Prim模板

通过代码探索Prim算法：最小生成树之旅在计算机科学领域，图算法占据了至关重要的位置，尤其是在设计高效的网络（无论是社交网络、计算机网络还是交通网）时。在这些算法中，寻找最小生成树（MST&am…...

编程日记 2024/2/10 1:57:07

CSS之盒子模型

盒子模型 01-选择器结构伪类选择器基本使用作用：根据元素的结构关系查找元素。 <!DOCTYPE html> <html lang"en"><head><meta charset"UTF-8"><meta http-equiv"X-UA-Compatible" content"IE…...

编程日记 2024/2/10 1:56:06

Linux系统安装（CentOS Vmware）

学习环境安装 VMware安装 VMware下载&安装访问官网：https://www.vmware.com 在此处可以选择语言点击China（简体中文） 点击产品，点击Workstation Pro 下滑，点击下载试用版下滑找到Workstation 17 Pro for Wi…...

编程日记 2024/2/10 1:46:54

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 …...

编程日记 2024/2/10 1:45:53

Window环境下使用go编译grpc最新教程

网上的grpc教程都或多或少有些老或者有些问题，导致最后执行生成文件时会报很多错。这里给出个人实践出可执行的编译命令与碰到的报错与解决方法。（ps:本文代码按照煎鱼的教程编写：4.2 gRPC Client and Server - 跟煎鱼学 Go (gitbook.io)&…...

编程日记 2024/2/10 1:41:48

STM32——FLASH（1）简单介绍、分类、读写流程及注意事项

文章目录 FLASH的特点Nor flash和nand flashflash的读写flash 的存储单位 flash的读写过程 FLASH的特点可擦写数据可修改可重写访问速度<ROM Nor flash和nand flash Nor flash 1、与SDRAM相似，用户可以直接运行装载到NORFLASH里面的代码，减少SRAM…...

编程日记 2024/2/10 1:36:44

MySQL的DML语言

DML：Data Manipulation Language（数据操作语言） DML语言用来对数据库中表的数据记录进行增、删、改操作。一、添加数据命令注意： 插入数据时，指定的字段顺序需要与值的顺序是一一对应的。字符串和日期型数据应该包…...

编程日记 2024/2/10 1:30:39

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…...

编程日记 2024/2/10 1:27:36

品牌如何营造生活感氛围？媒介盒子分享

「生活感」简而言之是指人们对生活的感受和意义，它往往没有充斥在各种重要的场合和事件中，而是更隐藏在细碎平凡的生活场景中。在营销越来越同质化的当下，品牌应该如何打破常规模式，洞察消费情绪，找到更能打动消费者心…...

编程日记 2024/2/10 1:25:34

Java 学习和实践笔记（2）

今天的学习进度： 注册并下载安装好了Java 8，之后进行以下配置。 1）path 是一个常见的环境变量，它告诉系统除了在当前的目标下妹寻找此程序外，还可以到path指定的目录下找。这句话是什么意思呢？以下举报例…...

编程日记 2024/2/10 1:24:33

Python：批量url链接保存为PDF

我的数据是先把url链接获取到存入excel中，后续对excel做的处理，各位也可以直接在程序中做处理，下面就是针对excel中的链接做批量处理 excel内容格式如下（涉及具体数据做了隐藏） 标题文件链接文件日期网页标题1http://…...

编程日记 2024/2/10 1:23:32

【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…...

编程日记 2024/2/10 1:21:30

形态学算法应用之连通分量提取的python实现——图像处理

原理连通分量提取是图像处理和计算机视觉中的一项基本任务，旨在识别图像中所有连通区域，并将它们作为独立对象处理。在二值图像中，连通分量通常指的是所有连接在一起的前景像素集合。这里的“连接”可以根据四连通或八连通的邻接关系来定义…...

编程日记 2024/2/10 1:20:29

Kafka系列之：Kafka集群同时设置基于时间和日志大小两种方式保存Topic的数据

Kafka系列之：Kafka集群同时设置基于时间和日志大小两种方式保存Topic的数据一、基于日志大小二、基于时间大小三、参数设置四、设置命令一、基于日志大小 "log.retention.bytes"是Apache Kafka中的一项配置参数，用于指定每个日志段文件的最大大小。当日志段文件的…...

编程日记 2024/2/10 1:18:27

pytest+allure批量执行测试用例

在 Pytest 中，可以使用装饰器 `@pytest.fixture` 来定义用例级别的前置和后置操作。下面是一个示例代码，演示了如何使用 Pytest 的前置和后置操作： ```python import pytest @pytest.fixture(scope="function") def setup_function(): print("Setup fu…...

编程日记 2024/2/10 1:15:24

SpringBoot和SpringMVC

目录一、springboot项目 （1）创建springboot项目 （2）目录介绍 （3）项目启动 （4）运行一个程序 （5）通过其他方式创建和运行springboot项目二、SpringMVC…...

编程日记 2024/2/10 1:11:20

ESP32读取DHT11温湿度数据

芯片：ESP32 环境：Arduino 一、安装DHT11传感器库红框的库，别安装错了二、代码注意，DATA口要连接在D15上 #include "DHT.h" // 包含DHT库#define DHTPIN 15 // 定义DHT11数据引脚连接到ESP32的GPIO15 #define D…...

编程新知 2025/12/16 5:07:06

P3 QT项目----记事本（3.8）

3.8 记事本项目总结项目源码 1.main.cpp #include "widget.h" #include <QApplication> int main(int argc, char *argv[]) {QApplication a(argc, argv);Widget w;w.show();return a.exec(); } 2.widget.cpp #include "widget.h" #include &q…...

编程新知 2025/12/12 10:35:10

SpringBoot+uniapp 的 Champion 俱乐部微信小程序设计与实现，论文初版实现

摘要本论文旨在设计并实现基于 SpringBoot 和 uniapp 的 Champion 俱乐部微信小程序，以满足俱乐部线上活动推广、会员管理、社交互动等需求。通过 SpringBoot 搭建后端服务，提供稳定高效的数据处理与业务逻辑支持；利用 uniapp 实现跨平台前…...

编程新知 2025/12/16 7:09:20

Mac软件卸载指南，简单易懂！

刚和Adobe分手，它却总在Library里给你写"回忆录"？卸载的Final Cut Pro像电子幽灵般阴魂不散？总是会有残留文件，别慌！这份Mac软件卸载指南，将用最硬核的方式教你"数字分手术"&#xff0…...

编程新知 2025/12/14 10:18:18

鸿蒙中用HarmonyOS SDK应用服务 HarmonyOS5开发一个生活电费的缴纳和查询小程序

一、项目初始化与配置 1. 创建项目 ohpm init harmony/utility-payment-app 2. 配置权限 // module.json5 {"requestPermissions": [{"name": "ohos.permission.INTERNET"},{"name": "ohos.permission.GET_NETWORK_INFO"…...

编程新知 2025/11/26 13:15:57

select、poll、epoll 与 Reactor 模式

在高并发网络编程领域，高效处理大量连接和 I/O 事件是系统性能的关键。select、poll、epoll 作为 I/O 多路复用技术的代表，以及基于它们实现的 Reactor 模式，为开发者提供了强大的工具。本文将深入探讨这些技术的底层原理、优缺点。一、I…...

编程新知 2025/12/12 20:44:23

OpenLayers 分屏对比(地图联动)

注：当前使用的是 ol 5.3.0 版本，天地图使用的key请到天地图官网申请，并替换为自己的key 地图分屏对比在WebGIS开发中是很常见的功能，和卷帘图层不一样的是，分屏对比是在各个地图中添加相同或者不同的图层进行对比查看。…...

编程新知 2025/9/30 5:23:30

CSS设置元素的宽度根据其内容自动调整

width: fit-content 是 CSS 中的一个属性值，用于设置元素的宽度根据其内容自动调整，确保宽度刚好容纳内容而不会超出。效果对比默认情况（width: auto）： 块级元素（如 <div>）会占满父容器…...

编程新知 2025/10/16 16:03:09

保姆级教程：在无网络无显卡的Windows电脑的vscode本地部署deepseek

文章目录 1 前言2 部署流程2.1 准备工作2.2 Ollama2.2.1 使用有网络的电脑下载Ollama2.2.2 安装Ollama（有网络的电脑）2.2.3 安装Ollama（无网络的电脑）2.2.4 安装验证2.2.5 修改大模型安装位置2.2.6 下载Deepseek模型 2.3 将deepse…...

编程新知 2025/12/16 15:06:54

AirSim/Cosys-AirSim 游戏开发（四）外部固定位置监控相机

这个博客介绍了如何通过 settings.json 文件添加一个无人机外的固定位置监控相机，因为在使用过程中发现 Airsim 对外部监控相机的描述模糊，而 Cosys-Airsim 在官方文档中没有提供外部监控相机设置，最后在源码示例中找到了，所以感…...

编程新知 2025/12/11 8:45:01

相关文章：