In [5]:

    

from PIL import Image
from numpy import *
from pylab import *
from scipy.ndimage import filters
import harris


    

In [6]:

    

im = array(Image.open('../data/empire.jpg').convert('L'))


    

In [36]:

    

im2 = zeros(im.shape)
figure(figsize=(8,8))
gray()
for sigma in range(6):
    im2 = filters.gaussian_filter(im, sigma*2)
    harrisim = harris.compute_harris_response(im2)
    filtered_coords = harris.get_harris_points(harrisim, 10, 0.1)
    print len(filtered_coords)
    harris.plot_harris_points(im2, filtered_coords)
show()


    

    




380






    






<matplotlib.figure.Figure at 0x7f0c28ad5150>






    













    




390






    













    




428






    













    




446






    













    




521






    













    




539






    

In [ ]:

    

#The number of points increased. Weird...