In [2]:

    

from numpy import *
from numpy import random
from scipy.ndimage import filters
from PIL import *
from pylab import *
from scipy.ndimage import measurements, morphology


    

In [3]:

    

im2 = array(Image.open('data/houses.png').convert('L'))
im2 = 1*(im2<64)
im_open = morphology.binary_opening(im2, ones((6, 3)), iterations=4)
figure(figsize=(8, 12))
subplot(1, 2, 1)
imshow(im2)
subplot(1, 2, 2)
imshow(im_open)
show()


    

In [4]:

    

labels2, nbr_objects2 = measurements.label(im_open)
print "Number of objects", nbr_objects2


    

    




Number of objects 47

In [5]:

    

center = measurements.center_of_mass(im_open, labels2, range(nbr_objects2))


    

    




/usr/local/lib/python2.7/dist-packages/scipy/ndimage/measurements.py:1301: RuntimeWarning: invalid value encountered in true_divide
  for dir in range(input.ndim)]

In [17]:

    

figure(figsize=(8, 8))
imshow(im_open)
for i in range(nbr_objects2):
    plot(center[i][1], center[i][0], "*")
show()