//**
// WissVis-Aufgabe 3 - Praxis WS 12/13 TUD Lehrstuhl CGV
//**
// mspehr

#include <vtkSLCReader.h>
#include <vtkMarchingCubes.h>
#include <vtkImageGradientMagnitude.h>


#include <vtkDataSetMapper.h>

/*
#include <vtkContourFilter.h>


#include <vtkThreshold.h>
#include <vtkLineSource.h>
#include <vtkGlyph3D.h>*/

#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkInteractorStyleTrackballCamera.h>
#include <vtkPointSource.h>
#include <vtkPolyDataMapper.h>
//---
#include <vtkSmartPointer.h>
// to make the code easier readable:
#define PTR vtkSmartPointer
// note: the usage of this smart pointer is optional; it just lowers the chances for inadvertent heap overflows
//       i.e.              { PTR<vtkClass> var = PTR<vtkClass>::New(); ...; }
//       is equivalent to  { vtkClass* var = vtkClass::New();          ...; var->Delete(); }
//---

#include "SamplePointSource.h"
#include <vtkDelaunay2D.h>
#include <vtkProperty.h>
#include <vtkImageData.h>
#include <vtkPointData.h>
#include <vtkImageMapper.h>
#include "MeshToRGBFilter.h"
#include <vtkVector.h>
#include <vtkCellArray.h>

vtkRenderWindow* createRenderWindowFromMapper(vtkMapper* mapper, bool wireFrame)
{
	PTR<vtkActor> actor = PTR<vtkActor>::New();
    actor->SetMapper(mapper);
	if(wireFrame)
		actor->GetProperty()->SetRepresentationToWireframe();

	PTR<vtkRenderer> renderer = PTR<vtkRenderer>::New();
    renderer->AddActor(actor);

	vtkRenderWindow* window = vtkRenderWindow::New();
    window->AddRenderer(renderer);
	return window; 
}

void doRenderingAndInteraction(vtkRenderWindow* window)
{
    // a) create interactor (with 'trackball' interactor style)
    PTR<vtkRenderWindowInteractor> interactor = PTR<vtkRenderWindowInteractor>::New();
    interactor->SetRenderWindow(window);
    interactor->SetInteractorStyle(PTR<vtkInteractorStyleTrackballCamera>::New());  //< (not necessary)

    // b) execute render/interaction loop
    interactor->Initialize();
    interactor->Start();
}

struct Triangle
{
	vtkVector2<double> p1;
	vtkVector2<double> p2;
	vtkVector2<double> p3;
};

void calcBarycentricCoords(vtkVector2<double> p, Triangle t, double coords[3])
{
	coords[0] = (p.X() - t.p2.X() - ((p.Y()-t.p2.Y())*(t.p3.X()-t.p2.X()))/(t.p3.Y()-t.p2.Y()))/
		(t.p1.X()-t.p2.X()-((t.p1.Y()-t.p2.Y())*(t.p3.X()-t.p2.X())/(t.p3.Y()-t.p2.Y())));
	coords[2] = (p.Y() - t.p2.Y() - coords[0] * (t.p1.Y() - t.p2.Y())) / (t.p3.Y() - t.p2.Y());
	coords[1] = 1 - coords[0] - coords[2];
}

int main()
{
    // -- begin of basic visualization network definition --
    // 1. create source
	PTR<SamplePointSource> source = PTR<SamplePointSource>::New();

	PTR<vtkDelaunay2D> delaunay = PTR<vtkDelaunay2D>::New();
	delaunay->SetInputConnection(source->GetOutputPort());
	delaunay->Update();
	
	MeshToRGBFilter* rgbFilter = new MeshToRGBFilter();
	rgbFilter->SetInputConnection(delaunay->GetOutputPort());
	rgbFilter->Update();

	
		
	PTR<vtkPolyDataMapper> mapper = PTR<vtkPolyDataMapper>::New();
	mapper->SetInputConnection(delaunay->GetOutputPort());
	vtkRenderWindow* window = createRenderWindowFromMapper(mapper, true);
	//doRenderingAndInteraction(window);
	window->Delete();

	PTR<vtkDataSetMapper> imMapper = PTR<vtkDataSetMapper>::New();
	imMapper->SetInputConnection(rgbFilter->GetOutputPort());
	//imMapper->SetInput(image);
	vtkRenderWindow* imWindow = createRenderWindowFromMapper(imMapper, false);
	doRenderingAndInteraction(imWindow);
	imWindow->Delete();
 	
    return 0;
}