用于三维数字减影血管造影图像的快速配准和拼接方法
【专利摘要】本发明公开了一种用于三维数字减影血管造影图像的快速配准和拼接方法，包括如下步骤：a)选取三维数字减影血管造影图像的两个体数据图像,提取重叠区域记为子体数据图像；b)将子体数据图像分别在冠状面、矢状面和横断面进行最大密度投影形成MIP1和MIP2图像；c)分别对MIP1和MIP2图像在上述二维平面上进行第一次配准；d)根据第一次配准结果对所述两个体数据图像在三维空间进行第二次配准，然后将重叠区域进行融合，最后拼接形成全景的三维减影体数据图像。本发明提供的用于三维数字减影血管造影图像的快速配准和拼接方法，能够避免图像数据质量对配准精度的影响，并且大大缩短配准和拼接处理时间。
【专利说明】用于三维数字减影血管造影图像的快速配准和拼接方法
【技术领域】
[0001]本发明涉及一种医学图像处理方法，尤其涉及一种用于三维数字减影血管造影图像的快速配准和拼接方法。
【背景技术】
[0002]血管造影检查目前是诊断血管性疾病的主要方法。在医学成像过程中，往往存在被检查部位远大于图像检测器面积的矛盾，需要将该部位分几次拍摄，再按照一定的规则拼接起来。以磁共振成像为例，由于MRI系统扫描板尺寸的限制，在进行下肢血管或全身血管扫描时，MRI会进行多段扫描、重建与减影，为了将多个分段的血管减影体数据合成一个全景血管减影体数据，拼接已经成为必要的图像后处理应用。如专利文献1:1anClarke, Graham Sellers, Zain Adam Yusuf.System and method for creating apanoramaimages from a plurality of source images[P].Epson research and development inc.US2005/0063608A1.Mar.24，2005给出了医学图像拼接方法的一般流程图。
[0003]对医学图像拼接方法的研究一部分是关注于扫描系统的设计，如专利文献2:Vernon Thomas Jensen，Richard Leparmentierj Vianney Pierre Battle.Systemand method for imagecomposition using position sensors [P].General ElectricCompany.US 7522701B2.Apr.21，2009，提出使用一种位置传感器，各子图的方位信息都可以通过位置传感器获取，在计算机后处理中凭借各子图的方位信息进行拼接，专利文献3:Shigehide Kuhara.Magnetic resonance imaginginvolving movement of patient’ scouch[P], Kabushiki Kaisha Toshiba.US 7423428B2.S印.9，2008 则是通过记录可移动床板的位置，实现多幅图像的拼接。
[0004]除此之外，对拼接方法的研究另一部分是集中在对拼接中使用的图像配准算法的改进 ° 专利文献 4:Benjamin Glocker, Nassir Navab, Christian ffachinger, JochenZeltner.Method forcombining images and magnetic resonance scanner[P].US2010/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 corporationintellectual propertydepartment.US 2006/0052686A1.Mar.9, 2006 针对三维 MR 血管造影法图像提出基于特征的配准方法，并给出实例，三个352*512*96的三维血管图像，最大可能的重叠区域是351*67*73，拼接耗时13秒。
[0005]由上可见，在各子图扫描期间病人床板的位移、每次扫描时技师对扫描范围的更改以及病人自身的运动等因素，都会导致扫描得到的各子图在病人坐标系中没有对准，所以必须经过坐标配准后，才能将各子图的重叠区域融合，进而拼成更大的全景图。对于多段的三维血管造影减影体数据图像而言，由于血管减影体数据中各片层的前景图像(即打了造影剂的高亮血管)信息量很少，更多的是噪声污染的黑色背景，而黑色背景中常包含未被消减彻底的骨头和软组织的残留信息，使得图像信噪比很低，这样的图像数据质量对三维图像配准算法的精度提出了极高的要求，并且往往血管减影体数据的片层数目较多，故三维图像的配准是比较耗时的，这样最终导致拼接的处理时间较长。此外，由于重叠区域都在图像的上下边缘，磁场梯度不均匀也常常会导致重叠区域的信息差异较大，这些因素会影响三维图像的配准算法精度。
【发明内容】
[0006]本发明所要解决的技术问题是提供一种用于三维数字减影血管造影图像的快速配准和拼接方法，能够避免图像数据质量对配准精度的影响，并且大大缩短配准和拼接处理时间。
[0007]本发明为解决上述技术问题而采用的技术方案是提供一种用于三维数字减影血管造影图像的快速配准和拼接方法，包括如下步骤:a)选取三维数字减影血管造影图像的两个体数据图像，提取所述两个体数据图像的重叠区域，记为子体数据图像Il和12 ；b)将子体数据图像Il和12分别在冠状面、矢状面和横断面进行最大密度投影，分别形成MIPl和MIP2图像；c)分别在冠状面、矢状面和横断面上对MIPl和MIP2图像在二维平面上进行第一次配准；记录冠状面、矢状面和横断面上的第一次配准结果的变换参数分别为Oaffi x，TCQR—Z)、(Tsag y, Tsag z)和(TAXI—x，Taxi y) ；d)根据所述第一次配准结果对所述两个体数据图像在三维空间进行第二次配准。
[0008]上述的用于三维数字减影血管造影图像的快速配准方法，其中，所述两个体数据图像为相邻扫描的图像。
[0009]上述的用于三维数字减影血管造影图像的快速配准方法，其中，所述第二次配准得到的两个体数据图像在三维空间的相对偏移量(tx，ty，tz)计算如下:
[0010]
【权利要求】
1.一种用于三维数字减影血管造影图像的快速配准方法，其特征在于，包括如下步骤: a)选取三维数字减影血管造影图像的两个体数据图像，提取所述两个体数据图像的重叠区域，记为子体数据图像Il和12 ; b)将子体数据图像Il和12分别在冠状面、矢状面和横断面进行最大密度投影，分别形成MIPl和MIP2图像；c)分别在冠状面、矢状面和横断面上对MIPl和MIP2图像在二维平面上进行第一次配准；记录冠状面、矢状面和横断面上的第一次配准结果的变换参数分别为(Taffi x, Tcoe z),(TsAG—Y，Tsag z)和(TAXI—X，Taxi y)； d)根据所述第一次配准结果对所述两个体数据图像在三维空间进行第二次配准。
2.如权利要求1所述的用于三维数字减影血管造影图像的快速配准方法，其特征在于，所述两个体数据图像为相邻扫描的图像。
3.如权利要求1所述的用于三维数字减影血管造影图像的快速配准方法，其特征在于，所述第二次配准得到的两个体数据图像在三维空间的相对偏移量(tx，ty，tz)计算如下:
4.如权利要求1所述的用于三维数字减影血管造影图像的快速配准方法，其特征在于，所述三维数字减影血管造影图像为MR血管造影法或CT血管造影法生成的图像。
5.一种用于三维数字减影血管造影图像的拼接方法，其特征在于，包括如下步骤: a)将如权利要求1~4任一项所述的配准方法获得的图像与所述两个体数据图像的重叠区进行融合； b)将所述两个体数据图像进行拼接，形成全景的三维数字减影血管造影图像。
【文档编号】G06T3/40GK103871036SQ201210533584
【公开日】2014年6月18日 申请日期:2012年12月12日 优先权日:2012年12月12日
【发明者】于文君, 其他发明人请求不公开姓名 申请人:上海联影医疗科技有限公司