python opencv 图像相减（原始图像与梯度图像）

 

import cv2
import numpy as np
import matplotlib.pyplot as plt
import sys

img1 = cv2.imread(file1)
img2 = cv2.imread(file2)
size = 512
img1 = cv2.fastNlMeansDenoisingColored(img1,None,10,10,7,21)
img2 = cv2.fastNlMeansDenoisingColored(img2,None,10,10,7,21)

img1 = cv2.resize(img1, (size, size)) 
img2 = cv2.resize(img2, (size, size)) 

#img1=cv2.Laplacian(img1,cv2.CV_64F)#CV_64F为图像深度
#img2=cv2.Laplacian(img2,cv2.CV_64F)#CV_64F为图像深度

#img1 = cv2.Canny(img1,100,200)#其他的默认
#img2 = cv2.Canny(img2,100,200)#其他的默认

is_grad = sys.argv[1] == '1'

if is_grad :
    x = cv2.Sobel(img1,cv2.CV_16S,1,0)
    y = cv2.Sobel(img1,cv2.CV_16S,0,1)
    absX = cv2.convertScaleAbs(x)   # 转回uint8
    absY = cv2.convertScaleAbs(y)
    img1 = cv2.addWeighted(absX,0.5,absY,0.5,0)

    x = cv2.Sobel(img2,cv2.CV_16S,1,0)
    y = cv2.Sobel(img2,cv2.CV_16S,0,1)
    absX = cv2.convertScaleAbs(x)   # 转回uint8
    absY = cv2.convertScaleAbs(y)
    img2 = cv2.addWeighted(absX,0.5,absY,0.5,0)

print('img1', img1.shape[1],img1.shape[0])
print('img2', img2.shape[1],img2.shape[0])

err = cv2.absdiff(img1,img2) 

fig = plt.figure('result')
plt.axis('off')  #关闭坐标轴

plt.subplot(2,2,1)     #将窗口分为两行两列四个子图，则可显示四幅图片
plt.title('img1')   #第一幅图片标题
plt.imshow(img1)  

plt.subplot(2,2,2)     
plt.title('img2')   
plt.imshow(img2)  

plt.subplot(2,2,4)     
plt.title('diff')   
plt.imshow(err)

fig.tight_layout()#调整整体空白
plt.subplots_adjust(wspace =0)#调整子图间距
plt.show()   #显示