At this point you undoubtedly have more than you need in answer to your original question, but here is one more suggestion. Instead of comparing the deformation fields directly you can apply both to the voxels of the image and evaluate the distance between the resulting positions. Note that for "dense" deformation fields this will be the same as comparing the deformation fields directly.<div>
<br></div><div>- Wes<br><br><div class="gmail_quote">On Tue, Oct 27, 2009 at 12:43 AM, Karthik Krishnan <span dir="ltr"><<a href="mailto:karthik.krishnan@kitware.com">karthik.krishnan@kitware.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">You could also try metrics to compare deformation fields.<br>
<br>
McCane, B., Novins, K., Crannitch, D. and Galvin, B., On Benchmarking<br>
Optical Flow, Computer Vision and Image Understanding, 84(1), pp<br>
126-143, 2001<br>
<br>
See Section 5 titled Error Metrics in [1]. It provides 2 metrics to<br>
compare deformation fields. It also surveys other metrics used to<br>
compare deformation fields in literature and illustrates their<br>
shortcomings.<br>
<br>
Their methods are openly available at<br>
<a href="http://sourceforge.net/projects/of-eval/files/" target="_blank">http://sourceforge.net/projects/of-eval/files/</a> .<br>
<br>
Motes : if you port them to ITK, and contribute it to IJ, I'd he happy<br>
to include them in the toolkit.<br>
<br>
Thanks<br>
--<br>
<font color="#888888">karthik<br>
</font><div><div></div><div class="h5"><br>
On Mon, Oct 26, 2009 at 1:32 PM, Luis Ibanez <<a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a>> wrote:<br>
> Hi Motes,<br>
><br>
> I would suggest that you subtract the Known deformation field<br>
> from the deformation field that results from the registration. This<br>
> will give you a "residual" deformation field. Then you could<br>
> compute the Average Magnitude of all the vectors in these<br>
> residual deformation field.<br>
><br>
> It may also be worth to report the "largest magnitude" of that<br>
> residual deformation field.<br>
><br>
> For subtracting two deformation fields, the itkSubtractImageFilter<br>
> should do the trick (although I have not tried this filter with images<br>
> of vectors).<br>
><br>
> The computation of the average and maximum magnitude may<br>
> require you to write a customized while-loop with Iterators.<br>
><br>
> It should take about 30~40 lines of code.... :-)<br>
><br>
><br>
> Please let us know if you find any problem,<br>
><br>
><br>
> Thanks<br>
><br>
><br>
> Luis<br>
><br>
><br>
> -----------------------------------------------------------------------------------------<br>
> On Mon, Oct 26, 2009 at 12:11 PM, motes motes <<a href="mailto:mort.motes@gmail.com">mort.motes@gmail.com</a>> wrote:<br>
>> I would like to measure the quality of my modfied BSplineDeformable<br>
>> transform. I have used one of the ITK examples to create a known<br>
>> deformation of the fixed image that I use as the moving image.<br>
>><br>
>> But how do I use this known deformation in the quality measure? Since<br>
>> the deformation field is applied to the fixed image and the result is<br>
>> the moving image I don't see what benefit it does when measuring the<br>
>> quality of the registration?<br>
>><br>
>> How are deformation fields used in this context?<br>
>> _____________________________________<br>
>> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
>><br>
>> Visit other Kitware open-source projects at<br>
>> <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
>><br>
>> Please keep messages on-topic and check the ITK FAQ at: <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
>><br>
>> Follow this link to subscribe/unsubscribe:<br>
>> <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
>><br>
> _____________________________________<br>
> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at<br>
> <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Kitware offers ITK Training Courses, for more information visit:<br>
> <a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
><br>
> Please keep messages on-topic and check the ITK FAQ at:<br>
> <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
><br>
_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Wesley D. Turner, Ph.D.<br>Kitware, Inc.<br>Technical Leader<br>28 Corporate Drive<br>Clifton Park, NY 12065-8662<br>Phone: 518-881-4920<br>
</div>