From 889df76ea7c38ad20031b678a582f77506e6c520 Mon Sep 17 00:00:00 2001
From: Alexander Mordvintsev <no@email>
Date: Sun, 10 Jun 2012 08:38:25 +0000
Subject: [PATCH] added lappyr.py sample

---
 samples/python2/common.py |  4 +++
 samples/python2/lappyr.py | 64 +++++++++++++++++++++++++++++++++++++++
 2 files changed, 68 insertions(+)
 create mode 100644 samples/python2/lappyr.py

diff --git a/samples/python2/common.py b/samples/python2/common.py
index 0eb49b8ea6..085accd2be 100644
--- a/samples/python2/common.py
+++ b/samples/python2/common.py
@@ -190,3 +190,7 @@ def mosaic(w, imgs):
     imgs = it.chain([img0], imgs)
     rows = grouper(w, imgs, pad)
     return np.vstack(map(np.hstack, rows))
+
+def getsize(img):
+    h, w = img.shape[:2]
+    return w, h
diff --git a/samples/python2/lappyr.py b/samples/python2/lappyr.py
new file mode 100644
index 0000000000..6199c22e7b
--- /dev/null
+++ b/samples/python2/lappyr.py
@@ -0,0 +1,64 @@
+''' An example of Laplacian Pyramid construction and merging.
+
+Level : Intermediate
+
+Usage : python lappyr.py [<video source>]
+
+References:
+  http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.54.299
+
+Alexander Mordvintsev 6/10/12
+'''
+
+import numpy as np
+import cv2
+import video
+from common import nothing, getsize
+
+def build_lappyr(img, leveln=6, dtype=np.int16):
+    img = dtype(img)
+    levels = []
+    for i in xrange(leveln-1):
+        next_img = cv2.pyrDown(img)
+        img1 = cv2.pyrUp(next_img, dstsize=getsize(img))
+        levels.append(img-img1)
+        img = next_img
+    levels.append(img)
+    return levels
+
+def merge_lappyr(levels):
+    img = levels[-1]
+    for lev_img in levels[-2::-1]:
+        img = cv2.pyrUp(img, dstsize=getsize(lev_img))
+        img += lev_img
+    return np.uint8(np.clip(img, 0, 255))
+
+
+if __name__ == '__main__':
+    import sys
+    print __doc__
+
+    try: fn = sys.argv[1]
+    except: fn = 0
+    cap = video.create_capture(fn)
+
+    leveln = 6
+    cv2.namedWindow('level control')
+    for i in xrange(leveln):
+        cv2.createTrackbar('%d'%i, 'level control', 5, 50, nothing)
+
+    while True:
+        ret, frame = cap.read()
+
+        pyr = build_lappyr(frame, leveln)
+        for i in xrange(leveln):
+            v = cv2.getTrackbarPos('%d'%i, 'level control') / 5
+            pyr[i] *= v
+        res = merge_lappyr(pyr)
+
+        cv2.imshow('laplacian pyramid filter', res)
+        
+        if cv2.waitKey(1) == 27:
+            break
+
+