Ejemplo n.º 1
0
        /// <summary>
        /// Constructor with data to init variables and the name of the passband
        /// </summary>
        /// <param name="row">Field info from query (DataRow)</param>
        /// <param name="band">Passband name (string)</param>
        public Band(DataRow row, string band)
        {
            if ("ugriz".IndexOf(band.ToLower()) < 0)
            {
                throw new ApplicationException("Parameter 'band' is invalid! Should be one of u,g,r,i or z...");
            }
            this.filter = band;
            this.node   = (double)row["node"];
            this.incl   = (double)row["incl"];
            string postfix = "_" + band;

            //this.a = (double)row["a"+postfix];
            //this.b = (double)row["b"+postfix];
            //this.c = (double)row["c"+postfix];
            //this.d = (double)row["d"+postfix];
            //this.e = (double)row["e"+postfix];
            //this.f = (double)row["f"+postfix];
            // this Convert is used because only dr8 gives problem it has float values
            // otherwise simple (double) cast works for all others
            this.a         = Convert.ToDouble(row["a" + postfix]);
            this.b         = Convert.ToDouble(row["b" + postfix]);
            this.c         = Convert.ToDouble(row["c" + postfix]);
            this.d         = Convert.ToDouble(row["d" + postfix]);
            this.e         = Convert.ToDouble(row["e" + postfix]);
            this.f         = Convert.ToDouble(row["f" + postfix]);
            this.pixperdeg = 9088;
            this.wcs       = this.Wcs();
        }
Ejemplo n.º 2
0
        /// <summary>
        /// Get WCS at any pixel in the image
        /// </summary>
        /// <param name="crpix1">pixel coord 1 (double)</param>
        /// <param name="crpix2">pixel coord 2 (double)</param>
        /// <returns>WCS</returns>
        public WCS Wcs(double crpix1, double crpix2)
        {
            SdssCoord[] corner;
            center = CoordsAtPix(crpix1, crpix2);
            corner = CornerCoords();
            double mu, nu, ra, dec;

            // build linear eqns and solve for affine transform
            double mumu = 0, nunu = 0, munu = 0, ramu = 0, ranu = 0, decmu = 0, decnu = 0;

            for (int i = 0; i < corner.Length; i++)
            {
                mu  = corner[i].mu - center.mu;
                nu  = corner[i].nu - center.nu;
                ra  = corner[i].ra - center.ra;
                dec = corner[i].dec - center.dec;
                // increment
                mumu  += mu * mu;
                munu  += mu * nu;
                nunu  += nu * nu;
                ramu  += ra * mu;
                ranu  += ra * nu;
                decmu += dec * mu;
                decnu += dec * nu;
            }

            double p, q, s, r;

            SolveSym2D(mumu, munu, nunu, ramu, ranu, out p, out q);
            SolveSym2D(mumu, munu, nunu, decmu, decnu, out s, out r);

            WCS wcs = new WCS();

            wcs.EPOCH  = "J2000";
            wcs.CTYPE1 = "DEC--TAN";
            wcs.CTYPE2 = "RA--TAN";
            wcs.CUNIT1 = "deg";
            wcs.CUNIT2 = "deg";

            wcs.NAXIS1 = this.npix1;
            wcs.NAXIS2 = this.npix2;

            wcs.CRPIX1 = crpix1 + 0.5;           // different SDSS vs FITS convention (Steve Kent)
            wcs.CRPIX2 = crpix2 + 0.5;           // -"-
            wcs.CRVAL1 = center.dec;
            wcs.CRVAL2 = center.ra;

            wcs.CD1_1 = s * c + r * f;
            wcs.CD1_2 = s * b + r * e;
            wcs.CD2_1 = p * c + q * f;
            wcs.CD2_2 = p * b + q * e;

            return(wcs);
        }