/usr/share/doc/python-pywt/examples/wavedec.py is in python-pywt 0.2.0-5.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 | #!/usr/bin/env python
# -*- coding: utf-8 -*-
import pywt
data = range(16)
wavelet = 'db4'
level = 2
mode = 'cpd'
print "original data:"
print data
print
# dec = [cA(n-1) cD(n-1) cD(n-2) ... cD(2) cD(1)]
dec = pywt.wavedec(data, wavelet, mode, level)
print "decomposition:"
print "cA%d:" % (len(dec)-1)
print ' '.join([("%.3f" % val) for val in dec[0]])
for i,d in enumerate(dec[1:]):
print "cD%d:" % (len(dec)-1-i)
print ' '.join([("%.3f" % val) for val in d])
print
print "reconstruction:"
print ' '.join([("%.3f" % val) for val in pywt.waverec(dec, wavelet, mode)])
|