X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src%2FgdcmOrientation.cxx;h=b6a33bb8e440caf0fe0bc5cf2e140bfa8b2f0cbd;hb=77cbc8232f9308aafb17bde18ef9549be212b741;hp=a8bb31d9e736deaf956505f8499136c6c2ca4002;hpb=8fce3ef7e4121ea7a4e5c940b6e38394abcacb07;p=gdcm.git diff --git a/src/gdcmOrientation.cxx b/src/gdcmOrientation.cxx index a8bb31d9..b6a33bb8 100644 --- a/src/gdcmOrientation.cxx +++ b/src/gdcmOrientation.cxx @@ -3,8 +3,8 @@ Program: gdcm Module: $RCSfile: gdcmOrientation.cxx,v $ Language: C++ - Date: $Date: 2005/09/29 14:23:58 $ - Version: $Revision: 1.13 $ + Date: $Date: 2005/10/03 14:54:16 $ + Version: $Revision: 1.16 $ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). All rights reserved. See Doc/License.txt or @@ -51,21 +51,42 @@ namespace gdcm * # 6 : Heart Sagital * # -6 : Heart Sagital invert */ -double Orientation::TypeOrientation( File *f ) + +static const char *OrientationTypeStrings[] = { + "Not Applicable", + "Axial", + "Coronal", + "Sagital", + "Heart Axial", + "Heart Coronal", + "Heart Sagital", + "Axial invert", + "Coronal invert", + "Sagital invert", + "Heart Axial invert", + "Heart Coronal invert", + "Heart Sagital invert", + NULL +}; + +const char* Orientation::GetOrientationTypeString(OrientationType const o) +{ + int k = (int)o; + if (k < 0) + k = -k + 6; + + return OrientationTypeStrings[k]; +} + +OrientationType Orientation::GetOrientationType( File *f ) { float iop[6]; bool succ = f->GetImageOrientationPatient( iop ); if ( !succ ) { - gdcmErrorMacro( "No Image Orientation (0020,0037) found in the file, cannot proceed." ) - return 0; + gdcmErrorMacro( "No Image Orientation (0020,0037)/(0020,0032) found in the file, cannot proceed." ) + return NotApplicable; } -/* -std::cout << " iop : "; -for(int i=0;i<6;i++) - std::cout << iop[i] << " "; -std::cout << std::endl; -*/ vector3D ori1; vector3D ori2; @@ -103,18 +124,11 @@ std::cout << std::endl; res=VerfCriterion( i, CalculLikelyhood2Vec(refA,refB,ori1,ori2), res ); res=VerfCriterion( -i, CalculLikelyhood2Vec(refB,refA,ori1,ori2), res ); } - return res.first; -/* -// i=0 -// res=[0,99999] ## [ , ] -// for plane in dicPlane: -// i=i+1 -// refA=plane[0] -// refB=plane[1] -// res=self.VerfCriterion( i , self.CalculLikelyhood2Vec(refA,refB,ori1,ori2) , res ) -// res=self.VerfCriterion( -i , self.CalculLikelyhood2Vec(refB,refA,ori1,ori2) , res ) -// return res[0] -*/ + // res thought looks like is a float value, but is indeed an int + // try casting it to int first then enum value to please VS7: + int int_res = (int)res.first; + gdcmAssertMacro( int_res <= 6 && int_res >= -6); + return (OrientationType)int_res; } Res @@ -131,16 +145,6 @@ Orientation::VerfCriterion(int typeCriterion, double criterionNew, Res const &in res.first = type; res.second = criterion; return res; -/* -// type = res[0] -// criterion = res[1] -// # if criterionNew<0.1 and criterionNew Personal remark from JPRx : +--> patient_position (0x0018,0x5100) can be "HFS ", "FFS ", "HFP ", "FFP " +--> the cosines may ahave any value -1< <+1, for MR images ! + +enum patient_slice_orientation_type + { + patient_slice_orientation_unknown = 0, + supine_headfirst_transaxial, + supine_headfirst_sagittal, + supine_headfirst_coronal, + supine_feetfirst_transaxial, + supine_feetfirst_sagittal, + supine_feetfirst_coronal, + prone_headfirst_transaxial, + prone_headfirst_sagittal, + prone_headfirst_coronal, + prone_feetfirst_transaxial, + prone_feetfirst_sagittal, + prone_feetfirst_coronal + }; + +void GetImageOrientationPatient(gdcm::File &h,F32 image_orientation_patient[6]) +{ + h.GetImageOrientationPatient(image_orientation_patient); +} + +#if 0 +// +// this is all completely cribbed from the xmedcon library, since +// we're trying to do what it does, mostly. +patient_slice_orientation_type +GetPatSliceOrient(gdcm::File &h) +{ + F32 image_orientation_patient[6]; + + // protected, do it the hard way + // h.GetImageOrientationPatient(image_orientation_patient); + GetImageOrientationPatient(h,image_orientation_patient); + + enum { headfirst, feetfirst } patient_orientation; + enum { supine, prone } patient_rotation; + enum { transaxial, sagittal, coronal } slice_orientation; + + std::string patient_position = h.GetEntryByNumber(0x0018,0x5100); + if(patient_position == "gdcm::Unfound") + { + patient_position = "HF"; + } + if(patient_position.find("HF") != std::string::npos) + { + patient_orientation = headfirst; + } + else if(patient_position.find("FF") != std::string::npos) + { + patient_orientation = feetfirst; + } + + if(patient_position.find("S") != std::string::npos) + { + patient_rotation = supine; + } + else if(patient_position.find("P") != std::string::npos) + { + patient_rotation = prone; + } + + if((image_orientation_patient[0] == 1 || image_orientation_patient[0] == -1) && + (image_orientation_patient[4] == +1 || image_orientation_patient[4] == -1)) + { + slice_orientation = transaxial; + } + else if((image_orientation_patient[1] == 1 || image_orientation_patient[1] == -1) && + (image_orientation_patient[5] == +1 || image_orientation_patient[5] == -1)) + { + slice_orientation = sagittal; + } + else if((image_orientation_patient[0] == 1 || image_orientation_patient[0] == -1) && + (image_orientation_patient[5] == +1 || image_orientation_patient[5] == -1)) + { + slice_orientation = coronal; + } + // + // combine results + patient_slice_orientation_type patient_slice_orientation = + patient_slice_orientation_unknown; + switch (patient_rotation) + { + case supine: + switch (patient_orientation) + { + case headfirst: + switch (slice_orientation) + { + case transaxial: + patient_slice_orientation = supine_headfirst_transaxial; + break; + case sagittal: + patient_slice_orientation = supine_headfirst_sagittal; + break; + case coronal: + patient_slice_orientation = supine_headfirst_coronal; + break; + } + break; + case feetfirst: + switch (slice_orientation) + { + case transaxial: + patient_slice_orientation = supine_feetfirst_transaxial; + break; + case sagittal: + patient_slice_orientation = supine_feetfirst_sagittal; + break; + case coronal: + patient_slice_orientation = supine_feetfirst_coronal; + break; + } + break; + } + break; + case prone: + switch (patient_orientation) + { + case headfirst: + switch (slice_orientation) + { + case transaxial: + patient_slice_orientation = prone_headfirst_transaxial; + break; + case sagittal: + patient_slice_orientation = prone_headfirst_sagittal; + break; + case coronal: + patient_slice_orientation = prone_headfirst_coronal; + break; + } + break; + case feetfirst: + switch (slice_orientation) + { + case transaxial: + patient_slice_orientation = prone_feetfirst_transaxial; + break; + case sagittal: + patient_slice_orientation = prone_feetfirst_sagittal; + break; + case coronal: + patient_slice_orientation = prone_feetfirst_coronal; + break; + } + break; + } + break; + } + if(patient_slice_orientation != patient_slice_orientation_unknown) + return patient_slice_orientation; + // + // this is what xmedcon does + if ((image_orientation_patient[0] == +1) && + (image_orientation_patient[4] == +1)) + patient_slice_orientation = supine_headfirst_transaxial; + else if ((image_orientation_patient[0] == -1) && + (image_orientation_patient[4] == +1)) + patient_slice_orientation = supine_feetfirst_transaxial; + else if ((image_orientation_patient[0] == -1) && + (image_orientation_patient[4] == -1)) + patient_slice_orientation = prone_headfirst_transaxial; + else if ((image_orientation_patient[0] == +1) && + (image_orientation_patient[4] == -1)) + patient_slice_orientation = prone_feetfirst_transaxial; + + else if ((image_orientation_patient[1] == +1) && + (image_orientation_patient[5] == -1)) + patient_slice_orientation = supine_headfirst_sagittal; + else if ((image_orientation_patient[1] == +1) && + (image_orientation_patient[5] == +1)) + patient_slice_orientation = supine_feetfirst_sagittal; + else if ((image_orientation_patient[1] == -1) && + (image_orientation_patient[5] == -1)) + patient_slice_orientation = prone_headfirst_sagittal; + else if ((image_orientation_patient[1] == -1) && + (image_orientation_patient[5] == +1)) + patient_slice_orientation = prone_feetfirst_sagittal; + + else if ((image_orientation_patient[0] == +1) && + (image_orientation_patient[5] == -1)) + patient_slice_orientation = supine_headfirst_coronal; + else if ((image_orientation_patient[0] == -1) && + (image_orientation_patient[5] == +1)) + patient_slice_orientation = supine_feetfirst_coronal; + else if ((image_orientation_patient[0] == -1) && + (image_orientation_patient[5] == -1)) + patient_slice_orientation = prone_headfirst_coronal; + else if ((image_orientation_patient[0] == +1) && + (image_orientation_patient[5] == +1)) + patient_slice_orientation = prone_feetfirst_coronal; + return patient_slice_orientation; +} +#else + +-------------------------------------------------------------------------*/ } // end namespace gdcm