1 2 #define DIR "#9/sky" 3 /* 4 * This code reflects many years of changes. There remain residues 5 * of prior implementations. 6 * 7 * Keys: 8 * 32 bits long. High 26 bits are encoded as described below. 9 * Low 6 bits are types: 10 * 11 * Patch is ~ one square degree of sky. It points to an otherwise 12 * anonymous list of Catalog keys. The 0th key is special: 13 * it contains up to 4 constellation identifiers. 14 * Catalogs (SAO,NGC,M,...) are: 15 * 31.........8|76|543210 16 * catalog # |BB|catalog name 17 * BB is two bits of brightness: 18 * 00 -inf < m <= 7 19 * 01 7 < m <= 10 20 * 10 10 < m <= 13 21 * 11 13 < m < inf 22 * The BB field is a dreg, and correct only for SAO and NGC. 23 * IC(n) is just NGC(n+7840) 24 * Others should be self-explanatory. 25 * 26 * Records: 27 * 28 * Star is an SAOrec 29 * Galaxy, PlanetaryN, OpenCl, GlobularCl, DiffuseN, etc., are NGCrecs. 30 * Abell is an Abellrec 31 * The Namedxxx records hold a name and a catalog entry; they result from 32 * name lookups. 33 */ 34 35 typedef enum 36 { 37 Planet, 38 Patch, 39 SAO, 40 NGC, 41 M, 42 Constel_deprecated, 43 Nonstar_deprecated, 44 NamedSAO, 45 NamedNGC, 46 NamedAbell, 47 Abell, 48 /* NGC types */ 49 Galaxy, 50 PlanetaryN, 51 OpenCl, 52 GlobularCl, 53 DiffuseN, 54 NebularCl, 55 Asterism, 56 Knot, 57 Triple, 58 Double, 59 Single, 60 Uncertain, 61 Nonexistent, 62 Unknown, 63 PlateDefect, 64 /* internal */ 65 NGCN, 66 PatchC, 67 NONGC 68 }Type; 69 70 enum 71 { 72 /* 73 * parameters for plate 74 */ 75 Pppo1 = 0, 76 Pppo2, 77 Pppo3, 78 Pppo4, 79 Pppo5, 80 Pppo6, 81 Pamdx1, 82 Pamdx2, 83 Pamdx3, 84 Pamdx4, 85 Pamdx5, 86 Pamdx6, 87 Pamdx7, 88 Pamdx8, 89 Pamdx9, 90 Pamdx10, 91 Pamdx11, 92 Pamdx12, 93 Pamdx13, 94 Pamdx14, 95 Pamdx15, 96 Pamdx16, 97 Pamdx17, 98 Pamdx18, 99 Pamdx19, 100 Pamdx20, 101 Pamdy1, 102 Pamdy2, 103 Pamdy3, 104 Pamdy4, 105 Pamdy5, 106 Pamdy6, 107 Pamdy7, 108 Pamdy8, 109 Pamdy9, 110 Pamdy10, 111 Pamdy11, 112 Pamdy12, 113 Pamdy13, 114 Pamdy14, 115 Pamdy15, 116 Pamdy16, 117 Pamdy17, 118 Pamdy18, 119 Pamdy19, 120 Pamdy20, 121 Ppltscale, 122 Pxpixelsz, 123 Pypixelsz, 124 Ppltra, 125 Ppltrah, 126 Ppltram, 127 Ppltras, 128 Ppltdec, 129 Ppltdecd, 130 Ppltdecm, 131 Ppltdecs, 132 Pnparam 133 }; 134 135 #define UNKNOWNMAG 32767 136 #define NPlanet 20 137 138 typedef float Angle; /* in radians */ 139 typedef int32 DAngle; /* on disk: in units of milliarcsec */ 140 typedef short Mag; /* multiplied by 10 */ 141 typedef int32 Key; /* known to be 4 bytes, unfortunately */ 142 143 /* 144 * All integers are stored in little-endian order. 145 */ 146 typedef struct NGCrec NGCrec; 147 struct NGCrec{ 148 DAngle ra; 149 DAngle dec; 150 DAngle dummy1; /* compatibility with old RNGC version */ 151 DAngle diam; 152 Mag mag; 153 short ngc; /* if >NNGC, IC number is ngc-NNGC */ 154 char diamlim; 155 char type; 156 char magtype; 157 char dummy2; 158 char desc[52]; /* 0-terminated Dreyer description */ 159 }; 160 161 typedef struct Abellrec Abellrec; 162 struct Abellrec{ 163 DAngle ra; 164 DAngle dec; 165 DAngle glat; 166 DAngle glong; 167 Mag mag10; /* mag of 10th brightest cluster member; in same place as ngc.mag*/ 168 short abell; 169 DAngle rad; 170 short pop; 171 short dist; 172 char distgrp; 173 char richgrp; 174 char flag; 175 char pad; 176 }; 177 178 typedef struct Planetrec Planetrec; 179 struct Planetrec{ 180 DAngle ra; 181 DAngle dec; 182 DAngle az; 183 DAngle alt; 184 DAngle semidiam; 185 double phase; 186 char name[16]; 187 }; 188 189 /* 190 * Star names: 0,0==unused. Numbers are name[0]=1,..,99. 191 * Greek letters are alpha=101, etc. 192 * Constellations are alphabetical order by abbreviation, and=1, etc. 193 */ 194 typedef struct SAOrec SAOrec; 195 struct SAOrec{ 196 DAngle ra; 197 DAngle dec; 198 DAngle dra; 199 DAngle ddec; 200 Mag mag; /* visual */ 201 Mag mpg; 202 char spec[3]; 203 char code; 204 char compid[2]; 205 char hdcode; 206 char pad1; 207 int32 hd; /* HD catalog number */ 208 char name[3]; /* name[0]=alpha name[1]=2 name[3]=ori */ 209 char nname; /* number of prose names */ 210 /* 36 bytes to here */ 211 }; 212 213 typedef struct Mindexrec Mindexrec; 214 struct Mindexrec{ /* code knows the bit patterns in here; this is a long */ 215 char m; /* M number */ 216 char dummy; 217 short ngc; 218 }; 219 220 typedef struct Bayerec Bayerec; 221 struct Bayerec{ 222 int32 sao; 223 char name[3]; 224 char pad; 225 }; 226 227 /* 228 * Internal form 229 */ 230 231 typedef struct Namedrec Namedrec; 232 struct Namedrec{ 233 char name[36]; 234 }; 235 236 typedef struct Namerec Namerec; 237 struct Namerec{ 238 int32 sao; 239 int32 ngc; 240 int32 abell; 241 char name[36]; /* null terminated */ 242 }; 243 244 typedef struct Patchrec Patchrec; 245 struct Patchrec{ 246 int nkey; 247 int32 key[60]; 248 }; 249 250 typedef struct Record Record; 251 struct Record{ 252 Type type; 253 int32 index; 254 union{ 255 SAOrec sao; 256 NGCrec ngc; 257 Abellrec abell; 258 Namedrec named; 259 Patchrec patch; 260 Planetrec planet; 261 /* PatchCrec is empty */ 262 } u; 263 }; 264 265 typedef struct Name Name; 266 struct Name{ 267 char *name; 268 int type; 269 }; 270 271 typedef struct Plate Plate; 272 struct Plate 273 { 274 char rgn[7]; 275 char disk; 276 Angle ra; 277 Angle dec; 278 }; 279 280 typedef struct Header Header; 281 struct Header 282 { 283 float param[Pnparam]; 284 int amdflag; 285 286 float x; 287 float y; 288 float xi; 289 float eta; 290 }; 291 typedef int32 Pix; 292 293 typedef struct Img Img; 294 struct Img 295 { 296 int nx; 297 int ny; /* ny is the fast-varying dimension */ 298 Pix a[1]; 299 }; 300 301 #define RAD(x) ((x)*PI_180) 302 #define DEG(x) ((x)/PI_180) 303 #define ARCSECONDS_PER_RADIAN (DEG(1)*3600) 304 #define MILLIARCSEC (1000*60*60) 305 306 int nplate; 307 Plate plate[2000]; /* needs to go to 2000 when the north comes */ 308 double PI_180; 309 double TWOPI; 310 double LN2; 311 int debug; 312 struct 313 { 314 float min; 315 float max; 316 float gamma; 317 float absgamma; 318 float mult1; 319 float mult2; 320 int neg; 321 } gam; 322 323 typedef struct Picture Picture; 324 struct Picture 325 { 326 int minx; 327 int miny; 328 int maxx; 329 int maxy; 330 char name[16]; 331 uchar *data; 332 }; 333 334 #ifndef _DRAW_H_ 335 typedef struct Image Image; 336 #endif 337 338 extern double PI_180; 339 extern double TWOPI; 340 extern char *progname; 341 extern char *desctab[][2]; 342 extern Name names[]; 343 extern Record *rec; 344 extern int32 nrec; 345 extern Planetrec *planet; 346 /* for bbox: */ 347 extern int folded; 348 extern DAngle ramin; 349 extern DAngle ramax; 350 extern DAngle decmin; 351 extern DAngle decmax; 352 extern Biobuf bout; 353 354 extern void saoopen(void); 355 extern void ngcopen(void); 356 extern void patchopen(void); 357 extern void mopen(void); 358 extern void constelopen(void); 359 extern void lowercase(char*); 360 extern void lookup(char*, int); 361 extern int typetab(int); 362 extern char*ngcstring(int); 363 extern char*skip(int, char*); 364 extern void prrec(Record*); 365 extern int equal(char*, char*); 366 extern int parsename(char*); 367 extern void radec(int, int*, int*, int*); 368 extern int btag(short); 369 extern int32 patcha(Angle, Angle); 370 extern int32 patch(int, int, int); 371 extern char*hms(Angle); 372 extern char*dms(Angle); 373 extern char*ms(Angle); 374 extern char*hm(Angle); 375 extern char*dm(Angle); 376 extern char*deg(Angle); 377 extern char*hm5(Angle); 378 extern int32 dangle(Angle); 379 extern Angle angle(DAngle); 380 extern void prdesc(char*, char*(*)[2], short*); 381 extern double xsqrt(double); 382 extern Angle dist(Angle, Angle, Angle, Angle); 383 extern Header* getheader(char*); 384 extern char* getword(char*, char*); 385 extern void amdinv(Header*, Angle, Angle, float, float); 386 extern void ppoinv(Header*, Angle, Angle); 387 extern void xypos(Header*, Angle, Angle, float, float); 388 extern void traneqstd(Header*, Angle, Angle); 389 extern Angle getra(char*); 390 extern Angle getdec(char*); 391 extern void getplates(void); 392 extern Img* dssread(char*); 393 extern void hinv(Pix*, int, int); 394 extern int input_bit(Biobuf*); 395 extern int input_nbits(Biobuf*, int); 396 extern int input_huffman(Biobuf*); 397 extern int input_nybble(Biobuf*); 398 extern void qtree_decode(Biobuf*, Pix*, int, int, int, int); 399 extern void start_inputing_bits(void); 400 extern Picture* image(Angle, Angle, Angle, Angle); 401 extern char* dssmount(int); 402 extern int dogamma(Pix); 403 extern void displaypic(Picture*); 404 extern void displayimage(Image*); 405 extern void plot(char*); 406 extern void astro(char*, int); 407 extern char* alpha(char*, char*); 408 extern char* skipbl(char*); 409 extern void flatten(void); 410 extern int bbox(int32, int32, int); 411 extern int inbbox(DAngle, DAngle); 412 extern char* nameof(Record*); 413 414 #define NINDEX 400 415