work on motempl.py sample
authorAlexander Mordvintsev <no@email>
Sun, 10 Jul 2011 21:23:32 +0000 (21:23 +0000)
committerAlexander Mordvintsev <no@email>
Sun, 10 Jul 2011 21:23:32 +0000 (21:23 +0000)
samples/python2/common.py
samples/python2/lk_track.py
samples/python2/motempl.py [new file with mode: 0644]

index e5eefd4..636ef69 100644 (file)
@@ -109,3 +109,9 @@ def make_cmap(name, n=256):
         ch = np.interp(xs, xp, yp)\r
         channels.append(ch)\r
     return np.uint8(np.array(channels).T*255)\r
+\r
+def nothing(*arg, **kw):\r
+    pass\r
+\r
+def clock():\r
+    return cv2.getTickCount() / cv2.getTickFrequency()
\ No newline at end of file
index 65346de..9987f26 100644 (file)
@@ -20,7 +20,7 @@ lk_params = dict( winSize  = (3, 3),
                   criteria = (cv2.TERM_CRITERIA_EPS | cv2.TERM_CRITERIA_COUNT, 10, 0.03),\r
                   derivLambda = 0.0 )    \r
 \r
-feature_params = dict( maxCorners = 500, \r
+feature_params = dict( maxCorners = 1000, \r
                        qualityLevel = 0.1,\r
                        minDistance = 5,\r
                        blockSize = 5 )\r
@@ -42,6 +42,8 @@ def main():
     try: video_src = sys.argv[1]\r
     except: video_src = video.presets['chess']\r
 \r
+    cv2.namedWindow('img', 0)\r
+\r
     track_len = 4\r
     tracks = []\r
     cam = video.create_capture(video_src)\r
@@ -63,7 +65,7 @@ def main():
                 tr.append((x, y))\r
                 if len(tr) > 10:\r
                     del tr[0]\r
-                cv2.circle(vis, (x, y), 2, (0, 255, 0), -2)\r
+                cv2.circle(vis, (x, y), 2, (0, 255, 0), -1)\r
             cv2.polylines(vis, [np.int32(tr) for tr in tracks], False, (0, 255, 0))\r
             draw_str(vis, (20, 20), ['new', 'old'][old_mode]+' mode')\r
             draw_str(vis, (20, 40), 'time: %.02f ms' % (dt*1000))\r
diff --git a/samples/python2/motempl.py b/samples/python2/motempl.py
new file mode 100644 (file)
index 0000000..f867b2e
--- /dev/null
@@ -0,0 +1,56 @@
+import numpy as np\r
+import cv2, cv\r
+import video\r
+from common import nothing, clock, draw_str\r
+\r
+MHI_DURATION = 1.0\r
+DEFAULT_THRESHOLD = 16\r
+MAX_TIME_DELTA = 0.5\r
+MIN_TIME_DELTA = 0.05\r
+\r
+\r
+if __name__ == '__main__':\r
+    import sys\r
+    try: video_src = sys.argv[1]\r
+    except: video_src = 'synth:class=chess:bg=../cpp/lena.jpg:noise=0.01'\r
+\r
+    cv2.namedWindow('motempl')\r
+    visuals = ['input', 'frame_diff', 'motion_hist', 'grad_orient']\r
+    cv2.createTrackbar('visual', 'motempl', 2, len(visuals)-1, nothing)\r
+    cv2.createTrackbar('threshold', 'motempl', DEFAULT_THRESHOLD, 255, nothing)\r
+\r
+    cam = video.create_capture(video_src)\r
+    ret, frame = cam.read()\r
+    h, w = frame.shape[:2]\r
+    prev_frame = frame.copy()\r
+    motion_history = np.zeros((h, w), np.float32)\r
+    hsv = np.zeros((h, w, 3), np.uint8)\r
+    hsv[:,:,1] = 255\r
+    while True:\r
+        ret, frame = cam.read()\r
+        frame_diff = cv2.absdiff(frame, prev_frame)\r
+        gray_diff = cv2.cvtColor(frame_diff, cv.CV_BGR2GRAY)\r
+        thrs = cv2.getTrackbarPos('threshold', 'motempl')\r
+        ret, motion_mask = cv2.threshold(gray_diff, thrs, 255, cv2.THRESH_BINARY)\r
+        timestamp = clock()\r
+        cv2.updateMotionHistory(motion_mask, motion_history, timestamp, MHI_DURATION)\r
+        mg_mask, mg_orient = cv2.calcMotionGradient( motion_history, MAX_TIME_DELTA, MIN_TIME_DELTA, apertureSize=5 );\r
+\r
+        visual_name = visuals[cv2.getTrackbarPos('visual', 'motempl')]\r
+        if visual_name == 'input':\r
+            vis = frame.copy()\r
+        elif visual_name == 'frame_diff':\r
+            vis = frame_diff.copy()\r
+        elif visual_name == 'motion_hist':\r
+            vis = np.uint8(np.clip((motion_history-(timestamp-MHI_DURATION)) / MHI_DURATION, 0, 1)*255)\r
+            vis = cv2.cvtColor(vis, cv.CV_GRAY2BGR)\r
+        elif visual_name == 'grad_orient':\r
+            hsv[:,:,0] = mg_orient/2\r
+            hsv[:,:,2] = mg_mask*255\r
+            vis = cv2.cvtColor(hsv, cv.CV_HSV2BGR)\r
+        draw_str(vis, (20, 20), visual_name)\r
+        cv2.imshow('motempl', vis)\r
+\r
+        prev_frame = frame.copy()\r
+        if cv2.waitKey(5) == 27:\r
+            break
\ No newline at end of file