<div dir="ltr">Hi all,<div><br></div><div>I would like to have some advise to do the following task.</div><div><br></div><div>I have a set of annotated DICOM (binary 0 - 255 intensity) with a given spacing, orientation and position. This manual segmentation is corresponding to a prostate from an initial T2W sequence.<br clear="all">
<div><br></div><div>A part of DICOM, I have an MRSI sequence from Siemens (*.rda) in which the header is containing information such as position, orientation and spacing of the csi grid matrix used during the acquisition.</div>
<div><br></div><div>My aim would be to resample the annotated DICOM serie into a new DICOM serie using the MRSI header. I attached the code that I tried to implement to do this task but I don't succeed to obtain a correct serie when saving into a DICOM serie. When opening the stored DICOM serie, it seems that the field regarding the z spacing and the position of the first voxel are not equal to the one of the rda manually affected during the resampling.</div>
<div><br></div><div>I was wondering to know the best way to do such manipulation. For instance, it could be better to create a manual image using the field of the rda file and just use the intensities</div><div><br></div>
<div>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(128,0,128)">ResampleFilterType</span><span style="color:rgb(0,0,0)">::</span><span style="color:rgb(128,0,128)">Pointer</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">resampler</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,0,128)">ResampleFilterType</span><span style="color:rgb(0,0,0)">::</span><span style="color:rgb(0,0,0)">New</span><span style="color:rgb(0,0,0)">();</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">resampler</span><span style="color:rgb(0,0,0)">-></span><span style="font-style:italic;color:rgb(0,0,0)">SetInput</span><span style="color:rgb(0,0,0)">(</span><span style="color:rgb(0,0,0)">inputAnnotatedVolume</span><span style="color:rgb(0,0,0)">-></span><span style="color:rgb(0,0,0)">GetOutput</span><span style="color:rgb(0,0,0)">());</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">resampler</span><span style="color:rgb(0,0,0)">-></span><span style="font-style:italic;color:rgb(0,0,0)">SetTransform</span><span style="color:rgb(0,0,0)">(</span><span style="color:rgb(0,0,0)">transform</span><span style="color:rgb(0,0,0)">);</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">resampler</span><span style="color:rgb(0,0,0)">-></span><span style="font-style:italic;color:rgb(0,0,0)">SetInterpolator</span><span style="color:rgb(0,0,0)">(</span><span style="color:rgb(0,0,0)">interpolator</span><span style="color:rgb(0,0,0)">);</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">Create</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">a</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">type</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">compatible</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">to</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">ITK</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">Define</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">the</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">new</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">origin</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">double</span><span style="color:rgb(0,0,0)">*</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputOrigin</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">new</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">double</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getPositionPatient</span><span style="color:rgb(0,0,0)">().</span><span style="color:rgb(0,0,0)">size</span><span style="color:rgb(0,0,0)">()];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">for</span><span style="color:rgb(0,0,0)">(</span><span style="color:rgb(128,128,0)">int</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">;</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)"><</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getPositionPatient</span><span style="color:rgb(0,0,0)">().</span><span style="color:rgb(0,0,0)">size</span><span style="color:rgb(0,0,0)">();</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(0,0,0)">++)</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputOrigin</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getPositionPatient</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">resampler</span><span style="color:rgb(0,0,0)">-></span><span style="font-style:italic;color:rgb(0,0,0)">SetOutputOrigin</span><span style="color:rgb(0,0,0)">(</span><span style="color:rgb(0,0,0)">outputOrigin</span><span style="color:rgb(0,0,0)">);</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">Define</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">the</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">new</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">spacing</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">double</span><span style="color:rgb(0,0,0)">*</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputSpacing</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">new</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">double</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getPixelSpacing</span><span style="color:rgb(0,0,0)">().</span><span style="color:rgb(0,0,0)">size</span><span style="color:rgb(0,0,0)">()];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">for</span><span style="color:rgb(0,0,0)">(</span><span style="color:rgb(128,128,0)">int</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">;</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)"><</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getPixelSpacing</span><span style="color:rgb(0,0,0)">().</span><span style="color:rgb(0,0,0)">size</span><span style="color:rgb(0,0,0)">();</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(0,0,0)">++)</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputSpacing</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getPixelSpacing</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">resampler</span><span style="color:rgb(0,0,0)">-></span><span style="font-style:italic;color:rgb(0,0,0)">SetOutputSpacing</span><span style="color:rgb(0,0,0)">(</span><span style="color:rgb(0,0,0)">outputSpacing</span><span style="color:rgb(0,0,0)">);</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">Define</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">the</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">new</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">orientation</span><span style="color:rgb(192,192,192)"> </span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">typename</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,0,128)">ResampleFilterType</span><span style="color:rgb(0,0,0)">::</span><span style="color:rgb(128,0,128)">DirectionType</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputDirection</span><span style="color:rgb(0,0,0)">;</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//double</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">**</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">outputDirection</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">new</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">double*</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">[3];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//for(int</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">i</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">0;</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">i</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)"><</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">3</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">;</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">i++)</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">outputDirection[i]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">new</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">double</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">[3];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">Affect</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">the</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">new</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">values</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputDirection</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getOrientationPatient</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputDirection</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getOrientationPatient</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputDirection</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getVoiNormalVector</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputDirection</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getOrientationPatient</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputDirection</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getOrientationPatient</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputDirection</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getVoiNormalVector</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputDirection</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getOrientationPatient</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputDirection</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getOrientationPatient</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputDirection</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">][</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getVoiNormalVector</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">];</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">resampler</span><span style="color:rgb(0,0,0)">-></span><span style="font-style:italic;color:rgb(0,0,0)">SetOutputDirection</span><span style="color:rgb(0,0,0)">(</span><span style="color:rgb(0,0,0)">inputAnnotatedVolume</span><span style="color:rgb(0,0,0)">-></span><span style="color:rgb(0,0,0)">GetOutput</span><span style="color:rgb(0,0,0)">()-></span>GetDirection<span style="color:rgb(0,0,0)">());</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">Define</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">the</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">size</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">of</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">the</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">ouput</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">Let's</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">comput</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">the</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">new</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">size</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,0,128)">InputImageType</span><span style="color:rgb(0,0,0)">::</span><span style="color:rgb(128,0,128)">SizeType</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputSize</span><span style="color:rgb(0,0,0)">;</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">typedef</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,0,128)">InputImageType</span><span style="color:rgb(0,0,0)">::</span><span style="color:rgb(128,0,128)">SizeType</span><span style="color:rgb(0,0,0)">::</span><span style="color:rgb(128,0,128)">SizeValueType</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,0,128)">SizeValueType</span><span style="color:rgb(0,0,0)">;</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputSize</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">static_cast</span><span style="color:rgb(0,0,0)"><</span><span style="color:rgb(128,0,128)">SizeValueType</span><span style="color:rgb(0,0,0)">>(</span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getCsiMatrixSize</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">0</span><span style="color:rgb(0,0,0)">]);</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputSize</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">static_cast</span><span style="color:rgb(0,0,0)"><</span><span style="color:rgb(128,0,128)">SizeValueType</span><span style="color:rgb(0,0,0)">>(</span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getCsiMatrixSize</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">1</span><span style="color:rgb(0,0,0)">]);</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">outputSize</span><span style="color:rgb(0,0,0)">[</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">]</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(128,128,0)">static_cast</span><span style="color:rgb(0,0,0)"><</span><span style="color:rgb(128,0,128)">SizeValueType</span><span style="color:rgb(0,0,0)">>(</span><span style="color:rgb(0,0,0)">rdaPatient</span><span style="color:rgb(0,0,0)">.</span><span style="color:rgb(0,0,0)">getCsiMatrixSize</span><span style="color:rgb(0,0,0)">()[</span><span style="color:rgb(0,0,128)">2</span><span style="color:rgb(0,0,0)">]);</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">resampler</span><span style="color:rgb(0,0,0)">-></span><span style="font-style:italic;color:rgb(0,0,0)">SetSize</span><span style="color:rgb(0,0,0)">(</span><span style="color:rgb(0,0,0)">outputSize</span><span style="color:rgb(0,0,0)">);</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">//</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">Resampling</span><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,128,0)">now</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(192,192,192)"> </span><span style="color:rgb(0,0,0)">resampler</span><span style="color:rgb(0,0,0)">-></span><span style="font-style:italic;color:rgb(0,0,0)">Update</span><span style="color:rgb(0,0,0)">();</span></pre>
<pre style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(0,0,0)"><br></span></pre><pre style="margin-top:0px;margin-bottom:0px">Thanks in advance for the advices,</pre></div>-- <br><div dir="ltr"><div><span style="font-family:'Times New Roman';font-size:medium"><table style="width:638px" border="0" cellpadding="0">
<tbody><tr style="height:37.5pt"><td style="width:228px;height:37.5pt;padding:0.75pt" align="left" valign="top"><b><span style="font-family:Arial,Helvetica;font-size:8pt">LEMAÎTRE Guillaume</span><span style="font-size:8pt;font-family:Helvetica"><br>
</span><font size="1"><span style="font-family:Arial">PhD Candiate<br>MSc Erasmus Mundus ViBOT (Vision-roBOTic)<br></span></font><span style="font-size:7.5pt;font-family:Arial" lang="EN-GB"><font color="#000000">MSc Business Innovation and Technology Management (in progress)<br>
</font></span></b><span style="font-size:7.5pt;font-family:Arial" lang="EN-GB"><b><br></b><font color="#0DB02B"><a href="mailto:g.lemaitre58@gmail.com" target="_blank">g.lemaitre58@gmail.com</a><br><br></font></span><p></p>
</td><td style="width:400px;height:37.5pt;padding:0.75pt" align="left" valign="top"><span style="font-size:8pt"><b><span style="font-family:Arial,Helvetica"><font color="#000000">ViCOROB - Computer Vision and Robotic Team</font></span></b></span><span style="font-size:7.5pt;font-family:Arial"><font color="#000000"><br>
Universitat de Girona, Campus Montilivi, Edifici P-IV 17071 Girona<br>Tel. +34 972 41 98 12 - Fax. +34 972 41 82 59 </font></span><br><span lang="EN-GB"><a href="http://vicorob.udg.es/" style="font-family:Arial;font-size:7.5pt" target="_blank">http://vicorob.udg.es/</a><br>
<font color="#000000"><b><font face="Arial, Helvetica"><span style="font-size:11px">LE2I - Le Creusot</span></font><br></b><font face="Arial" size="1">IUT Le Creusot, Laboratoire LE2I, 12 rue de la Fonderie, 71200 Le Creusot<br>
Tel. +33 3 85 73 10 90</font></font><span style="font-family:Arial;font-size:x-small"> - Fax. +33 3 85 73 10 97</span><span style="font-family:Arial;font-size:x-small"> </span><span style="font-family:Arial;font-size:x-small"><br>
</span><span lang="EN-GB" style="font-family:Arial;font-size:7.5pt"><a href="http://le2i.cnrs.fr" target="_blank">http://le2i.cnrs.fr</a></span></span></td></tr></tbody></table></span></div><div style="text-align:left"><img src="https://sites.google.com/site/glemaitre58/_/rsrc/1340103962984/config/le2i.png"> <img src="https://sites.google.com/site/glemaitre58/_/rsrc/1340103388011/config/vico.png"> <img src="https://sites.google.com/site/glemaitre58/_/rsrc/1340103887254/config/ub.png"> <img src="https://sites.google.com/site/glemaitre58/_/rsrc/1340103809482/config/udg.png"> <img src="https://sites.google.com/site/glemaitre58/_/rsrc/1340104225210/config/vibot.png"></div>
<a href="https://sites.google.com/site/glemaitre58/" target="_blank">https://sites.google.com/site/glemaitre58/</a><br>Vice - Chairman of A.S.C. Fours UFOLEP<br>Chairman of A.S.C. Fours FFC<br>Webmaster of <a href="http://ascfours.free.fr" target="_blank">http://ascfours.free.fr</a></div>
</div></div>