itk matlab,一种itk和matlab形变场双向通道建立方法

一种itk和matlab形变场双向通道建立方法

【技术领域】

[0001] 本发明涉及图像配准领域，尤其是涉及一种ITK和MATLAB形变场双向通道建立方 法。

【背景技术】

[0002] 医学图像配准，即通过寻找一种空间变换，使两幅医学图像上对应的点达到空间 位置及解剖位置的完全一致。配准的结果应使两幅医学图像上所有的解剖点，或至少是所 有具有诊断意义上的点都对应匹配。给定一幅参考图像和一幅目标图像，非参数图像配准 被看作是一个优化问题，目标是去找到每个像素点的位移场，去找到一个图像的最佳变换。 这个变换可以用形变场来表示，表示所有图像上的点从目标图像空间到参考图像空间上的 空间变换。大部分非参数图像配准可以看成，基于相似性准则，通过不同的计算方式不断计 算更新的形变场，并以此来优化形变场。因此，非参数配准的目的也可以说是获得最优的形 变场，而对形变场的分析也必然成为配准工作中的重点。

[0003] 图像配准是通过将模板图像变形为目标图像实现两图像中相同结构间的良好对 应，这种对应性用配准得到的空间形变场来表达。形变场在ITK下可以用nii.gz格式的 文件表达，在MATLAB下可以用矩阵的形式表达。针对ITK下nii. gz格式的形变场文件不 能再MATLAB读取，而且MATLAB矩阵形式的形变场不能在ITK下读取，利用CMake跨平台自 动编译工具，本发明提出了一种ITK和MATLAB形变场双向通道建立方法。

【发明内容】

[0004] 本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种便于读取和分 析、双向通道的ITK和MATLAB形变场双向通道建立方法。

[0005] 本发明的目的可以通过以下技术方案来实现：

[0006] 一种ITK和MATLAB形变场双向通道建立方法，包括以下步骤：

[0007] 1)获取ITK到MATLAB通道CPP文件；

[0008] 2)获取MATLAB到ITK通道CPP文件；

[0009] 3)获取 CMkeList. txt ;

[0010] 4)根据CMkeList. txt中限定的编译条件，对ITK到MATLAB通道CPP文件和 MATLAB到ITK通道CPP文件通过CMake进行交叉编译，得到ITK和MATLAB形变场双向通 道；

[0011] 5)检验所建立的ITK和MATLAB形变场双向通道的准确性。

[0012] 所述的步骤1)具体包括以下步骤：

[0013] 11)获取ITK下形变场文件的参数：

[0014] 读取ITK下形变场的区域变量、尺寸大小和索引参数信息，根据参数信息定位并 获取形变场数据；

[0015] 12)初始化MATLAB形变场矩阵；

[0016] 13)根据ITK下形变场的参数信息为初始化的MATLAB形变场矩阵赋值。

[0017] 所述的步骤2)具体包括以下步骤：

[0018] 21)定义和初始化形变场以及缓冲图像：

[0019] 定义形变场在ITK下的形变场文件的参数，包括形变场区域变量、尺寸大小和索 弓丨，初始化形变场，并初始化缓冲图像；

[0020] 22)形变场数据交换，在ITK下初始化形变场文件之后，将MATLAB下矩阵形式的数 据赋值给形变场文件；

[0021] 23)获取形变场文件。

[0022] 所述的步骤3)中的CMakeList. txt文件包括对生成目标、编译的标示和选项声 明。

[0023] 所述的步骤5)具体包括以下步骤：

[0024] 51)利用基于ITK配准程序配准目标图像和参考图像，产生原始形变场；

[0025] 52)根据建立的ITK到MATLAB通道读取原始形变场，得到转换后的形变场；

[0026] 53)分别在MATLAB平台采用转换后的形变场将测试图像去形变，以及在ITK框架 下采用原始形变场去形变，判断两种方法形变后的图像是否相同，若是，则通道建立正确， 若否，则通道建立不正确，返回步骤1)，重新建立通道。

[0027] 所述的步骤5)包括以下步骤：

[0028] 54)基于MATLAB配准程序配准目标图像和参考图像，产生原始形变场；

[0029] 55)根据建立的MATLAB到ITK通道写出原始形变场，得到转换后的形变场；

[0030] 56)分别在MATLAB平台采用原始形变场将测试图像去形变，以及在ITK框架下采 用转换后的形变场去形变，判断两种方法形变后的图像是否相同，若是，则通道建立正确， 若否，则通道建立不正确，返回步骤1)，重新建立通道。

[0031] 与现有技术相比，本发明具有以下优点：

[0032] 1)便于读取和分析：本发明利用CMake交叉编译，建立了 ITK到MATLAB形变场通 道，使得ITK框架下的nii. gz格式的形变场能够在MATLAB中读取并进行分析，比如形变场 任意区域的叠加。

[0033] 2)双向通道：本发明利用CMake交叉编译，建立了 MATLAB到ITK形变场通道，使 得MATLAB下产生的形变场数据能写为ITK下能读取的nii. gz格式的文件，该通道能够将 ITK下和MATLAB下产生的不同形式的形变场相互转换，对MATLAB和ITK平台配准算法形变 场分析有重要的意义。

