<div dir="ltr">Marie,<div><br></div><div>Your observation is correct.</div><div><br></div><div>There is a bit of history on this.</div><div><br></div><div>The short story is that when some of the ImageToImage metrics were refactored</div>
<div>to support multi-threading, the itkNormalizedCorrelationImageToImageMetric,</div><div>was not refactored along with them.</div><div><br></div><div>Therefore, it doesn't uses some of the API methods that were added at the time</div>
<div>to the base class, the itkImageToImageMetric. In particular, it does not take into</div><div>account the collection of pixel samples that are generated in the base class.</div><div>(which includes the options for ignoring pixels based on a threshold).</div>
<div><br></div><div><br></div><div>  HTH</div><div><br></div><div><br></div><div>     Luis</div><div><br></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Jan 13, 2014 at 12:17 PM, Cherel, Marie C <span dir="ltr"><<a href="mailto:cherel@email.unc.edu" target="_blank">cherel@email.unc.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">Hi everyone,
<br>
<br>
I am writing because I have some issues with the class itkNormalizedCorrelationImageToImageMetric.<br>
I use this class in order to compute locally the normalized cross correlation between two brain MRIs that have been registered. 
<br>
<br>
So for each voxel I define a small neighborhood around this voxel (with a radius of 2 or 3 pixels for example) and then this region is set to the metric via the method SetFixedImageRegion. Then I get the value of the metric for a transform set to identity.
<br>
<br>
But after I tried to use the method SetFixedImageSamplesIntensityThreshold in order to avoid the background and I did not get any difference. So I tried to do the intensity check by my self and to use the method SetFixedImageIndexes but the metric was always
 equal to 0. <br>
<br>
So I looked into the source code of itkImageToImageMetric and itkNormalizedCorrelationImageToImageMetric. It seems that the parent class creates and fills a FixedImageIndexContainer which contains only the pixels that satisfy the threshold limitation, but this
 variable is not use in the itkNormalizedCorrelationImageToImageMetric class. itkNormalizedCorrelationImageToImageMetric only uses the fixedImageRegion and the masks, so it is what I used to figured out this problem. I also skip the first index selection which
 is useless by setting SetNumberOfFixedImageSamples to 0.<br>
<br>
Now the results seems to be fine, but please let me know if I misunderstood something.
<br>
Many thanks for considering my request, <br>
<br>
Marie<br>
</div>
</div>

<br>_______________________________________________<br>
Community mailing list<br>
<a href="mailto:Community@itk.org">Community@itk.org</a><br>
<a href="http://public.kitware.com/cgi-bin/mailman/listinfo/community" target="_blank">http://public.kitware.com/cgi-bin/mailman/listinfo/community</a><br>
<br></blockquote></div><br></div>