title: How to make ndimage.filters.maximum_filter work like MATLAB's imregionalmax function?
body: <p>After reading <a href="https://stackoverflow.com/questions/27598103/what-is-the-difference-between-imregionalmax-of-matlab-and-scipy-ndimage-filte">this post</a>, and play also with SciKit-image I found a difference in Python compared to MATLAB's function <code>imregionalmax</code>.</p>
  
  <p>I have these lines of code:</p>
  
  <pre><code>from skimage.feature import peak_local_max
  
  manos = np.ones([5,5])
  manos[2,2] = 0.
  manos[2,4] = 2.
  
  giannis = peak_local_max(manos,min_distance=1, indices=False, exclude_border=False)
  
  giorgos = ndimage.filters.maximum_filter(manos, footprint=np.ones([3,3]))
  giorgos = (giorgos == manos)
  </code></pre>
  
  <p>I would expect a 2D array with only one <code>True</code> value ([2,4]) for variables <code>giannis</code> or <code>giorgos</code> as I get in MATLAB. Instead I take more than one maximum. </p>
  
  <p>Any idea why this works this way and how to make it work like in MATLAB?</p>