使用 Meanshift 和 Camshift 算法在视频中找到并跟踪
Meanshift 算法
Meanshift 算法的基本原理是和很简单的。假设我们有一堆点,和一个小的圆形窗口,Meanshift 算法就是不断移动小圆形窗口,直到找到圆形区域内最大灰度密度处为止。
# -*- coding: utf-8 -*-
"""
Created on Mon Jan 27 08:02:04 2017
@author: zhang
"""
import numpy as np
import cv2
#打开视频文件
cap = cv2.VideoCapture('slow.flv')
#获取视频的第一帧
ret,frame = cap.read()
# 设置窗口的初始位置
r,h,c,w = 250,90,400,125 # 写死的值
track_window = (c,r,w,h)
# 建立跟踪的ROI
roi = frame[r:r+h, c:c+w]
#图片从RGB空间转换为HSV空间
hsv_roi = cv2.cvtColor(frame, cv2.COLOR_BGR2HSV)
#将低亮度的值忽略掉
mask = cv2.inRange(hsv_roi, np.array((0., 60.,32.)), np.array((180.,255.,255.)))
#计算图像直方图
roi_hist = cv2.calcHist([hsv_roi],[0],mask,[180],[0,180])
#归一化处理
cv2.normalize(roi_hist,roi_hist,0,255,cv2.NORM_MINMAX)
# 设置终止标准,10次迭代或至少移动1 pt
term_crit = ( cv2.TERM_CRITERIA_EPS | cv2.TERM_CRITERIA_COUNT, 10, 1 )
while(1):
ret ,frame = cap.read()
if ret == True:
hsv = cv2.cvtColor(frame, cv2.COLOR_BGR2HSV)
dst = cv2.calcBackProject([hsv],[0],roi_hist,[0,180],1)
# apply meanshift to get the new location
ret, track_window = cv2.meanShift(dst, track_window, term_crit)
# Draw it on image
x,y,w,h = track_window
img2 = cv2.rectangle(frame, (x,y), (x+w,y+h), 255,2)
cv2.imshow('img2',img2)
k = cv2.waitKey(60) & 0xff
if k == 27:
break
else:
cv2.imwrite(chr(k)+".jpg",img2)
else:
break
cv2.destroyAllWindows()
cap.release()