ref: af23bb343a067d5346696b788f4cce5b744cd51d
dir: /sys/src/cmd/scat/sky.h/
#define DIR "/lib/sky" /* * This code reflects many years of changes. There remain residues * of prior implementations. * * Keys: * 32 bits long. High 26 bits are encoded as described below. * Low 6 bits are types: * * Patch is ~ one square degree of sky. It points to an otherwise * anonymous list of Catalog keys. The 0th key is special: * it contains up to 4 constellation identifiers. * Catalogs (SAO,NGC,M,...) are: * 31.........8|76|543210 * catalog # |BB|catalog name * BB is two bits of brightness: * 00 -inf < m <= 7 * 01 7 < m <= 10 * 10 10 < m <= 13 * 11 13 < m < inf * The BB field is a dreg, and correct only for SAO and NGC. * IC(n) is just NGC(n+7840) * Others should be self-explanatory. * * Records: * * Star is an SAOrec * Galaxy, PlanetaryN, OpenCl, GlobularCl, DiffuseN, etc., are NGCrecs. * Abell is an Abellrec * The Namedxxx records hold a name and a catalog entry; they result from * name lookups. */ typedef enum { Planet, Patch, SAO, NGC, M, Constel_deprecated, Nonstar_deprecated, NamedSAO, NamedNGC, NamedAbell, Abell, /* NGC types */ Galaxy, PlanetaryN, OpenCl, GlobularCl, DiffuseN, NebularCl, Asterism, Knot, Triple, Double, Single, Uncertain, Nonexistent, Unknown, PlateDefect, /* internal */ NGCN, PatchC, NONGC, }Type; 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, }; #define UNKNOWNMAG 32767 #define NPlanet 20 typedef float Angle; /* in radians */ typedef long DAngle; /* on disk: in units of milliarcsec */ typedef short Mag; /* multiplied by 10 */ typedef long Key; /* known to be 4 bytes, unfortunately */ #pragma pack on /* * All integers are stored in little-endian order. */ typedef struct NGCrec NGCrec; struct NGCrec{ DAngle ra; DAngle dec; DAngle dummy1; /* compatibility with old RNGC version */ DAngle diam; Mag mag; short ngc; /* if >NNGC, IC number is ngc-NNGC */ char diamlim; char type; char magtype; char dummy2; char desc[52]; /* 0-terminated Dreyer description */ }; typedef struct Abellrec Abellrec; struct Abellrec{ DAngle ra; DAngle dec; DAngle glat; DAngle glong; Mag mag10; /* mag of 10th brightest cluster member; in same place as ngc.mag*/ short abell; DAngle rad; short pop; short dist; char distgrp; char richgrp; char flag; char pad; }; typedef struct Planetrec Planetrec; struct Planetrec{ DAngle ra; DAngle dec; DAngle az; DAngle alt; DAngle semidiam; double phase; char name[16]; }; /* * Star names: 0,0==unused. Numbers are name[0]=1,..,99. * Greek letters are alpha=101, etc. * Constellations are alphabetical order by abbreviation, and=1, etc. */ typedef struct SAOrec SAOrec; struct SAOrec{ DAngle ra; DAngle dec; DAngle dra; DAngle ddec; Mag mag; /* visual */ Mag mpg; char spec[3]; char code; char compid[2]; char hdcode; char pad1; long hd; /* HD catalog number */ char name[3]; /* name[0]=alpha name[1]=2 name[3]=ori */ char nname; /* number of prose names */ /* 36 bytes to here */ }; typedef struct Mindexrec Mindexrec; struct Mindexrec{ /* code knows the bit patterns in here; this is a long */ char m; /* M number */ char dummy; short ngc; }; typedef struct Bayerec Bayerec; struct Bayerec{ long sao; char name[3]; char pad; }; #pragma pack off /* * Internal form */ typedef struct Namedrec Namedrec; struct Namedrec{ char name[36]; }; typedef struct Namerec Namerec; struct Namerec{ long sao; long ngc; long abell; char name[36]; /* null terminated */ }; typedef struct Patchrec Patchrec; struct Patchrec{ int nkey; long key[60]; }; typedef struct Record Record; struct Record{ Type type; long index; union{ SAOrec sao; NGCrec ngc; Abellrec abell; Namedrec named; Patchrec patch; Planetrec planet; /* PatchCrec is empty */ }; }; typedef struct Name Name; struct Name{ char *name; int type; }; 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 Pix; typedef struct Img Img; struct Img { int nx; int ny; /* ny is the fast-varying dimension */ Pix a[1]; }; #define RAD(x) ((x)*PI_180) #define DEG(x) ((x)/PI_180) #define ARCSECONDS_PER_RADIAN (DEG(1)*3600) #define MILLIARCSEC (1000*60*60) int nplate; Plate plate[2000]; /* needs to go to 2000 when the north comes */ double PI_180; double TWOPI; double LN2; int debug; struct { float min; float max; float gamma; float absgamma; float mult1; float mult2; int neg; } gam; typedef struct Picture Picture; struct Picture { int minx; int miny; int maxx; int maxy; char name[16]; uchar *data; }; typedef struct Image Image; extern double PI_180; extern double TWOPI; extern char *progname; extern char *desctab[][2]; extern Name names[]; extern Record *rec; extern long nrec; extern Planetrec *planet; /* for bbox: */ extern int folded; extern DAngle ramin; extern DAngle ramax; extern DAngle decmin; extern DAngle decmax; extern Biobuf bout; extern void saoopen(void); extern void ngcopen(void); extern void patchopen(void); extern void mopen(void); extern void constelopen(void); extern void lowercase(char*); extern void lookup(char*, int); extern int typetab(int); extern char*ngcstring(int); extern char*skip(int, char*); extern void prrec(Record*); extern int equal(char*, char*); extern int parsename(char*); extern void radec(int, int*, int*, int*); extern int btag(short); extern long patcha(Angle, Angle); extern long patch(int, int, int); extern char*hms(Angle); extern char*dms(Angle); extern char*ms(Angle); extern char*hm(Angle); extern char*dm(Angle); extern char*deg(Angle); extern char*hm5(Angle); extern long dangle(Angle); extern Angle angle(DAngle); extern void prdesc(char*, char*(*)[2], short*); extern double xsqrt(double); extern Angle dist(Angle, Angle, Angle, Angle); extern Header* getheader(char*); extern char* getword(char*, char*); extern void amdinv(Header*, Angle, Angle, float, float); extern void ppoinv(Header*, Angle, Angle); extern void xypos(Header*, Angle, Angle, float, float); extern void traneqstd(Header*, Angle, Angle); extern Angle getra(char*); extern Angle getdec(char*); extern void getplates(void); extern Img* dssread(char*); extern void hinv(Pix*, int, int); extern int input_bit(Biobuf*); extern int input_nbits(Biobuf*, int); extern int input_huffman(Biobuf*); extern int input_nybble(Biobuf*); extern void qtree_decode(Biobuf*, Pix*, int, int, int, int); extern void start_inputing_bits(void); extern Picture* image(Angle, Angle, Angle, Angle); extern char* dssmount(int); extern int dogamma(Pix); extern void displaypic(Picture*); extern void displayimage(Image*); extern void plot(char*); extern void astro(char*, int); extern char* alpha(char*, char*); extern char* skipbl(char*); extern void flatten(void); extern int bbox(long, long, int); extern int inbbox(DAngle, DAngle); extern char* nameof(Record*); #define NINDEX 400