用于三维数字减影血管造影图像的快速配准和拼接方法

摘要

本发明公开了一种用于三维数字减影血管造影图像的快速配准和拼接方法，包括如下步骤：a)选取三维数字减影血管造影图像的两个体数据图像,提取重叠区域记为子体数据图像；b)将子体数据图像分别在冠状面、矢状面和横断面进行最大密度投影形成MIP1和MIP2图像；c)分别对MIP1和MIP2图像在上述二维平面上进行第一次配准；d)根据第一次配准结果对所述两个体数据图像在三维空间进行第二次配准，然后将重叠区域进行融合，最后拼接形成全景的三维减影体数据图像。本发明提供的用于三维数字减影血管造影图像的快速配准和拼接方法，能够避免图像数据质量对配准精度的影响，并且大大缩短配准和拼接处理时间。

说明书

技术领域

本发明涉及一种医学图像处理方法，尤其涉及一种用于三维数字减影血管造影图像的快速配准和拼接方法。

背景技术

血管造影检查目前是诊断血管性疾病的主要方法。在医学成像过程中，往往存在被检查部位远大于图像检测器面积的矛盾，需要将该部位分几次拍摄，再按照一定的规则拼接起来。以磁共振成像为例，由于MRI系统扫描板尺寸的限制，在进行下肢血管或全身血管扫描时，MRI会进行多段扫描、重建与减影，为了将多个分段的血管减影体数据合成一个全景血管减影体数据，拼接已经成为必要的图像后处理应用。如专利文献1：Ian Clarke,Graham Sellers,Zain Adam Yusuf.System and method for creating apanorama images from a plurality of source images[P].Epson research and development inc.US2005/0063608A1.Mar.24,2005给出了医学图像拼接方法的一般流程图。

对医学图像拼接方法的研究一部分是关注于扫描系统的设计，如专利文献2：Vernon Thomas Jensen,Richard Leparmentier,Vianney Pierre Battle.System and method for imagecomposition using position sensors[P].General Electric Company.US 7522701B2.Apr.21,2009，提出使用一种位置传感器，各子图的方位信息都可以通过位置传感器获取，在计算机后处理中凭借各子图的方位信息进行拼接，专利文献3：Shigehide Kuhara.Magnetic resonance imaginginvolving movement of patient’s couch[P].Kabushiki Kaisha Toshiba.US 7423428B2.Sep.9,2008则是通过记录可移动床板的位置，实现多幅图像的拼接。

除此之外，对拼接方法的研究另一部分是集中在对拼接中使用的图像配准算法的改进。专利文献4：Benjamin Glocker,Nassir Navab,Christian Wachinger,Jochen Zeltner.Method forcombining images and magnetic resonance scanner[P].US 2010/0067762A 1.Mar.18,2010提出了一种配准策略，通过迭代的方式不断缩小上下两幅图与加权叠加的重叠区域之间的差异，当迭代收敛时，更新的加权叠加的重叠区域被视为拼接而成的重叠区域；专利文献5：Ram Nathaniel,Dan Rappaport,Ishay Goldin.Method and system for stitching multiple images into apanoramic image[P].US 2011/0188726A1.Aug.4,2011在两张图像的重叠区域上设置标志点，通过配准标志点来实现两张图像的位置配准。专利文献6：Li Zhang,Michal Sofka,Ulf Schafer.Feature-based composing for 3D MR angiography images[P].Siemens corporation intellectual propertydepartment.US 2006/0052686A1.Mar.9,2006针对三维MR血管造影法图像提出基于特征的配准方法，并给出实例，三个352*512*96的三维血管图像，最大可能的重叠区域是351*67*73，拼接耗时13秒。

