[Insight-users] possible problem with warper in ITK 3.16

[Luis Ibanez]

Hi Stefan,

Thanks for pointing this out.

The source of the problem is that ITK 3.16 is more formal
about managing the image orientation.

What you have found is that the WarpImageFilter, now
requires you to set the Direction of the output image.


This is done by adding the call:

Index: DeformableRegistration2.cxx
===================================================================
RCS file: /cvsroot/Insight/Insight/Examples/Registration/DeformableRegistration2.cxx,v
retrieving revision 1.35
diff -r1.35 DeformableRegistration2.cxx
342a343
>   warper->SetOutputDirection( fixedImage->GetDirection() );

To the configuration of the WarpImageFilter.



We just committed a fix to the CVS repository,
for the files in the directory:

            Examples/Registration.

They are:

http://public.kitware.com/cgi-bin/viewcvs.cgi/Examples/Registration/DeformableRegistration5.cxx?root=Insight&r1=1.11&r2=1.12&sortby=date
http://public.kitware.com/cgi-bin/viewcvs.cgi/Examples/Registration/LandmarkWarping2.cxx?root=Insight&r1=1.13&r2=1.14&sortby=date
http://public.kitware.com/cgi-bin/viewcvs.cgi/Examples/Registration/DeformableRegistration16.cxx?root=Insight&r1=1.2&r2=1.3&sortby=date
http://public.kitware.com/cgi-bin/viewcvs.cgi/Examples/Registration/DeformableRegistration10.cxx?root=Insight&r1=1.6&r2=1.7&sortby=date
http://public.kitware.com/cgi-bin/viewcvs.cgi/Examples/Registration/DeformableRegistration17.cxx?root=Insight&r1=1.2&r2=1.3&sortby=date
http://public.kitware.com/cgi-bin/viewcvs.cgi/Examples/Registration/DeformableRegistration2.cxx?root=Insight&r1=1.35&r2=1.36&sortby=date
http://public.kitware.com/cgi-bin/viewcvs.cgi/Examples/Registration/DeformableRegistration9.cxx?root=Insight&r1=1.7&r2=1.8&sortby=date
http://public.kitware.com/cgi-bin/viewcvs.cgi/Examples/Registration/DeformationFieldInitialization.cxx?root=Insight&r1=1.10&r2=1.11&sortby=date
http://public.kitware.com/cgi-bin/viewcvs.cgi/Examples/Registration/DeformableRegistration3.cxx?root=Insight&r1=1.18&r2=1.19&sortby=date


If you have a chance,
please give it a try and let us know
if you see any problems.


     Thanks


           Luis


--------------------------------------------
On Wed, Nov 11, 2009 at 3:25 AM, Stefan Bauer <cactux at web.de> wrote:
> Hello,
>
> I m wondering if there s a problem with the warper in ITK 3.16.
> I changed the deformableregistration2.cxx example to work with 3d data. (I changed line 67, 68, 69 and 119 to 3 instead of 2).
> If I run it on ITK 3.16 with the SRI24 anatomy spgr.hdr dataset (from http://www.nitrc.org/frs/?group_id=214&release_id=453  -> sri24_anatomy_analyze.zip) to register the data to itself, the output after warping is completely black and all the pixels are assigned zero values. However the deformation field looks okay (zero, as it should be). So the problem must be due to the warping.
>
> I compiled the exactly same code with ITK 3.8 and the deformed output image was alright.
>
> Thanks for your help
> Stefan
>
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:
> http://www.kitware.com/products/protraining.html
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users
>