#!/bin/sh
# Beispiel für VTK, Carsten Zerbst zerbst@tu-harburg.de
# \
exec vtk $0 $@

source vtkInt.tcl

# read in some structured points
#
vtkPNMReader reader
  reader SetFileName elev.ppm

vtkImageLuminance luminance
  luminance SetInput [reader GetOutput]

vtkStructuredPointsGeometryFilter geometry
  geometry SetInput [luminance GetOutput]

  
vtkWarpScalar warp
  warp SetInput [geometry  GetOutput]
  warp SetScaleFactor 0.1


vtkTriangleFilter tri
tri SetInput [warp GetPolyDataOutput]

vtkDecimate deci; 
    deci SetInput [tri GetOutput]
    deci SetTargetReduction 0.5
    deci SetAspectRatio 20
    deci SetInitialError 0.0002
    deci SetErrorIncrement 0.0005
    deci SetMaximumIterations 6
    deci SetInitialFeatureAngle 45


#Eine Farbtabelle mit den Farben von rot nach blau
vtkLookupTable lut
    lut SetNumberOfColors 50
    lut SetHueRange 0.3333 1

#
# use merge to put back scalars from image file
#
vtkMergeFilter merge
  merge SetGeometry [deci GetOutput]
  merge SetScalars  [reader GetOutput]

vtkDataSetMapper mapper
  mapper SetInput [merge GetOutput]
  mapper SetScalarRange 0 1000
#  mapper ImmediateModeRenderingOff
  mapper SetLookupTable lut
vtkActor actor
  actor SetMapper mapper
  
#[actor GetProperty] SetRepresentationToWireframe
  

# Create renderer stuff
#
vtkRenderer ren1
vtkRenderWindow renWin
    renWin AddRenderer ren1
vtkRenderWindowInteractor iren
    iren SetRenderWindow renWin

# Add the actors to the renderer, set the background and size
#
ren1 AddActor actor
[ren1 GetActiveCamera] Azimuth 20
[ren1 GetActiveCamera] Elevation 30
ren1 SetBackground 1 1 1 
ren1 ResetCameraClippingRange

renWin SetSize 450 450

# render the image
#
iren SetUserMethod {wm deiconify .vtkInteract}
set cam1 [ren1 GetActiveCamera]
$cam1 Zoom 10
renWin Render

renWin SetFileName "tux.ppm"
renWin SaveImageAsPPM

# prevent the tk window from showing up then start the event loop
wm withdraw .