def preproc_dog(bg):
"""test low level image filter
"""
#special dependencies
from scipy.ndimage import gaussian_filter
from pyrankfilter import rankfilter,__version__
bg = bg.astype(float)
f1 = gaussian_filter(bg,sigma=14.5)
f2 = gaussian_filter(bg,sigma=15)
f = ((f1-f2+3)*10).astype('uint8')
loc_max = rankfilter(f,'highest',20,infSup=[0,0])
r = bg.copy()
r[loc_max>0]=0
plt.figure()
plt.subplot(2,2,1)
plt.imshow(bg)
plt.subplot(2,2,2)
plt.imshow(f)
plt.colorbar()
plt.subplot(2,2,3)
plt.imshow(loc_max)
plt.subplot(2,2,4)
plt.imshow(r)
plt.show()
def preproc(bg):
"""test low level image filter
"""
#special dependencies
from pyrankfilter import rankfilter,__version__
print __version__
#custom structuring element
xx,yy = npy.meshgrid(npy.linspace(-50,50,100),npy.linspace(-50,50,100))
zz = npy.sqrt(xx**2+yy**2)
struct_elem = npy.logical_and(zz>10,zz<20)
plt.figure()
f = rankfilter(bg,'median',struct_elem=struct_elem.astype('uint8'))
plt.imshow(f)
plt.show()