由上可见，在各子图扫描期间病人床板的位移、每次扫描时技师对扫描范围的更改以及病人自身的运动等因素，都会导致扫描得到的各子图在病人坐标系中没有对准，所以必须经过坐标配准后，才能将各子图的重叠区域融合，进而拼成更大的全景图。对于多段的三维血管造影减影体数据图像而言，由于血管减影体数据中各片层的前景图像（即打了造影剂的高亮血管）信息量很少，更多的是噪声污染的黑色背景，而黑色背景中常包含未被消减彻底的骨头和软组织的残留信息，使得图像信噪比很低，这样的图像数据质量对三维图像配准算法的精度提出了极高的要求，并且往往血管减影体数据的片层数目较多，故三维图像的配准是比较耗时的，这样最终导致拼接的处理时间较长。此外，由于重叠区域都在图像的上下边缘，磁场梯度不均匀也常常会导致重叠区域的信息差异较大，这些因素会影响三维图像的配准算法精度。

发明内容

本发明所要解决的技术问题是提供一种用于三维数字减影血管造影图像的快速配准和拼接方法，能够避免图像数据质量对配准精度的影响，并且大大缩短配准和拼接处理时间。

本发明为解决上述技术问题而采用的技术方案是提供一种用于三维数字减影血管造影图像的快速配准和拼接方法，包括如下步骤：a)选取三维数字减影血管造影图像的两个体数据图像,提取所述两个体数据图像的重叠区域，记为子体数据图像I1和I2；b)将子体数据图像I1和I2分别在冠状面、矢状面和横断面进行最大密度投影，分别形成MIP1和MIP2图像；c)分别在冠状面、矢状面和横断面上对MIP1和MIP2图像在二维平面上进行第一次配准；记录冠状面、矢状面和横断面上的第一次配准结果的变换参数分别为(T_{COR_X}，T_{COR_Z})、(T_{SAG_Y}，T_{SAG_Z})和(T_{AXI_X}，T_{AXI_Y})；d)根据所述第一次配准结果对所述两个体数据图像在三维空间进行第二次配准。

上述的用于三维数字减影血管造影图像的快速配准方法，其中，所述两个体数据图像为相邻扫描的图像。

上述的用于三维数字减影血管造影图像的快速配准方法，其中，所述第二次配准得到的两个体数据图像在三维空间的相对偏移量(t_x，t_y，t_z)计算如下：

$>t_x=\frac{T_{\mathrm{COR}\_X}+T_{\mathrm{AXI}\_X}}{2}$>

$>t_y=\frac{T_{\mathrm{SAG}\_Y}+T_{\mathrm{AXI}\_Y}}{2}$>

$>t_z=\frac{T_{\mathrm{COR}\_Z}+T_{\mathrm{SAG}\_Z}}{2}$>

上述的用于三维数字减影血管造影图像的快速配准方法，其中，所述三维血管减影体数据图像为MR血管造影法或CT血管造影法生成的图像。

本发明为解决上述技术问题还提供一种用于三维数字减影血管造影图像的快速拼接方法，包括如下步骤：a）将上述配准方法获得的图像与所述两个体数据图像的重叠区进行融合；b)将所述两个体数据图像进行拼接，形成全景的三维数字减影血管造影图像。

本发明对比现有技术有如下的有益效果：本发明提供的用于三维数字减影血管造影图像的快速配准和拼接方法，利用最大密度投影将数字减影血管造影体数据的片层间原本分散的信息集中起来且最大化，为基于像素灰度统计的配准算法提供了更多的有效信息，并将三维图像的配准简化为二维图像的配准，从而避免了图像数据质量对配准精度的影响，并且大大缩短了处理时间。

附图说明

图1为本发明用于三维数字减影血管造影图像的快速配准流程示意图；

图2为本发明实施例拼接后的血管减影体数据在各角度的MIP图像。

具体实施方式

下面结合附图和实施例对本发明作进一步的描述。

图1为本发明用于三维数字减影血管造影图像的快速配准流程示意图。

请参见图1，本发明提供的用于三维数字减影血管造影图像的快速配准方法包括如下步骤：

S101：选取两个相邻血管的三维血管减影体数据图像沿重叠区域裁剪分割，提取所述两个体数据图像的重叠区域，获得两个子体数据图像，记为子体数据图像I1和I2；所述三维血管减影体数据图像为MR血管造影法、CT血管造影法或其它成像方式生成的图像。以MR三维血管造影法减影体数据图像的处理为例，MR血管造影法采集得到的三维减影体数据图像位于冠状面，则首先根据DICOM（Digital Imagingand Communication of Medicine，即数字影像和通信标准）头信息的（0020 0032）图像位置信息将相邻两个减影体数据的重叠区域裁剪出来制作成两个子体数据图像。

