<div dir="ltr">Hi Gavin,<div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br>
<div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px">1. I have checked out the </span><span style="font-family:arial,sans-serif;font-size:13px">CompositeTransform code, however, I can not find where the problem is treated. I am trying to understand how deformation at out of bound locations are treated. My current solution is:</span></div>
<div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px"> T2(T1(x)) = T2(T1(x)), if T1(x) within dimension of T2</span></div>
<div><span style="font-family:arial,sans-serif;font-size:13px"> T2(T1(x)) = T1(x), if T1(x) is out of dimension of T2</span></div><div><br></div><div>Assuming the dimension of T2 is the image size. My current problem is that discontinuity happens at locations where <span style="font-family:arial,sans-serif;font-size:13px">T1(x) out of dimension in the final composited deformation field.</span></div>
</div></blockquote><div><br></div><div>This is a reasonable solution. Note that T1(x) is not actually being computed in continuous space, but a discrete sampled space, and an interpolator is used on T2. A slight improvement is then to use an extrapolator for T2 when T1(x) is outside T2. Depending on the situation, throwing an exception and not counting the point may be most appropriate.</div>
<div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">
<div><br></div><div>2. Yes, I am comparing the displacement field directly. My problem here is also with the out of bound displacements. |T1-T2| will be large if there is some locations in T2 has large out of bound displacements. My question is whether these out of bound displacement should be reduced in their contribution to the displacement error. Because the out of bound displacement does not affect the re-sampled image.</div>
</div></blockquote><div><br></div><div>If comparing two displacement fields that are defined on the same input domain, it should not matter if any of the displacements end at location that is out of bounds of the input domain.</div>
<div><br></div><div>Hope this helps,</div><div>Matt</div><div><br></div><div><br></div></div></div></div>