Hello,<div><br></div><div>I have a little problem that I think is just down to some small detail in my understanding of DICOM and ITK. I want avoid resampling an image when applying a rigid transform by applying the transform to the directional cosines (C) and the position vector (V) of the image instead. I do this by first getting the rotation matrix (M) and the offset (O) from the transform (an Euler3DTransform) and setting the new directional cosines as the matrix multiplication of the rotation matrix and the directional cosines, </div>
<div><br></div><div>C2=M*C, </div><div><br></div><div>and the new position vector as </div><div><br></div><div>V2 = M*V + O.</div><div><br></div><div>I noticed that in order to get correct results I must subtract the offset O, so </div>
<div><br></div><div>V2 = M*V - O.</div><div><br></div><div>I check the result by resampling my new image using an identity transform and my fixed image as reference and then take the difference between them. The result is slightly worse than if I just resample my original moving image with the resulting transform of the registration. I use bsplines = 3 when resampling in both cases. It seems that my way of applying the transform is off by a voxel or so. </div>
<div><br></div><div>What am I doing wrong? Since I wanted to add the offset instead of subtracting it I guess I don't fully understand the parameters I'm working with so I guess I'm making more mistakes. Can someone point them out to me? Maybe there is already a filter that does what I want to do?</div>
<div><br></div><div>Thanks in advance!</div><div><br></div><div>Regards</div><div>===================================
<div>| <span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px">Patrik Brynolfsson</span></div>
<div><span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px">| Ph.D. student at Dept of Radiation Physics</span></div>
<div><span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px"><span style="font-size:small">| </span></span><span style="border-collapse:collapse;font-family:arial, sans-serif">Umeå University<br>
| SE-901 87 Umeå</span><br></div>
<div>
<div><span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px">| Sweden</span></div>
<div><span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px">| Phone: +46 (0)90 785 2296</span></div></div>
<div><font face="arial, sans-serif"><span style="border-collapse:collapse"><span style="border-collapse:separate;font-family:arial">===================================</span></span></font></div><br>
</div>