Program: bbtk
Module: $RCSfile: bbvtkConcatImages.cxx,v $
Language: C++
- Date: $Date: 2010/04/08 14:39:45 $
- Version: $Revision: 1.1 $
+ Date: $Date: 2010/04/28 22:10:07 $
+ Version: $Revision: 1.2 $
=========================================================================*/
/* ---------------------------------------------------------------------
void ConcatImages::Process()
{
int dim[3];
+ int nbComponents;
+ int scalarType;
int nb = (unsigned int)bbGetInputIn().size();
- if (nb == 0) {} // ??? JPR
+ int consistentNb(nb);
+
+ if (nb == 0) {
+ // ??? JPR
+ }
bbGetInputIn()[0]->GetDimensions(dim);
if (dim[2] > 1 || nb == 1) // dim[2] > 1 : Hopeless for vtk : the first file contains already a 3D object
//std::cout << "--------PrintSelf firstImage " << std::endl;
// firstImage->PrintSelf(std::cout, vtkIndent(2));
- int nbComponents = firstImage->GetNumberOfScalarComponents();
+ nbComponents = firstImage->GetNumberOfScalarComponents();
mConcat->SetNumberOfScalarComponents(nbComponents);
- mConcat->SetScalarType(firstImage->GetScalarType( ));
- firstImage->GetDimensions(dim);
- dim[2]=nb;
- mConcat->SetDimensions(dim);
- mConcat->SetExtent(0, dim[0]-1, 0, dim[1]-1, 0, dim[2]-1);
+
+ mConcat->SetScalarType(firstImage->GetScalarType( ));
+ scalarType = firstImage->GetScalarType( );
+ mConcat->SetScalarType(scalarType);
+
mConcat->SetSpacing(firstImage->GetSpacing());
+ firstImage->GetDimensions(dim);
+ // brute way to perform an ultra-mini consistency check :
+ // First image is supposed to be the reference image,
+ // any unconsistent image is just discarted...
+ for(int i=0; i<nb; i++)
+ {
+ if ( !CheckConsistency(firstImage, dim, nbComponents, scalarType))
+ consistentNb--;
+ }
+
+ dim[2]=consistentNb;
+ mConcat->SetDimensions(dim);
+ mConcat->SetExtent(0, dim[0]-1, 0, dim[1]-1, 0, dim[2]-1);
+
mConcat->AllocateScalars();
mConcat->Update();
for(index_image=0; index_image<nb; index_image++)
{
+ if (!CheckConsistency(bbGetInputIn()[index_image], dim, nbComponents, scalarType))
+ continue; // discard unconsistent image !
for(index_lignes=0; index_lignes<dim[0];index_lignes++)
{
for(index_colonnes=0; index_colonnes<dim[1];index_colonnes++)
bbSetOutputOut(mConcat);
}
+bool ConcatImages::CheckConsistency( vtkImageData *curImage, int dim[], int nbComponents, int scalarType )
+{
+ assert (curImage);
+
+ if (scalarType != curImage->GetScalarType( ) )
+ return false;
+
+ if (nbComponents != curImage->GetNumberOfScalarComponents() )
+ return false;
+
+ int curDim[3];
+ curImage->GetDimensions(curDim);
+ if ( (curDim[0] =! dim[0]) || (curDim[1] =! dim[1]) )
+ return false;
+
+ return true;
+}
+
+
}//namespace bbtk
#endif // _USE_VTK_