[发明专利]图像处理方法及装置

权利要求书

1.一种图像处理方法，多次反复进行以下工序以进行第1及第2图像的对位：

第1工序，按照包含同一拍摄对象的所述第1及第2图像中的每个坐标，决定该坐标上的像素值的梯度所涉及的特征向量，其中当该坐标上的像素值的梯度向量的大小为既定阈值以上时，将以该梯度向量的大小对所述梯度向量进行归一化后的值决定为该坐标上的特征向量，当所述梯度向量的大小小于所述既定阈值时，将0（零）向量决定为该坐标上的特征向量；

第2工序，按照所述第1及第2图像中互相对应的每个坐标，算出与计算该坐标上的特征向量彼此的内积的绝对值的N次方而得的值相当的相关值，求出包含按每个所述坐标算出的相关值的累计值的评价值，其中N为自然数；以及

第3工序，以使所述评价值变得更大的方式变更所述第2图像。

2.一种图像处理装置，其中包括：

决定单元，按照包含同一拍摄对象的第1及第2图像中的每个坐标，当该坐标上的像素值的梯度向量的大小为既定阈值以上时，将以该梯度向量的大小对所述梯度向量进行归一化后的值决定为该坐标上的特征向量，当所述梯度向量的大小小于所述既定阈值时，将0（零）向量决定为该坐标上的特征向量；

算出单元，按照所述第1及第2图像中互相对应的每个坐标，算出与计算该坐标上的特征向量彼此的内积的绝对值的N次方而得的值相当的相关值，求出包含按每个所述坐标算出的相关值的累计值的评价值，其中N为自然数；以及

变更单元，以使所述评价值变得更大的方式变更所述第2图像，

多次反复进行所述决定单元、算出单元及变更单元所进行的处理，进行所述第1及第2图像的对位。

3.如权利要求2所述的图像处理装置，其中所述梯度向量为图像中各坐标轴向的1次偏微分。

4.如权利要求2所述的图像处理装置，其中所述相关值是计算所述特征向量彼此的内积的平方而得的值。

5.如权利要求2所述的图像处理装置，其中所述评价值是以所述第1或第2图像的大小对所述相关值的累计值进行归一化而得的值。

6.如权利要求2所述的图像处理装置，其中所述变更包含平行移动、旋转及变形之中的至少一种。

7.如权利要求2所述的图像处理装置，其中所述决定单元基于对所述第1及第2图像分别进行构造物强调滤波处理和/或噪声降低处理而得到的图像，求出所述梯度向量。

8.如权利要求2所述的图像处理装置，其中所述第1及第2图像为医用图像。

9.如权利要求8所述的图像处理装置，其中所述第1及第2图像中摄像模式互相不同。

10.如权利要求9所述的图像处理装置，其中所述第1及第2图像的一个为超声波图像、即US图像。

11.如权利要求2所述的图像处理装置，其中反复进行所述决定单元、算出单元及变更单元所进行的处理，直至进行该处理的次数达到既定数，或直至所述评价值实质上收敛。

12.如权利要求2所述的图像处理装置，其中还包括推定所述第1及第2图像的至少一个的对象区域中的噪声量，并基于该噪声量决定所述阈值的单元。

13.如权利要求12所述的图像处理装置，其中所述对象区域为所述第1图像与所述第2图像之间共同的解剖学特征点的周边区域。

14.如权利要求12所述的图像处理装置，其中所述对象区域为所述至少一个图像的中央区域。

15.如权利要求12所述的图像处理装置，其中所述对象区域为病变部的周边区域。

16.如权利要求12所述的图像处理装置，其中所述对象区域为全部图像区域。

17.如权利要求12所述的图像处理装置，其中所述决定阈值的单元求出所述对象区域中的各像素上的梯度强度，并基于该梯度强度的直方图来决定所述阈值。

18.如权利要求17所述的图像处理装置，其中所述决定阈值的单元根据以下的（a）～（f）式的任一个或多个的组合的平均值来算出噪声量Nqt：

Nqt＝M_fmax＋σ×HWHM_L（a）

Nqt＝M_fmax＋σ×HWHM_R（b）

Nqt＝M_fmax＋σ×(HWHM_L＋HWHM_R)/2（c）

Nqt＝σ×M_fmax（d）

Nqt＝σ×M_mom1（e）

Nqt＝σ×M_mom2（f）

其中，M_fmax是在所述直方图分布上出现的峰值之中、梯度强度最低的峰值给出众数值的梯度强度；

σ为任意常数，

HWHM_L是从所述直方图分布上的M_fmax观看低值侧的半值半宽度；

HWHM_R是从所述直方图分布上的M_fmax观看高值侧的半值半宽度；

M_mom1是与所述直方图分布上的梯度强度从0到HWHM_R的范围的重心相当的梯度强度；

M_mom2是与所述直方图分布上的梯度强度从HWHM_L到HWHM_R的范围的重心相当的梯度强度。