int outputdims[3];
imageIn->GetDimensions (inputdims);
unsigned short *input = (unsigned short *)(imageIn->GetScalarPointer());
int outputdims[3];
imageIn->GetDimensions (inputdims);
unsigned short *input = (unsigned short *)(imageIn->GetScalarPointer());
unsigned int div = (unsigned int)ceil(sqrt( (double)numberOfImagesInMosaic ) );
outputdims[0] = inputdims[0] / div;
unsigned int div = (unsigned int)ceil(sqrt( (double)numberOfImagesInMosaic ) );
outputdims[0] = inputdims[0] / div;
vtkImageOut = vtkImageData::New();
vtkImageOut->SetDimensions( outputdims );
vtkImageOut->SetExtent(0,outputdims[0]-1,0,outputdims[1]-1,0,outputdims[2]-1);
vtkImageOut = vtkImageData::New();
vtkImageOut->SetDimensions( outputdims );
vtkImageOut->SetExtent(0,outputdims[0]-1,0,outputdims[1]-1,0,outputdims[2]-1);
vtkImageOut->SetWholeExtent(0,outputdims[0]-1,0,outputdims[1]-1,0,outputdims[2]-1);
vtkImageOut->SetNumberOfScalarComponents(1);
//vtkImageOut->SetSpacing( blabla );
vtkImageOut->SetScalarType( VTK_UNSIGNED_SHORT );
vtkImageOut->AllocateScalars();
vtkImageOut->Update();
vtkImageOut->SetWholeExtent(0,outputdims[0]-1,0,outputdims[1]-1,0,outputdims[2]-1);
vtkImageOut->SetNumberOfScalarComponents(1);
//vtkImageOut->SetSpacing( blabla );
vtkImageOut->SetScalarType( VTK_UNSIGNED_SHORT );
vtkImageOut->AllocateScalars();
vtkImageOut->Update();