三帧差分法python代码

以下是一个简单的三帧差分法的Python代码示例：

import cv2

# 读取视频文件
cap = cv2.VideoCapture('test.mp4')

# 读取第一帧
ret, frame1 = cap.read()

# 读取第二帧
ret, frame2 = cap.read()

# 循环处理每一帧
while cap.isOpened():
    # 读取第三帧
    ret, frame3 = cap.read()
    if not ret:
        break

    # 将图像转换为灰度图像
    gray1 = cv2.cvtColor(frame1, cv2.COLOR_BGR2GRAY)
    gray2 = cv2.cvtColor(frame2, cv2.COLOR_BGR2GRAY)
    gray3 = cv2.cvtColor(frame3, cv2.COLOR_BGR2GRAY)

    # 计算第一帧和第二帧的差异
    diff1 = cv2.absdiff(gray1, gray2)

    # 计算第二帧和第三帧的差异
    diff2 = cv2.absdiff(gray2, gray3)

    # 将两个差异图像进行逐像素比较，得到最终的差异图像
    diff = cv2.bitwise_and(diff1, diff2)

    # 对差异图像进行二值化处理
    thresh = cv2.threshold(diff, 25, 255, cv2.THRESH_BINARY)[1]

    # 对二值化后的图像进行膨胀操作
    kernel = cv2.getStructuringElement(cv2.MORPH_RECT, (5, 5))
    thresh = cv2.dilate(thresh, kernel, iterations=2)

    # 在原始图像上绘制差异区域的边框
    contours, hierarchy = cv2.findContours(thresh, cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_SIMPLE)
    for c in contours:
        if cv2.contourArea(c) < 5000:
            continue
        (x, y, w, h) = cv2.boundingRect(c)
        cv2.rectangle(frame2, (x, y), (x+w, y+h), (0, 255, 0), 2)

    # 显示处理后的图像
    cv2.imshow('frame', frame2)

    # 更新帧
    frame1 = frame2
    frame2 = frame3

    # 按下 q 键退出循环
    if cv2.waitKey(1) & 0xFF == ord('q'):
        break

# 释放视频文件句柄并关闭所有窗口
cap.release()
cv2.destroyAllWindows()

这个代码示例中，我们使用了OpenCV库来读取视频文件、进行图像处理和显示图像。我们首先读取前两帧，并将它们转换为灰度图像。然后，我们读取第三帧，并计算前两帧和后两帧之间的差异。我们通过逐像素比较前两帧和后两帧的差异，得到最终的差异图像。我们将差异图像进行二值化和膨胀处理，以便于识别差异区域。最后，我们在原始图像上绘制差异区域的边框，并显示处理后的图像。我们还可以通过按下“q”键来退出循环。