]> Creatis software - creaContours.git/blobdiff - lib/kernel_ManagerContour_NDimensions/KernelManagerContour.cxx
no message
[creaContours.git] / lib / kernel_ManagerContour_NDimensions / KernelManagerContour.cxx
index 30cff3614c071f943e2a35faf29e2dc23a61d2c1..1fb5e3f9416463d5a98bc6b7a83885e8789bb151 100644 (file)
@@ -12,7 +12,7 @@ KernelManagerContour::KernelManagerContour(){
 **
 **/
        inredo = 0;
-       inundo = 0;     
+       inundo = 0;
 
        time_t seconds;
        seconds = time (NULL);
@@ -42,7 +42,7 @@ KernelManagerContour::KernelManagerContour(std::vector<vtkImageData*> images,std
 **
 **/
        inredo = 0;
-       inundo = 0;     
+       inundo = 0;
 
        time_t seconds;
        seconds = time (NULL);
@@ -51,7 +51,7 @@ KernelManagerContour::KernelManagerContour(std::vector<vtkImageData*> images,std
 
        stundoredo = tmpdir+"/temp"+intToString(time);
        _currentIndex = 0;
-       
+
        _contourPropagation = NULL;
 #if(WIN32)
        mkdir(tmpdir.c_str());
@@ -64,35 +64,58 @@ KernelManagerContour::KernelManagerContour(std::vector<vtkImageData*> images,std
 /**
 **  FINISH PERSISTANCE
 **/
-       setVectImages(images);
+       setVectImages(images);        
        initializeEnvironment(datadir);
 
 }
-KernelManagerContour::~KernelManagerContour(){
-}      
 
-std::vector<vtkImageData*> KernelManagerContour::getVectImages(){
+KernelManagerContour::~KernelManagerContour()
+{
+}
+
+std::vector<vtkImageData*> KernelManagerContour::getVectImages()
+{
        return vectimages;
 }
-void KernelManagerContour::setVectImages(std::vector<vtkImageData*> vectimg){
 
-       for(int i = 0; i < vectimg.size(); i++){
+
+void KernelManagerContour::setVectImages(std::vector<vtkImageData*> vectimg)
+{
+       double spc[3];
+       std::vector<double> vectspc;
+       
+       for(int i = 0; i < (int)(vectimg.size()); i++){
                vtkImageData* img = vectimg[i];
                vtkImageChangeInformation* change = vtkImageChangeInformation::New();
                change->SetInformationInput(img);
                change->SetInputConnection(img->GetProducerPort());
 
+               img->GetSpacing(spc);
                change->SetOutputSpacing(1,1,1);
-               
                change->Update();
 
+               vectspc.clear();
+               vectspc.push_back(spc[0]);
+               vectspc.push_back(spc[1]);
+               vectspc.push_back(spc[2]);
+               vectimagesSpacing.push_back(vectspc);
                vectimg[i] = change->GetOutput();
        }
 
        vectimages = vectimg;
 }
 
+
+void KernelManagerContour::GetSpacing(double *vecspc, int iImage)
+{
+       vecspc[0]=vectimagesSpacing[iImage][0];
+       vecspc[1]=vectimagesSpacing[iImage][1];
+       vecspc[2]=vectimagesSpacing[iImage][2];
+}
+
+
 void KernelManagerContour::initializeEnvironment(std::string datadir){
+       _datadir = datadir;
        std::string conceptsFN                  =  datadir+"holaConceptsFile.cf";
        std::string imageSourcesFN              =  datadir+"holaImagesInstantsFile.of";
        std::string imageSectionsFN             = "";
@@ -101,19 +124,19 @@ void KernelManagerContour::initializeEnvironment(std::string datadir){
        std::map<std::string, ImageSourceThing *> * sourcesMap          = new std::map<std::string, ImageSourceThing *>();
        std::map<std::string, ImageSectionThing *>* sectionsMap         = new std::map<std::string, ImageSectionThing *>();
        std::map<std::string, AxeThing *>* axesMap                                      = new std::map<std::string, AxeThing *>();
-       std::map<std::string, ContourThing *>* outlinesMap                      = new std::map<std::string, ContourThing *>();  
+       std::map<std::string, ContourThing *>* outlinesMap                      = new std::map<std::string, ContourThing *>();
 
-       for(int i = 0; i < vectimages.size(); i++){
+       for(int i = 0; i < (int)(vectimages.size()); i++){
                vtkImageData* selectedimage = vectimages[i];
                ImageSourceThing * thing                                                                        = new ImageSourceThing(selectedimage);
                std::string imgstring = "Source Image "+intToString(i+1);
                sourcesMap->insert(std::pair<std::string, ImageSourceThing *>( imgstring, thing));
        }
-       
+
 
        OutlineModelBuilder * _builder                                                          = new OutlineModelBuilder( conceptsFN ,datadir);
        _builder->buildImageSource_Envornment( imageSourcesFN, sourcesMap );
-       _builder->buildImageSection_Envornment( imageSectionsFN, sectionsMap );         
+       _builder->buildImageSection_Envornment( imageSectionsFN, sectionsMap );
        _builder->buildAxe_Envornment(axeThingsFN, axesMap );
        _builder->buildCountour_Envornment( imageSectionsFN, outlinesMap );
 
@@ -130,7 +153,7 @@ OutlineModelManager* KernelManagerContour::getOutlineModelManager(){
        return modelManager;
 }
 
-std::string KernelManagerContour::createOutline(manualContourModel * manModelContour,std::vector<int> instantVector){
+std::string KernelManagerContour::createOutline(manualBaseModel * manModelContour,std::vector<int> instantVector){
        return modelManager->createOutline( manModelContour, instantVector );
 }
 
@@ -143,7 +166,7 @@ std::string KernelManagerContour::intToString(int num){
                while (k > 0){
                        char temp = k % 10 + 48;
                        k = k / 10;
-                       result = temp + result; 
+                       result = temp + result;
                }
        }
        return result;
@@ -153,8 +176,8 @@ std::vector<std::string> KernelManagerContour::GetLstNameThingsStatic(){
        return modelManager->GetLstNameThingsStatic();
 }
 
-void KernelManagerContour::SaveThingName(FILE* pFile, std::string name ){
-       modelManager->SaveThingName(pFile, name);
+void KernelManagerContour::SaveThingName(FILE* pFile, FILE *pFileData, std::string name ){
+       modelManager->SaveThingName(pFile, pFileData, name);
 }
 std::vector<std::string> KernelManagerContour::GetLstNameThings(){
        return modelManager->GetLstNameThings();
@@ -168,7 +191,7 @@ bool KernelManagerContour::IsPartOfStaticList(std::string theKeyName){
 }
 
 void KernelManagerContour::deleteCModel(std::string theKeyName){
-       manualContourModel* cModel              = modelManager->getOutlineByKeyName(theKeyName)->getModel();
+       manualBaseModel* cModel         = modelManager->getOutlineByKeyName(theKeyName)->getModel();
        modelManager->removeOutline( theKeyName );
 
        delete cModel;
@@ -207,7 +230,9 @@ std::string KernelManagerContour::createCopyContourOf ( std::string anExistingKN
 {
        return modelManager->createCopyContourOf(anExistingKName, instantNoTouchData);
 }
-manualContourModel* KernelManagerContour::getOutlineByKeyName(std::string cloneName){
+
+manualBaseModel* KernelManagerContour::getOutlineByKeyName(std::string cloneName)
+{
        return modelManager->getOutlineByKeyName (cloneName )->getModel();
 }
 
@@ -215,16 +240,16 @@ bool KernelManagerContour::onRedo(std::string& filename){
        if(inredo > 0){
                inredo--;
                inundo++;
-               std::string str = intToString(inundo);        
-               filename = stundoredo + str + ".roi";   
-               //loadState(temp);              
+               std::string str = intToString(inundo);
+               filename = stundoredo + str + ".roi";
+               //loadState(temp);
                return true;
        }
        return false;
 }
 
 bool KernelManagerContour::onUndo(std::string& filename){
-       if(inundo>0){           
+       if(inundo>0){
                inredo++;
                inundo--;
 
@@ -232,9 +257,9 @@ bool KernelManagerContour::onUndo(std::string& filename){
                //itoa(inundo, str, 10);
                std::string str = intToString(inundo);
 
-        
+
                filename = stundoredo + str + ".roi";
-       
+
                return true;//loadState(temp);
        }
        return false;
@@ -258,7 +283,7 @@ bool KernelManagerContour::onUndoSaveFile(std::string& filename){
 
                        filename = stundoredo + str + ".roi";
                        return true;
-                       //saveFileWithContours(temp);                           
+                       //saveFileWithContours(temp);
 
                }
        }
@@ -270,55 +295,57 @@ void KernelManagerContour :: changeContourOfManager(std::string keyName, Instant
        modelManager->ChangeContourOfList(keyName, instant);
 }
 
-void KernelManagerContour ::resetAppend(){
-       if(isInitContourPropagation()){
+void KernelManagerContour ::resetAppend()
+{
+    if (_contourPropagation!=NULL)
+       {
                _contourPropagation->resetAppend();
-       }       
-}
-
-bool KernelManagerContour::isInitContourPropagation(){
-       if(_contourPropagation == NULL){
-               _contourPropagation = new ContourPropagation();
-       }else{
-               return true;
        }
-       return false;
-
 }
 
-std::string KernelManagerContour::onSpreadAdd( std::vector<double> *vecX, std::vector<double> *vecY, std::vector<double> *vecZ, std::vector<int> instants){
-       if(isInitContourPropagation()){
-               if (vecX->size()!=0){
-                       
-                       int i,size=vecZ->size();
-                       int actualSlice = instants[1];
-                       for ( i=0 ; i<size ; i++ )
-                       {
-                               (*vecZ)[i] = actualSlice;
-                       } // for
 
-                       _contourPropagation->appendContour(vecX , vecY , vecZ);
-                       return intToString(actualSlice);                
-               }
-       }
+std::string KernelManagerContour::onSpreadAdd( std::vector<double> *vecX, std::vector<double> *vecY, std::vector<double> *vecZ, std::vector<int> instants)
+{
+    if (_contourPropagation==NULL)
+    {
+               _contourPropagation = new ContourPropagation();
+    }
+
+    if (vecX->size()!=0){
+
+        int i,size=vecZ->size();
+        int actualSlice = instants[1];
+        for ( i=0 ; i<size ; i++ )
+        {
+            (*vecZ)[i] = actualSlice;
+        } // for
+
+        _contourPropagation->appendContour(vecX , vecY , vecZ);
+        return intToString(actualSlice);
+    }
        return "";
 }
 
-void KernelManagerContour::getMaxMinZ(double *minZ,double *maxZ){
-       if(isInitContourPropagation()){
+
+void KernelManagerContour::getMaxMinZ(double *minZ,double *maxZ)
+{
+    if (_contourPropagation!=NULL)
+       {
                _contourPropagation->getMaxMinZ(minZ, maxZ);
        }
 }
-void KernelManagerContour::CalculeSplinePropagation(){
+
+
+void KernelManagerContour::CalculeSplinePropagation()
+{
        _contourPropagation->setInterpolationNumber(100);
        _contourPropagation->CalculeSplinePropagation();
 }
 
-manualContourModel* KernelManagerContour::GetPoints(int z,int type, std::vector<double>* vecCtrlPointX,std::vector<double>* vecCtrlPointY,std::vector<double>* vecCtrlPointZ, std::string& theName,int typeofcontour, std::vector<int> tempVector){
-    
-       
+manualBaseModel* KernelManagerContour::GetPoints(int z,int type, std::vector<double>* vecCtrlPointX,std::vector<double>* vecCtrlPointY,std::vector<double>* vecCtrlPointZ, std::string& theName,int typeofcontour, std::vector<int> tempVector)
+{
        bool addedModel = false;
-       manualContourModel* manModelContour=NULL;
+       manualBaseModel* manModelContour=NULL;
 
        if (_contourPropagation->ifSliceKeyContourExist(z)==false){
 
@@ -327,9 +354,9 @@ manualContourModel* KernelManagerContour::GetPoints(int z,int type, std::vector<
                int idTmp = _contourPropagation->FindIdWithZ(z);
 
                if (type==0) // Initial Points
-               {   
+               {
                        _contourPropagation->GetInitialControlPoints( idTmp , vecCtrlPointX,vecCtrlPointY,vecCtrlPointZ);
-               } 
+               }
                if (type==1)  // Automatique Method
                {
                        _contourPropagation->GetControlPoints( idTmp  ,vecCtrlPointX,vecCtrlPointY,vecCtrlPointZ);
@@ -343,7 +370,7 @@ manualContourModel* KernelManagerContour::GetPoints(int z,int type, std::vector<
                for (int j=0 ; j<sizeCtrPt ; j++)
                {
 //JSTG_16-07-08_----------------------------------------------------------------
-                       manModelContour->AddPoint( (*vecCtrlPointX)[j] , (*vecCtrlPointY)[j] , -900  );
+                       manModelContour->AddPoint( (*vecCtrlPointX)[j] , (*vecCtrlPointY)[j] , 900  );
 //--------------------------------------------------------------------
                } // for j
 
@@ -353,15 +380,17 @@ manualContourModel* KernelManagerContour::GetPoints(int z,int type, std::vector<
                if(!addedModel){
                        manModelContour = NULL;
                }
-               
+
        }// ifSliceKeyContourExist
        return manModelContour;
-
 }
 
-manualContourModel * KernelManagerContour::factoryManualContourModel(int typeContour)
-{      
-       manualContourModel *manModelContour=NULL;
+manualBaseModel * KernelManagerContour::factoryManualContourModel(int typeContour)
+{
+       manualBaseModel *manModelContour=NULL;
+
+       // Creating the model
+       // NOTE: The view and the controler are created in the wxVtkBaseView_SceneManager class, configureViewControlTo method
 
        // spline
        if (typeContour==0)
@@ -393,6 +422,18 @@ manualContourModel * KernelManagerContour::factoryManualContourModel(int typeCon
                manModelContour = new manualContourModelLine();
        }
 
+       // points
+       if (typeContour==7)
+       {
+               manModelContour = new manualBaseModel();
+       }
+
+       // polygon
+       if (typeContour==10)
+       {
+               manModelContour = new manualContourModelPolygon();
+       }
+
        return manModelContour;
 }
 
@@ -400,20 +441,20 @@ std::vector<std::string> KernelManagerContour::getOutlinesNameAtInstant(std::vec
        Instant instant(&tempvector);
        std::vector<ContourThing**> vectcont = modelManager->getOutlinesAtInstant( &instant );
        std::vector<std::string> vectname;
-       for(int i = 0; i < vectcont.size(); i++){
+       for(int i = 0; i < (int)(vectcont.size()); i++){
                ContourThing **contourthing = vectcont[i];
                vectname.push_back((*contourthing)->getName());
        }
-    return vectname;   
+    return vectname;
 }
 
 std::vector<ContourThing**> KernelManagerContour::getOutlinesAtInstant(Instant* instant ){
        return modelManager->getOutlinesAtInstant(instant);
 }
 
-std::vector<manualContourModel*> KernelManagerContour::ExploseEachModel( std::vector<manualContourModel*> lstManConMod ){
-       std::vector<manualContourModel*> lstTmp;
-       std::vector<manualContourModel*> lstResult;
+std::vector<manualBaseModel*> KernelManagerContour::ExploseEachModel( std::vector<manualBaseModel*> lstManConMod ){
+       std::vector<manualBaseModel*> lstTmp;
+       std::vector<manualBaseModel*> lstResult;
        int j,jSize;
        int i,iSize=lstManConMod.size();
        for (i=0;i<iSize;i++)
@@ -436,7 +477,7 @@ void KernelManagerContour::getConceptsInformation(std::vector<std::string>& conc
 vtkImageData* KernelManagerContour::getImageAtInstant(std::vector<int> inst){
        int index = inst[5]-1;
 
-       if(index < vectimages.size()&&index!=_currentIndex){
+       if(index < (int)(vectimages.size())&&index!=_currentIndex){
                _currentIndex=index;
                return vectimages[index];
        }
@@ -451,4 +492,30 @@ void KernelManagerContour::setCurrentFileName(std::string filenam){
        this->filename = filenam;
 }
 
+std::string KernelManagerContour::parseOsirixFile(std::string filename){
+
+
+#ifdef ParserOsirix_BUILD
+       vtkImageData* sourceimage;
+       std::string xsdfile;
+
+       xsdfile = _datadir;
+
+       xsdfile.append("\\XML\\osirixschema.xsd");
+
+       sourceimage = getSourceImage();
+       OsirixParser p(xsdfile.c_str(), sourceimage->GetSpacing(), sourceimage->GetExtent());
+
+       if(p.ParseFile(filename.c_str())!= 0){
+
+       }
+
+       return p.getContoursFileName();
+#else
+       return "";
+#endif
+
+
+}
+