<div dir="ltr"><div>Hi Zhang,</div><div><br></div><div>The links in David Cole's email get to the heart of the issue, so please read them.  I think that we are using different definitions of "spacing" and "thickness".  The reason I think this, is because of your statement "<span style="color:rgb(80,0,80);font-size:12.8px">obtain the image whose thickness is 0.4 from the image whose </span><span style="color:rgb(80,0,80);font-size:12.8px">thickness is 0.8", which makes no sense to me.  Let us start with my definitions of "spacing" and "thickness":</span></div><div><br></div><div>The "spacing" is sample spacing, e.g. the distance from the center of one slice to the center of the next slice.</div><div><br></div><div>The "thickness" refers to the collimated X-ray beam width in CT, or to the selection slice thickness in 2D MRI.  Usually "thickness" is only approximate, because the slice profile is not uniform.  Essentially, the "slice profile" is the point-spread function in the slice direction, and "slice thickness" is a measure of the width of the slice profile.</div><div><br></div><div>Do you agree with these definitions?</div><div><br></div><div>So if a CT is acquired with a slice thickness of 0.8 mm, it is impossible to convert it to 0.4 mm.  that would break the laws of physics, so to speak.  It is possible to increase the slice thickness by averaging slices together or by applying other filtering operations, but decreasing the thickness is not possible.</div><div><br></div><div>The operations that vtkImageReslice can do are: 1) increasing or decreasing the spacing (by resampling with interpolation), and 2) increasing the thickness (by averaging, or by maximum intensity projection).</div><div><br></div><div> - David</div><div><br></div><br><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Aug 2, 2017 at 8:48 AM, David Cole via vtkusers <span dir="ltr"><<a href="mailto:vtkusers@vtk.org" target="_blank">vtkusers@vtk.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">vtkImageData is a regularly spaced collection of *point* samples,<br>
there is no "thickness"...<br>
<br>
Please refer to the classic "A Pixel Is Not A Little Square" :<br>
<a href="http://alvyray.com/Memos/CG/Microsoft/6_pixel.pdf" rel="noreferrer" target="_blank">http://alvyray.com/Memos/CG/<wbr>Microsoft/6_pixel.pdf</a><br>
<br>
Also, a shorter summary here:<br>
<a href="http://graphics.cs.wisc.edu/WP/cs559-fall2014/2014/08/29/what-is-a-pixel-and-what-is-a-point-sample/" rel="noreferrer" target="_blank">http://graphics.cs.wisc.edu/<wbr>WP/cs559-fall2014/2014/08/29/<wbr>what-is-a-pixel-and-what-is-a-<wbr>point-sample/</a><br>
<br>
<br>
HTH,<br>
David C.<br>
<div class="gmail-HOEnZb"><div class="gmail-h5"><br>
<br>
<br>
On Wed, Aug 2, 2017 at 10:39 AM, zhang qiang <<a href="mailto:zhang07101@126.com">zhang07101@126.com</a>> wrote:<br>
> Dear David:<br>
><br>
>      Thanks a lot. And the vtkImageData also do not record the thickness<br>
> information. I wonder why VTK ignore the thickness information? I think the<br>
> thickness information is quite important for DICOM image. My guess is that<br>
> we can not obtain the image whose thickness is 0.4 from the image whose<br>
> thickness is 0.8. So,  the thickness information is useless in vtkImageData?<br>
> Am I right?<br>
><br>
> zhang qiang<br></div></div></blockquote></div></div></div>