[Insight-users] ResamplingFilter problem

[Ariel Hernán Curiale]

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/4a534751/attachment-0001.htm>