Source Code for Module lib.hpfilter

 1  #!/usr/bin/env python 
 2   
 3   
 4  import gpkimgclass 
 5   
 6  import sys 
 7  import os 
 8  import voice_misc 
 9   
10   
11  if __name__ == '__main__': 
12          cutoff_freq = float(sys.argv[1]) 
13          assert cutoff_freq > 0 
14   
15          d = gpkimgclass.read(sys.argv[2]) 
16          print 'shape=', d.d.shape 
17          filtered = voice_misc.hipass_sym_butterworth(d.d[:,0], cutoff_freq*d.dt(), order=4) 
18          gpkimgclass.gpk_img(d.hdr, filtered).write(sys.argv[3]) 
19