/*=========================================================================

  */

#ifndef __axisExtractor_h
#define __axisExtractor_h



#include <list>
#include <stack>
#include <time.h>
#include <stdio.h>
#include <sys/types.h>
#include <sys/timeb.h>
#include <string.h>



#include <vtkPolyDataSource.h>
#include <vtkImageThreshold.h>
#include <vtkImageCast.h>
#include <vtkImageSeedConnectivity.h>
#include <vtkImageData.h>
#include <vtkMarchingCubes.h>
#include <vtkDoubleArray.h>
#include <vtkPointData.h>
#include <vtkExtractVOI.h>
#include <vtkPoints.h>
#include <vtkCellArray.h>
#include <vtkImageResample.h>





class  axisExtractor : public vtkPolyDataSource
{
public:
	static axisExtractor *New();
	void PrintSelf(ostream& os, vtkIndent indent);

	virtual void SetInput(vtkImageData *input);
	vtkImageData *GetInput();

	void SetMaxPropRadio(double value);
	double GetMaxPropRadio();

	void SetHumbral(double value);
	double GetHumbral();

	void SetMaxPropMasa(double value);
	double GetMaxPropMasa();

	void SetMinPropMasa(double value);
	double GetMinPropMasa();

	void SetPoint(double puntoactualprov[3] );

	vtkImageData *GetVolumen();

  
  
protected:  
	axisExtractor();
	
	~axisExtractor() {};

	void Execute();
 
  

private:
	axisExtractor(const axisExtractor&);  // Not implemented.

	void operator=(const axisExtractor&);  // Not implemented.
  
	void searc(int i, int j, int k );

	void find_components( );

	unsigned short maximo();

	void blanquear();

	void blanquear2();

	
	
	void corte(double punto1[3], double punto2[3], double punto3[3], double centro[3],  double radio );

	void avanzar();

	
	void redondear(vtkImageData *data );

	
	
	void copiar(vtkImageData *data, vtkImageData *data2 );

	double distancia(double a[3], double b[3] );	
	
	int envolumen(int a[3], vtkImageData *datae );

	void realtoreal(double a[3], double b[3] );


	void realtoreal2(double a[3], double b[3] );

	void realtoindex(double a[3], int b[3] );



	void indextoreal(int a[3], double b[3] );


	void indextoreal(double a[3], double b[3] );

	



	vtkImageData *dataprov;
	vtkImageData *datatotal;
	
	unsigned char label;
	unsigned char label2;

	unsigned long vector[50][4] ;
	unsigned long vectorb[50][4] ;

	vtkExtractVOI *extrac;
	vtkImageThreshold *thresh;
	vtkImageCast *cast;
	vtkImageSeedConnectivity *connect;

	vtkImageResample *resample;
	
	int iter;

	vtkPoints *points;

	vtkCellArray *lineas; 

	double humbral;
	double maxpropradio;
	double maxpropmasa;
	double minpropmasa; 

	int buenos;

	std::stack< double  >						m_Stack0;
	std::stack< double  >						m_Stack1;
	std::stack< double  >						m_Stack2;
	std::stack< double  >						m_Stack3;
	std::stack< double  >						m_Stack4;
	std::stack< double  >						m_Stack5;
	std::stack< int  >						m_Stack;

	int flagg;
	
};

#endif