[vtkusers] problem with isosurface and colors
Julia Faurie
julia at lrde.epita.fr
Mon Jul 22 10:40:26 EDT 2002
hi!
i have an isosurface and i want to assign a color to each point of the
isosurface according to his coordinates (the three coordinates are the
red green blue components)
So i used a look up table (in my vtk file i associated scalars value to
each point, and the data is an colors histogram so every point has a
different scalar value)
I get an isosurface with rights colors but there are some blobs on the
surface which have wrong colors and i can't find why...
So if anyone knows what i'm doing wrong or if anyone knows a easier way
to do that...
thanks
my code:
//reader for the isosurface
vtkStructuredPointsReader *reader =
vtkStructuredPointsReader::New();
reader->SetFileName(nomfichier);
reader->SetScalarsName("scalars");
reader->Update();
//Reader fot the colors (scalars values)
vtkStructuredPointsReader *coul =
vtkStructuredPointsReader::New();
coul->SetFileName(nomfichier);
coul->SetScalarsName("scaltab");
coul->Update();
//creation of the look up table
vtkLookupTable *clut=vtkLookupTable::New();
clut->SetNumberOfColors(hsize*hsize*hsize);
clut->Build();
int numerop=0;
for (int r = 0; r < hsize; ++r)
for (int g = 0; g < hsize; ++g)
{
for (int b = 0; b < hsize; ++b)
{
//division to have datas between 0 and 1
float re=((float)r)/248;
float gr=((float)g)/248;
float bl=((float)b)/248;
clut->SetTableValue(numerop,re,gr,bl,1.0);
numerop++;
}
}
...
and later:
//colors+ datas
vtkProbeFilter *probe2=vtkProbeFilter::New();
probe2->SetInput (iso->GetOutput());
probe2->SetSource (coul->GetOutput());
vtkCastToConcrete *cast2=vtkCastToConcrete::New();
cast2->SetInput(probe2->GetOutput());
//normals
vtkPolyDataNormals *normals=vtkPolyDataNormals::New();
normals->SetInput(cast2->GetPolyDataOutput());
normals->SetFeatureAngle(45);
//mapper
isoMap->SetInput(normals->GetOutput());
isoMap->ScalarVisibilityOn();
isoMap->SetScalarRange (scalrang[0],scalrang[1]);
isoMap->SetColorModeToMapScalars();
isoMap->SetLookupTable(lut);
More information about the vtkusers
mailing list