OpenCV图像减法与Numpy减法

时间:2017-08-22 11:58:42

标签: python python-3.x opencv numpy image-processing

我试图跟踪下一张图片与前一张图片的差异,假设场景中有一些移动。决定在两个jpg图像之间应用相应像素值的减法,然后计算所得矩阵的平均值,以便检查它是否低于或低于某个阈值水平(用于进一步分析)。

通过cv2.subtract和np.subtract方法进行减法。我注意到结果差异很大。似乎numpy以某种方式拉伸直方图并将结果值归一化,但为什么呢?

图像是通过cv2.open加载的。我知道这种方法使用频道的BGR顺序,但它没有解释发生了什么。加载的图像是numpy nd.array与np.uint值。使用Python 3.7开发Spyder。

编辑:cv2.imread中的参数0告诉以灰度加载图像

OpenCV subtraction result

Numpy subtraction result

#loading images

img_cam0 = cv2.imread(r'C:\Users\Krzysztof\Desktop\1.jpg',0)
img_cam1 = cv2.imread(r'C:\Users\Krzysztof\Desktop\2.jpg', 0)
print('img0 type:',type(img_cam0), 'and shape:', img_cam0.shape)
print('img1 type:',type(img_cam1),'and shape:', np.shape(img_cam1))
print('\n')

#opencv subtraction

cv2_subt = cv2.subtract(img_cam0,img_cam1)
cv2_mean = cv2.mean(cv2_subt)

print('open cv mean is:', cv2_mean)
f.show_im(cv2_subt, 'cv2_subtr')

#np subtraction and mean

np_subtr = np.subtract(img_cam0, img_cam1)
np_mean = np.mean(np_subtr)

print('numpy mean is:', np_mean)
f.show_im(np_subtr, 'np_subtr')

1 个答案:

答案 0 :(得分:11)

差异很简单 - 饱和度与无饱和度。

cv2.subtract执行饱和度。根据文件:

dst(I) = saturate(src1(I) - src2(I))

numpy.subtract只执行常规减法,因此结果以integer overflow为准(即值包围)。

  

饱和度表示当输入值v超出目标类型的范围时,结果不仅仅是通过输入低位来形成,而是形成值被修剪。例如:

uchar a = saturate_cast<uchar>(-100); // a = 0 (UCHAR_MIN)
short b = saturate_cast<short>(33333.33333); // b = 32767 (SHRT_MAX)
     

当目标类型为unsigned charsigned charunsigned shortsigned short时,会执行此类裁剪。对于32位整数,不进行削波。

实施例

>>> import cv2
>>> import numpy as np

>>> a = np.arange(9, dtype=np.uint8).reshape(3,3)
>>> a
array([[0, 1, 2],
       [3, 4, 5],
       [6, 7, 8]], dtype=uint8)
>>> b = np.full((3,3), 4, np.uint8)
>>> b
array([[4, 4, 4],
       [4, 4, 4],
       [4, 4, 4]], dtype=uint8)

>>> np.subtract(b,a)
array([[  4,   3,   2],
       [  1,   0, 255],
       [254, 253, 252]], dtype=uint8)

>>> cv2.subtract(b,a)
array([[4, 3, 2],
       [1, 0, 0],
       [0, 0, 0]], dtype=uint8)