/*  * This file is part of the UCB release of Plan 9. It is subject to the license * terms in the LICENSE file found in the top-level directory of this * distribution and at http://akaros.cs.berkeley.edu/files/Plan9License. No * part of the UCB release of Plan 9, including this file, may be copied, * modified, propagated, or distributed except according to the terms contained * in the LICENSE file. */ #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 */ /* * 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;}; /* * 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