肺部 ct 的三维配准方法
【专利摘要】一种肺部CT的三维配准方法，其包括以下步骤：载入第一肺部CT体数据和第二肺部CT体数据；以互信息为配准测度，对第一肺部CT体数据和第二肺部CT体数据进行三维刚性配准，使二者中的肺部区域整体对齐，记录X轴、Y轴和Z轴三个方向的位移参数和旋转参数，作为第一配准参数；在刚性配准的基础上进行非刚性配准，记录非刚性配准的最优位移场，作为第二配准参数；用第一配准参数和第二配准参数对第一肺部CT体数据和第二肺部CT体数据中的结节进行匹配，建立对应关系。本发明方法能够实现精确配准，并利用配准结果将同一位置的结节病灶在不同CT数据中的成像位置进行自动对应，从而能够方便医生对结节病情发展进行对比。
【专利说明】肺部CT的三维配准方法
【技术领域】
[0001]本发明属于计算机视觉及图像处理【技术领域】，具体涉及一种肺部CT的三维配准方法。
【背景技术】
[0002]根据WH02008年的全球肿瘤流行病统计数据(GL0B0CAN2008)，仅2008年全球肿瘤病例为1270万，死亡760万例[I]，其中，中国所有的肿瘤新病例为280万，是美国新发肿瘤病例的2倍(美国为140万)；同年中国肿瘤死亡病例为196万，而美国的死亡病例为57万，二者的差距为3.4倍。恶性肿瘤占我国居民死亡原因的20%以上。据WHO国际肿瘤研究理事会负责人Bernard Stewart报道,到2020年全球肿瘤发病率还将上升50%。
[0003]提高肿瘤病人的五年存活率的关键在于如何尽早的发现并确诊肿瘤病情，并进行治疗。随着CT技术的发展，高分辨率CT因为其良好的成像质量，成为肿瘤筛查的重要工具，被越来越广泛地在临床上采用，但是随之而来的问题是医生阅片工作量的增加，尤其是在对比阅片时在大量CT数据中手动进行相同解剖位置图像匹配。
[0004]因此利用计算机辅助诊断工具帮助医生高效准确地对肿瘤病情进行诊断有着巨大的现实意义。目前相关的计算机辅助诊断工具的研究都集中如何从一套CT数据中进行肺结节的识别和分割，以及根据特征对结节的良恶性进行判断。在肿瘤的诊断和治疗中，对比同一位置病灶特征的变化有着重要的临床意义，而在如何提高医生在进行病情对比时的阅片效率方面的研究很少。
[0005]在一些PACS (Picture Archiving and Communication Systems)系统中，虽然提供了自动定位相近断层图像的功能，但是这种功能精确度不高，并且只是进行了断层与断层之间的匹配，如果一个断层中出现多个结节，这种方法并不能提供有效的对应结果。
【发明内容】
[0006]本发明的目的是提供一种肺部CT的三维配准方法，该三维配准方法能够将同一位置的结节在不同CT数据中的成像位置进行自动对应。
[0007]本发明中提出的配准技术应用的对象是高分辨率肺部CT图像，其技术思路为:首先进行刚性配准，通过互信息将CT体数据中的肺部区域进行整体对齐；然后在刚性配准的基础上利用提取的匹配特征点对结合局部单点互信息(SMI)最大化的方式进行非刚性配准，得到最终的精确配准结果；然后利用配准结果提供的空间变换信息，进行相同结节成像位置的对应，从而方便结节的特征变化进行对比。
[0008]本发明的具体技术方案为:
[0009]一种肺部CT的三维配准方法，包括以下步骤:
[0010]载入第一肺部CT体数据和第二肺部CT体数据；
[0011]以互信息为配准测度，对第一肺部CT体数据和第二肺部CT体数据进行三维刚性配准，使二者中的肺部区域整体对齐，记录X轴、Y轴和Z轴三个方向的位移参数和旋转参数，作为第一配准参数；
[0012]在刚性配准的基础上进行非刚性配准，记录非刚性配准的最优位移场，作为第二配准参数；
[0013]用第一配准参数和第二配准参数对第一肺部CT体数据和第二肺部CT体数据中的结节进行匹配，建立对应关系。
[0014]在上述的肺部CT的三维配准方法中，优选地，所述非刚性配准包括特征点提取、特征点匹配和位移场搜索。
[0015]在上述的肺部CT的三维配准方法中，优选地,采用三维的Fdrstner;算子提取特征点，特征点匹配中被匹配为特征点对的两个特征点r、f满足以下两个条件:一、参考图像中的特征点r在浮动图像中最邻近的点是特征点f，且特征点f在参考图像中的最邻近的点也是特征点r ;二、最邻近特征点与次临近特征点间的距离比小于预设的第一阈值。
[0016]在上述的肺部CT的三维配准方法中，优选地，采用BBF算法确定最邻近的特征点。
[0017]在上述的肺部CT的三维配准方法中，优选地，所述非刚性配准用互信息结合特征点限制作为相似性测度。
[0018]在上述的肺部CT的三维配准方法中，优选地，所述位移场搜索包括以下步骤:计算位移场变化量；将位移场变化量与位移场相加，更新位移场；利用特征点对位移场进行限制，包括计算特征点X处的位移场Un+1 (X)，计算该位移场Un+1 (X)与V(X)的差值，该差值小于预设的第二阈值时，使位移场Un+1 (x)) V (X)，否则Un+1 (x)保持不变，V (x)为构成特征点对的两个特征点之间的位移，该特征点对包括特征点X ;计算浮动图像插值结果，计算该插值结果与参考图像之间的互信息，当互信息停止增大或迭代次数达到限定次数时，将当前位移场作为最优位移场，否则重复上述步骤。
[0019]在上述的肺部CT的三维配准方法中，优选地，更新位移场前还包括平滑位移场变化量的步骤，在利用特征点对位移场进行限制后还包括对位移场进行加权高斯平滑的步骤，所述加权高斯平滑在特征点处采用第一权值，在其它位置采用第二权值。
[0020]在上述的肺部CT的三维配准方法中，优选地，采用对称搜索的方式计算位移场变化量。
[0021]在上述的肺部CT的三维配准方法中，优选地，该三维配准方法还包括在进行配准前将第一肺部CT体数据和第二肺部CT体数据每个体素代表的物理尺寸降采样到相同尺度上。
[0022]在上述的肺部CT的三维配准方法中，优选地，采用浮动图像与参考图像之间的联合灰度分布直方图计算互信息，其中采用PV插值的方式来计算联合灰度分布直方图。
[0023]本发明方法能够实现精确配准，并利用配准结果将同一位置的结节病灶在不同CT数据中的成像位置进行自动对应，能够方便医生对结节病情发展进行对比，从而能够大大提闻医生的诊断效率。
【专利附图】
【附图说明】
[0024]图1为一些实施例肺部CT的三维配准方法的流程图；
[0025]图2为采用P&rstnei算子提取特征点匹配后的结果；
[0026]图3为浮动图像配准后的插值图像与参考图像相减图。【具体实施方式】
[0027]下面结合附图和实施例对本发明进一步说明。这些更详细的描述旨在帮助理解本发明，而不应被用于限制本发明。根据本发明公开的内容，本领域技术人员明白，可以不需要一些或者所有这些特定细节即可实施本发明。而在其它情况下，为了避免将发明创造淡化，未详细描述众所周知的操作过程。
[0028]如图1所示，本肺部CT的三维配准方法包括以下步骤:
[0029]步骤S1、载入第一肺部CT体数据和第二肺部CT体数据。第一肺部CT体数据和第二肺部CT体数据是同一被测者在不同时间或不同检测条件下的肺部CT体数据。
[0030]步骤S2、以互信息为配准测度，对第一肺部CT体数据和第二肺部CT体数据进行三维刚性配准，使二者中的肺部区域整体对齐，记录X轴、Y轴和Z轴三个方向的位移参数和旋转参数，作为第一配准参数。
[0031]三维刚性配准采用互信息作为相似性测度，因此不要对肺实质进行预分割处理。互信息的定义式如下:
【权利要求】
1.一种肺部CT的三维配准方法，其特征在于，该三维配准方法包括以下步骤: 载入第一肺部CT体数据和第二肺部CT体数据； 以互信息为配准测度，对第一肺部CT体数据和第二肺部CT体数据进行三维刚性配准，使二者中的肺部区域整体对齐，记录X轴、Y轴和Z轴三个方向的位移参数和旋转参数，作为第一配准参数； 在刚性配准的基础上进行非刚性配准，记录非刚性配准的最优位移场，作为第二配准参数； 用第一配准参数和第二配准参数对第一肺部CT体数据和第二肺部CT体数据中的结节进行匹配，建立对应关系。
2.根据权利要求1所述的肺部CT的三维配准方法，其特征在于，所述非刚性配准包括特征点提取、特征点匹配和位移场搜索。
3.根据权利要求2所述的肺部CT的三维配准方法，其特征在于,采用三维的Far伽w算子提取特征点，特征点匹配中被匹配为特征点对的两个特征点r、f满足以下两个条件:一、参考图像中的特征点r在浮动图像中最邻近的点是特征点f，且特征点f在参考图像中的最邻近的点也是特征点r ;二、最邻近特征点与次临近特征点间的距离比小于预设的第一阈值。
4.根据权利要 求3所述的肺部CT的三维配准方法，其特征在于，采用BBF算法确定最邻近的特征点。
5.根据权利要求2所述的肺部CT的三维配准方法，其特征在于，所述非刚性配准用互信息结合特征点限制作为相似性测度。
6.根据权利要求2所述的肺部CT的三维配准方法，其特征在于，所述位移场搜索包括以下步骤: 计算位移场变化量； 将位移场变化量与位移场相加，更新位移场； 利用特征点对位移场进行限制，包括计算特征点X处的位移场Un+1 (x)，计算该位移场Un+1(x)与V(X)的差值，该差值小于预设的第二阈值时，使位移场Un+1 (X) =V(X)否则Un+1(x)保持不变，V(X)为构成特征点对的两个特征点之间的位移,该特征点对包括特征点X ； 计算浮动图像插值结果，计算该插值结果与参考图像之间的互信息，当互信息停止增大或迭代次数达到限定次数时，将当前位移场作为最优位移场，否则重复上述步骤。
7.根据权利要求6所述的肺部CT的三维配准方法，其特征在于，更新位移场前还包括平滑位移场变化量的步骤，在利用特征点对位移场进行限制后还包括对位移场进行加权高斯平滑的步骤，所述加权高斯平滑在特征点处采用第一权值，在其它位置采用第二权值。
8.根据权利要求6所述的肺部CT的三维配准方法，其特征在于，采用对称搜索的方式计算位移场变化量。
9.根据权利要求1所述的肺部CT的三维配准方法，其特征在于，该三维配准方法还包括在进行配准前将第一肺部CT体数据和第二肺部CT体数据每个体素代表的物理尺寸降采样到相同尺度上。
10.根据权利要求1所述的肺部CT的三维配准方法，其特征在于，采用浮动图像与参考图像之间的联合灰度分布直方图计算互信息，其中采用PV插值的方式来计算联合灰度分布直方 图。
【文档编号】G06T7/00GK104021547SQ201410209850
【公开日】2014年9月3日 申请日期:2014年5月17日 优先权日:2014年5月17日
【发明者】袁克虹, 于川汇 申请人:清华大学深圳研究生院