In [13]: a Out[13]: array([[ 0., 0., 0., 0., 0.], [ 0., 1., 1., 1., 0.], [ 0., 1., 2., 1., 0.], [ 0., 1., 1., 1., 0.], [ 0., 0., 0., 0., 0.]]) In [14]: ndimage.grey_erosion(a,size=(3,3)) Out[14]: array([[ 0., 0., 0., 0., 0.], [ 0., 0., 0., 0., 0.], [ 0., 0., 1., 0., 0.], [ 0., 0., 0., 0., 0.], [ 0., 0., 0., 0., 0.]]) In [15]: ndimage.grey_dilation(a,size=(3,3)) Out[15]: array([[ 1., 1., 1., 1., 1.], [ 1., 2., 2., 2., 1.], [ 1., 2., 2., 2., 1.], [ 1., 2., 2., 2., 1.], [ 1., 1., 1., 1., 1.]])