一种基于目标物像素投影判断的断层切片图像配准方法
【专利摘要】本发明公开了一种基于目标物像素投影判断的断层切片图像配准方法，首先将参考图像和浮动图像进行二值化处理，将其目标物和背景区分开来，通过将浮动图像按照指定的微小的角度增量旋转，每旋转一次计算浮动图像在每一像素列上目标像素的数目，并计算求其与参考图像每一像素列上目标像素数目之差，将所有差值的绝对值求和，当其值达到最小时判断此时浮动图像的角度配准达到最优，其次对其平移配准，在平面的X和Y方向内，浮动图像在给定的范围内以单个像素为增量进行平动，每平动一次计算浮动图像和参考图像在相同坐标点处都是目标物像素的数目，当其值达到最大时候判断浮动图像已经与参考图像达到比较完美配准。
【专利说明】一种基于目标物像素投影判断的断层切片图像配准方法
【技术领域】
[0001] 本发明公开了一种基于目标物像素投影判断的断层切片图像配准方法，涉及图像 处理【技术领域】。
【背景技术】
[0002] 图像配准是图像处理领域内的一个非常常见的研究内容，从20世纪90年代初期 此项研究才逐渐得到重视。其作为一种基础的图像处理方面的技术，其在诸多领域有着广 泛的应用，从图像配准的提出到目前为止的几十年时间里，随着计算机科学技术的飞速发 展和普及，图像配准被用于许多应用领域，诸如医学应用、计算机辅助制造、机器人等。
[0003] 图像配准在医学应用领域尤其有用。例如，在诸如病人随访和手术或疗法计划等 临床场景中。为了使医学专家能够以精确的方式比较医学图像，这些图像需要被配准以便 去除病人的姿势变化。例如，所述配准可能导致两个医学图像中的相应解剖结构的对齐。
[0004] 一直存在提高图像配准系统的准确度和稳健性以及减少图像配准所花的时间和 所需的计算资源的需要。而以往图像进行匹配时，存在计算复杂，计算量大，所耗时间长，效 率上提1?不上去等弊端。
【发明内容】
[0005] 本发明所要解决的技术问题是：针对现有技术的缺陷，提供一种基于目标物像素 投影判断的断层切片图像配准方法，在保证配准精度的基础上，减少配准计算时间，从而大 大提1?其配准效率。
[0006] 本发明为解决上述技术问题采用以下技术方案：
[0007] -种基于目标物像素投影判断的断层切片图像配准方法，具体步骤包括：
[0008] 步骤一、对参考图像和浮动图像进行二值化操作，将目标物和背景进行区分；
[0009] 步骤二、以参考图像为基准，对浮动图像进行角度配准操作，得到二值化浮动图像 的旋转角度；
[0010] 步骤三、以参考图像为基准，对浮动图像进行平移配准操作，得到二值化浮动图像 在平面内横向和纵向移动的距离；
[0011] 步骤四、将所得到的二值化浮动图像旋转的角度以及在平面内横向和纵向移动的 距离分别运用于实际原图像中，使得其按照指定的角度和距离进行旋转和平动，矫正原浮 动图像，使其与参考图像现实配准。
[0012] 作为本发明的进一步优选方案，所述步骤二具体包括：
[0013] (201):对经过二值化处理的浮动图像设置其对应的角度旋转范围，当所述浮动图 像在所述角度旋转范围内旋转时，浮动图像将有机会与参考图像的角度一直；
[0014] (202):在（201)设定的角度旋转范围内，对经过二值化处理的浮动图像设置其对 应的旋转角度增量，浮动图像以旋转角度增量不停的旋转，每旋转一次分别计算二值化浮 动图像和二值化参考图像上每一像素列上目标像素的数目；
[0015] (203):根据（202)得出的二值化参考图像和二值化浮动图像每一列的目标物像 素数，列出两个向量表达式：
[0016] R = [R(1),R(2)……R(P)],F = [F(1),F(2)……F(P)]；
[0017] 其中，R表示二值化参考图像中每一列的目标物像素数，F表示二值化浮动图像中 每一列的目标物像素数，P表示列数；
[0018] (204):将步骤（203)中，将两个向量表达式中前面的0元素全部去除，使得第一个 非零元素作为向量的第一个元素；
[0019] (205):在经过（204)处理的向量中截取其前N个元素，使得这N个元素均能将目 标物包含在内，最终得出如下向量表达式：
[0020] R，= [R，（1)，R，（2)，R，（3)......R，（N-2)，R，（N_l)，R，（N)]，
[0021] F，= [F，（1)，F，（2)，F，（3)......F，（N_2)，F，（N_l)，F，（N)];
[0022] 其中，R'和F'分别为R和F对应的变化后表达形式；
[0023] (206):每旋转一次浮动图像分别计算向量R'与F'各对应元素之差的绝对值之 和，即表达式 F'(z')l的值，最终得到关于Μ的一个向量：
[0024] Μ= [Μ(1),Μ(2),Μ(3), ......M(j-2), M(j-l), M(j)]；
[0025] 求出其中的最小值 M (k) = min [Μ (1)，Μ (2)，Μ (3)，……Μ (j_2)，Μ (j_l)，Μ (j)]，当 Μ 的取值为M(k)时，浮动图像矫正完毕。
[0026] 作为本发明的进一步优选方案，所述步骤三具体包括：
[0027] (301):给定二值化浮动图像在X轴和Y轴方向移动的范围，在此范围内浮动图像 将有机会与参考图像重置；
[0028] (302):在浮动图像的移动过程中，计算浮动图像和参考图像在相同的坐标点都是 目标物像素点的像素数S ;
[0029] (303):经过步骤（302)得出一个关于S值分布的向量：
[0030] [S ⑴，S ⑵，S (3)，……S (j-2)，S (j-1)，S (j)];
[0031] 求得其最大值：
[0032] S (m) = max [S (1)，S (2)，S (3)，...... S (j-2), S (j-1), S (j)]；
[0033] 在二值化浮动图像平动的过程中，当S的取值为S(m)时，二值化浮动图像与参考 图像实现配准。
[0034] 作为本发明的进一步优选方案，所述参考图像和浮动图像的尺寸大小一致。
[0035] 作为本发明的进一步优选方案，步骤（201)中，所述角度旋转范围为0°至360°。
[0036] 本发明采用以上技术方案与现有技术相比，具有以下技术效果：本发明所公开的 图像配准方法，显著提高了图像配准系统的准确度和稳健性，减少图像配准所花的时间和 所需的计算资源，计算量小、速度快，明显的提升了系统的效率。
【专利附图】
【附图说明】
[0037] 图1为二值化后的参考图像，其中白色为目标物，黑色为背景；
[0038] 图2为二值化后的浮动图像，其中白色为目标物，黑色为背景；
[0039] 图3为二值化浮动图像每一像素列上目标像素的数目坐标图；
[0040] 图4为二值化参考图像每一像素列上目标像素数目坐标图；
[0041] 图5.为角度得到矫正后的二值化浮动图像；
[0042] 图6为四幅简单二值化图像示意图；
[0043] 图7为最终得到配准的二值化浮动图像；
[0044] 图8为图像配准流程图。
【具体实施方式】
[0045] 下面结合附图对本发明的技术方案做进一步的详细说明：
[0046] 对浮动图像的配准操作过程中，其主要包括角度配准以及平移配准两方面。
[0047] 本发明所公开的方法中，图像配准流程如图8所示。所述基于目标物像素投影判 断的断层切片图像配准方法，具体步骤包括：
[0048] 步骤一、对参考图像和浮动图像进行二值化操作，将目标物和背景进行区分；
[0049] 步骤二、以参考图像为基准，对浮动图像进行角度配准操作，得到二值化浮动图像 的旋转角度；
[0050] 步骤三、以参考图像为基准，对浮动图像进行平移配准操作，得到二值化浮动图像 在平面内横向和纵向移动的距离；
[0051] 步骤四、将所得到的二值化浮动图像旋转的角度以及在平面内横向和纵向移动的 距离分别运用于实际原图像中，使得其按照指定的角度和距离进行旋转和平动，矫正原浮 动图像，使其与参考图像现实配准。
[0052] 出于方便操作的考虑，首先将参考图像和浮动图像均进行二值化的操作，将目标 物和背景进行区分。二值化后的参考图像和浮动图像分别如图1、图2所示。首先对浮动图 像进行角度配准的操作。
[0053] 给二值化浮动图像一个角度旋转的范围，使得二值化浮动图像的角度在这个旋转 范围内能有机会能与参考图像一致。在给定的角度范围内，二值化浮动图像以一个很小的 角度增量不停地旋转，每旋转一次分别计算二值化浮动图像和二值化参考图像上每一像素 列上目标像素的数目，分别得到如图3、图4所示的坐标图。
[0054] 根据坐标图可以明显看出，在刚开始一段，目标像素数目是0,这是由于目标物并 非是从最左边开始，由于参考图像和浮动图像的目标物的角度存在一定的差异，因此两坐 标图的形状存在一定的差异。
[0055] 根据二值化参考图像和二值化浮动图像每一列的目标物像素数，这里可以列出两 个向量R= [R(1)，R(2)……R(700)]以及F = [F(1)，F(2)……F(700)]，其向量的每个元 素的值表示对应力的列上的像素数目。两个向量前面的若干元素值都是〇,到后面的第一个 非零值时说明图像的该列已经有目标物的像素，到最后所有元素值又统一为〇,说明在图像 的最后数列内并没有目标物的像素，所以R和F应该具有的形式为：
[0056] R= [0,0,0......R(i-l), R(i), R(i+1)......0,0,0];
[0057] F = [0, 0, 0……F (i-1)，F (i)，F (i+1)……0,0,0]；
[0058] 将两个向量的前面的所有0元素全部去除，使得第一个非零元素作为向量的第一 个元素，得到两个新的向量之后同时截取前面的N个元素，使得这N个元素列都能将目标物 包含在内，最终得到的两个向量为：
[0059] R，= [R，（1)，R，（2)，R，⑶......R，（N_2)，R，（N_l)，R，（N)];
[0060] F，= [F，（1)，F，（2)，F，（3)......F，（N_2)，F，（N_l)，F，（N)];
[0061] 每旋转一次浮动图像后，分别计算
【权利要求】
1. 一种基于目标物像素投影判断的断层切片图像配准方法，其特征在于，具体步骤包 括： 步骤一、对参考图像和浮动图像进行二值化操作，将目标物和背景进行区分； 步骤二、以参考图像为基准，对浮动图像进行角度配准操作，得到二值化浮动图像的旋 转角度； 步骤三、以参考图像为基准，对浮动图像进行平移配准操作，得到二值化浮动图像在平 面内横向和纵向移动的距离； 步骤四、将所得到的二值化浮动图像旋转的角度以及在平面内横向和纵向移动的距离 分别运用于实际原图像中，使得其按照指定的角度和距离进行旋转和平动，矫正原浮动图 像，使其与参考图像现实配准。
2. 如权利要求1所述的一种基于目标物像素投影判断的断层切片图像配准方法，其特 征在于，所述步骤二具体包括： (201) :对经过二值化处理的浮动图像设置其对应的角度旋转范围，当所述浮动图像在 所述角度旋转范围内旋转时，浮动图像将有机会与参考图像的角度一致； (202) :在（201)设定的角度旋转范围内，对经过二值化处理的浮动图像设置其对应的 旋转角度增量，浮动图像以旋转角度增量不停的旋转，每旋转一次分别计算二值化浮动图 像和二值化参考图像上每一像素列上目标像素的数目； (203) :根据（202)得出的二值化参考图像和二值化浮动图像每一列的目标物像素数， 列出两个向量表达式： R= [R(1)，R(2)……R(P)],F= [F(I),F(2)……F(P)]； 其中，向量R表示二值化参考图像中每一列的目标物像素数，向量F表示二值化浮动图 像中每一列的目标物像素数，P表示列数； (204) :将步骤（203)中，将两个向量表达式中前面的0元素全部去除，使得第一个非零 元素作为向量的第一个元素； (205) :在经过（204)处理的向量中截取其前N个元素，使得这N个元素均能将目标物 包含在内，最终得出如下向量表达式： R，= [R，（1)，R，（2)，R，⑶......R，（N-2)，R，（N-I)，R，（N)]; F' = [F，（1)，F，（2)，F，（3)......F'（N-2)，F'（N-I)，F'（N)]; 其中，R'和F'分别为R和F对应的变化后表达形式； (206) :每旋转一次浮动图像分别计算向量R'与F'各对应元素之差的绝对值之和，即 表达式从=Zl/?'(/)_，'(/)I的值，最终得到关于M的一个向量： M= [M(1)，M(2)，M(3)，......M(j-2)，M(j-l)，M(j)]; 向量M的元素的个数即浮动图像所旋转的次数。找出向量M所有元素中值最小的元素M(k) =min[M⑴，M⑵，M(3)，……M(j_2)，M(j_l)，M(j)]，当M的取值为M(k)时，浮动图像 矫正完毕。
3. 如权利要求1或2所述的一种基于目标物像素投影判断的断层切片图像配准方法， 其特征在于，所述步骤三具体包括： (301):给定二值化浮动图像在X轴和Y轴方向移动的范围，在此范围内浮动图像将有 机会与参考图像重叠； (302) :在浮动图像的移动过程中，计算浮动图像和参考图像在相同的坐标点都是目标 物像素点的像素数S; (303) :经过步骤（302)得出一个关于S值分布的向量： [S(I)，S(2)，S(3)，……S(j-2)，S(j-1)，S(j)]; 求得其最大值： S(m) =max[S(I)，S(2)，S(3)，……S(j-2),S(j-1),S(j)]； 在二值化浮动图像平动的过程中，当S的取值为S(m)时，二值化浮动图像与参考图像 实现配准。
4. 如权利要求1所述的一种基于目标物像素投影判断的断层切片图像配准方法，其特 征在于：所述参考图像和浮动图像的尺寸大小一致。
5. 如权利要求2所述的一种基于目标物像素投影判断的断层切片图像配准方法，其特 征在于：步骤（201)中，所述角度旋转范围为0°至360°。
【文档编号】G06T7/00GK104268894SQ201410555187
【公开日】2015年1月7日 申请日期:2014年10月17日 优先权日:2014年10月17日
【发明者】袁铁军, 孟龙晖, 张秀丽 申请人:盐城工学院