<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">3D Slicer uses RAS coordinate system, unlike DICOM which uses LPS. Thus x and y coordinates are inverted. Slicer takes care of this inversion for volume images, but not for polygonal meshes. In other words, you cannot have a mesh which appears correctly in both Slicer and ParaView.</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">HTH</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 17, 2015 at 7:17 AM, David Gobbi <span dir="ltr"><<a href="mailto:david.gobbi@gmail.com" target="_blank">david.gobbi@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">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><span class="HOEnZb"><font color="#888888"><div><br></div><div> - David</div></font></span><div><div class="h5"><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></div></div>
<br>_______________________________________________<br>
Powered by <a href="http://www.kitware.com" rel="noreferrer" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the VTK FAQ at: <a href="http://www.vtk.org/Wiki/VTK_FAQ" rel="noreferrer" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=vtkusers" rel="noreferrer" target="_blank">http://markmail.org/search/?q=vtkusers</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/vtkusers" rel="noreferrer" target="_blank">http://public.kitware.com/mailman/listinfo/vtkusers</a><br>
<br></blockquote></div><br></div>