[Insight-users] NCC in itk

[Dan Mueller]

I grepped the testing folder for ImageKernelOperator and had only one hit:
  D:\InsightToolkit-3.16.0\Testing\Code\Review\itkReviewPrintTest.cxx(40):#include
"itkImageKernelOperator.h"
  D:\InsightToolkit-3.16.0\Testing\Code\Review\itkReviewPrintTest.cxx(147):
 itk::ImageKernelOperator<float> kernelOperator;
  D:\InsightToolkit-3.16.0\Testing\Code\Review\itkReviewPrintTest.cxx(148):
 std::cout << "--------------ImageKernelOperatorObj ";

"If it's not tested, it's broken".

Cheers, Dan

On 10 November 2010 18:10, David Doria <daviddoria at gmail.com> wrote:
> On Wed, Nov 10, 2010 at 10:49 AM, somi <seesomi at gmail.com> wrote:
>> Hi,
>> I am trying to use itk::NormalizedCorrelationImageFilter , with template as
>> another image.
>> To set up another image as as templete, I use   itk::ImageKernelOperator
>>   itk::ImageKernelOperator<float, Dimension> kernelOperator;
>>     kernelOperator.SetImageKernel(readerTemplate->GetOutput());
>>     kernelOperator.CreateToRadius(radius);
>>
>> typedef itk::NormalizedCorrelationImageFilter<InputImageType,
>> InputImageType, InputImageType> NCCFilterType;
>>     NCCFilterType::Pointer filter = NCCFilterType::New();
>>     filter->SetInput(reader->GetOutput());
>>     filter->SetTemplate(kernelOperator);
>>
>> However, ther output of NCC is all infinity. I found all the coefficients
>> of   itk::ImageKernelOperator were zero.
>> Is there any step I am missing ? My test code and inputs are attached.
>> Usage:
>> NCC main.png templete.png out.nii
>>
>> Thanks,
>> Somi
>
> I also have had no luck with this operation:
> http://www.itk.org/Wiki/ITK/Examples/Broken/Images/NormalizedCorrelationImageFilter
>
> Has anyone done this successfully?
>
> David
> _____________________________________
> 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
>