[Numpy-discussion] 2D Hamming window

Henrik Ronellenfitsch henrik-web@phoboid....
Sun Jul 27 04:10:37 CDT 2008


Hi,
Gary Ruben wrote:
> import numpy as np
> import scipy.signal as ss
> 
> # read heightmap here - in my case it's a square numpy float array
> 
> # build 2d window
> hm_len = heightmap.shape[0]
> bw2d = np.outer(ss.hamming(hm_len), np.ones(hm_len))
> bw2d = np.sqrt(bw2d * bw2d.T) # I don't know whether the sqrt is correct
> 
> # window the heightmap
> heightmap *= bw2d
> 
> --
> Gary R.

Thanks very much for your solution, this is exactly what I needed!
If I'm not mistaken, though, you can achieve the same result with

h = hamming(n)
ham2d = sqrt(outer(h,h))

which is a bit more compact.

Regards,
Henrik


More information about the Numpy-discussion mailing list