#include "vtkImageReader.h"
#include "vtkImageGradient.h"
#include "vtkImageViewer.h"
#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkActor.h"


int main()
{
	vtkImageReader *reader = vtkImageReader::New();
	reader->SetDataByteOrderToLittleEndian();
	reader->SetDataExtent(0, 63, 0, 63, 0, 93);
	reader->SetFilePrefix("E:/Data/headsq/quarter");
	reader->SetDataMask(0x7fff);

	vtkImageGradient *gradient = vtkImageGradient::New();
	gradient->SetInput(reader->GetOutput());
	gradient->SetDimensionality(3);

	
	vtkImageViewer *viewer = vtkImageViewer::New();
	viewer->SetInput(gradient->GetOutput());
	viewer->SetZSlice(22);
	viewer->SetColorWindow(400);
	viewer->SetColorLevel(0);

	vtkRenderer *ren = vtkRenderer::New();
	vtkRenderWindow *renWin = vtkRenderWindow::New();
	renWin->AddRenderer(ren);
	vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
	iren->SetRenderWindow(renWin);

	reader->Delete();
	gradient->Delete();
	viewer->Delete();
	
	return 0;
}