Ejemplo n.º 1
0
        /// <summary>
        /// Creates a new FITS image.
        /// </summary>
        /// <param name="File">File backing the image.</param>
        /// <param name="Width">Image width.</param>
        /// <param name="Height">Image height.</param>
        /// <param name="Transform">WCS transformation.</param>
        /// <param name="BitPix">BITPIX value.</param>
        /// <param name="ExtraProperties">Extra image properties to write in the header.</param>
        /// <param name="ReverseAxis">Reverses the order of the axis in the header.</param>
        public FitsImage(FitsFile File, uint Width, uint Height, WCSViaProjection Transform, int BitPix, List <ImageProperties> ExtraProperties = null, bool ReverseAxis = false) : this()
        {
            if (Width > MaxSize || Height > MaxSize)
            {
                throw new FITSFormatException("Image too large for Umbrella2.");
            }
            this.File      = File;
            this.Width     = Width;
            this.Height    = Height;
            this.Transform = Transform;
            BytesPerPixel  = (byte)Math.Abs(BitPix / 8);
            var RW = GetRW(BitPix);

            Reader  = RW.Item1;
            Writer  = RW.Item2;
            RAFirst = !ReverseAxis;
            Header  = GetHeader(BitPix);
            if (ExtraProperties != null)
            {
                foreach (ImageProperties prop in ExtraProperties)
                {
                    foreach (ElevatedRecord er in prop.GetRecords())
                    {
                        Header.Add(er.Name, er);
                    }
                }
            }
            File.SetPrimaryHeaders(Header);
        }
Ejemplo n.º 2
0
        /// <summary>
        /// Retrieves an image from a FITS file.
        /// </summary>
        /// <param name="File">Input file.</param>
        /// <param name="Number">Image number in multi-image (MEF) FITS files.</param>
        /// <param name="SkipWCS">Whether to parse the WCS headers.</param>
        public FitsImage(FitsFile File, int Number = 0, bool SkipWCS = false) :
            this(ParseHeaderTable(Number, Number == 0 ? File.PrimaryTable : File.MEFHeaderTable[Number - 1], SkipWCS))
        {
            try
            {
                /* Computes BytesPerPixel and selects reading/writing functions */
                BytesPerPixel = (byte)Math.Abs((Header["BITPIX"].Int / 8));
                var RW = GetRW(Header["BITPIX"].Int);

                Reader = RW.Item1;
                Writer = RW.Item2;

                /* Loads SWarp scaling to the properties dictionary */

                try { GetProperty <KnownKeywords.SWarpScaling>(); }
                catch { }
            }
            catch (Exception ex) { throw new FITSFormatException("Cannot understand FITS file.", ex); }
            this.File = File;
        }
Ejemplo n.º 3
0
        /// <summary>
        /// Retrieves an image from a FITS file.
        /// </summary>
        /// <param name="File">Input file.</param>
        /// <param name="Number">Image number in multi-image (MEF) FITS files.</param>
        public FitsImage(FitsFile File, int Number = 0) : this()
        {
            ImageNumber = Number;
            if (Number == 0)
            {
                Header = File.PrimaryTable;
            }
            else
            {
                Header = File.MEFHeaderTable[Number - 1];
            }
            try
            {
                /* Parse image size */
                Width  = (uint)Header["NAXIS1"].Int;
                Height = (uint)Header["NAXIS2"].Int;
                if (Width > MaxSize || Height > MaxSize)
                {
                    throw new FITSFormatException("Image too large for Umbrella2.");
                }

                /* Parse axis types and projection algorithm */
                string Axis1     = Header["CTYPE1"].GetFixedString;
                string Axis2     = Header["CTYPE2"].GetFixedString;
                string Algorithm = Axis1.Substring(5, 3);
                string Nm1       = Axis1.Substring(0, 4);
                string Nm2       = Axis2.Substring(0, 4);

                /* Parses the order of the axes and checks consistency of projection algorithm between the axes. */
                if (Nm1.ToUpper() == "RA--" && Nm2.ToUpper() == "DEC-")
                {
                    RAFirst = true;
                }
                else if (Nm1.ToUpper() == "DEC-" && Nm2.ToUpper() == "RA--")
                {
                    RAFirst = false;
                }
                else
                {
                    throw new FITSFormatException("Cannot understand axis format");
                }
                if (Axis2.Substring(5, 3) != Algorithm)
                {
                    throw new Exception("Projection Algorithm Mismatch.");
                }
                /* Computes the linear transformation part of the WCS projection */
                double     RA0  = (RAFirst ? Header["CRVAL1"] : Header["CRVAL2"]).FloatingPoint;
                double     Dec0 = (RAFirst ? Header["CRVAL2"] : Header["CRVAL1"]).FloatingPoint;
                double     X0   = Header["CRPIX1"].FloatingPoint;
                double     Y0   = Header["CRPIX2"].FloatingPoint;
                WCSLinPart linpart;
                if (RAFirst)
                {
                    linpart = new WCSLinPart(Header["CD1_1"].FloatingPoint, Header["CD1_2"].FloatingPoint, Header["CD2_1"].FloatingPoint, Header["CD2_2"].FloatingPoint, X0, Y0);
                }
                else
                {
                    linpart = new WCSLinPart(Header["CD2_1"].FloatingPoint, Header["CD2_2"].FloatingPoint, Header["CD1_1"].FloatingPoint, Header["CD1_2"].FloatingPoint, X0, Y0);
                }


                if (Header["CUNIT1"].GetFixedString != "deg     " || Header["CUNIT2"].GetFixedString != "deg     ")
                {
                    throw new FITSFormatException("Wrong unit types for axes");
                }

                /* Retrieves the projection algorithm */
                WCSProjectionTransform ipt;
                try { ipt = Umbrella2.WCS.Projections.WCSProjections.GetProjectionTransform(Algorithm, RA0 * Math.PI / 180, Dec0 * Math.PI / 180); }
                catch (KeyNotFoundException ex) { throw new FITSFormatException("Cannot understand projection algorithm", ex); }

                Transform = new WCSViaProjection(ipt, linpart);

                /* Computes BytesPerPixel and selects reading/writing functions */
                BytesPerPixel = (byte)Math.Abs((Header["BITPIX"].Int / 8));
                var RW = GetRW(Header["BITPIX"].Int);

                Reader = RW.Item1;
                Writer = RW.Item2;

                /* Loads SWarp scaling to the properties dictionary */
                try { GetProperty <KnownKeywords.SWarpScaling>(); }
                catch { }
            }
            catch (Exception ex) { throw new FITSFormatException("Cannot understand FITS file.", ex); }
            this.File = File;
        }