图像配准概述

article/2025/8/5 1:30:39

图像配准在医学图像领域是一项重要的技术,在许多的临床诊断中,为了分析患者的状况,常常需要采集患者的扫描影像,例如,X线、MRI、CT和超声,这些扫描影像可以对患者的诊断提供依据,然而,在独立的影像中,临床信息并不是十分明显,所以采用单模或者多模的方法用以获得额外的临床信息,出于此目的,发现这些影像间的空间关系就变得很重要。图像配准就是找到一幅图像像素到另一幅图像像素间的空间映射关系。

配准过程需要两幅图像,一幅图像(moving image)$ I_M(x) $ 形变到另一幅图像(fixed image)$ I_F(x) $ ,并且这两幅图像都是在各自的空间域上定义且具有相同的维度$ d , 配 准 就 是 要 找 到 位 移 ,配准就是要找到位移 u(x) 使 得 使得 使 I_M(x+u(x)) $ 空间对其到$ I_F(x) , 用 一 个 等 价 的 公 式 表 达 就 是 找 到 一 个 变 换 ,用一个等价的公式表达就是找到一个变换 T(x)=x+u(x) $ 使得$ I_M(T(x))$ 空间对其到 $ I_F(x) , 这 个 变 换 定 义 为 从 ,这个变换定义为从 fixed image $ 到 $ moving image $ 的映射
左边是fixed image 右边是moving image

这里写图片描述
$ T:\Omega_F \subset \Bbb R^d \to \Omega_M \subset \Bbb R^d $ 对齐的质量通过距离或者相似性测度 $ S 定 义 , 例 如 平 方 差 和 ( 定义,例如平方差和( SSD ) 、 相 关 比 率 、 互 信 息 ( )、相关比率、互信息( MI ) , 由 于 这 个 问 题 不 适 用 于 非 刚 性 变 换 ),由于这个问题不适用于非刚性变换 T$,通常引入正则化或惩罚系数 $ \mathcal P$ 来约束$ T$
一般来讲,配准问题是一个使得代价函数 $ \mathcal C $ 最小的优化问题,用公式表达$ T$:
T ^ = arg ⁡ min ⁡ T C ( T ; T F , I M ) \hat T = \arg \min_{T} \mathcal C (\boldsymbol{T};T_F,I_M) T^=argTminC(T;TF,IM)
C ( T ; T F , I M ) = − S ( T ; I F , I M ) + γ P ( T ) \mathcal C (\boldsymbol{T};T_F,I_M) = -\mathcal S (\boldsymbol{T};I_F,I_M)+ \gamma \mathcal P(\boldsymbol{T}) C(T;TF,IM)=S(T;IF,IM)+γP(T)

$ \gamma $ 权重因子 为了解决上述最小化问题,有两个基本的方法:参数化和非参数化,这里不讨论非参数化方法,在参数化方法中,引入变换的参数化模型来约束可能变换的数量,原来的优化问题就变成了:
T ^ = arg ⁡ min ⁡ T u C ( T u ; T F , I M ) \hat T = \arg \min_{T_u} \mathcal C (\boldsymbol{T_u};T_F,I_M) T^=argTuminC(Tu;TF,IM)

下标 u u u 表示变换已被参数化,向量 u u u 包含变换参数值,例如:当该变换为2D刚性变换,参数向量 u u u 包含一个旋转角度和 x x x y y y方向上的位移,所以,等式也可以等价于:
u ^ = arg ⁡ min ⁡ u C ( u ; T F , I M ) \hat u = \arg \min_{u} \mathcal C (\boldsymbol{u};T_F,I_M) u^=arguminC(u;TF,IM)
通过这个等式,原来的问题就变得简化了,通过优化元素 u u u取代了优化变换函数T。

图像
相似性测度

在文献中可以发现相似性测度的几种共同选择,接下来描述几点共同的选择

均方差($ MSD$):定义

