[发明专利]用图像质量评估作为相似性测度的血管减影图像配准方法

权利要求书

1.一种用图像质量评估作为相似性侧度的的血管减影图像配准方法，其特征在于具体步骤如下：

首先，进行特征点提取；特征点采用手动选取，或者通过特征点选择器选取；在每个m×m像素的窗口中只采一个特征点，m为20-100的整数；

然后，进行特征点对应关系计算；即对提取出的特征点在另一幅图像中寻找对应点；对于一个特征点，另一幅图像中的局部图像和这个特征点最相似的点被选择为这个特征点的对应点；也就是，对于每个特征点(x，y)，在另一幅图像中找到满足以下条件的对应点(x_p，y_p)：

(xp,yp)=argmaxxp,ypS(Il(x,y),Im(xp,yp))---(1)]]>

其中，S(.，.)是SSIM测度，I_m和I_l是在蒙片和活片上对应的点的局部图像区域；

限制|x_p-x|和|y-y_p|分别不超过n；达到最大相似度的特征点的对应关系是通过穷举这个区域内的所有点的相似度，并且找出其中的相似度最大的点来达到的；n取10-50的整数；

最后，根据得到的特征点的对应关系计算出两幅图片之间的非线性变换，并且通过这个非线性变换配准两幅图像；其中，使用薄板样条函数来实现配准；得到非线性变换的过程如下：

根据得到的特征点的对应关系U＝{u_a：a＝1，2，...，n}和V＝{v_a：a＝1，2，...，n}；薄板样条函数表示成以下的形式：

ETPS(f)=Σa=1n||ua-f(va)||2+λ∫∫[(∂2f∂x2)+2(∂2f∂x∂y)+(∂2f∂y2)]dxdy---(2)]]>

其中f是在v_a和u_a之间的变换，λ是一个用来调节薄板样条函数平滑性的约束函数；λ取值为0.1-10；对于固定的λ，存在着函数f(v)可以表示成：

f(v)＝v.d+φ(v).w                    (3)

其中v是需要变换的点集，d是一个3×3的仿射变换，w是一个n×3的非仿射变换矩阵；φ(v)是一个由薄板样条核生成的1×n向量；对于v中的每一个点，存在φ_a(v)，表示为φ_a(v)＝c||v-v_a||²log||v-v_a||，其中c是一个常量；

将(3)式代入(2)式，可以得到

E_TPS(d，w)＝‖U-Vd-φw||²+λtrace(w^TΦw)         (4)

其中U和V分别是由u_a和v_a构成的矩阵；Φ是一个由φ(v_a)构成的n×n的矩阵；

然后通过QR分解分离出仿射和非仿射的变换：

V=[Q1Q2]R0---(5)]]>

其中Q₁和Q₂是正交矩阵；R是一个上三角阵；最后，w和d的解表示为：

w=Q2(Q2TΦQ2+λIN-3)-1Q2TU---(6)]]>

d=R-1Q1T(U-Φw)---(7)]]>

由此得非线性变换函数f(v)，然后函数对图像进行非线性变换，即得到配准后的图像。