image_info.h

/*
 * image_info.h
 *
 *  Created on: Jun 26, 2012
 *      Author: mpetkova
 */
 
#ifndef IMAGE_INFO_H_
#define IMAGE_INFO_H_
 
#include "Kist.h"
#include "pointlist.h"
#include "Tree.h"
//#include <grid_maintenance.h>
 
enum GATE {unchecked,yes,no};
 
 
struct ImageInfo{
 
 ImageInfo();
  ImageInfo(const ImageInfo & rhs){
    imagekist = new Kist<Point>;
    innerborder = new Kist<Point>;
    outerborder = new Kist<Point>;
 
    copy(rhs,true);
  }
 ~ImageInfo();
  
  ImageInfo & operator=(const ImageInfo & rhs){
    if(this == &rhs) return *this;
    copy(rhs,true);
    return *this;
  }
  
  void Empty();
  
  ImageInfo & operator+=(ImageInfo & rhs);
  
  //Point *points;
  //unsigned long Npoints;
  Kist<Point> * imagekist;
  PosType gridrange[3];
  PosType centroid[2];
  PosType area;
  PosType area_error;
  Kist<Point> * innerborder;
  Kist<Point> * outerborder;
  //short Nencircled;
 
  short ShouldNotRefine;
  GATE uniform_mag;
 
  unsigned long getNimagePoints(){return imagekist->Nunits();}
  KappaType aveTimeDelay();
  KappaType aveInvMag();
  RAY closestRay(const Point_2d &y);
  Point * closestPoint(const Point_2d &y);
  RAY highestSurfaceBrightnessRay();
  void PrintImageInfo();
  void copy(const ImageInfo & image,bool copykists = true);
  
  bool IsMergedImages();
 
  void ArcInfo(PosType *area,PosType *area_circ,PosType theta);
  void FindArc(PosType &radius,PosType *xc,PosType *arc_c,PosType &arclength,PosType &width
            ,PosType resolution,PosType threshold);
  bool constant(LensingVariable lenvar,PosType tol);
  
  PosType ConcaveHullImageArea(bool useborder);
  PosType ConcaveHullSourceArea(bool useborder);
};
 
typedef struct OldImageInfo{
 
 OldImageInfo();
 ~OldImageInfo();
 
  Point *points;
  unsigned long Npoints;
  //Kist<Point> * imagekist;
  PosType gridrange[3];
  PosType centroid[2];
  PosType area;
  PosType area_error;
  Kist<Point> * innerborder;
  Kist<Point> * outerborder;
  short ShouldNotRefine;
 
} OldImageInfo;
 
 
#endif /* IMAGE_INFO_H_ */