[Insight-users] ResamplingFilter problem

[alex Dowson]

Hi

Try to produce small example that we can debug it and check where is problem

for small example see this

http://www.itk.org/Wiki/ITK/Examples/ImageProcessing/ResampleImageFilter





From: Ariel Hernán Curiale 
Sent: Tuesday, October 02, 2012 5:04 PM
To: alex Dowson 
Cc: Alessandro ; Insight User 
Subject: Re: [Insight-users] ResamplingFilter problem

I forgot to say that I use the interpolator only in the registration and the registration modified my transform which is a BSplineTransform. 

Saludos
__________________________________
| Ariel Hernán Curiale Ph.D Student
| ETSI Telecomunicación
| Universidad de Valladolid
| Campus Miguel Delibes
| 47011 Valladolid, Spain
| Phone: 983-423000 ext. 5590
| Web: www.curiale.com.ar
|_________________________________

El 02/10/2012, a las 13:25, Ariel Hernán Curiale escribió:


  Hi,  
  I'm using this code and works ok for me:

  typedef itk::InterpolateImageFunction< MovingImageType,double> InterpolatorType;
  typename InterpolatorType::Pointer   interpolator;

  .....
  .....

  interpolator = BsplineInterpolatorType::New(); //Generic interpolator type
  typedef   BsplineInterpolatorType *           InterpolatorPointer;
  InterpolatorPointer int_aux = dynamic_cast< InterpolatorPointer>(interpolator.GetPointer());
  int_aux->SetSplineOrder(3);


  typedef itk::ResampleImageFilter<FixedImageType,MovingImageType >ResampleFilterType;

  typename ResampleFilterType::Pointer resample_grid = ResampleFilterType::New();

  resampleGridImages.push_back(resample_grid);

  resample_grid->SetTransform( transform );
  resample_grid->SetInput( gridImage );

  resample_grid->SetSize(    fixedImage->GetLargestPossibleRegion().GetSize() );
  resample_grid->SetOutputOrigin(  fixedImage->GetOrigin() );
  resample_grid->SetOutputSpacing( fixedImage->GetSpacing() );
  resample_grid->SetOutputDirection( fixedImage->GetDirection() );
  resample_grid->SetDefaultPixelValue( 100 );
  resample_grid->Update();


  Saludos
  __________________________________
  | Ariel Hernán Curiale Ph.D Student
  | ETSI Telecomunicación
  | Universidad de Valladolid
  | Campus Miguel Delibes
  | 47011 Valladolid, Spain
  | Phone: 983-423000 ext. 5590
  | Web: www.curiale.com.ar
  |_________________________________

  El 02/10/2012, a las 12:47, alex Dowson escribió:




    Hi

    You didn’t set the Transformation anywhere. Can you show me code again with correction .


    -----Original Message----- From: Alessandro
    Sent: Tuesday, October 02, 2012 4:15 PM
    To: Insight User
    Subject: Re: [Insight-users] ResamplingFilter problem

    Hi,

    I tried but still nothing seems to come out...

    Ale

    Il giorno 02/ott/2012, alle ore 12.39, alex Dowson ha scritto:




      Hi





      Try calling the resample->Update() before the resample->GetOutput()









      -----Original Message----- From: Alessandro

      Sent: Tuesday, October 02, 2012 4:07 PM

      To: Insight User

      Subject: [Insight-users] ResamplingFilter problem



      Hello,



      I tried to Resample a 3d Slice Volume using a reference image. I wrote the following function (inside a class) to perform the task:



      ImageVolumeType::Pointer VolumeRetriever::applyResamplingFromImage (ImageVolumeType::Pointer img, ImageVolumeType::Pointer ref) {



      typedef itk::IdentityTransform<double, 3> TransformType;

      typedef itk::ResampleImageFilter<ImageVolumeType, ImageVolumeType> ResampleImageFilterType;



      ResampleImageFilterType::Pointer resample = ResampleImageFilterType::New();



      resample->SetInput(img);

      resample->UseReferenceImageOn();

      resample->SetReferenceImage(ref);

      resample->UpdateLargestPossibleRegion();



      return resample->GetOutput();

      }



      where ImageVolumeType is  typedef itk::Image< InputPixelType, 3 > ImageVolumeType;



      Problem is: when I run it the return is a NULL pointer...



      What am I missing?



      Thanks,

      Ale



      _____________________________________

      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.php



      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


    _____________________________________
    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.php

    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 
    _____________________________________
    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.php

    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


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20121002/20a5d528/attachment.htm>