<div dir="ltr"><div class="gmail_quote"><div dir="ltr"><div class="gmail_quote"><div dir="ltr"><div><span style="font-family:monospace,monospace">Dear all:</span><br></div><font face="monospace, monospace"><br>I would like to transforming the tracks to the ROI space which image valve is integer from 0 to 490,<br>but not sure how to apply the transformation matrix.</font><div><font face="monospace, monospace">So I use MR::Image::Interp to access the voxel intensity of ROI image.<br><br>Here is how I use,<br>===<br></font><div><span style="font-family:monospace,monospace">std::vector<MR::Point> tck; <span style="background-color:rgb(255,255,255)"><font color="#6aa84f">// one track got from tck file</font></span></span></div><div><font face="monospace, monospace">MR::Image::Object image_object;<br>image_object.open(ROI.nii); <font color="#6aa84f">// which is coregistered to DWI</font></font></div><div><div><font face="monospace, monospace">MR::Image::Interp interp (image_object);</font></div><div><span style="font-family:monospace,monospace">MR::Point track_point, </span><span style="font-family:monospace,monospace">voxel_position;</span><font face="monospace, monospace"><br></font><div><font face="monospace, monospace"><div>for (std::vector< MR::Point>::iterator i = tck.begin(); i != tck.end(); ++i)</div><div> {</div><div>  float *p = i->get();</div><div>  track_point.set(p);</div><div>  voxel_position = interp.R2P(track_point);</div><div>  interp.P (voxel_position);</div><div>  value = interp.value();</div><div>  if(value != 0)</div><div>  {</div><div>   cout << value << endl ; </div><div>  }</div><div> }</div><div>===</div>Is this a correct method to using Interp? </font><span style="font-family:monospace,monospace">Cause I get a float value.</span></div><div><font face="monospace, monospace">I hope I can get a integer value show that fiber is passing through the specific region.</font></div><div><span style="font-family:monospace,monospace"><br></span></div><div><font face="monospace, monospace">Should I round down, up or nearest to the point of track(track_point) to access right voxel intensity?</font></div><div><font face="monospace, monospace"><br></font></div><div><span style="font-family:monospace,monospace">Thank you for your help!</span><br></div></div></div></div><span class="HOEnZb"><font color="#888888"><span><div><span style="color:rgb(0,0,0);font-family:monospace,monospace">Kuo, Kuan-Tsen</span><br></div></span></font></span></div>
</div></div></div><br></div>