图像配准装置及其图像配准的方法
【专利摘要】一种图像配准装置，包括：获取模块获取人脑计算机断层(CT)图像；分割模块对人脑CT图像进行脑分割以得到分割图像；插值模块对分割图像进行插值获得插值图像，其中获取模块还获取磁共振(MRI)脑图谱；配准模块将所述MRT脑图谱变形到插值图像空间以得到初始配准图像；提取模块提取初始配准图像的中矢状面，其中，配准模块还对中矢状面进行配准，得到中矢配准结果，根据中矢配准结果以及灰度信息对初始配准图像进行再次配准以获得再次配准结果，插值模块还对再次配准结果进行插值以获得配准图像。本发明提高了图像配准的精准度。
【专利说明】图像配准装置及其图像配准的方法
【技术领域】
[0001]本发明涉及医学图像处理技术，尤其涉及一种图像配准装置及其图像配准的方法。
【背景技术】
[0002]医学影像分析作为一种计算机辅助诊断方式已经在临床上得到了广泛的应用。对于那些病情复杂，容易误诊的重大疾病，医学图像分析方法与系统可以为诊断提供更为客观的参考信息和技术支持，一定程度的避免了由于主观因素造成误诊的情况。
[0003]X射线计算机断层(CT)图像和磁共振图像(MRI)作为当今临床上两种重要的影像手段具有很强的互补性。CT图像的优点是空间分辨率高且成像速度快，能清晰地反映出骨骼特征，但是CT图像不能有效地区分软组织信息。MRI的优势在于能提供清晰的软组织信息，特别适合于软组织信息丰富的人脑的成像。CT成像设备相对较低廉，有更多的医院拥有该设备；MRI设备价格昂贵，只有较少的医院拥有该设备；一些疾病的诊断(如脑中风)的首选影像是CT而不是MRI。因此，如何利用现有的技术来解读软组织信息不太丰富的头颅CT图像具有重大的意义。
[0004]CT (X射线计算机断层成像)和MRI (磁共振成像)是临床上用于诊断疾病的两种重要的手段。由于成像方式的差异，这两种方式反映了不同的病理信息。临床上经常需要综合CT图像的MR(磁共振)图像的信息来进行诊断，特别是对于人脑这样结构复杂的器官时，经常需要将CT图像和MRI图像进行配准融合来获取更多的信息。
[0005]图像配准技术是图像处理领域的一项十分重要而且非常困难的技术，其中多模态的医学图像配准算法更是研究的热点。目前的医学图像配准方法大体上可以分为两类:一类是基于特征的配准方法，这里的特征可以是边缘和表面特征，也可以是一些特征点。迭代最近点(Iterative Closest Point)算法是基于特征点配准的经典的算法。基于特征的配准方法利用图像的结构信息，速度较快，但是忽略了图像的灰度信息。另一类配准方法是基于图像灰度信息进行配准，但是没有考虑到图像的结构信息。这两种方法各有优劣，而且具有很大的互补性，因此通过研究混合算法使得在配准过程中兼顾这两类信息从而获得高精度的配准具有重要意义，也是当前研究的热点。
[0006]2001Hsu等提出了一种基于特征的多模态配准方法。在算法中应用了图像的表面和边缘作为特征，然后应用迭代最近点算法将提取出的特征进行配准。
[0007]Vercauteren等在2009年提出了一种基于灰度信息的配准算法。这种算法是在原始的demons算法的基础上发展而来的。它使demons算法的优化过程在微分同胚空间中进行，解决了大形变的问题。
[0008]Hellier等在2003年提出了一种基于特征点约束和灰度信息的配准算法，该算法是在分割的基础上进行的。首先提取出脑沟，分割出来的脑沟上的点作为特征点，对应关系可以不严格，因此不是作为硬约束作为下一步配准算法的约束。
[0009]Mattes等在2003年提出了一种基于B样条的自由形变模型的多模态配准方法。在这个算法中以互信息为相似度准则，以近似牛顿法作为优化方法并且在优化过程中采用了多分辨率策略。
[0010]Andronache等在2008年提出了一种互信息和相关比结合多模态配准的方法。此算法中通过不断地判断在分层过程中互信息数值的连续性来实现从互信息到相关比的转换，在多模态如脑部CT和MRI的配准中得到了很好的效果。
[0011]然而，Hsu等提出的方法只考虑到了图像的结构信息而忽略了图像的灰度信息，在远离特征点的区域配准的结果会很差。
[0012]Vercauteren提出的基于微分同胚的demons算法没有摆脱demons算法的固有缺陷，该算法只适用于单模态配准中。
[0013]Hellier等提出的算法利用的形变模型是一种基于光流场的模型，这个模型属于非参数模型。该算法同demons算法一样适用于单一模态的配准中。而且把脑沟这个解剖结构作为约束也是有局限性的，因为脑沟在一些成像方式如CT图像上是很难分割出来的。
[0014]Mattes等提出的算法是一种基于灰度信息的方法，该算法忽略了图像的结构信息。算法中的形变模型使得形变是自由的，这就可能导致某些刚性组织如脑的中矢状面会发生非刚性形变。
[0015]Andronache等提出的算法也没有考虑图像的结构信息，这样可能会导致一些图像上的一些关键点(如一些有解剖意义的点)配准不好，从而影响配准精度。
【发明内容】
[0016]有鉴于此，有必要提供一种图像配准装置以及及其图像配准的方法。
[0017]本发明提供的图像配准装置，包括:获取模块、分割模块、插值模块、配准模块以及提取模块，其中获取模块用于获取人脑计算机断层(CT)图像；分割模块用于对所述人脑CT图像进行脑分割以得到分割图像；插值模块用于对所述分割图像进行插值获得插值图像，其中所述获取模块还用于获取磁共振(MRI)脑图谱；配准模块用于将所述MRT脑图谱变形到所述插值图像空间以得到初始配准图像；提取模块用于提取所述初始配准图像的中矢状面，其中，所述配准模块还用于对所述中矢状面进行配准，得到中矢配准结果，根据所述中矢配准结果以及灰度信息对所述初始配准图像进行再次配准以获得再次配准结果，所述插值模块还用于对所述再次配准结果进行插值以获得配准图像。
[0018]本发明提供的图像配准的方法，包括:获取人脑计算机断层(CT)图像；对所述人脑CT图像进行脑分割以得到分割图像；对所述分割图像进行插值获得插值图像；获取磁共振(MRI)脑图谱，并将所述MRT脑图谱变形到所述插值图像空间以得到初始配准图像；提取所述初始配准图像的中矢状面；对所述中矢状面进行配准，得到中矢配准结果；根据所述中矢配准结果以及灰度信息对所述初始配准图像进行再次配准以获得再次配准结果；以及对所述再次配准结果进行插值以获得配准图像。
[0019]本发明所提供的图像配准装置10及其图像配准的方法改善了中矢状面和表面的配准结果，提高了图像配准的精准度。
【专利附图】
【附图说明】
[0020]图1为本发明一实施方式中图像配准装置的模块图；[0021]图2为本发明一实施方式中基于仿射变换的初始配准框架；
[0022]图3为本发明一实施方式中图像配准的方法的流程图。
【具体实施方式】
[0023]下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
[0024]在本发明的描述中，术语“内”、“外”、“纵向”、“横向”、“上”、“下”、“顶”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明而不是要求本发明必须以特定的方位构造和操作，因此不能理解为对本发明的限制。
[0025]请参阅图1，图1所示为本发明一实施方式中图像配准装置10的模块图。
[0026]在本实施方式中，图像配准装置10包括:获取模块110、分割模块120、插值模块130、配准模块140以及提取模块150。
[0027]在本实施方式中，获取模块110用于获取人脑计算机断层(CT)图像。在本实施方式中获取模块110为CT。
[0028]在本实施方式中，分割模块120用于对所述人脑CT图像进行脑分割以得到分割图像。在本实施方式中，所述脑分割是指去除颅骨，得到脑(包括大脑和小脑)内的灰质、白质与脑脊液。
[0029]在本实施方式中，所述分割模块120利用模糊c均值算法求取灰度阈值以及利用所述灰度阈值对所述CT图像进行分割。
[0030]在本实施方式中，插值模块130用于对所述分割图像进行插值获得插值图像。
[0031]在本实施方式中，CT图像采用轴向位置扫描，原始的CT图像的尺寸是512x512x切片数，切片数介于10-20之间，对应的X与Y轴方向的像素大小是相同的(约0.5毫米)，Z轴方向的像素大小为5-10毫米。为便于后续的空间变换(如旋转、非线性形变等)，在预处理阶段需要对CT图像进行插值操作。可选的插值方法有很多，最简单的是最邻近内插法:插值后的图像在非整数位置的灰度值是等于与其位置最接近像素的灰度值。为了方便，插值图像的像素尺寸和图像大小与MRI图谱一致，图像大小为149x188x148 (重采样去掉了一些不包含脑组织的像素使得所有脑组织的像素都被包含在插值后的图像中)。
[0032]在本实施方式中，所述获取模块110还用于获取磁共振(MRI)脑图谱。在本实施方式中，MRT脑图谱预先存储在图像配准装置10中，是LONI Atlases中的ICBM452TlAtlas。
[0033]在本实施方式中，配准模块140用于将所述MRT脑图谱变形到所述插值图像空间以得到初始配准图像。
[0034]在本实施方式中，所述初始配准通过仿射变换模型来实现。
[0035]在本实施方式中，仿射变换模型是一种全局形变模型，由12个参数控制。其中，像素P点(X，y，Ζ)τ处经过仿射变换后的坐标变为(x’，y’，z’)τ,变换公式为:
【权利要求】
1.一种图像配准装置，包括: 获取模块，用于获取人脑计算机断层(CT)图像； 分割模块，用于对所述人脑CT图像进行脑分割以得到分割图像； 插值模块，用于对所述分割图像进行插值获得插值图像，其中所述获取模块还用于获取磁共振(MRI)脑图谱； 配准模块，用于将所述MRT脑图谱变形到所述插值图像空间以得到初始配准图像； 提取模块，用于提取所述初始配准图像的中矢状面，其中，所述配准模块还用于对所述中矢状面进行配准，得到中矢配准结果，根据所述中矢配准结果以及灰度信息对所述初始配准图像进行再次配准以获得再次配准结果，所述插值模块还用于对所述再次配准结果进行插值以获得配准图像。
2.如权利要求1所述的图像配准装置，其特征在于，所述脑分割是指去除颅骨，得到脑(包括大脑和小脑)内的灰质、白质与脑脊液。
3.如权利要求2所述的图像配准装置，其特征在于，所述分割模块利用模糊c均值算法求取灰度阈值以及利用所述灰度阈值对所述CT图像进行分割。
4.如权利要求1所述的图像配准装置，其特征在于，所述初始配准通过仿射变换模型来实现。
5.如权利要求1所述的图像配准装置，其特征在于，所述中矢状面是将脑分成左右半脑的曲面的平面逼近，由脑纵裂组成，基于灰度局部对称性在轴向切片上定位所述脑纵裂。
6.如权利要求5所述的图像配准装置，其特征在于，所述配准模块通过将所述MRI图谱配准到所述CT图像的形变场，并对所述形变场进行采样，得到所述CT图像与所述初始配准图像的对应关系，构成中矢状面约束，其中，所述中矢状面约束为所述中矢配准结果。
7.一种图像配准的方法，包括: 获取人脑计算机断层(CT)图像； 对所述人脑CT图像进行脑分割以得到分割图像； 对所述分割图像进行插值获得插值图像； 获取磁共振(MRI)脑图谱，并将所述MRT脑图谱变形到所述插值图像空间以得到初始配准图像； 提取所述初始配准图像的中矢状面； 对所述中矢状面进行配准，得到中矢配准结果； 根据所述中矢配准结果以及灰度信息对所述初始配准图像进行再次配准以获得再次配准结果；以及 对所述再次配准结果进行插值以获得配准图像。
8.如权利要求7所述的图像配准的方法，其特征在于，所述脑分割是指去除颅骨，得到脑(包括大脑和小脑)内的灰质、白质与脑脊液。
9.如权利要求8所述的图像配准的方法，其特征在于，所述脑分割包括以下步骤: 利用模糊c均值算法求取灰度阈值； 利用所述灰度阈值对所述CT图像进行分割。
10.如权利要求7所述的图像配准的方法，其特征在于，所述初始配准通过仿射变换模型来实现。
11.如权利要求7所述的图像配准的方法，其特征在于，所述中矢状面是将脑分成左右半脑的曲面的平面逼近，由脑纵裂组成，基于灰度局部对称性在轴向切片上定位所述脑纵裂。
12.如权利要求11所述的图像配准的方法，其特征在于，“对所述中矢状面进行配准”的步骤包括: 将所述MRI图谱配准到所述CT图像的形变场； 对所述形变场进行采样，得到所述CT图像与所述初始配准图像的对应关系，构成中矢状面约束，其中，所述中矢状面约束为所述中矢配准结果。
【文档编号】G06T7/00GK103793904SQ201210424003
【公开日】2014年5月14日 申请日期:2012年10月29日 优先权日:2012年10月29日
【发明者】胡庆茂, 高爱丽, 陈明扬, 张晓东, 罗火灵, 贾富仓 申请人:深圳先进技术研究院