Subversion Repositories planix.SVN

#define RAD(x)  ((x)*PI_180)
#define DEG(x)  ((x)/PI_180)
#define ARCSECONDS_PER_RADIAN   (DEG(1)*3600)
#define input_nybble(infile)    input_nbits(infile,4)

typedef float   Angle;  /* in radians */

enum
{
        /*
         * parameters for plate
         */
        Pppo1   = 0,
        Pppo2,
        Pppo3,
        Pppo4,
        Pppo5,
        Pppo6,
        Pamdx1,
        Pamdx2,
        Pamdx3,
        Pamdx4,
        Pamdx5,
        Pamdx6,
        Pamdx7,
        Pamdx8,
        Pamdx9,
        Pamdx10,
        Pamdx11,
        Pamdx12,
        Pamdx13,
        Pamdx14,
        Pamdx15,
        Pamdx16,
        Pamdx17,
        Pamdx18,
        Pamdx19,
        Pamdx20,
        Pamdy1,
        Pamdy2,
        Pamdy3,
        Pamdy4,
        Pamdy5,
        Pamdy6,
        Pamdy7,
        Pamdy8,
        Pamdy9,
        Pamdy10,
        Pamdy11,
        Pamdy12,
        Pamdy13,
        Pamdy14,
        Pamdy15,
        Pamdy16,
        Pamdy17,
        Pamdy18,
        Pamdy19,
        Pamdy20,
        Ppltscale,
        Pxpixelsz,
        Pypixelsz,
        Ppltra,
        Ppltrah,
        Ppltram,
        Ppltras,
        Ppltdec,
        Ppltdecd,
        Ppltdecm,
        Ppltdecs,
        Pnparam,
};

typedef struct  Plate   Plate;
struct  Plate
{
        char    rgn[7];
        char    disk;
        Angle   ra;
        Angle   dec;
};

typedef struct  Header  Header;
struct  Header
{
        float   param[Pnparam];
        int     amdflag;

        float   x;
        float   y;
        float   xi;
        float   eta;
};
typedef long    Type;

typedef struct  Image   Image;
struct  Image
{
        int     nx;
        int     ny;     /* ny is the fast-varying dimension */
        Type    a[1];
};

int     nplate;
Plate   plate[2000];            /* needs to go to 2000 when the north comes */
double  PI_180;
double  TWOPI;
int     debug;
struct
{
        float   min;
        float   max;
        float   del;
        double  gamma;
        int     neg;
} gam;

char*   hms(Angle);
char*   dms(Angle);
double  xsqrt(double);
Angle   dist(Angle, Angle, Angle, Angle);
Header* getheader(char*);
char*   getword(char*, char*);
void    amdinv(Header*, Angle, Angle, float, float);
void    ppoinv(Header*, Angle, Angle);
void    xypos(Header*, Angle, Angle, float, float);
void    traneqstd(Header*, Angle, Angle);
Angle   getra(char*);
Angle   getdec(char*);
void    getplates(void);

Image*  dssread(char*);
void    hinv(Type*, int, int);
int     input_bit(Biobuf*);
int     input_nbits(Biobuf*, int);
void    qtree_decode(Biobuf*, Type*, int, int, int, int);
void    start_inputing_bits(void);
Bitmap* image(Angle, Angle, Angle, Angle);
int     dogamma(int);