X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=lib%2FcreaVtk%2FvtkVectorsTensorsVisuBase.cpp;h=9c53d977ce873a530fcb5abf4e2f8794beef7d9d;hb=refs%2Fheads%2Fvtk9itk5wx3-macos;hp=8da1d1a0cf5854268dc2b5e91ceb2d84de7b2200;hpb=7b056a20a4768265dac820ba75338779e04c787d;p=creaVtk.git

diff --git a/lib/creaVtk/vtkVectorsTensorsVisuBase.cpp b/lib/creaVtk/vtkVectorsTensorsVisuBase.cpp
index 8da1d1a..9c53d97 100644
--- a/lib/creaVtk/vtkVectorsTensorsVisuBase.cpp
+++ b/lib/creaVtk/vtkVectorsTensorsVisuBase.cpp
@@ -105,41 +105,24 @@ void vtkLookupTableMapDirVectorEED(vtkLookupTableDirectionVector *self, T *input
                           int inIncr, int outFormat)
 {
 	double tmp,sum;
-//	double *mag;
 	int i, j;
 	double dirx,diry,dirz;
-	
-	
-inIncr=3;
-//	printf("EED length %d  %p\n", length,input);
-//	mag = new double[length];
+	double angle;
 	for (i = 0; i < length; ++i)
     {
-		dirx	= 0;
-		diry	= 0;
-		dirz	= 0;
-		sum	= 0;
-		for (j = 0; j < inIncr; ++j)
-		{
-			if (j==0) dirx= static_cast<T>(*input);  
-			if (j==1) diry= static_cast<T>(*input);  
-			if (j==2) dirz= static_cast<T>(*input);  
-			tmp = static_cast<T>(*input);  
-			sum += (tmp * tmp);
-			++input;
-		}
-		sum=sqrt(sum);
-		*output++ = (unsigned char) abs( (255*dirx/sum) );
-		*output++ = (unsigned char) abs( (255*diry/sum) );
-		*output++ = (unsigned char) abs( (255*dirz/sum) );
-		*output++ = 255;
-//		printf("%d %d %d   ",(int)(255*dirx/sum),(int)(255*diry/sum),(int)(255*dirz/sum));
-//		printf(" C     %d        %f %f %f   \n",inIncr,dirx,diry,dirz);
-    }
-	
-//	vtkLookupTableMapData(self, mag, output, length, 1, outFormat);
-	
-//	delete [] mag;
+			dirx	= static_cast<T>(input[0]);
+			diry	= static_cast<T>(input[1]);
+			dirz	= static_cast<T>(input[2]);
+			input	= input+inIncr;
+			sum		= sqrt( dirx*dirx + diry*diry + dirz*dirz );
+			*output++ = (unsigned char) abs( (255*dirx/sum) );
+			*output++ = (unsigned char) abs( (255*diry/sum) );
+			*output++ = (unsigned char) abs( (255*dirz/sum) );
+			*output++ = 255;
+	//		printf("%d %d %d   ",(int)(255*dirx/sum),(int)(255*diry/sum),(int)(255*dirz/sum));
+	//		printf(" C     %d        %f %f %f   \n",inIncr,dirx,diry,dirz);
+
+	} // for
 }
 
 
@@ -151,7 +134,6 @@ void vtkLookupTableDirectionVector::MapScalarsThroughTable2(void *input,
                                              int inputIncrement,
                                              int outputFormat)
 {
-//  printf("vtkLookupTableEED::MapScalarsThroughTable2 inputIncrement=%d   inputDataType=%d\n",inputIncrement,inputDataType);
 // if (this->UseMagnitude && inputIncrement > 1)
 // {
      switch (inputDataType)
@@ -196,7 +178,8 @@ vtkVectorsTensorsVisuBase::vtkVectorsTensorsVisuBase()
 	_typeForm		= 0;
 	_orientation	= 0;
 
-	_LutEED 	= vtkLookupTableDirectionVector::New();
+	_LutEED 		= vtkLookupTableDirectionVector::New();
+	_externalLut 	= NULL;
 	//_LutEED->SetVectorMode(0);
 	//_LutEED->SetVectorModeToMagnitude();
 	//_LutEED->SetVectorModeToComponent();
@@ -300,7 +283,6 @@ void vtkVectorsTensorsVisuBase::Process()    // virtual
 //------------------------------------------------------------------------------
 void vtkVectorsTensorsVisuBase::VisibilityActor()
 {
-
 	if ( (_active==true) && (_actorAdded==false) ){
 		if (GetRenderer()!=NULL)
 		{
@@ -351,4 +333,9 @@ int vtkVectorsTensorsVisuBase::GetOrientation()
 	return _orientation;
 }
 
+//------------------------------------------------------------------------------
+void vtkVectorsTensorsVisuBase::SetExternalLut(vtkScalarsToColors* lut)
+{
+	_externalLut=lut;
+}