Hi,<br><br>I&#39;m rendering elevation data. Originally data is a structured grid and very dense. But looks fine. When I try to reduce number of polygons using vtkQuadricDecimation I loose all colors and get a gray scale image. I wonder why? What does decimation do to the data so that mapper cannot color it any more? <br>
<br>Btw. Using vtkDecimatePro solves the coloring problem. But the question remains: why colors disappear when I use vtkQuadricDecimation? Why? <br><br>The pipeline I use is below. (No smart pointers, I know. It&#39;s a sample and easier to read, when not using smart pointers). It works just fine if I skip quadric decimation.<br>
<br>Thanks,<br><br>Olli<br><br><br>    vtkImageImport* reader = vtkImageImport::New();<br>        // import settings<br>        // ...<br>    reader-&gt;Update();<br><br>    vtkDataSetSurfaceFilter* filter = vtkDataSetSurfaceFilter::New();<br>
    filter-&gt;SetInputConnection(reader-&gt;GetOutputPort());<br>    reader-&gt;Delete();<br>    filter-&gt;ReleaseDataFlagOn();<br>    filter-&gt;Update();<br><br>    vtkWarpScalar* warp = vtkWarpScalar::New();<br>    warp-&gt;SetInputConnection(filter-&gt;GetOutputPort());<br>
    filter-&gt;Delete();<br>    warp-&gt;SetScaleFactor(3);<br>    warp-&gt;Update();<br><br>    vtkQuadricClustering* cluster = vtkQuadricClustering::New();<br>    cluster-&gt;SetInputConnection(warp-&gt;GetOutputPort());<br>
    warp-&gt;Delete();<br>    cluster-&gt;AutoAdjustNumberOfDivisionsOff();<br>    cluster-&gt;SetNumberOfDivisions(240, 240, 10);<br>    cluster-&gt;PreventDuplicateCellsOn();<br>    cluster-&gt;UseInputPointsOn();<br>    cluster-&gt;Update();<br>
<br>    vtkQuadricDecimation *dec=vtkQuadricDecimation ::New();<br>    dec-&gt;SetInputConnection(cluster-&gt;GetOutputPort());<br>    cluster-&gt;Delete();<br>    dec-&gt;SetTargetReduction(.80);<br>    dec-&gt;Update();<br>
<br>    // color map for terrain<br>    vtkColorTransferFunction* lut = vtkColorTransferFunction::New();<br>    lut-&gt;AddHSVSegment(1, 0.3, 0.7, 0.7, 2000, 0.1, 0.7, 0.7); // terrain<br>    lut-&gt;AddHSVPoint(0, 0.667, 0.7, 1); // sea level<br>
<br>    vtkPolyDataMapper* mapper = vtkPolyDataMapper::New();<br>    mapper-&gt;SetInputConnection(dec-&gt;GetOutputPort());<br>    dec-&gt;Delete();<br>    mapper-&gt;ScalarVisibilityOn();<br>    mapper-&gt;SetScalarRange(0, 1000);<br>
    mapper-&gt;SetLookupTable(lut);<br>    lut-&gt;Delete();<br>    mapper-&gt;Update();<br>   ...<br>