Hi Kishore,<br>This is indeed very helpful.  Thanks for your time.<br><br>Regards,<br>Seth<br><br><div class="gmail_quote">On Fri, Sep 24, 2010 at 10:00 AM, Kishore Mosaliganti <span dir="ltr">&lt;<a href="mailto:kishoreraom@gmail.com">kishoreraom@gmail.com</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Is this class relevant to your problem?<br>
<br>
<a href="http://www.itk.org/Doxygen318/html/classitk_1_1ZeroFluxNeumannBoundaryCondition.html" target="_blank">http://www.itk.org/Doxygen318/html/classitk_1_1ZeroFluxNeumannBoundaryCondition.html</a><br>
<font color="#888888"><br>
Kishore<br>
</font><div class="im"><br>
On Thu, Sep 23, 2010 at 5:42 PM, Seth Gilchrist &lt;<a href="mailto:seth@mech.ubc.ca">seth@mech.ubc.ca</a>&gt; wrote:<br>
</div><div><div></div><div class="h5">&gt; Hello,<br>
&gt; I know what you are thinking...why are you up-sampling an image?  Here is<br>
&gt; the situation and problem and I hope that you can help.<br>
&gt;<br>
&gt; I have an image of vectors (3D image, 3D vectors) that is of a course<br>
&gt; resolution (35,20,20).  These vectors represent the displacement of a<br>
&gt; grouping of pixels in a fixed image. Each voxel in the vector image<br>
&gt; represents the average displacement of 21x21x21 voxels in the fixed image.<br>
&gt; What I want to do is linearly interpolate the data from the course vector<br>
&gt; image into the grid (they occupy the same space) of the fixed image.  The<br>
&gt; new vector image will then be used to morph the fixed image into the moving<br>
&gt; image space, giving a deformable-like transform, and then the fixed and<br>
&gt; moving image can be operated on in the same space.<br>
&gt;<br>
&gt; The problem is that when interpolating the values at the edge of the vector<br>
&gt; image the interpolator uses a single value for the out-of-bounds pixels.<br>
&gt; This means that the vectors at the edge of my new vector image converge to<br>
&gt; the displacement represented by the out-of-bounds value (a constant).  What<br>
&gt; I would like, is to have the edge given a zero derivative (Neumann)<br>
&gt; condition (or even better, be able to try a few different boundary<br>
&gt; conditions).  I have looked through the interpolate and resample filters and<br>
&gt; can&#39;t seem to find this already implemented.  Has this been done, or should<br>
&gt; I just get on writing it myself?<br>
&gt;<br>
&gt; Thanks for your time and help.<br>
&gt;<br>
&gt; Regards,<br>
&gt; Seth<br>
&gt;<br>
</div></div><div><div></div><div class="h5">&gt; _____________________________________<br>
&gt; Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
&gt;<br>
&gt; Visit other Kitware open-source projects at<br>
&gt; <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
&gt;<br>
&gt; Kitware offers ITK Training Courses, for more information visit:<br>
&gt; <a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
&gt;<br>
&gt; Please keep messages on-topic and check the ITK FAQ at:<br>
&gt; <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
&gt;<br>
&gt; Follow this link to subscribe/unsubscribe:<br>
&gt; <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
&gt;<br>
&gt;<br>
</div></div></blockquote></div><br>