<div dir="ltr">Hi Jacky,<div><br></div><div>Can you attach the code that you use to read the image, to create the transform, and to apply the transform?</div><div><br></div><div>Also note that, for vtk-dicom, there is a method vtkDICOMReader::SetMemoryRowOrderToFileNative() that will ensure that the layout of pixels in memory is identical to the layout of pixels in the file.  And the class vtkDICOMToRAS allows conversion of the image to the NIfTI coordinate system, though I'm guessing that you want to work in DICOM coordinates.</div><div><br></div><div> - David</div><div> <br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 17, 2015 at 3:44 AM, Jacky Ko <span dir="ltr"><<a href="mailto:jackkykokoko@gmail.com" target="_blank">jackkykokoko@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello,<div><br></div><div>I am using  vtkDICOMImageReader  to  read  dicom  series.  I  need to  export  the  rendered  surface  in stl format  for  further  analysis.  However,  I  found  that the exported  surface 's  orientation  is  different  from the  dicom  one  (e.g. loading both  dicom  series and  STL surfaces in 3D slicer).  But  the  case  doesn't  happen in Paraview. </div><div><br></div><div>I noticed that  there are two options for  extracting scanning  positions and  orientations by  GetImagePositionPatient() and  GetImageOrientationPatient() under  vtkDICOMImageReader class.</div><div><br></div><div>I also  read  the  code  in  vtk-dicom  as an extension.  However,  some problem occurs  when I am  using the package (I am developing in python).  I need to write  my  own transformation matrix and apply  to the image volume to get a spatial corrected surfaces. The problem is ,  how to get  the correct transformation matrix ?</div><div><br></div><div>Thanks.</div></div>
<br></blockquote></div><br></div></div></div>