<DIV>Hi all!</DIV>
<DIV> </DIV>
<DIV>I have to do an aplication , that must be able to register volumes having different orientations. </DIV>
<DIV> </DIV>
<DIV>For reading the series of DICOM images I tried both libraries : DICOMParser, gdcm <BR></DIV>
<DIV>Between the 2 libraries I observed the following difference regarding the computation of normal used to sort the files by Image Position Patient </DIV>
<DIV> </DIV>
<DIV>- in DICOMParser: file DICOMAppHelper.cxx ,functionGetImagePositionPatientFilenamePairs, line 2244 :<BR> normal[1] = ((*sn_iter).second.ImageOrientationPatient[0] <BR>*(*sn_iter).second.ImageOrientationPatient[5])<BR> - <BR>((*sn_iter).second.ImageOrientationPatient[2]* <BR>(*sn_iter).second.ImageOrientationPatient[3]);</DIV>
<DIV> </DIV>
<DIV>- in gdcm : file gdcmHeaderHelper.cxx, function ImagePositionPatientOrdering(), line 184 :<BR> normal[1] = cosines[2]*cosines[3] - cosines[0]*cosines[5];</DIV>
<DIV> </DIV>
<DIV>This difference influence the coordinate system of the resulting image volume.<BR>How is correct?</DIV>
<DIV> </DIV>
<DIV>I want to use itkOrientImageFilter to bring all the volume to a single coordinate system.This filter has to know the coordinate system of the input image. </DIV>
<DIV><BR>The question I have is : How can I found out the coordinate system of the image I get from the reader in both cases : when I use DICOMParser and when I use gdcm ? (in both cases I want the files to be read ascending , by Image Position Patient) </DIV>
<DIV> </DIV>
<DIV>I use itk 2.0</DIV>
<DIV> </DIV>
<DIV>Thank you in advance! <BR>Roxana</DIV><p>__________________________________________________<br>Do You Yahoo!?<br>Tired of spam? Yahoo! Mail has the best spam protection around <br>http://mail.yahoo.com