X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;ds=sidebyside;f=packages%2Fvtk%2Fsrc%2FbbvtkCreateImage.cxx;fp=packages%2Fvtk%2Fsrc%2FbbvtkCreateImage.cxx;h=8940726e583c6a4db8f7bd65b424e890b2fc5413;hb=9155aa82faa81612fbfe56671f80c10103e69595;hp=074e9de8a37d013d7a1f10b8ae6963035ab529ff;hpb=ddc63e39c2ec439a7b383fd0be55aff6b1d46a84;p=bbtk.git

diff --git a/packages/vtk/src/bbvtkCreateImage.cxx b/packages/vtk/src/bbvtkCreateImage.cxx
index 074e9de..8940726 100644
--- a/packages/vtk/src/bbvtkCreateImage.cxx
+++ b/packages/vtk/src/bbvtkCreateImage.cxx
@@ -53,13 +53,9 @@ void CreateImage::Process()
 //    * TYPE is the C++ type of the input/output
 //      (the one provided in the attribute 'type' of the tag 'input')
 
-	
-
-	double	spc[6];
+	double	spc[3];
 	int		dim[3];	
 	int		outputformat = VTK_UNSIGNED_SHORT;
-	
-
 	if (bbGetInputOutputFormat()=="VTK_BIT")					outputformat = VTK_BIT;                 // 1
 	else if (bbGetInputOutputFormat()=="VTK_CHAR")				outputformat = VTK_CHAR;				// 2
 	else if (bbGetInputOutputFormat()=="MET_CHAR")				outputformat = VTK_CHAR;				// 2
@@ -82,49 +78,53 @@ void CreateImage::Process()
 	else if (bbGetInputOutputFormat()=="MET_FLOAT")				outputformat = VTK_FLOAT;               // 10
 	else if (bbGetInputOutputFormat()=="VTK_DOUBLE")			outputformat = VTK_DOUBLE;              // 11  		
 	else if (bbGetInputOutputFormat()=="MET_DOUBLE")			outputformat = VTK_DOUBLE;              // 11  
-
-
-
 	spc[0] = bbGetInputSpacing()[0];
 	spc[1] = bbGetInputSpacing()[1];
 	spc[2] = bbGetInputSpacing()[2];
-	
 	dim[0] =  bbGetInputDimensions()[0];
 	dim[1] =  bbGetInputDimensions()[1];
 	dim[2] =  bbGetInputDimensions()[2];
-	
+	if ((dim[0]<=0)||(dim[1]<=0)||(dim[2]<=0))
+	{
+		dim[0]=1;
+		dim[1]=1;
+		dim[2]=1;
+	}
+	if ((spc[0]<=0)||(spc[1]<=0)||(spc[2]<=0))
+	{
+		spc[0]=1;
+		spc[1]=1;
+		spc[2]=1;
+	}
+//EED 2017-12-7
+	if (imageoutput==NULL) 
+	{
+		imageoutput = vtkImageData::New();
+	}
 	imageoutput->Initialize();
 	imageoutput->SetScalarType( outputformat );	
 	imageoutput->SetSpacing( spc );
 	imageoutput->SetDimensions(  dim[0], dim[1], dim[2] );
 	imageoutput->AllocateScalars();
-	
-
-		int i,j,k;
-		long sizeBlock;	
-
-		for (i=0; i<dim[0]; i++)
-		{
-			imageoutput->SetScalarComponentFromDouble(i,0,0, 0, bbGetInputInitialValue() );
-		}
-		
-		sizeBlock= dim[0] * imageoutput->GetScalarSize();
-                for (j=1; j<dim[1]; j++)
-		{
-			memcpy(  imageoutput->GetScalarPointer(0,j,0), imageoutput->GetScalarPointer()  , sizeBlock );
-		}		
-
-		sizeBlock= dim[0] * dim[1] * imageoutput->GetScalarSize();
-		for (k=0; k<dim[2]; k++)
-		{
-			memcpy(  imageoutput->GetScalarPointer(0,0,k) ,imageoutput->GetScalarPointer()  , sizeBlock );
-		} // for k
-
-		imageoutput->Update();
-		imageoutput->Modified();
-		bbSetOutputOut(imageoutput);
-
-
+	int i,j,k;
+	long sizeBlock;	
+	for (i=0; i<dim[0]; i++)
+	{
+		imageoutput->SetScalarComponentFromDouble(i,0,0, 0, bbGetInputInitialValue() );
+	}	
+	sizeBlock= dim[0] * imageoutput->GetScalarSize();
+            for (j=1; j<dim[1]; j++)
+	{
+		memcpy(  imageoutput->GetScalarPointer(0,j,0), imageoutput->GetScalarPointer()  , sizeBlock );
+	}		
+	sizeBlock= dim[0] * dim[1] * imageoutput->GetScalarSize();
+	for (k=0; k<dim[2]; k++)
+	{
+		memcpy(  imageoutput->GetScalarPointer(0,0,k) ,imageoutput->GetScalarPointer()  , sizeBlock );
+	} // for k
+	imageoutput->Update();
+	imageoutput->Modified();
+	bbSetOutputOut(imageoutput);
 }	
 	
 //===== 
@@ -161,7 +161,7 @@ void CreateImage::bbUserInitializeProcessing()
 //    but this is where you should allocate the internal/output pointers 
 //    if any 
 
-  	imageoutput = vtkImageData::New();
+//EED 2017-12-7  	imageoutput = vtkImageData::New();
 
 }