#ifndef __VTK__VTK__DIJKSTRA__IMAGE__DATA__H
 #define __VTK__VTK__DIJKSTRA__IMAGE__DATA__H
 
 #include "vtkPolyDataSource.h"
 #include "vtkImageData.h"
 #include "vtkPoints.h"
 #include "vtkCellArray.h"
 #include "vtkPolyData.h"

 
 class vtkIntArray;
 class vtkDataArray;
 class vtkPriorityQueue;
 class vtkIdList;
 class vtkDijkstraImageData : public vtkPolyDataSource
 {
   public:
   static vtkDijkstraImageData *New();
   vtkTypeMacro(vtkDijkstraImageData,vtkPolyDataSource);
   void PrintSelf(ostream& os, vtkIndent indent);
   
   vtkSetMacro(SourceID,int);
   vtkGetMacro(SourceID,int);
 
 
   vtkSetMacro(SinkID,int);
   vtkGetMacro(SinkID,int);
 
   vtkGetObjectMacro(Parent,vtkIntArray);
 
   
   
   virtual void SetBoundaryScalars(vtkDataArray*);
   vtkGetObjectMacro(BoundaryScalars,vtkDataArray);
 
   
   vtkSetMacro(NumberOfInputPoints,int);
   vtkGetMacro(NumberOfInputPoints,int);
   
   
   vtkSetMacro(NumberOfGraphNodes,int);
   vtkGetMacro(NumberOfGraphNodes,int);
   
   
   vtkIdList* GetShortestPathIdList();
   
   
   vtkSetMacro(UseInverseDistance,int);
   vtkGetMacro(UseInverseDistance,int);
   
  
   vtkSetMacro(UseInverseSquaredDistance,int);
   vtkGetMacro(UseInverseSquaredDistance,int);
 
  
   vtkSetMacro(UseInverseExponentialDistance,int);
   vtkGetMacro(UseInverseExponentialDistance,int);
 
   vtkSetMacro(UseSquaredDistance,int);
   vtkGetMacro(UseSquaredDistance,int);
 
       
   // Description:
   unsigned long GetMTime();
   
   // path operations
   void InitTraversePath();
   int GetNextPathNode();
   int GetNumberOfPathNodes();
   
 
   void SetInput(vtkImageData*);
   vtkImageData* GetInput();
   
   
   void CreateGraph(vtkImageData *inData);
   void DeleteGraph();
   void FindNeighbors(vtkIdList *list,int id, vtkDataArray *scalars); 
   
   int findClosestPointInGraph(vtkDataArray *scalars,int id,int dim0,int dim1, int dim2); 
   
 
  
   
   void InitSingleSource(int startv);
  
   void RunDijkstra(vtkDataArray *scalars,int source, int sink);
   
   float EdgeCost(vtkDataArray *scalars, int u, int v);
   void BuildShortestPath(int start,int end);
  
         
   
   
   int SourceID;
   
   
   int SinkID;  
   
   
   vtkDataArray *BoundaryScalars;
 
  
   int NumberOfInputPoints;
   int NumberOfGraphNodes;
 
   int UseInverseDistance;
   int UseInverseSquaredDistance;
   int UseInverseExponentialDistance;
   int UseSquaredDistance;
 
  
   
   vtkPriorityQueue *PQ;
 
   
   vtkIdList *ShortestPathIdList;
   
   vtkIntArray *Parent;
   
   vtkIntArray *Visited;  
  
  
   int PathPointer;
  
   int StopWhenEndReached;


   vtkPoints *puntos; 
   vtkCellArray *lineas;
   vtkImageData *input;



  	void init(vtkImageData *inData);

	void printPointData(int pointID);

	float fuerzaAtraccion(int u, float w);
	
 
 protected:
   vtkDijkstraImageData();
   ~vtkDijkstraImageData();
    void Execute();

   
   FILE *logger;
 
 private:
//   vtkDijstraImageData(const vtkDijkstraImageData&);
   void operator=(const vtkDijkstraImageData&);
 };
 
 #endif //__VTK__VTK__DIJKSTRA__IMAGE__DATA__H