add local test

This commit is contained in:
Olivier Debeir
2012-10-10 17:32:35 +02:00
parent fee2df8d42
commit 29d133b290
2 changed files with 29 additions and 1 deletions
+1 -1
View File
@@ -26,7 +26,7 @@ See Wikipedia_ for more details on the algorithm.
"""
import numpy as np
from scipy import ndimage
e
import matplotlib.pyplot as plt
from skimage.morphology import watershed, is_local_maximum
@@ -0,0 +1,28 @@
import numpy as np
import matplotlib.pyplot as plt
import gdal
from skimage.morphology import disk
import skimage.rank as rank
filename = 'iko_pan_Ja1.tif'
im16 = gdal.Open(filename).ReadAsArray().astype(np.uint16)
plt.figure()
plt.imshow(im16,cmap=plt.cm.gray)
plt.colorbar()
f0 = rank.median(im16,disk(1))
f1 = rank.bilateral_mean(im16,disk(20),s0=200,s1=200)
f2 = rank.equalize(f1,disk(10))
f3 = rank.bottomhat(f1,disk(1))
plt.figure()
plt.imshow(f2,cmap=plt.cm.gray,interpolation='nearest')
plt.colorbar()
plt.show()