+ VectorMapInfoDicom vectormapinfodicom;
+ uint16_t group;
+ uint16_t elem;
+ std::string key;
+ std::string strTagValue; // read on disc
+ int iTag,sizeDicomTagsVector=bbGetInputDicomTags().size();
+ GDCM_NAME_SPACE::FileList::const_iterator iitt = l->begin();
+ float px,py,pz;
+
+ std::vector<std::string> lstFileNames;
+
+ for ( ; iitt != l->end(); ++iitt)
+ {
+ lstFileNames.push_back( (*iitt)->GetFileName() );
+ // EED 2022-12-15
+ if ( iitt == l->begin() ) // The first ordered element IPP
+ {
+ group = 0;
+ elem = 0;
+ strTagValue = std::string("void");
+ sscanf("D0020_0032","D%04hx_%04hx",&group,&elem);
+ strTagValue = (*iitt)->GetEntryString(group,elem);
+ sscanf(strTagValue.c_str(),"%f\\%f\\%f\\%f", &px , &py , &pz);
+ std::vector<double> v_ipp;
+ v_ipp.push_back( (double)px );
+ v_ipp.push_back( (double)py );
+ v_ipp.push_back( (double)pz );
+ bbSetOutputIPP(v_ipp );
+ } // if first elemetn get IPP
+
+ MapInfoDicom mapinfodicom;
+ for (iTag=0; iTag<sizeDicomTagsVector; iTag++)
+ {
+ key = bbGetInputDicomTags()[iTag];
+ group = elem = 0;
+ if ( (key.size()==10) && (key[0] == 'D') && (key[5] == '_') )
+ {
+ sscanf(key.c_str(),"D%04hx_%04hx ",&group,&elem);
+ }// if key
+ strTagValue = (*iitt)->GetEntryString(group,elem);
+ mapinfodicom.insert ( std::pair<std::string,std::string>(key,strTagValue) );
+ } // for iTag
+ vectormapinfodicom.push_back( mapinfodicom );
+ } // for iitt
+ bbSetOutputDicomInfo( vectormapinfodicom );
+ bbSetOutputOutFileNames( lstFileNames );
+
+ if (bbGetInputReadRaw()==true)
+ {
+ if (reader!=NULL)
+ {
+ reader->Delete();
+ reader = NULL;
+ }
+ reader = vtkGdcmReader::New();
+ //EED 21 mars 2012 FLIP probleme ..PLOP..
+ reader->SetFlipY(false);
+ //reader->SetFileName( bbGetInputIn().c_str() );
+ reader->SetCoherentFileList(l);
+ reader->Update();
+ reader->GetOutput();
+ bbSetOutputOut( reader->GetOutput() );
+ } // if bGetInputReadRaw
+}
+#endif
+// endif USE_GDCM
+
+
+#if defined USE_GDCM2
+void GetXCoherentInfoGdcmReader::Process()
+{
+// printf("GetXCoherentInfoGdcmReader %d files.\n", bbGetInputIn().size() );
+ // Read the *first* image file (a SET of file names is given as input)
+ gdcm::Reader *read = new gdcm::Reader();
+ reader->SetFileName( bbGetInputIn()[0].c_str() );
+ bool res = read->Read();
+ if ( !res )
+ {
+ delete read;
+ bbSetOutputOut(NULL);
+ return;
+ }
+
+ // Get info from the first image file
+ const gdcm::File &f = read->GetFile();
+ int i;
+ std::vector<double> v_iop;
+ const gdcm::DataElement &deIop = f.GetDataSet().GetDataElement(gdcm::Tag(0x0020, 0x0037));
+ std::stringstream ss;
+ deIop.GetValue().Print(ss);
+ gdcm::Element<gdcm::VR::DS,gdcm::VM::VM6> iop;
+ iop.Read( ss );
+ for(i=0; i< 6; i++)
+ {
+ v_iop.push_back((float)(iop[i]));
+ }
+ bbSetOutputIOP(v_iop );
+
+ std::vector<double> v_ipp;
+ const gdcm::DataElement &deIpp = f.GetDataSet().GetDataElement(gdcm::Tag((0x0020,0x0032)));
+ deIpp.GetValue().Print(ss);
+ gdcm::Element<gdcm::VR::DS,gdcm::VM::VM3> ipp;
+ ipp.Read( ss );
+ for(i=0; i< 3; i++)
+ {
+ v_ipp.push_back((float)(ipp[i]));
+ } // for 3
+ bbSetOutputIPP(v_ipp );
+
+// Add *all the files* to the IPPsorter
+ gdcm::IPPSorter s;
+ s.SetComputeZSpacing( true );
+ s.SetZSpacingTolerance( 1e-3 );
+ std::vector<std::string> gii = bbGetInputIn();
+ s.Sort(gii);
+
+ //l = sh->GetFirstSingleSerieUIDFileSet();
+ // no Test if we have multiple series
+ //s.GetFilenames();
+ //int nbFiles;
+
+ double zspacing = 0.;
+ zspacing = s.GetZSpacing();
+ std::vector<double> v_pixelspacing = gdcm::ImageHelper::GetSpacingValue(f);
+ v_pixelspacing.push_back( v_pixelspacing[0] );
+ v_pixelspacing.push_back( v_pixelspacing[1] );
+ v_pixelspacing.push_back( zspacing );
+ bbSetOutputPixelSpacing(v_pixelspacing);
+