#3320 ReadImages box

[creaImageIO.git] / src / creaImageIOSimpleView.cpp

diff --git a/src/creaImageIOSimpleView.cpp b/src/creaImageIOSimpleView.cpp
index 43bbdf088e346d3a88ed435fcd0199d7eec7af71..0219927b6bdbb4a654a222abfa66aea9dda2c5a7 100644 (file)
--- a/src/creaImageIOSimpleView.cpp
+++ b/src/creaImageIOSimpleView.cpp
@@ -146,4 +146,104 @@ namespace creaImageIO
                        }
                        return bresult;
        }
+
+       
+       vtkImageData* SimpleView::getVolumeSelected( std::vector<vtkImageData*> lstImages)
+       {
+               vtkImageData *result;
+
+//        creaImageIO::WxSimpleDlg dlg(0,crea::std2wx(bbGetInputTitle()),"localdatabase_Descriptor.dscp","Local Database"); 
+//             dlg.SetAttrDicomTags( bbGetInputDicomTags() );
+//             dlg.ShowModal(); 
+//             bbSetOutputDicomInfo( dlg.getDicomInfoImagesSelected() );
+
+               if (lstImages.size()!=0){ 
+//               bbSetOutputFileName( std::string( dlg.getInfoImage().mb_str() ) );
+                 if(lstImages.size() ==1)
+                 {
+                       result = lstImages[0] ;  
+                 } else {
+                               // FCY: it will be a big problem if we have several kind of data in the same folder.
+                               // creation of a huge vtkImageData!!!!
+                               vtkImageData* first = lstImages[0];
+                               int dim[3];
+                               first->GetDimensions(dim);
+                               if (dim[2]==1)
+                               {
+                                       vtkImageData *out = vtkImageData::New();
+                                       int ext[6];
+//EED 2017-01-01 Migration VTK7
+#if VTK_MAJOR_VERSION <= 5
+                                       first->GetWholeExtent(ext); 
+#else
+                                       first->GetExtent(ext); 
+#endif
+                                       if(ext[5] == 0)
+                                       {
+                                          ext[5] = lstImages.size()-1;
+                                        } else {
+                                               ext[5] = ext[5] * lstImages.size()-1; // to deal with multiframes 
+                                       } // ext
+                                       out->SetExtent(ext);
+                                       double spac[3];
+                                       first->GetDimensions(dim);
+                                       first->GetSpacing(spac);
+                                       if (spac[0]==spac[1])
+                                       {
+                                               spac[2]=spac[0];  
+                               }
+                                       out->SetSpacing(spac);
+                                       out->SetDimensions(dim[0], dim[1], lstImages.size() );
+//EED 2017-01-01 Migration VTK7
+#if VTK_MAJOR_VERSION <= 5
+
+                                       out->SetScalarType(first->GetScalarType());
+                                       out->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
+                                       out->AllocateScalars();
+                                       out->Update();
+#else
+                                       out->AllocateScalars(first->GetScalarType(),first->GetNumberOfScalarComponents());
+#endif
+
+
+                                       unsigned long imsize = dim[0] * dim[1];
+                                       imsize = imsize * dim[2];  // deal with multiframes here
+                                       // differents formats char , short, etc...
+                                       // differents components 1..3  ex. jpg ->RGB 3
+                                       imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();
+                               int slice,sizeImageVector=lstImages.size();
+                                       for (slice=0 ; slice<sizeImageVector ; slice++)
+                                       {
+                                               vtkImageData *img = lstImages[slice];
+                                               memcpy(out->GetScalarPointer(0,0,slice), img->GetScalarPointer(0,0,0), imsize);
+                                               //img->Delete();
+                                       } // for slice  
+                                       out->Modified();
+//EED This is NOT working. The iteration over the pointer vtkImageData
+//                             int slice=0;
+//                             std::vector<vtkImageData*>::iterator it;
+//                             for (it=dlg.getImagesSelected().begin(); it!=dlg.getImagesSelected().end(); ++it) 
+//                                     {
+//                                               memcpy(out->GetScalarPointer(0,0,slice), (*it)->GetScalarPointer(0,0,0), imsize);
+//                                               slice++;
+//                                     }       
+//                             //for (it=dlg.getImagesSelected().begin(); it!=dlg.getImagesSelected().end(); ++it) 
+//                                     //{
+//                                     //              (*it)->Delete();
+//                                     //}
+
+                                       result =out;
+
+                               } else {
+                                       result = first;
+                               } // dim
+                       } // if getImagesSelected size
+               } else { 
+               result = NULL;  
+               } // getImagesSelected().size()!=0
+                       result->Modified();
+               return result;
+       }
+       
+       
 }