{"id":68,"date":"2006-07-07T14:11:04","date_gmt":"2006-07-07T18:11:04","guid":{"rendered":"https:\/\/blogs.mathworks.com\/steve\/?p=68"},"modified":"2019-10-22T13:59:26","modified_gmt":"2019-10-22T17:59:26","slug":"spatial-transformations-translation-confusion","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/steve\/2006\/07\/07\/spatial-transformations-translation-confusion\/","title":{"rendered":"Spatial transformations: Translation confusion"},"content":{"rendered":"
\r\n

The last time I wrote about spatial transformations<\/a>, I explained that imtransform<\/tt> uses the function findbounds<\/tt> to locate the transformed image in output space. By default, imtransform<\/tt> computes an output image grid that is just big enough to capture the output image, wherever it is located.\r\n <\/p>\r\n

We designed imtransform<\/tt> this way because we thought it was what users would expect most of the time. Generally, this default behavior works out\r\n just fine. The big exception, though, is when users call imtransform<\/tt> to apply a pure translation. Let's try it and see.\r\n <\/p>

I = imread('pout.tif'<\/span>);\r\nimshow(I)<\/pre> 
Now make an affine tform struct representing a pure translation:<\/p>
T = [1 0 0; 0 1 0; 50 100 1];\r\ntform = maketform('affine'<\/span>, T);<\/pre>
Apply the translation:<\/p>
I2 = imtransform(I, tform);<\/pre>
And compare the two images:<\/p>
subplot(1,2,1)\r\nimshow(I)\r\ntitle('original image'<\/span>)\r\n\r\nsubplot(1,2,2)\r\nimshow(I2)\r\ntitle('translated image'<\/span>)<\/pre> 
They're the same!  Hence, the confusion.<\/p>\r\n   
Here imtransform<\/tt> has been a little \"smarter\" than users want.  By automatically adjusting the output-space grid to capture the output image\r\n      wherever it moves, the translation gets removed.\r\n   <\/p>\r\n   
There are a couple of things you can do to make this work better.  First, you can ask imtransform<\/tt> to give you more information about the output space coordinates of the output image, and then you can use this information\r\n      when displaying the image.  Here's how that might work:\r\n   <\/p>
[I2,xdata,ydata] = imtransform(I,tform);<\/pre>
xdata<\/tt> is a two-element vector specifying where the first and last columns of I2<\/tt> are located in output space.  Similarly, ydata<\/tt> specifies where the first and last rows are located.  You can pass this information along to the display function, imshow<\/tt>.\r\n   <\/p>
subplot(1,2,1)\r\nimshow(I)\r\naxis on<\/span>\r\naxis([0 400 0 400])\r\n\r\nsubplot(1,2,2)\r\nimshow(I2,'XData'<\/span>,xdata,'YData'<\/span>,ydata)\r\naxis on<\/span>\r\naxis([0 400 0 400])<\/pre> 
A second technique is to tell imtransform<\/tt> what output-space rectangle to use.\r\n   <\/p>
I3 = imtransform(I,tform,'XData'<\/span>,[1 290],'YData'<\/span>,[1 391]);\r\nclf\r\nimshow(I3)<\/pre> 
If you understand the XData and YData parameters, you can use these techniques separately or in combination to implement whatever\r\n      translation effect you'd like.\r\n   <\/p>\r\n   
Since translation is a regular source of calls to tech support, we are exploring ways to alleviate the confusion.  If you\r\n      have suggestions, please leave a blog comment.  Thanks!\r\n   <\/p>\r\n