//////////////////////////////////////////////////////////////////////////////////
/////////////////////////          Skull Test          ///////////////////////////
//////////////////////////////////////////////////////////////////////////////////



// Include head files
#include "vtkPolyDataReader.h"
#include <string>
#include <iostream>

#include "vtkRadialThicknessCalculate.h"


int main( int argc, char *argv[] )
{
	// class object initialization
	vtkRadialThicknessCalculate radialThicknessCal;

	// read VTK data 1, which is skull inner surface
	vtkPolyDataReader* reader1 = vtkPolyDataReader::New();
	std::string fileName = "skull_inner.vtk";
	std::cout << "Reading " << fileName.c_str() << std::endl; 
	reader1->SetFileName(fileName.c_str());
	reader1->Update();
	vtkPolyData* inner = reader1->GetOutput();

	// read VTK data 2, which is skull outer surface
	vtkPolyDataReader* reader2 = vtkPolyDataReader::New();
	fileName = "skull_outer.vtk";
	std::cout << "Reading " << fileName.c_str() << std::endl; 
	reader2->SetFileName(fileName.c_str());
	reader2->Update();
	vtkPolyData* outer = reader2->GetOutput();


	//////////////////////////////////////////////////////////////////////////////////
	// Example 1: set skull inner surface as the master mesh
	//////////////////////////////////////////////////////////////////////////////////

	printf("\n\n Example 1 \n\n");
	radialThicknessCal.SetMasterMesh(inner);
	radialThicknessCal.SetSupplementaryMesh(outer);
	// set thickness file name as the output
	radialThicknessCal.SetThicknessFileName("skull_inner_thickness");
	radialThicknessCal.StartThicknessCalculate();
	printf("Done!\n\n");


	//////////////////////////////////////////////////////////////////////////////////
	// Example 2: set skull outer surface as the master mesh
	//////////////////////////////////////////////////////////////////////////////////

	printf("\n\n Example 2 \n\n");
	radialThicknessCal.SetMasterMesh(outer);
	radialThicknessCal.SetSupplementaryMesh(inner);
	// set thickness file name as the output
	radialThicknessCal.SetThicknessFileName("skull_outer_thickness");
	radialThicknessCal.StartThicknessCalculate();
	printf("Done!\n\n");

	return 0;
}