第6章医学图象配准61-62.pptVIP

应用逆矩阵运算或组合矩阵变换 应用逆矩阵运算或组合矩阵变换可以很容易实现多种形式的图像配准。 例如，有一个刚体变换矩阵M将第一幅图像A配准到第二幅图像B，则简单的逆矩阵会使第二幅图像B配准到原来的第一幅图像A。 类似地，若有一个刚体变换矩阵将第一幅图像A配准到第二幅图像B，又有另一个刚体变换矩阵将第二幅图像B配准到第三幅图像C，则简单的矩阵乘法得到的刚体变换矩阵会使第一幅图像A直接与第三幅图像C配准。 2. 三维刚体变换 在3D情况，需要6个参数描述一个刚体变换。即沿x轴、y轴及z轴的平移，与绕三个轴的旋转。 变换公式如下： 沿y轴的平移 沿z轴的平移 沿x轴的平移 其中，p，q，r是沿x轴、y轴及z轴的平移刚体变换参数 绕x轴的旋转 绕y轴的旋转 绕z轴的旋转 其中,θ,?,及?是绕x轴、y轴及z轴的旋转刚体变换参数。 沿x轴的平移 沿y轴的平移 沿z轴的平移 绕x轴的旋转 绕y轴的旋转 绕z轴的旋转 (6.26) 3D 变换的矩阵表示 与二维刚体变换类似，旋转与平移的顺序会影响变换的结果，绕不同轴旋转顺序不同，结果也不一样。可能的旋转变换参数的顺序有6种： 结合先旋转、后平移，及先平移、后旋转，共有12种不同的组合顺序。 先旋转、后平移 同样的参数，做逆顺序做变换就会得到完全不同的结果： 3D刚体变换数值举例：一个图像先绕x轴旋转7?，绕y轴旋转11?，再绕z轴旋转13?，接着沿x轴平移2单位，沿y轴平移3单位，最后再沿z轴平移5单位。 旋转的过程图示: 从初始位置到终止位置的变换序列不是唯一的。若沿x轴平移2.2045单位，沿y轴平移2.7793单位，再沿z轴平移5.0414单位,旋转顺序也与原来相反,即先绕z轴旋转14.4198?，绕y轴旋转9.0397?，最后再绕x轴旋转9.4030?。这样变换后得到相同的结果。 6.2.2 全局尺度变换 在3D刚体变换的基础上增加一个新参数m（尺度变换系数）描述不同条件下获取的图像尺度的差异。如果将一幅图像全局缩放m倍的话可用下面的变换： 在顺序的基本变换中，尺度变换与平移变换的顺序对变换结果有影响，但与尺度变换与旋转变换顺序无关。 7个参数的变换矩阵 : 一个图像先绕x轴旋转7?，绕y轴旋转11?，再绕z轴旋转13?，然后全局放大2倍；接着沿x轴平移2单位，沿y轴平移3单位，最后再沿z轴平移5单位。 6.2.3 九参数仿射变换 全局尺度变换是三个方向的缩放倍数相同。但有时待配准的两幅图像也可能各方向上尺度变换系数不同，即 使用九个参数的仿射变换意味待配准的两幅图像存在固有的不对称性。 第一种仿射变换： 对应第二种情况的仿射变换： 6.2.4 一般仿射变换 一般的形式的二维仿射变换有6个独立的参数，三维仿射变换有12个独立的参数。变换的几何约束是变换前后保持直线的平行性。 二维仿射变换： ? 三维仿射变换： ? 由于仿射变换无须计算正弦和余弦，因此对所有变换参数 不再有小于1的约束。 6.2.5 透视变换 主要特点是，直线经变换后仍然是直线，但平行的直线经变换后却可能相交。 在二维情况下，有两个控制透视的参数（f和g）。变换矩阵为 6.2.6 非线性空间变换 1. 二阶多项式变形模型 涉及30个变换参数 2. 基于薄板样条函数 薄板样条函数变换可以表示为仿射变换与径向基函数的线性组合： 3. 基于B样条的变换 基于三次B样条的自由形变变换可以表示为 第6章 医学图像配准-1 基于多种原因，临床上通常需要对同一个病人进行多种模式或同一种模式的多次成像。即同时从几幅图像获得信息，进行综合分析。 使用同种成像设备在不同时间成像，可以观察病灶生长，对比手术前后的治疗效果等。 当一种成像设备所提供的信息不能满足需要时，可以采用多种模式成像。例如，在CT上观察骨组织，而从MRI上得到软组织信息；或将来源于PET、SPECT的功能信息与来源于CT、MRI的解剖信息结合起来分析。 人脑多种模式成像 6.1 图像配准的概述 6.1.1 图象配准的概念 对几幅不同的图像作定量分析，首先要解决这几幅图像的严格对齐问题，这就是我们所说的图像的配准（Image Registration）。 医学图像配准是指对于一幅医学图像寻求一种（或一系列）空间变换，使它与另一幅医学图像上的对应点达到空间上的一致。这种一致是指人体上的同一解剖点在两张