# get the interactor ui and colors
source vtkInt.tcl
source "colors.tcl"

# First create the render master
vtkRenderMaster rm;

# Now create the RenderWindow, Renderer and both Actors
#
set renWin [rm MakeRenderWindow];
set ren1   [$renWin MakeRenderer];
set iren [$renWin MakeRenderWindowInteractor];

# create pipeline
#
vtkMCubesReader reader;
    reader SetFilename "../../data/pineRoot/pine_root.tri";
    reader FlipNormalsOff;
    reader DebugOn;
vtkDecimate deci;
    deci SetInput [reader GetOutput];
    deci SetTargetReduction 0.9;
    deci SetAspectRatio 20;
    deci SetInitialError 0.0005;
    deci SetErrorIncrement 0.001;
    deci SetMaximumIterations 6;
    deci SetInitialFeatureAngle 30;
    deci DebugOn;
vtkConnectivityFilter connect;
    connect SetInput [deci GetOutput];
    connect ExtractLargestRegion;
    connect DebugOn;
vtkDataSetMapper isoMapper;
    isoMapper SetInput [connect GetOutput];
    isoMapper ScalarsVisibleOff;
vtkActor isoActor;
    isoActor SetMapper isoMapper;
    eval [isoActor GetProperty] SetColor $raw_sienna;

vtkOutlineFilter outline;
    outline SetInput [reader GetOutput];
vtkPolyMapper outlineMapper;
    outlineMapper SetInput [outline GetOutput];
vtkActor outlineActor;
    outlineActor SetMapper outlineMapper;
    [outlineActor GetProperty] SetColor 0 0 0;

# Add the actors to the renderer, set the background and size
#
$ren1 AddActors outlineActor;
$ren1 AddActors isoActor;
$ren1 SetBackground 1 1 1;
$renWin SetSize 750 750;
eval $ren1 SetBackground $slate_grey;

# render the image
#
$iren SetUserMethod {wm deiconify .vtkInteract};
set cam [$ren1 GetActiveCamera];
  $cam SetFocalPoint 40.6018 37.2813 50.1953
  $cam SetPosition 40.6018 -280.533 47.0172;
  $cam CalcViewPlaneNormal;
  $cam SetClippingRange 26.1073 1305.36;
  $cam SetViewAngle 20.9219;
  $cam SetViewUp 0.0 0.0 1.0
$renWin Render;

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

$iren Start;
