<div dir="ltr">Hey,<br><br>I need the linearly interpolated intensity values of some certain points in a 3d volume. I am trying to use LinearInterpolateImageFunction. The code is as follows.<br>There are no errors while compiling or linking, but it concts off before the Evaluate. Will appreciate any sort of help.<br>
<br>Thank you.<br><br>Manisha<br><br><br>&nbsp;&nbsp;&nbsp; typedef float InputPixelType;<br>&nbsp;&nbsp;&nbsp; const&nbsp;&nbsp; unsigned int&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Dimension = 3;<br>&nbsp;&nbsp;&nbsp; typedef itk::Image&lt; InputPixelType,&nbsp; Dimension &gt;&nbsp;&nbsp;&nbsp; InputImageType;<br>&nbsp;&nbsp;&nbsp; typedef itk::ImageFileReader&lt; InputImageType&nbsp; &gt;&nbsp; ReaderType;<br>
&nbsp;&nbsp;&nbsp; ReaderType::Pointer reader = ReaderType::New();<br>&nbsp;&nbsp;&nbsp; const char * inputFilename&nbsp; = argv[1];<br>&nbsp;&nbsp;&nbsp; reader-&gt;SetFileName( inputFilename&nbsp; );<br>&nbsp;&nbsp;&nbsp; typedef itk::LinearInterpolateImageFunction&lt;InputImageType, double&gt; InterpolationFunction;<br>
&nbsp;&nbsp;&nbsp; InterpolationFunction::Pointer interpolator = InterpolationFunction::New();<br><br>&nbsp;&nbsp;&nbsp; typedef itk::Point &lt;double, 3&gt; PointType;<br>&nbsp;&nbsp; PointType point;<br>&nbsp;&nbsp; point[0]=25;<br>&nbsp;&nbsp; point[1]=25;<br>&nbsp;&nbsp; point[2]=25;<br>
&nbsp;&nbsp; float I = interpolator-&gt;Evaluate(point);<br>&nbsp;&nbsp; cout &lt;&lt; I &lt;&lt; &quot;&nbsp; &quot;&nbsp; &lt;&lt; endl;<br></div>