Hello, >From the following script: import numpy from scipy.ndimage.interpolation import zoom n = 5 Z = numpy.ones((n)) print Z print zoom(Z,1,order=0) print zoom(Z.reshape((1,n)),(1,1),order=0) I get the following result: [ 1. 1. 1. 1. 1.] [ 1. 1. 1. 1. 1.] [[ 0. 0. 0. 0. 0.]] Is that the expected behavior ? Nicolas