[发明专利]补偿机械运动误差的图像配准方法

摘要

本发明公开了一种补偿机械运动误差的图像配准方法，包括步骤：对旋转减影得到的蒙片图像和造影图像进行高斯滤波；对高斯滤波后的造影图像进行分块；使造影图像块在蒙片图像中平移一定距离，并计算对应图像块之间的相关性系数，多次平移不同的距离得到一系列相关性系数，选择最小的相关性系数所对应的平移距离，使造影图像块在蒙片图像中平移该距离后进行减影，从而得到减影图像。本发明有效地消除了减影图像中的伪影。

主权项

1.一种补偿机械运动误差的图像配准方法，其特征在于，包括以下步骤：S1：分别对旋转减影得到的蒙片图像和造影图像进行高斯滤波，高斯函数如下：G(uj)=12πσexp(-uj22σ2);]]>其中u_j是坐标值，j是整数，j∈[0，N)，N是高斯滤波系数的个数，是大于3的奇数，σ是标准差，其中，u_j＝[j/(N-1)×2.0-1.0]×h，h是坐标的放大或缩小系数，用于调整u_j坐标的范围，使u_j∈[-h，h]；S2：造影图像和蒙片图像都具有H行W列像素，对高斯滤波后的造影图像进行分块，每个图像块有K×L个像素，其中L是行数，K是列数，K和L取值为2的倍数，且L≤H，K≤W；S3：取出第(i，j)个造影图像块，其中造影图像像素在整幅造影图像中的位置是(x，y)，其中，造影图像块左上角像素的位置作为造影图像块进行相关性比较的基准点，所述基准点在整幅造影图像中的位置为(x₀，y₀)，其中x₀＝i×K，y₀＝j×L；S4：将第(i，j)个造影图像块放在对应的蒙片图像(x，y)位置上，对应的K×L大小的蒙片图像块向周围延伸行列；S5：让K×L大小的造影图像块在2K×2L大小的蒙片图像块内平移(r，s)，其中平移后，对每个K×L造影图像块与对应的K×L蒙片图像块进行相关性比较，相关性判别用如下公式：C(r,s)=|A(r,s)-B(r,s)|(KL)2]]>其中，A(r,s)=KLΣx=x0x0+K-1Σy=y0y0+L-1(gc(x,y)-gm(x+r,y+s))2]]>B(r,s)=(Σx=x0x0+K-1Σy=y0y0+L-1(gc(x,y)-gm(x+r,y+s)))2]]>其中，g_c(x，y)是造影图像在坐标(x，y)处的像素值，g_m(x，y)是蒙片图像在坐标(x，y)的像素值，根据C(r，s)的值来判断相关性，C(r，s)越小相关性越大；S6：重复步骤S5若干次，每次取不同的r和s值，得到若干相关性系数C(r，s)，对这些系数排序，取其中最小的C(r，s)，并保存r和s的值；S7：将造影图像块在蒙片图像块中平移(r，s)，其中，r和s为步骤S6中保存的r和s，然后按以下公式减影，得到减影图像：g_d(x，y)＝g_c(x，y)-g_m(x+r，y+s)+d_mid其中，d_mid是中间灰度，g_d(x，y)为减影图像在坐标(x，y)处的像素值。