X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src2%2FcreaImageIOGimmickView.cpp;h=30ca0d5d094b9bd430e6e1824da392bb64057bf7;hb=a13b9d55cc9b8ca9b886e99a73946dd8eccefe65;hp=c7e0cf15c0f1aedb86f6ee447a76b2b70e24abae;hpb=da80530068246081505e9649523d10e8da8e6acd;p=creaImageIO.git diff --git a/src2/creaImageIOGimmickView.cpp b/src2/creaImageIOGimmickView.cpp index c7e0cf1..30ca0d5 100644 --- a/src2/creaImageIOGimmickView.cpp +++ b/src2/creaImageIOGimmickView.cpp @@ -2,6 +2,11 @@ #include #include "boost/filesystem.hpp" +#include // JPR +#include // JPR +#include // JPR +#include + namespace fs = boost::filesystem; namespace creaImageIO @@ -13,22 +18,21 @@ namespace creaImageIO public: ImageExtent(const std::string& x, const std::string& y, const std::string& z, const std::string& t) { - sscanf(x.c_str(),"%d",&mExtent[0]); - sscanf(y.c_str(),"%d",&mExtent[1]); - sscanf(z.c_str(),"%d",&mExtent[2]); - sscanf(t.c_str(),"%d",&mExtent[3]); - if(x==""){mExtent[0]=1;} - if(y==""){mExtent[1]=1;} - if(z==""){mExtent[2]=1;} - if(t==""){mExtent[3]=1;} - - if (mExtent[3]>1) mDim=4; - else if (mExtent[2]>1) mDim=3; - else if (mExtent[1]>1) mDim=2; - else if (mExtent[0]>1) mDim=1; - else mDim=0; + sscanf(x.c_str(),"%d",&mExtent[0]); + sscanf(y.c_str(),"%d",&mExtent[1]); + sscanf(z.c_str(),"%d",&mExtent[2]); + sscanf(t.c_str(),"%d",&mExtent[3]); + if(x==""){mExtent[0]=1;} + if(y==""){mExtent[1]=1;} + if(z==""){mExtent[2]=1;} + if(t==""){mExtent[3]=1;} + + if (mExtent[3]>1) mDim=4; + else if (mExtent[2]>1) mDim=3; + else if (mExtent[1]>1) mDim=2; + else if (mExtent[0]>1) mDim=1; + else mDim=0; } - ///Clears the extent void Clear() { mExtent[0] = mExtent[1] = mExtent[2] = mExtent[3] = 1; } @@ -66,7 +70,6 @@ namespace creaImageIO // Anciently started the threads ... // Threads now automatically start at first image request //mReader.Start(); - } //====================================================================== @@ -88,7 +91,7 @@ namespace creaImageIO mReaderStarted=false; } //====================================================================== - + //====================================================================== /// Finalize void GimmickView::Finalize() @@ -118,7 +121,6 @@ namespace creaImageIO } - //====================================================================== //====================================================================== @@ -133,7 +135,7 @@ namespace creaImageIO GimmickError("INTERNAL ERROR : GimmickView::UpdateTreeView : '" <second->UpdateLevel(l); + i->second->UpdateLevel(l); } //====================================================================== @@ -148,7 +150,7 @@ namespace creaImageIO } //====================================================================== - + //====================================================================== bool ImageExtent::IsCompatible(const ImageExtent& ie) { @@ -163,7 +165,7 @@ namespace creaImageIO } //====================================================================== - + //====================================================================== void ImageExtent::Add(const ImageExtent& ie) { @@ -283,22 +285,31 @@ namespace creaImageIO ///Reads Images (Non Threaded) void GimmickView::ReadImagesNotThreaded(std::vector& s, std::vector im, int dimension) { + +/* remember! + +#define GIMMICK_NO_IMAGE_SELECTION 0 +#define GIMMICK_2D_IMAGE_SELECTION 2 +#define GIMMICK_3D_IMAGE_SELECTION 3 +#define GIMMICK_4D_IMAGE_SELECTION 4 + +#define NATIVE 0 +#define _2D 2 +#define _3D 3 + +*/ // Create the output data if (im.size()==1) - { - // Only one image : give it - vtkImageData* out = vtkImageData::New(); - GimmickDebugMessage(3, "State Check: Full Filename: " - <ShallowCopy(mReader.GetImage(im.front())); - s.push_back( out ); - } - else if (im.size()>1) // Test inutile ? JPR { + // O0. + } + else if (im.size()>1) // Test inutile ? JPR + { + // Read the FIRST image to get info ... about all the other ones. ?!? // JPR vtkImageData* first = mReader.GetImage( im.front()); + if (dimension == 2) - { + { // n3D std::vector::iterator it; for (it=im.begin(); it!=im.end(); ++it) @@ -308,20 +319,20 @@ void GimmickView::ReadImagesNotThreaded(std::vector& s, std::vect s.push_back(out); } } - else + else if (dimension == 3) { - // n2D to 3D + // n*2D to 3D :: --> NON ! Ca va tres bien pour du 2D+T (le ZSpacing, on s'en tape) + // // ne convient pas pour du n*2D a  lire comme un 3D + // // Il faut un gdcmSerieHelper ! vtkImageData* out = vtkImageData::New(); out->CopyStructure(first); out->SetScalarType(first->GetScalarType()); + out->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents()); int ext[6]; - //first->GetExtent(ext); // JPR - first->GetWholeExtent(ext); -std::cout <<"in GimmickView::ReadImagesNotThreaded GetWholeExtent ext ="; -for (int jjj=0;jjj<6;jjj++) - std:cout << " [" << jjj << "]=" << ext[jjj]; -std::cout << endl; - + //first->GetExtent(ext); // JPR + first->GetWholeExtent(ext); // renvoie egalement 0,0 en Z // JPR + // Normal : + if(ext[5] == 0) { ext[5] = im.size()-1; @@ -331,37 +342,99 @@ std::cout << endl; ext[5] = ext[5] * im.size()-1; // to deal with multiframes - JPR } out->SetExtent(ext); - // LG : TODO : Z Spacing ? - - out->AllocateScalars(); - //first->Print(std::cout); - // out->Print(std::cout); + // LG : TODO : Z Spacing ? // Voir + loin // JPR int dim[3]; first->GetDimensions(dim); + + out->SetDimensions(dim[0], dim[1], im.size() ); + out->AllocateScalars(); + out->Update(); + unsigned long imsize = dim[0] * dim[1]; - imsize = imsize * dim[2] ; // deal with multiframes// JPR + imsize = imsize * dim[2] ; // deal with multiframes // JPR + + //EED 03-11-2009 - imsize = imsize * first->GetScalarSize(); + // differents formats char , short, etc... + // differents components 1..3 ex. jpg ->RGB 3 + imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents(); + + // Order the file name vector + + double spc[3]; + out->GetSpacing(spc); + spc[2]=OrderTheFileNameVector(im); + out->SetSpacing(spc); int slice = 0; std::vector::iterator it; + for (it=im.begin(); it!=im.end(); ++it) { - - vtkImageData* cur = mReader.GetImage( (*it)); - // void* src = cur->GetScalarPointer(0,0,0); - // void* dst = out->GetScalarPointer(0,0,slice); - memcpy(out->GetScalarPointer(0,0,slice),cur->GetScalarPointer(0,0,0),imsize); - - slice++; + vtkImageData* cur = mReader.GetImage( (*it) ); + memcpy(out->GetScalarPointer(0,0,slice), cur->GetScalarPointer(0,0,0), imsize); + slice++; } s.push_back(out); + + } // dimension == 3 + + } // size >1 + +} + //====================================================================== + + +double GimmickView::OrderTheFileNameVector(std::vector &im) +{ + double spacing=1; + typedef std::vector FileList; + FileList fileVector; + //GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New(); + GDCM_NAME_SPACE::SerieHelper *sh = GDCM_NAME_SPACE::SerieHelper::New(); + std::vector lstAux; + std::vector::iterator it; + for (it=im.begin(); it!=im.end(); ++it) + { + ///\TODO liberer les GDCM_NAME_SPACE::File a la fin! // JPR + GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New(); + f->SetFileName(*it); + f->Load(); + if (f->IsReadable()) + { + fileVector.push_back(f); + } else { + lstAux.push_back(*it); } + } // for + + + if ((fileVector.size()>1) && (sh->IsCoherent( &fileVector ))) + { + sh->OrderFileList(&fileVector); + spacing= sh->GetZSpacing(); + im.clear(); + int i; + for (i=0; iGetFileName() ); + } + for (i=0; i& s, std::vector im, int dimension) { @@ -382,9 +455,10 @@ void GimmickView::ReadImagesNotThreadedInVector(std::vector& s, s if (dimension == 2) { // n3D - std::vector::iterator it; + std::vector::iterator it; for (it=im.begin(); it!=im.end(); ++it) { + std::cout << "----------------------------------------n3D" << std::endl; vtkImageData* out = vtkImageData::New(); out->ShallowCopy(mReader.GetImage(*it)); s.push_back(out); @@ -394,11 +468,15 @@ void GimmickView::ReadImagesNotThreadedInVector(std::vector& s, s { // n2D to 3D // NO! // n *2D + T in a vector : - + std::cout << "---------------------------------------- n *2D + T" << std::endl; std::vector::iterator it; for (it=im.begin(); it!=im.end(); ++it) { - vtkImageData* out = mReader.GetImage( (*it)); + std::cout << "---------------------------------------- [" << *it << "]" << std::endl; + //vtkImageData* out = mReader.GetImage( (*it)); + //s.push_back(out); + vtkImageData* out = vtkImageData::New(); + out->ShallowCopy(mReader.GetImage(*it)); s.push_back(out); } } @@ -411,15 +489,15 @@ void GimmickView::ReadImagesNotThreadedInVector(std::vector& s, s void GimmickView::RequestReading(tree::Node* n, int prio, int selection_index, boost::shared_ptr p) { - if(!mReaderStarted) - { - mReader.Start(); - mReaderStarted=true; - } - ImageEventType t(n,selection_index); - t.pointerHolder = p; - mImageEventMap[n->GetAttribute("FullFileName")] = t; - mReader.Request(this,n->GetAttribute("FullFileName"),prio); + if(!mReaderStarted) + { + mReader.Start(); + mReaderStarted=true; + } + ImageEventType t(n,selection_index); + t.pointerHolder = p; + mImageEventMap[n->GetAttribute("FullFileName")] = t; + mReader.Request(this,n->GetAttribute("FullFileName"),prio); } //======================================================================