M S D ( μ ; I F , I M ) = 1 ∣ Ω F ∣ ∑ x i ∈ Ω F ( I F ( x i ) − I M ( T μ ( x i ) ) ) 2 MSD(\mu; I_F,I_M)=\frac{1} {|\Omega_F|}\sum_{x_i\in\Omega_F}(I_F(x_i)-I_M(\boldsymbol T_\mu(x_i)))^2 MSD(μ;IF,IM)=ΩF1xiΩF(IF(xi)IM(Tμ(xi)))2

$ \Omega_F $ 为 fixed image 的定义域,$ |\Omega_F| $ 为像素的数量,给定变换$ \boldsymbol T$,这个度量可以通过循环迭代fixed image中的像素来实现

归一化互相关(NCC)

NCC定义为:
N C C ( μ ; I F , I M ) = ∑ x i ∈ Ω F ( I F ( x i ) ) − I F ‾ ) ( I M ( T μ ( x i ) ) − I M ‾ ) ∑ x i ∈ Ω F ( I F ( x i ) − I F ‾ ) 2 ∑ x i ∈ Ω F ( I M ( T μ ( x i ) ) − I M ‾ ) 2 NCC(\mu;I_F,I_M)=\frac{\sum_{x_i \in \Omega_F}(I_F(x_i))-\overline{I_F})(I_M(\boldsymbol T_\mu(x_i))-\overline{I_M})}{ \sqrt {\sum_{x_i \in \Omega_F}(I_F(x_i)-\overline{I_F})^2 \sum_{x_i \in \Omega_F}(I_M(\boldsymbol T_\mu(x_i))-\overline{I_M})^2}} NCC(μ;IF,IM)=xiΩF(IF(xi)IF)2xiΩF(IM(Tμ(xi))IM)2 xiΩF(IF(xi))IF)(IM(Tμ(xi))IM)


http://chatgpt.dhexx.cn/article/nPXq7xgZ.shtml

相关文章

基于MATLAB的医学图像配准算法仿真

目录 一、理论基础 二、案例背景 1.问题描述 2.思路流程 三、部分MATLAB程序 四、仿真结论分析 五、算法相关应用 六、参考文献 一、理论基础 其中h表示二维空间坐标变换,g表示灰度或辐射变换,描述因传感器类型的不同或辐射变形所引入的图像变换…

图像配准简介

图像配准在目标检测、模型重建、运动估计、特征匹配,肿瘤检测、病变定位、血管造影、地质勘探、航空侦察等领域都有广泛的应用。 每一种配准方法通常都针对某个具体问题而设计的,众多方法中,唯一的共性就是每个配准问题最终都要在变换空间中寻找一种最有的变换,这种变换能…

基于深度学习的单模医学图像配准综述(附VoxelMorph配准实例)

本文是基于深度学习的单模态医学图像配准的综述,除了介绍配准任务、配准过程之外,还会从实际操作出发,以经典的VoxelMorph为例做详细介绍。如果有什么讲的不清楚的地方欢迎大家留言讨论,如果有什么错误的地方,也恳请大…

图像配准

1、定义 维基百科上的定义:图像配准与相关是图像处理研究领域中的一个典型问题和技术难点,其目的在于比较或融合针对同一对象在不同条件下获取的图像,例如图像会来自不同的采集设备,取自不同的时间,不同的拍摄视角等等…

【Image Registration】图像配准综述

