title: `scipy.ndimage.zoom` vs `skimage.transform.rescale` with `order=0`
body: <p>I'm creating a colour-over-grayscale image to show classification of image blocks. I need to upsample a low-resolution classification image to match the original image size, using zero-order.</p>
  
  <p>When I use <code>skimage.transform.rescale(I,16,order=0)</code>, I get the desired result:
  <a href="https://i.stack.imgur.com/dChfX.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/dChfX.png" alt="rescale result"></a>
  However, when I use <code>scipy.ndimage.zoom(I,16,order=0)</code>, I get the following image with misalignment of the upsampled classification blocks and the original image patches (note the non-square blocks at the edges):
  <a href="https://i.stack.imgur.com/kcu3q.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/kcu3q.png" alt="zoom result"></a>
  My questions are:</p>
  
  <ol>
  <li>Is this expected behaviour?</li>
  <li>Why the difference?</li>
  </ol>
  
  <p>I could see the <code>zoom</code> behaviour causing unexpected grief, since I only noticed it here with the obvious misalignment of blocks.</p>
  
  <p>[edit 1]: <a href="https://stackoverflow.com/questions/33721988/scipy-ndimage-zoom-result-depends-on-image-size">perhaps related question</a>.</p>
  
  <p>[edit 2]: Should be able to pull <a href="http://github.com/jessexknight/bayes-texture" rel="nofollow noreferrer">this repo</a> and run python main.py with one of two imports at the top commented (flagged).</p>