{"id":63,"date":"2006-06-13T08:06:53","date_gmt":"2006-06-13T12:06:53","guid":{"rendered":"https:\/\/blogs.mathworks.com\/steve\/?p=63"},"modified":"2019-10-22T12:29:30","modified_gmt":"2019-10-22T16:29:30","slug":"spatial-transformations-findbounds","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/steve\/2006\/06\/13\/spatial-transformations-findbounds\/","title":{"rendered":"Spatial transformations: findbounds"},"content":{"rendered":"
\r\n

I've written previously<\/a> about how Image Processing Toolbox uses inverse mapping to implement spatial transforms. In this method, you set up a grid\r\n in output space. For each pixel in the output space grid, you use the inverse transform to determine the corresponding input\r\n space location, and you interpolate to get the input image pixel value at that point.\r\n <\/p>\r\n

I've also written<\/a> about the question of where to establish the grid in output space. If you don't somehow pay attention to the specific transformation,\r\n you might choose a grid that's too small:\r\n <\/p>\r\n

<\/p>\r\n

Or your grid might be too big:<\/p>\r\n

<\/p>\r\n

Or your grid might be in the wrong place completely:<\/p>\r\n

<\/p>\r\n

To avoid these problems (and to avoid the tech support calls!), the function imtransform<\/tt> uses the function findbounds<\/tt> to determine automatically where to place the output grid and how big it should be.\r\n <\/p>\r\n

findbounds<\/tt> first creates a grid of points. These points are located in input space at each corner of the image, as well as in between\r\n each corner and in the middle. It looks like this:\r\n <\/p>

I = imread('rice.png'<\/span>);\r\nh = imshow(I);\r\nset(h,'AlphaData'<\/span>,0.3);\r\naxis on<\/span>, grid on<\/span>\r\nin_points = [ ...<\/span>\r\n    0.5000    0.5000\r\n    0.5000  256.5000\r\n  256.5000    0.5000\r\n  256.5000  256.5000\r\n    0.5000  128.5000\r\n  128.5000    0.5000\r\n  128.5000  128.5000\r\n  128.5000  256.5000\r\n  256.5000  128.5000];\r\nhold on<\/span>\r\nplot(in_points(:,1), in_points(:,2), '.'<\/span>, 'MarkerSize'<\/span>, 18)\r\nhold off<\/span><\/pre> 
findbounds<\/tt> then calls tformfwd<\/tt> to transform these points into output space.  (For this example I'm going to construct an affine tform struct that combines\r\n      scaling, rotation, shear, and a large translation.)\r\n   <\/p>
tform = maketform('affine'<\/span>, ...<\/span>\r\n    [1.1067 -0.2341 0; 0.5872 1.1769 0; 1000 -300 1]);\r\nout_points = tformfwd(tform, in_points)<\/pre>
\r\nout_points =\r\n\r\n  1.0e+003 *\r\n\r\n    1.0008   -0.2995\r\n    1.1512    0.0018\r\n    1.2842   -0.3595\r\n    1.4345   -0.0582\r\n    1.0760   -0.1489\r\n    1.1425   -0.3295\r\n    1.2177   -0.1789\r\n    1.2928   -0.0282\r\n    1.3593   -0.2088\r\n\r\n<\/pre>
plot(out_points(:,1), out_points(:,2), '.'<\/span>, 'MarkerSize'<\/span>, 18)\r\naxis ij<\/span>, axis image<\/span>\r\ngrid on<\/span><\/pre> 
The bounding box of the points in output space tells us where to put the output space grid.<\/p>\r\n   
When you call imtransform<\/tt>, you can use optional output arguments to determine where the image is located in output space.\r\n   <\/p>
[J,XData,YData] = imtransform(I, tform);\r\nXData<\/pre>
\r\nXData =\r\n\r\n  1.0e+003 *\r\n\r\n    1.0017    1.4337\r\n\r\n<\/pre>
YData<\/pre>
\r\nYData =\r\n\r\n -358.7527    1.2473\r\n\r\n<\/pre>
You can use this information together with imshow<\/tt> to place the image in the right place on an axes that contains other X-Y data.\r\n   <\/p>
h = imshow(J,'XData'<\/span>,XData,'YData'<\/span>,YData);\r\nset(h,'AlphaData'<\/span>,0.3)\r\nhold on<\/span>\r\nplot(out_points(:,1), out_points(:,2), '.'<\/span>, 'MarkerSize'<\/span>, 18)\r\naxis on<\/span>\r\ngrid on<\/span>\r\nhold off<\/span>\r\naxis ij<\/span>, axis image<\/span><\/pre> 
So that's the method imtransform<\/tt> uses to determine where the output image should be in output space.\r\n   <\/p>\r\n   
In a future posting I'll cover two more topics related to findbounds<\/tt> and imtransform<\/tt>:\r\n   <\/p>\r\n   
\r\n      
    \r\n         
  • What to do when there's no forward transform in the tform structure<\/li>\r\n         
  • The special problem of pure translation<\/li>\r\n      <\/ul>\r\n   <\/div>\r\n