S102：分别在冠状面、矢状面和横断面上将这两个子体数据图像I1和I2进行最大密度投影（maximum intensity projection，MIP），各平面都得到两张MIP图像，分别形成MIP1和MIP2图像；最大密度投影有时又称为“最大亮度投影”，是在可视化平面之上投射三维空间数据的一种计算机可视化方法；其中，沿着从视点到投影平面的平行光线，各个体素密度值的所呈现的亮度将以某种方式加以衰减，并且最终在投影平面上呈现的是亮度最大的体素。

S103：接着分别在上述三个平面上，对得到的MIP图像进行第一次图像配准，记录冠状面、矢状面和横断面上的第一次配准结果的变换参数分别为(T_{COR_X}，T_{COR_Z})、(T_{SAG_Y}，T_{SAG_Z})和(T_{AXI_X}，T_{AXI_Y})；MIP图像上血管信息丰富，配准算法精度可以得到保障。

S104：根据三个平面的二维图像配准结果对所述两个体数据图像在三维空间进行第二次配准，计算出两个血管减影体数据图像在三维空间的相对偏移量后进行坐标配准，对三组配准结果中，基于像素灰度进行综合统计，方法如下：冠状面上的配准结果的变换参数：(T_{COR_X}，T_{COR_Z})；矢状面上的配准结果的变换参数：(T_{SAG_Y}，T_{SAG_Z})；横断面上的配准结果的变换参数：(T_{AXI_X}，T_{AXI_Y})；则第二个血管减影体数据相对于第一个血管减影体数据的偏移量(t_x，t_y，t_z)如下：

$>t_x=\frac{T_{\mathrm{COR}\_X}+T_{\mathrm{AXI}\_X}}{2}$>

$>t_y=\frac{T_{\mathrm{SAG}\_Y}+T_{\mathrm{AXI}\_Y}}{2}$>

$>t_z=\frac{T_{\mathrm{COR}\_Z}+T_{\mathrm{SAG}\_Z}}{2}$>

本发明对三个平面上的配准结果进行综合分析和处理，通过两次配准得到各个减影体数据在病人坐标系中相对偏移量(t_x，t_y，t_z)，等效于实现了使用三维图像配准的目的。凭借最大密度投影，将血管造影减影体数据的片层间原本分散的信息集中起来且最大化，为基于像素灰度统计的配准算法提供了更多的有效信息；同时传统的三维图像配准被简化为二维图像配准，极大地减少了处理时间。根据上述结果对各个减影体数据重新定位，利用(t_x，t_y，t_z)将相邻两个三维血管减影体数据图像进行坐标配准，然后将重叠区域进行融合，最后拼接形成全景的三维减影体数据图像。

例如对三个冠状面血管减影体数据进行拼接，每个减影体数据尺寸是384*512*88，根据DICOM头信息可知相邻两个血管减影体数据的最大重叠区域尺寸是384*72*88，在应用中对这三个减影体数据进行拼接，计算机系统计时处理时间约40s；使用本专利所提方法所需时间是2.075s，运行环境如下，Intel i5-2400处理器，主频是3.10GHz，4GB内存，64位操作系统。

医学影像工作站中，通常将血管减影体数据转为3D MIP显示，目的是可以从各角度观察血管走向。现将拼接而成的血管减影体数据沿轴向（头-脚方向）逆时针旋转，并在冠状面上进行最大密度投影，得到多角度的MIP，图2中展示拼接后的血管减影体数据在各角度的MIP图像，其中，图2中(a)、(b)、(c)、(d)、(e)、(f)、(g)、(h)和(i)分别对应0°、45°、90°、135°、180°、225°、270°、315°和360°的MIP图像，从血管在缝合处的连贯性可以验证本发明配准和拼接结果的准确性。

虽然本发明已以较佳实施例揭示如上，然其并非用以限定本发明，任何本领域技术人员，在不脱离本发明的精神和范围内，当可作些许的修改和完善，因此本发明的保护范围当以权利要求书所界定的为准。