[Scipy-tickets] [SciPy] #515: ndimage.zoom introduces artefacts at right and bottom edges

SciPy scipy-tickets@scipy....
Wed Oct 17 19:35:20 CDT 2007

#515: ndimage.zoom introduces artefacts at right and bottom edges
 Reporter:  0ion9                                  |        Owner:  stefan  
     Type:  defect                                 |       Status:  assigned
 Priority:  normal                                 |    Milestone:  0.7     
Component:  scipy.ndimage                          |      Version:  devel   
 Severity:  normal                                 |   Resolution:          
 Keywords:  ndimage, data, representation, spline  |  
Comment (by stefan):

 Thanks for reporting this problem.  It is a result of a more complicated
 issue: the way we represent images.  Given values


 does the '0' represent a little block/pixel of size 1x1, or does it just
 represent a data-point at (0,0).  The code in ndimage confuses the issue
 somewhat: it represents the data internally as data-points, but clearly
 algorithms like zoom see the data as pixels.

 Let me explain further.  Zoom currently functions like this: if the input
 size is (m,n) then the output size is (2m,2n).  Therefore, at each
 coordinate in the output, the corresponding coordinate in the input is
 expected to be (x/2,y/2).  This isn't quite right.  Say that we start with
 a (5,5) array, and map it to a (10,10).  Now we are saying that position
 (9,9) should map to (4.5,4.5).  But (4.5,4.5) is *outside* the input
 array, and that is why you see the zeroes.  If the image were represented
 as pixel blocks, (4.5,4.5) would have been valid.

 Then, it would seem like a simple matter to change the image
 representation to be pixels, instead of data-points.  Unfortunately, that
 would require hacking the spline engine.  I'm not sure it would achieve
 much, either.  Say you extend your data using a constant value and
 interpolate, then zoom would still end up producing a black border around
 the bottom edges.

 I would therefore propose the following:

 1) Input data represent data points.  The borders for


 therefore lie at row 0, column 0, row 2 and column 2 (on top of the
 values, not to the left, right, top or bottom of them).  Coordinates like
 (2.5,2.5) fall *outside* the image, and must be extended according to the
 given mode.

 2) All ndimage functions should respect this data-layout.  For example,
 zooming by a factor of two should be internally translated to zooming by a
 factor of (output_dim - 1)/(input_dim - 1).  This ensures that the input
 boundary maps onto the output boundary.

 3) Make 'mirror' and 'reflect' extension modes behave the same.  Due to
 the data-point representation, only 'reflect' makes any sense (I got this
 the wrong way around in the current version):

 Length: 3
 Index:         0, 1, 2, 3, 4, 5
 Reflect Index: 0, 1, 2, 1, 0, 1

Ticket URL: <http://scipy.org/scipy/scipy/ticket/515#comment:2>
SciPy <http://www.scipy.org/>
SciPy is open-source software for mathematics, science, and engineering.

More information about the Scipy-tickets mailing list