第五次作业整理（视频与图像处理2/20211015）

目录

利用均值哈希判断两张图片是否相同

利用判断图片是否相同的方法判断分镜

均值哈希

感知哈希

直方图

​ 利用ffmpeg截取一分钟视频

---

利用均值哈希判断两张图片是否相同

import cv2
import numpy as np
import matplotlib.pyplot as plt
import os
from PIL import Image
 
print(os.getcwd())
os.chdir(r"/Users/mac/Desktop/courses/ad/python/20211015")
 
#均值哈希算法
def aHash(img):
    #缩放为8*8
    plt.imshow(img)
    plt.axis('off')
    plt.show()
    img=cv2.resize(img,(8,8))
    plt.imshow(img)
    plt.axis('off')
    plt.show()
    
    #转换为灰度图
    gray=cv2.cvtColor(img,cv2.COLOR_BGR2GRAY)

    #s为像素和，初值为0；hash_str为hash值，初值为"
    s=0
    hash_str=""
    #遍历累加求像素和
    for i in range(8):
        for j in range(8):
            s=s+gray[i,j]
    #求平均灰度
    avg=s/64
    #灰度大于平均值为1相反为0生成图片的hash值
    for i in range(8):
        for j in range(8):
            if gray[i,j]>avg:
                hash_str=hash_str+'1'
            else:
                hash_str=hash_str+'0'
    return hash_str
 
#Hash值对比
def cmpHash(hash1,hash2):
    n=0
    print(hash1)
    print(hash2)
    #hash长度不同则返回-1代表传参出错
    if len(hash1)!=len(hash2):
        return-1
    #遍历判断
    for i in range(len(hash1)):
        #不相等则n计数+1，n最终为相似度
        if hash1[i]!=hash2[i]:
            n=n+1
    return n
 
img1 = cv2.imread('/Users/mac/Desktop/courses/ad/python/20211008/pic/image34.jpg')
img2 = cv2.imread('/Users/mac/Desktop/courses/ad/python/20211008/pic/image57.jpg')
 
 
hash1 = aHash(img1)
hash2 = aHash(img2)
n = cmpHash(hash1, hash2)
print('均值哈希算法相似度：', n)

利用判断图片是否相同的方法判断分镜

均值哈希

import cv2
import numpy as np
import matplotlib.pyplot as plt
import os
from PIL import Image
 
print(os.getcwd())
os.chdir(r"/Users/mac/Desktop/courses/ad/python/20211015")
 
#均值哈希算法
def aHash(img):
    #缩放为8*8
    plt.imshow(img)
    plt.axis('off')
    plt.show()
    img=cv2.resize(img,(8,8))
    plt.imshow(img)
    plt.axis('off')
    plt.show()
    
    #转换为灰度图
    gray=cv2.cvtColor(img,cv2.COLOR_BGR2GRAY)
    #s为像素和，初值为0；hash_str为hash值，初值为"
    s=0
    hash_str=""
    #遍历累加求像素和
    for i in range(8):
        for j in range(8):
            s=s+gray[i,j]
    #求平均灰度
    avg=s/64
    #灰度大于平均值为1相反为0生成图片的hash值
    for i in range(8):
        for j in range(8):
            if gray[i,j]>avg:
                hash_str=hash_str+'1'
            else:
                hash_str=hash_str+'0'
    return hash_str
 
#Hash值对比
def cmpHash(hash1,hash2):
    n=0
    print(hash1)
    print(hash2)
    #hash长度不同则返回-1代表传参出错
    if len(hash1)!=len(hash2):
        return-1
    #遍历判断
    for i in range(len(hash1)):
        #不相等则n计数+1，n最终为相似度
        if hash1[i]!=hash2[i]:
            n=n+1
    return n
 
for i in range(549):
    img1=cv2.imread('/Users/mac/Desktop/courses/ad/python/20211008/pic/image{}.jpg'.format(i))
    img2=cv2.imread('/Users/mac/Desktop/courses/ad/python/20211008/pic/image{}.jpg'.format(i+1))
    hash1=aHash(img1)
    hash2=aHash(img2)
    n=cmpHash(hash1,hash2)
    if (n>32):
        print('均值哈希算法相似度：',n/64) 
        cv2.imwrite('/Users/mac/Desktop/courses/ad/python/20211015/pic/image{}.jpg'.format(i+1),img2)

 运行结果

感知哈希

def pHash(img):
    # 缩放32*32
    img = cv2.resize(img, (32, 32))  # , interpolation=cv2.INTER_CUBIC
 
    # 转换为灰度图
    gray = cv2.cvtColor(img,cv2.COLOR_BGR2GRAY)
    # 将灰度图转为浮点型，再进行dct变换
    dct = cv2.dct(np.float32(gray))
    # opencv实现的掩码操作
    dct_roi = dct[0:8, 0:8]
 
    hash = []
    avreage = np.mean(dct_roi)
    for i in range(dct_roi.shape[0]):
        for j in range(dct_roi.shape[1]):
            if dct_roi[i, j] > avreage:
                hash.append(1)
            else:
                hash.append(0)
    return hash

直方图

import cv2
import numpy as np
import matplotlib.pyplot as plt
import os
os.chdir(r"/Users/mac/Desktop/courses/ad/python/20211015")
 
 
#通过得到RGB每个通道的直方图来计算相似度
def classify_hist_with_split(image1, image2, size=(256, 256)):
    # 将图像resize后，分离为RGB三个通道，再计算每个通道的相似值
    image1 = cv2.resize(image1, size)
    image2 = cv2.resize(image2, size)
    
    plt.imshow(image1)
    plt.axis('off')
    plt.show()
    
    plt.imshow(image2)
    plt.axis('off')
    plt.show()
    
    sub_image1 = cv2.split(image1)
    sub_image2 = cv2.split(image2)
    sub_data = 0
    
    for im1, im2 in zip(sub_image1, sub_image2):
        sub_data += calculate(im1, im2)
    sub_data = sub_data / 3
    return sub_data
 
# 计算单通道的直方图的相似值
def calculate(image1, image2):
    hist1 = cv2.calcHist([image1], [0], None, [256], [0.0, 255.0])
    hist2 = cv2.calcHist([image2], [0], None, [256], [0.0, 255.0])
    plt.plot(hist1, color="r")
    plt.plot(hist2, color="g")
    plt.show()
    # 计算直方图的重合度
    degree = 0
    for i in range(len(hist1)):
        if hist1[i] != hist2[i]:
            degree = degree + (1 - abs(hist1[i] - hist2[i]) / max(hist1[i], hist2[i]))
        else:
            degree = degree + 1    #统计相似
    degree = degree / len(hist1)
    return degree
 
for i in range(549):
    img1=cv2.imread('/Users/mac/Desktop/courses/ad/python/20211008/pic/image{}.jpg'.format(i))
    img2=cv2.imread('/Users/mac/Desktop/courses/ad/python/20211008/pic/image{}.jpg'.format(i+1))
    n = classify_hist_with_split(img1, img2)
    if (n<0.6):
        print('三直方图算法相似度：',n) 
        cv2.imwrite('/Users/mac/Desktop/courses/ad/python/20211015/pic2/image{}.jpg'.format(i+1),img2)

运行结果

 利用ffmpeg截取一分钟视频

ffmpeg -ss 00:00:00 -t 00:01:00 -accurate_seek -i /Users/mac/Desktop/courses/ad/python/20211015/fsj.mp4 -codec copy fsj-output.mp4