医学图像笔记(七)医学图像配准

1、度量标准

从Elastix手册中查到
Releases · SuperElastix/elastix (github.com)
在这里插入图片描述
在这里插入图片描述

总结
MSD度量是一种仅适用于具有相等强度分布的两个图像的度量,即来自相同模态的图像。
NCC相对宽松,它假设固定图像和运动图像的强度值之间存在线性关系,因此可以更频繁地使用。
MI度量更加普遍:仅假设固定图像和运动图像的强度的概率分布之间的关系。众所周知,对于多模态图像,它不仅适用于单模态图像,也适用于多模态图像对,这种方法通常是图像配准的好选择。
NMI度量就像MI一样,适用于单模态和多模态配准。
KS测量特别意味着记录二进制图像(分割)。它测量分段的“重叠”。

参考:图像配准常用评价方法

2、DTI配准之Slicer (Elastix or BRAI NS)配准篇

多模态影像的配准原则:

①优先以分辨率高、层厚薄T1为基准图像;

②若CT层薄,可以先以薄层CT为基准图像配准T1,然后以配准后的T1像为基准图像配准其他系列影像;
③优先相同扫描仪器、相同扫描系列的配准(如CTA与CT配准,T1与T1增强等)。

配准的原理:

①刚体变换(rigid body transformation):只需要经过空间的平移和旋转,核磁图像中,空间可沿着x、y、z轴平移或旋转,因此刚性变换只需要6个自由度。

②仿射变换(affine transformation):除了需要空间平移和旋转以外,还需要图像的拉伸(Scaling)如放大、缩小,和图像的倾斜(Skews/Shears),图像的拉伸和倾斜也有6个自由度。因此,仿射变换的自由度是12个。
③非线性变换:除了平移、旋转、拉伸、倾斜外,还需要图像局部的形变,故需要的自由度在12个以上。
DTI由包含多个方向的扫描序列,故配准比较特殊。可以通过Slicer自带的General Registration(BRAINS)及插件General Registration(Elastix)配准。还可以通过SPM联合BRAINS或SPM联合Elastix配准。

参考:DTI配准之Slicer (Elastix or BRAI NS)配准篇

### SPM12、NCCT 和 MNI 的转换与标化指南 #### 安装与置 SPM12 是一款广泛应用于神经影像学研究的软件工具包,主要用于处理和分析脑成像数据。要使用 SPM12 进行 NCCT 图像以及 MNI 标空间的操作,需先确保 MATLAB 已经正确安装并激活了必要的工具箱。 对于 SPM12 的安装,可以从官方网站下载最新版本,并按照官方文档说明完成设置过程[^1]。一旦成功安装,可以通过启动 MATLAB 并运行 `spm` 命令来加载 SPM 用户界面。 #### 数据预处理流程 针对 NCCT 影像,在导入至 SPM12 后通常会经历如下几个主要阶段: - **重采样 (Resampling)**:调整原始扫描的空间分辨率以匹目标模板的要求。 - **分割 (Segmentation)**:将 T1 加权 MRI 划分为灰质、白质和其他成分,这一步骤同样适用于 CT 扫描中的软组织对比度差异较大的情况[^2]。 - ** (Registration/Normalization)**:通过变形场映射个体解剖结构到标坐标系(如 MNI),从而实现不同受试者间的一致性比较。此过程中可以指定特定参数控制最终输出体素尺寸,比如 3×3×3 mm³ 来降低计算复杂度而不显著牺牲精度。 #### 实现具体操作实例 下面给出一段用于批量执行上述步骤的脚本示例,特别是如何读取本地存储路径下的 PET 或其他模态图像文件列表,并应用预先定义好的作业流(`jobfile`)来进行自动化处理: ```matlab % 设置工作目录及获取文件名列表 global sub filelist pathname; pathname = 'C:\Users\86191\Desktop\AD数据整理'; filelist = dir(fullfile(pathname, '*')); filelist = {filelist([filelist.isdir]' == false).name}'; for i = 1:length(filelist) sub = filelist{i}; % 置批处理任务 nrun = 1; % 设定运行次数 jobfile = {'C:\Users\86191\Desktop\2.0ruijinprocess\petmni_job.m'}; jobs = repmat(jobfile, 1, nrun); inputs = cell(0, nrun); spm('defaults', 'PET'); spm_jobman('run', jobs, inputs{:}); end ``` 这段代码展示了怎样遍历指定文件夹内的所有符合条件的医学影像文件,并调用已有的作业描述文件 (`petmni_job.m`) 对这些文件实施一系列预定的预处理措施,包括但不限于向 MNI 空间的变换和标化[^3]。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

落花逐流水

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值