X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;ds=sidebyside;f=packages%2Fvtk%2Fsrc%2FbbvtkRescaleSlopeIntercept.cxx;h=3f4494eb978a272ea7e8294d3088b14de565e7ac;hb=ce89f878310f2d109f9e71f5632377b7c14a10fe;hp=8e9ea4c06768d31278b6dc695611f917bb9c7e4d;hpb=5aadf601c2cac6a4b998840f66fef9ba26ae67cd;p=bbtk.git

diff --git a/packages/vtk/src/bbvtkRescaleSlopeIntercept.cxx b/packages/vtk/src/bbvtkRescaleSlopeIntercept.cxx
index 8e9ea4c..3f4494e 100644
--- a/packages/vtk/src/bbvtkRescaleSlopeIntercept.cxx
+++ b/packages/vtk/src/bbvtkRescaleSlopeIntercept.cxx
@@ -3,6 +3,8 @@
 //===== 
 #include "bbvtkRescaleSlopeIntercept.h"
 #include "bbvtkPackage.h"
+
+
 namespace bbvtk
 {
 
@@ -28,13 +30,13 @@ void RescaleSlopeIntercept::Process()
 //    bbSetOutputOut( bbGetInputIn() );
 //    std::cout << "Output value = " <<bbGetOutputOut() << std::endl;
   
+	
+	printf("EED  RescaleSlopeIntercept::Process start\n");		
+	
 
 	double range[2];
-	double spc[6];
-	int dim[3];
+	double spc[3];
 	int ext[6];
-	double value;
-	double newValue;
 	double A=bbGetInputA();
 	double B=bbGetInputB();
 	
@@ -53,58 +55,68 @@ void RescaleSlopeIntercept::Process()
 		else if (bbGetInputOutputFormat()=="VTK_FLOAT")			outputformat = VTK_FLOAT;
 		else if (bbGetInputOutputFormat()=="VTK_DOUBLE")		outputformat = VTK_DOUBLE;
 		
-		imageoutput->Initialize();
-		imageoutput->SetScalarType( outputformat );	
+		
 		bbGetInputIn()->GetScalarRange(range);
 		bbGetInputIn()->GetSpacing(spc);
-		imageoutput->SetSpacing( spc );
-		bbGetInputIn()->GetWholeExtent(ext);
-		dim[0]=ext[1]-ext[0]+1;
-		dim[1]=ext[3]-ext[2]+1;
-		dim[2]=ext[5]-ext[4]+1;
-		imageoutput->SetDimensions( dim[0] , dim[1], dim[2] );
-		imageoutput->AllocateScalars();
-
+		bbGetInputIn()->GetExtent( ext );				
+		
 		double range0plusrange1 = range[0]+range[1];
-		double difrange = range[1]-range[0];
-		double newdifrange = B-A;
+		double difrange			= range[1]-range[0];
+		double newdifrange		= B-A;
 		
-		int i,j,k;
-		for (i=0; i<dim[0]; i++)
-		{
-			for (j=0; j<dim[1]; j++)
-			{
-				for (k=0; k<dim[2]; k++)
-				{
-					value = bbGetInputIn()->GetScalarComponentAsDouble(i,j,k,0);
-
-					if (bbGetInputType()==0) // Rescale Slope/Intercept
-					{ 
-						newValue = value*A+B;
-					} // if type 0
-
-					if (bbGetInputType()==1) // Invert
-					{ 
-						newValue = value*(-1) + (range0plusrange1);
-					} // if type 1
-					
-					if (bbGetInputType()==2) // Redimension
-					{ 
-						newValue = ((value-range[0])/difrange)*newdifrange + A;
-					} // if type 2
-
-					if (bbGetInputType()==3) // Invert redimension
-					{ 
-						newValue = ((value*(-1)+range[1])/difrange)*newdifrange + A;
-					} // if type 3
-					
-					imageoutput->SetScalarComponentFromDouble (i,j,k, 0, newValue);
-				} // for k
-			} // for j
-		} // for i
-			
-		bbSetOutputOut(imageoutput);
-	} // if
+		double shift;
+		double scale;
+		
+		if (bbGetInputType()==0) // Rescale Slope/Intercept
+		{ 
+			//			newValue = value*A+B;
+			scale = A;
+			shift = B;
+		} // if type 0
+		
+		if (bbGetInputType()==1) // Invert
+		{ 
+			//			newValue = value*(-1) + (range0plusrange1);
+			scale=-1;
+			shift=range0plusrange1;
+		} // if type 1
+		
+		if (bbGetInputType()==2) // Redimension
+		{ 
+			//			newValue = ((value-range[0])/difrange)*newdifrange + A;
+			scale=newdifrange/difrange;
+			shift= A - (range[0]*newdifrange/difrange);
+		} // if type 2
+		
+		if (bbGetInputType()==3) // Invert redimension
+		{ 
+			//			newValue = ((value*(-1)+range[1])/difrange)*newdifrange + A;
+			scale=-newdifrange/difrange;
+			shift= A + (range[1]*newdifrange/difrange);
+		} // if type 3
+		
+		filter1->SetInput( bbGetInputIn() );
+		filter1->SetOutputScalarType( VTK_DOUBLE );
+		filter1->SetShift( 0 );
+		filter1->SetScale( scale );
+		
+		filter2->SetInput( filter1->GetOutput() );
+		filter2->SetOutputScalarType( outputformat );
+		filter2->SetShift( shift );
+		filter2->SetScale( 1 );
+		
+		mchange->SetInput( filter2->GetOutput() );
+		mchange->SetExtentTranslation( -ext[0], -ext[2], -ext[4] );
+		mchange->SetOutputSpacing (spc);
+		mchange->Update();    //important
+		bbSetOutputOut( mchange->GetOutput() );
+	 	
+	} else  {
+		bbSetOutputOut( NULL );
+	}		// if bbGetInputIn
+	
+	printf("EED  RescaleSlopeIntercept::Process end\n");		
+	
 }
 	
 //===== 
@@ -120,7 +132,10 @@ void RescaleSlopeIntercept::bbUserSetDefaultValues()
 	bbSetInputA(1);
 	bbSetInputB(0);
 	bbSetInputOutputFormat("SAME");
-	imageoutput=NULL;
+
+	filter1	= NULL;
+	filter2	= NULL;
+	mchange	= NULL;
 }
 //===== 
 // Don't edit this file. This file is generated from xml description.. 
@@ -133,7 +148,10 @@ void RescaleSlopeIntercept::bbUserInitializeProcessing()
 //    but this is where you should allocate the internal/output pointers 
 //    if any 
 
-	imageoutput = vtkImageData::New();
+	filter1	= vtkImageShiftScale::New();
+	filter2	= vtkImageShiftScale::New();
+	mchange = vtkImageChangeInformation::New();
+
   
 }
 //===== 
@@ -146,8 +164,11 @@ void RescaleSlopeIntercept::bbUserFinalizeProcessing()
 //    Here does nothing 
 //    but this is where you should desallocate the internal/output pointers 
 //    if any
-	imageoutput->Delete();
+	mchange->Delete();		
+	filter1->Delete();
+	filter2->Delete();
 }
+	
 }
 // EO namespace bbvtk