【附图说明】

[0034] 图1为本发明的流程示意图。

[0035] 图2为实施例的ITK到MATLAB形变场通道实验图。

[0036] 图3为实施例的MATLAB到ITK形变场通道实验图。

【具体实施方式】

[0037] 下面结合附图和具体实施例对本发明进行详细说明。

[0038] 实施例：

[0039] 一、运行环境：

[0040] 该方法基于MATLAB平台，针对医学图像配准算法产生的形变场，建立了 ITK和 MATLAB形变场双向通道一ITK到MATLAB形变场通道和MATLAB到ITK形变场通道。

[0041] 二、实施步骤：

[0042] 如图1所示，一种ITK和MATLAB形变场双向通道建立方法，该方法建立了 ITK和 MTLAB形变场双向通道，具体包括步骤：

[0043] 步骤1.编写ITK到MATLAB通道CPP文件：读取ITK下形变场文件参数；初始化 MTLAB下形变场矩阵；读出形变场文件数据并赋值。

[0044] a)读取ITK下形变场文件参数；

[0045] 首先读取ITK下形变场的区域变量、尺寸大小和索引，有了这些信息才能正确定 位形变场数据并提取出来。此过程主要代码如下：

[0046] //读取形变场区域变量、尺寸大小和索引

【主权项】

1. 一种ITK和MATLAB形变场双向通道建立方法，其特征在于，包括以下步骤： 1) 获取ITK到MATLAB通道CPP文件； 2) 获取MATLAB到ITK通道CPP文件； 3) 获取 CMkeList. txt ; 4) 根据CMkeList. txt中限定的编译条件，对ITK到MATLAB通道CPP文件和MATLAB到 ITK通道CPP文件通过CMake进行交叉编译，得到ITK和MATLAB形变场双向通道； 5) 检验所建立的ITK和MATLAB形变场双向通道的准确性。

2. 根据权利要求1所述的一种ITK和MATLAB形变场双向通道建立方法，其特征在于， 所述的步骤1)具体包括以下步骤： 11) 获取ITK下形变场文件的参数： 读取ITK下形变场的区域变量、尺寸大小和索引参数信息，根据参数信息定位并获取 形变场数据； 12) 初始化MATLAB形变场矩阵； 13) 根据ITK下形变场的参数信息为初始化的MTLAB形变场矩阵赋值。

3. 根据权利要求1所述的一种ITK和MATLAB形变场双向通道建立方法，其特征在于， 所述的步骤2)具体包括以下步骤： 21) 定义和初始化形变场以及缓冲图像： 定义形变场在ITK下的形变场文件的参数，包括形变场区域变量、尺寸大小和索引，初 始化形变场，并初始化缓冲图像； 22) 形变场数据交换，在ITK下初始化形变场文件之后，将MATLAB下矩阵形式的数据赋 值给形变场文件； 23) 获取形变场文件。

4. 根据权利要求1所述的一种ITK和MATLAB形变场双向通道建立方法，其特征在于， 所述的步骤3)中的CMakeList. txt文件包括对生成目标、编译的标示和选项声明。

5. 根据权利要求1所述的一种ITK和MATLAB形变场双向通道建立方法，其特征在于， 所述的步骤5)具体包括以下步骤： 51) 利用基于ITK配准程序配准目标图像和参考图像，产生原始形变场； 52) 根据建立的ITK到MATLAB通道读取原始形变场，得到转换后的形变场； 53) 分别在MATLAB平台采用转换后的形变场将测试图像去形变，以及在ITK框架下采 用原始形变场去形变，判断两种方法形变后的图像是否相同，若是，则通道建立正确，若否， 则通道建立不正确，返回步骤1)，重新建立通道。

6. 根据权利要求1所述的一种ITK和MATLAB形变场双向通道建立方法，其特征在于， 所述的步骤5)包括以下步骤： 54) 基于MATLAB配准程序配准目标图像和参考图像，产生原始形变场； 55) 根据建立的MATLAB到ITK通道写出原始形变场，得到转换后的形变场； 56) 分别在MATLAB平台采用原始形变场将测试图像去形变，以及在ITK框架下采用转 换后的形变场去形变，判断两种方法形变后的图像是否相同，若是，则通道建立正确，若否， 则通道建立不正确，返回步骤1)，重新建立通道。

【专利摘要】本发明涉及一种ITK和MATLAB形变场双向通道建立方法，包括以下步骤：1)获取ITK到MATLAB通道CPP文件；2)获取MATLAB到ITK通道CPP文件；3)获取CMkeList.txt；4)根据CMkeList.txt中限定的编译条件，对ITK到MATLAB通道CPP文件和MATLAB到ITK通道CPP文件通过CMake进行交叉编译，得到ITK和MATLAB形变场双向通道；5)检验所建立的ITK和MATLAB形变场双向通道的准确性。与现有技术相比，本发明具有便于读取和分析、双向通道等优点。

【IPC分类】G06T7-00

【公开号】CN104835171

【申请号】CN201510259692

【发明人】刘刚, 赵龙, 朱凯, 张万

【申请人】上海电力学院

【公开日】2015年8月12日

【申请日】2015年5月20日