文章目录 一、图像配准定义二、图像配准应用场景2.1 医学图像领域2.2 其他领域 三、图像配准分类四、图像配准过程4.1 特征检测(Feature detection)4.2 特征匹配(Feature matching)4.2.1 基于区域的方法(Area-based me…

计算机视觉 什么是图像配准?

一、图像配准概述 图像配准是叠加两个或多个来自不同来源、在不同时间和角度拍摄的图像的过程。图像配准过程是一种自动或手动操作,它试图发现两张照片之间的匹配点并在空间上对齐它们以最小化所需的误差,即两幅图像之间的统一邻近度测量。医学、遥感和计算机视觉都使用图像配…

1.图像配准简述

1.图像配准:作为图像融合的一个预处理步骤,图像配准是对两幅图像,通过寻找一种空间变换把一幅图像映射到另一幅图像,使得两图中对应于空间同一位置的点一一对应起来,从而达到信息融合的目的。 2.原理: 3.分…

机器视觉(九):图像配准

目录: 机器视觉(一):概述 机器视觉(二):机器视觉硬件技术 机器视觉(三):摄像机标定技术 机器视觉(四):空域图像增强 …

【图像处理】什么是图像配准?

一、图像配准概述 图像配准是叠加两个或多个来自不同来源、在不同时间和角度拍摄的图像的过程。图像配准过程是一种自动或手动操作,它试图发现两张照片之间的匹配点并在空间上对齐它们以最小化所需的误差,即两幅图像之间的统一邻近度测量。医学、遥感和计…

图解矩阵的秩

1.图解矩阵的秩 1.图解矩阵的秩1.1 满秩矩阵1.2 非满秩矩阵1.3 零矩阵1.4 秩的性质1.4.1 满秩矩阵复合的性质1.4.2 一般矩阵复合的性质 1.图解矩阵的秩 笔记来源于:《马同学图解线性代数》 详细过程: 1.1 满秩矩阵 1.2 非满秩矩阵 1.3 零矩阵 1.4 秩…

java 矩阵求秩_线性代数精华3——矩阵的初等变换与矩阵的秩

矩阵的初等变换这个概念可能在很多人听来有些陌生,但其实我们早在初中的解多元方程组的时候就用过它。只不过在课本当中,这种方法叫做消元法。我们先来看一个课本里的例子: 假设我们要解这个方程,怎么做呢? 首先&#…

矩阵秩的定义和相关结论汇总

秩的定义:对于矩阵,以下陈述为真。(如果,则用共轭转置替换下述转置) rank(A)矩阵A经过行初等变换,所得行阶梯形矩阵的非零行数rank(A)矩阵A经过行初等变换,所得行阶梯形矩阵的主元数rank(A)矩阵…

秩为1的矩阵的性质总结

网上东拼西凑找到的,只能用手抄下来,方便自己复习看👀。有什么不对的还希望大家指出!

秩一矩阵的优良性质

前言:仅个人小记 秩一矩阵非常漂亮的五个性质: (1)秩一矩阵一定能够拆解为两个列向量 a ⃗ \vec{a} a , b ⃗ \vec{b} b 矩阵乘积的形式,具体为 A a ⃗ b ⃗ T A\vec{a}{\vec{b}}^{T} Aa b T这种形式 &…

伴随矩阵秩的证明

伴随矩阵是引出n阶矩阵逆计算的一个重要矩阵工具,课本中关于伴随矩阵的涉及的讲解并不多,本文将从一下方面讲解: 基础重要计算公式(引入逆的求解);通过和逆的关系求解伴随矩阵以及证明相关运算律&#xff…

使用 Amazon Amplify快速创建简单的 Android 应用程序

背景: 亚马逊云科技提供了100余种产品免费套餐。其中,计算资源Amazon EC2首年12个月免费,750小时/月;存储资源 Amazon S3 首年12个月免费,5GB标准存储容量。 大家好,我是坚果,由于最近一直疫情…

基于Android的个人时间管理设计与开发

选题依据(包括目的、意义、国内外现状和发展趋势,主要参考文献): 课题的意义随着科学技术高速发展,手机普遍率越来越高,大学生甚至小学生几乎人人拥有一台智能手机,移动端应用发展也越来越快&am…

基于Android平台的酒店预订管理系统软件设计的论文

基于Android平台的酒店预订管理系统软件设计 摘要 随着智能手机的快速普及,智能手机操作系统市场风生水起。为了让智能手机用户能够随时随地查询互联网所提供的服务,一种高效的办法就是将应用系统的功能拓展到手机终端上,让手机能够通过移动网…

Android Studio项目结构详解

Android Studio 是谷歌推出的一个Android集成开发工具,基于IntelliJ IDEA. 类似 Eclipse ADT,Android Studio 提供了集成的 Android 开发工具用于开发和调试。 在IDEA的基础上,Android Studio 提供: 1.基于Gradle的构建支持。2.An…