Пример #1
0
        internal Grib1WaveSpectra2DDirFreq(System.IO.Stream raf)
        {
            SupportClass.Skip(raf, 4);

            byte[] ver = GribNumbers.ReadBytes(raf, 4);
            System.Text.Encoding enc = System.Text.Encoding.ASCII;
            version = enc.GetString(ver);

            SupportClass.Skip(raf, 2);

            directionNumber     = raf.ReadByte();
            frequencyNumber     = raf.ReadByte();
            numberOfDirections  = raf.ReadByte();
            numberOfFrequencies = raf.ReadByte();
            dirScaleFactor      = GribNumbers.int4(raf);
            freqScaleFactor     = GribNumbers.int4(raf);

            SupportClass.Skip(raf, 37);

            directions  = new double[numberOfDirections];
            frequencies = new double[numberOfFrequencies];

            for (int i = 0; i < numberOfDirections; i++)
            {
                int val = GribNumbers.int4(raf);
                directions[i] = (double)val / (double)dirScaleFactor;
            }

            for (int i = 0; i < numberOfFrequencies; i++)
            {
                int val = GribNumbers.int4(raf);
                frequencies[i] = (double)val / (double)freqScaleFactor;
            }
        }
Пример #2
0
        // *** constructors *******************************************************

        /// <summary> Constructs a <tt>Grib2LocalUseSection</tt> object from a raf.</summary>
        /// <param name="raf">
        /// </param>
        /// <throws>  IOException  if raf contains no valid GRIB product </throws>
        //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'"
        public Grib2LocalUseSection(System.IO.FileStream raf)
        {
            // octets 1-4 (Length of GDS)
            length = GribNumbers.int4(raf);
            //System.out.println( "LUS length=" + length );

            section = raf.ReadByte();             // This is section 2
            //System.out.println( "LUS Local Use is 2, section=" + section );

            if (section != 2)
            {
                // no local use section
                SupportClass.Skip(raf, -5);
                return;
            }
            else
            {
                //SupportClass.Skip(raf, length - 5);
                bytes = new byte[length - 5];
                int nb = raf.Read(bytes, 0, length - 5);
                if (nb != length - 5)
                {
                    throw new NoValidGribException("Failed to read Local Use Section data");
                }
            }
        }         // end of Grib2LocalUseSection
Пример #3
0
        private int Xlength;       // length of the x axis

        // *** constructors *******************************************************

        /// <summary> Constructor for a Grib2 Data Section.</summary>
        /// <param name="getData">
        /// </param>
        /// <param name="raf">
        /// </param>
        /// <param name="gds">
        /// </param>
        /// <param name="drs">
        /// </param>
        /// <param name="bms">
        /// </param>
        /// <throws>  IOException </throws>
        //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'"
        public Grib2DataSection(bool getData, System.IO.FileStream raf, Grib2GridDefinitionSection gds, Grib2DataRepresentationSection drs, Grib2BitMapSection bms)
        {
            //System.out.println( "raf.FilePointer=" + raf.FilePointer() );
            // octets 1-4 (Length of DS)
            length = GribNumbers.int4(raf);
            //System.out.println( "DS length=" + length );
            //System.out.println( "DS calculated end=" + ( raf.getFilePointer() + length -4 ));
            // octet 5  section 7
            section = raf.ReadByte();
            //System.out.println( "DS is 7, section=" + section );
            if (!getData)
            {
                // skip data read
                //System.out.println( "raf.position before reposition="+raf.getFilePointer());
                //System.out.println( "raf.length=" + raf.length() );
                // sanity check for erronous ds length
                if (length > 0 && length < raf.Length)
                {
                    SupportClass.Skip(raf, length - 5);
                    //System.out.println( "raf.skipBytes = " + (length -5) );
                }
                else
                {
                    length = 5;                     // only read length and section
                }
                //System.out.println( "raf.position after skip=" + raf.getFilePointer() );
                return;
            }
            int dtn = drs.DataTemplateNumber;

            //System.out.println( "DS dtn=" + dtn );
            if (dtn == 0 || dtn == 1)
            {
                // 0: Grid point data - simple packing
                // 1: Matrix values - simple packing
                simpleUnpacking(raf, gds, drs, bms);
            }
            else if (dtn == 2)
            {
                // 2:Grid point data - complex packing
                complexUnpacking(raf, gds, drs);
            }
            else if (dtn == 3)
            {
                // 3: complex packing with spatial differencing
                complexUnpackingWithSpatial(raf, gds, drs);
            }
            else if (dtn == 40 || dtn == 40000)
            {
                // JPEG 2000 Stream Format
                jpeg2000Unpacking(raf, gds, drs, bms);
            }
        }         // end Grib2DataSection
Пример #4
0
        /// <summary> Creates an ensemble object for the product PDS.</summary>
        /// <param name="raf">RandomAccessFile.
        /// </param>
        /// <param name="parameterNumber">
        /// </param>
        /// <throws>  IOException </throws>
        //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'"
        internal Grib1Ensemble(System.IO.Stream raf, int parameterNumber)
        {
            // skip 12 bytes to start of ensemble
            SupportClass.Skip(raf, 12);

            // octet 41 id's ensemble
            int app = raf.ReadByte();

            if (app != 1)
            {
                System.Console.Out.WriteLine("not ensemble product");
                return;
            }

            // octet 42 Type
            eType = raf.ReadByte();

            // octet 43 Identification number
            eId = raf.ReadByte();

            // octet 44 Product Identifier
            eProd = raf.ReadByte();

            // octet 45 Spatial Identifier
            eSpatial = raf.ReadByte();

            if (parameterNumber == 191 || parameterNumber == 192)
            {
                // octet 46 Probability product definition
                epd = raf.ReadByte();

                // octet 47 Probability type
                ept = raf.ReadByte();

                // octet 48-51 Probability lower limit
                epll = GribNumbers.int4(raf);

                // octet 52-55 Probability upper limit
                epul = GribNumbers.int4(raf);

                // octet 56-60 reserved
                int reserved = GribNumbers.int4(raf);

                if (eType == 4 || eType == 5)
                {
                    // octet 61 Ensemble size
                    eSize = raf.ReadByte();

                    // octet 62 Cluster size
                    eCSize = raf.ReadByte();

                    // octet 63 Number of clusters
                    eCNumber = raf.ReadByte();

                    // octet 64 Clustering Method
                    eCMethod = raf.ReadByte();

                    // octet 65-67 Northern latitude of clustering domain
                    nLat = GribNumbers.int3(raf);

                    // octet 68-70 Southern latitude of clustering domain
                    sLat = GribNumbers.int3(raf);

                    // octet 71-73 Eastern latitude of clustering domain
                    eLat = GribNumbers.int3(raf);

                    // octet 74-76 Western latitude of clustering domain
                    wLat = GribNumbers.int3(raf);

                    if (eType == 4)
                    {
                        // octets 77-86 Cluster Membership
                        cm = new sbyte[10];
                        SupportClass.ReadInput(raf, cm, 0, cm.Length);
                    }
                }
            }
        }
Пример #5
0
        // *** constructors *******************************************************

        /// <summary> Constructs a <tt>Grib1ProductDefinitionSection</tt> object from a raf.
        ///
        /// </summary>
        /// <param name="raf">with PDS content
        ///
        /// </param>
        /// <throws>  NotSupportedException  if raf contains no valid GRIB file </throws>
        //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'"
        public Grib1ProductDefinitionSection(System.IO.FileStream raf)
        {
            // octets 1-3 PDS length
            length = (int)GribNumbers.uint3(raf);
            //System.out.println( "PDS length = " + length );

            // Paramter table octet 4
            table_version = raf.ReadByte();

            // Center  octet 5
            center_id = raf.ReadByte();

            // octet 6 Generating Process - See Table A
            process_id = raf.ReadByte();

            // octet 7 (id of grid type) - not supported yet
            grid_id = raf.ReadByte();

            //octet 8 (flag for presence of GDS and BMS)
            int exists = raf.ReadByte();

            //BKSystem.IO.BitStream s = new BKSystem.IO.BitStream(8);
//            s.WriteByte((byte)raf.ReadByte());
//            s.Position = 0;
//            sbyte exists;
//            s.Read(out exists);
//            bms_exists = (exists & 64) == 64;

            gds_exists = (exists & 128) == 128;
            bms_exists = (exists & 64) == 64;

            // octet 9 (parameter and unit)
            parameterNumber = raf.ReadByte();

            // octets 10-12 (level)
            int levelType   = raf.ReadByte();
            int levelValue1 = raf.ReadByte();
            int levelValue2 = raf.ReadByte();

            level = new GribPDSLevel(levelType, levelValue1, levelValue2);

            // octets 13-17 (base time for reference time)
            int year   = raf.ReadByte();
            int month  = raf.ReadByte();
            int day    = raf.ReadByte();
            int hour   = raf.ReadByte();
            int minute = raf.ReadByte();

            // get info for forecast time
            // octet 18 Forecast time unit
            timeUnit = raf.ReadByte();

            switch (timeUnit)
            {
            case 0:                      // minute
                tUnit = "minute";
                break;

            case 1:                      // hours
                tUnit = "hour";
                break;

            case 2:                      // day
                tUnit = "day";
                break;

            case 3:                      // month
                tUnit = "month";
                break;

            case 4:                      //1 year
                tUnit = "1year";
                break;

            case 5:                      // decade
                tUnit = "decade";
                break;

            case 6:                      // normal
                tUnit = "day";
                break;

            case 7:                      // century
                tUnit = "century";
                break;

            case 10:                      //3 hours
                tUnit = "3hours";
                break;

            case 11:                      // 6 hours
                tUnit = "6hours";
                break;

            case 12:                      // 12 hours
                tUnit = "12hours";
                break;

            case 254:                      // second
                tUnit = "second";
                break;


            default:
                System.Console.Error.WriteLine("PDS: Time Unit " + timeUnit + " is not yet supported");
                break;
            }

            // octet 19 & 20 used to create Forecast time
            p1 = raf.ReadByte();
            p2 = raf.ReadByte();

            // octet 21 (time range indicator)
            timeRangeValue = raf.ReadByte();
            // forecast time is always at the end of the range
            //System.out.println( "PDS timeRangeValue =" + timeRangeValue );
            switch (timeRangeValue)
            {
            case 0:
                timeRange    = "product valid at RT + P1";
                forecastTime = p1;
                break;

            case 1:
                timeRange    = "product valid for RT, P1=0";
                forecastTime = 0;
                break;

            case 2:
                timeRange    = "product valid from (RT + P1) to (RT + P2)";
                forecastTime = p2;
                break;

            case 3:
                timeRange    = "product is an average between (RT + P1) to (RT + P2)";
                forecastTime = p2;
                break;

            case 4:
                timeRange    = "product is an accumulation between (RT + P1) to (RT + P2)";
                forecastTime = p2;
                break;

            case 5:
                timeRange    = "product is the difference (RT + P2) - (RT + P1)";
                forecastTime = p2;
                break;

            case 6:
                timeRange    = "product is an average from (RT - P1) to (RT - P2)";
                forecastTime = -p2;
                break;

            case 7:
                timeRange    = "product is an average from (RT - P1) to (RT + P2)";
                forecastTime = p2;
                break;

            case 10:
                timeRange = "product valid at RT + P1";
                // p1 really consists of 2 bytes p1 and p2
                forecastTime = p1 = GribNumbers.int2(p1, p2);
                p2           = 0;
                break;

            case 51:
                timeRange    = "mean value from RT to (RT + P2)";
                forecastTime = p2;
                break;

            default:
                System.Console.Error.WriteLine("PDS: Time Range Indicator " + timeRangeValue + " is not yet supported");
                break;
            }

            // octet 22 & 23
            int avgInclude = GribNumbers.int2(raf);

            // octet 24
            int avgMissing = raf.ReadByte();

            // octet 25
            int century = raf.ReadByte() - 1;

            // octet 26, sub center
            subcenter_id = raf.ReadByte();

            // octets 27-28 (decimal scale factor)
            decscale = GribNumbers.int2(raf);

            refTime       = new DateTime(century * 100 + year, month, day, hour, minute, 0, DateTimeKind.Utc);
            baseTime      = refTime;
            referenceTime = refTime.ToString(dateFormat);

            //System.out.println( "PDS referenceTime raw = " +
            //(century * 100 + year) +"-" +
            //month + "-" + day + "T" + hour +":" + minute +":00Z" );
            //System.out.println( "PDS referenceTime     = " + referenceTime );

            // TODO

            /*
             *          parameter_table = GribPDSParamTable.getParameterTable(center_id, subcenter_id, table_version);
             *          parameter = parameter_table.getParameter(parameterNumber);
             */
            parameter = new Parameter();


            if (center_id == 7 && subcenter_id == 2)
            {
                // ensemble product
                epds = new Grib1Ensemble(raf, parameterNumber);
            }
            // Special handling of 2D Wave Spectra (single)
            else if (table_version == 140 && parameterNumber == 251)
            {
                SupportClass.Skip(raf, 12);
                int extDef = raf.ReadByte(); // Extension definition

                // Read ECMWF extension for 2D wave spectra single
                if (extDef == 13)
                {
                    waveSpectra2DDirFreq = new Grib1WaveSpectra2DDirFreq(raf);
                }
            }
            else if (length != 28)
            {
                // check if all bytes read in section
                //lengthErr = true;
                int extra;
                for (int i = 29; i <= length; i++)
                {
                    extra = raf.ReadByte();
                }
            }
        }         // end Grib1ProductDefinitionSection
Пример #6
0
        // *** constructors *******************************************************

        /// <summary> Constructs a Grib2ProductDefinitionSection  object from a raf.
        ///
        /// </summary>
        /// <param name="raf">RandomAccessFile with PDS content
        ///
        /// </param>
        /// <throws>  IOException  if raf contains no valid GRIB file </throws>
        //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'"
        public Grib2ProductDefinitionSection(System.IO.FileStream raf)
        {
            // octets 1-4 (Length of PDS)
            length = GribNumbers.int4(raf);
            //System.out.println( "PDS length=" + length );

            // octet 5
            section = raf.ReadByte();
            //System.out.println( "PDS is 4, section=" + section );

            // octet 6-7
            coordinates = GribNumbers.int2(raf);
            //System.out.println( "PDS coordinates=" + coordinates );

            // octet 8-9
            productDefinition = GribNumbers.int2(raf);
            //System.out.println( "PDS productDefinition=" + productDefinition );

            switch (productDefinition)
            {
            // Analysis or forecast at a horizontal level or in a horizontal
            // layer at a point in time
            case 0:
            case 1:
            case 2:
            case 3:
            case 4:
            case 5:
            case 6:
            case 7:
            case 8:  {
                // octet 10
                parameterCategory = raf.ReadByte();
                //System.out.println( "PDS parameterCategory=" +
                //parameterCategory );

                // octet 11
                parameterNumber = raf.ReadByte();
                //System.out.println( "PDS parameterNumber=" + parameterNumber );

                // octet 12
                typeGenProcess = raf.ReadByte();
                //System.out.println( "PDS typeGenProcess=" + typeGenProcess );

                // octet 13
                backGenProcess = raf.ReadByte();
                //System.out.println( "PDS backGenProcess=" + backGenProcess );

                // octet 14
                analysisGenProcess = raf.ReadByte();
                //System.out.println( "PDS analysisGenProcess=" +
                //analysisGenProcess );

                // octet 15-16
                hoursAfter = GribNumbers.int2(raf);
                //System.out.println( "PDS hoursAfter=" + hoursAfter );

                // octet 17
                minutesAfter = raf.ReadByte();
                //System.out.println( "PDS minutesAfter=" + minutesAfter );

                // octet 18
                timeRangeUnit = raf.ReadByte();
                //System.out.println( "PDS timeRangeUnit=" + timeRangeUnit );

                // octet 19-22
                forecastTime = GribNumbers.int4(raf);
                //System.out.println( "PDS forecastTime=" + forecastTime );

                // octet 23
                typeFirstFixedSurface = raf.ReadByte();
                //System.out.println( "PDS typeFirstFixedSurface=" +
                //     typeFirstFixedSurface );

                // octet 24
                int scaleFirstFixedSurface = raf.ReadByte();
                //System.out.println( "PDS scaleFirstFixedSurface=" +
                //     scaleFirstFixedSurface );

                // octet 25-28
                int valueFirstFixedSurface = GribNumbers.int4(raf);
                //System.out.println( "PDS valueFirstFixedSurface=" +
                //     valueFirstFixedSurface );

                //UPGRADE_WARNING: Data types in Visual C# might be different.  Verify the accuracy of narrowing conversions. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1042'"
                FirstFixedSurfaceValue = (float)((scaleFirstFixedSurface == 0 || valueFirstFixedSurface == 0)?valueFirstFixedSurface:System.Math.Pow(valueFirstFixedSurface, -scaleFirstFixedSurface));

                // octet 29
                typeSecondFixedSurface = raf.ReadByte();
                //System.out.println( "PDS typeSecondFixedSurface=" +
                //typeSecondFixedSurface );

                // octet 30
                int scaleSecondFixedSurface = raf.ReadByte();
                //System.out.println( "PDS scaleSecondFixedSurface=" +
                //scaleSecondFixedSurface );

                // octet 31-34
                int valueSecondFixedSurface = GribNumbers.int4(raf);
                //System.out.println( "PDS valueSecondFixedSurface=" +
                //valueSecondFixedSurface );

                //UPGRADE_WARNING: Data types in Visual C# might be different.  Verify the accuracy of narrowing conversions. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1042'"
                SecondFixedSurfaceValue = (float)((scaleSecondFixedSurface == 0 || valueSecondFixedSurface == 0)?valueSecondFixedSurface:System.Math.Pow(valueSecondFixedSurface, -scaleSecondFixedSurface));


                // Individual ensemble forecast, control and perturbed, at a
                // horizontal level or in a horizontal layer at a point in time
                if (productDefinition == 1)
                {
                    System.Console.Out.WriteLine("PDS productDefinition == 1 not done");

                    //Derived forecast based on all ensemble members at a horizontal
                    // level or in a horizontal layer at a point in time
                }
                else if (productDefinition == 2)
                {
                    System.Console.Out.WriteLine("PDS productDefinition == 2 not done");

                    // Derived forecasts based on a cluster of ensemble members over
                    // a rectangular area at a horizontal level or in a horizontal layer
                    // at a point in time
                }
                else if (productDefinition == 3)
                {
                    System.Console.Out.WriteLine("PDS productDefinition == 3 not done");

                    // Derived forecasts based on a cluster of ensemble members
                    // over a circular area at a horizontal level or in a horizontal
                    // layer at a point in time
                }
                else if (productDefinition == 4)
                {
                    System.Console.Out.WriteLine("PDS productDefinition == 4 not done");

                    // Probability forecasts at a horizontal level or in a horizontal
                    //  layer at a point in time
                }
                else if (productDefinition == 5)
                {
                    System.Console.Out.WriteLine("PDS productDefinition == 5 not done");

                    // Percentile forecasts at a horizontal level or in a horizontal layer
                    // at a point in time
                }
                else if (productDefinition == 6)
                {
                    System.Console.Out.WriteLine("PDS productDefinition == 6 not done");

                    // Analysis or forecast error at a horizontal level or in a horizontal
                    // layer at a point in time
                }
                else if (productDefinition == 7)
                {
                    System.Console.Out.WriteLine("PDS productDefinition == 7 not done");

                    // Average, accumulation, and/or extreme values at a horizontal
                    // level or in a horizontal layer in a continuous or non-continuous
                    // time interval
                }
                else if (productDefinition == 8)
                {
                    //System.out.println( "PDS productDefinition == 8 " );
                    //  35-41 bytes
                    int year   = GribNumbers.int2(raf);
                    int month  = (raf.ReadByte()) - 1;
                    int day    = raf.ReadByte();
                    int hour   = raf.ReadByte();
                    int minute = raf.ReadByte();
                    int second = raf.ReadByte();
                    //System.out.println( "PDS date:" + year +":" + month +
                    //":" + day + ":" + hour +":" + minute +":" + second );

                    // 42 - 46
                    int timeRanges = raf.ReadByte();
                    //System.out.println( "PDS timeRanges=" + timeRanges ) ;
                    int missingDataValues = GribNumbers.int4(raf);
                    //System.out.println( "PDS missingDataValues=" + missingDataValues ) ;
                    // 47 - 48
                    int outmostTimeRange = raf.ReadByte();
                    //System.out.println( "PDS outmostTimeRange=" + outmostTimeRange )
                    ;
                    int missing = raf.ReadByte();
                    //System.out.println( "PDS missing=" + missing ) ;
                    // 49 - 53
                    int statisticalProcess = raf.ReadByte();
                    //System.out.println( "PDS statisticalProcess=" + statisticalProcess ) ;
                    int timeIncrement = GribNumbers.int4(raf);
                    //System.out.println( "PDS timeIncrement=" + timeIncrement ) ;

                    // 54 - 58
                    int indicatorTR = raf.ReadByte();
                    //System.out.println( "PDS indicatorTR=" + indicatorTR ) ;

                    int lengthTR = GribNumbers.int4(raf);

                    if (timeRanges == 1)
                    {
                        this.forecastTime += calculateIncrement(statisticalProcess, timeIncrement);
                    }
                    //System.out.println( "PDS lengthTR=" + lengthTR ) ;

                    //int indicatorSF = raf.read();
                    //System.out.println( "PDS indicatorSF=" + indicatorSF ) ;

                    //int incrementSF = GribNumbers.int4( raf );
                    //System.out.println( "PDS incrementSF=" + incrementSF ) ;
                }
                break;
            }                             // cases 0-8

            // Radar product

            case 20:  {
                parameterCategory = raf.ReadByte();
                //System.out.println( "PDS parameterCategory=" +
                //parameterCategory );

                parameterNumber = raf.ReadByte();
                //System.out.println( "PDS parameterNumber=" + parameterNumber );

                typeGenProcess = raf.ReadByte();
                //System.out.println( "PDS typeGenProcess=" + typeGenProcess );

                backGenProcess = raf.ReadByte();
                //System.out.println( "PDS backGenProcess=" + backGenProcess );

                hoursAfter = GribNumbers.int2(raf);
                //System.out.println( "PDS hoursAfter=" + hoursAfter );

                minutesAfter = raf.ReadByte();
                //System.out.println( "PDS minutesAfter=" + minutesAfter );

                timeRangeUnit = raf.ReadByte();
                //System.out.println( "PDS timeRangeUnit=" + timeRangeUnit );

                forecastTime = GribNumbers.int4(raf);
                //System.out.println( "PDS forecastTime=" + forecastTime );

                break;
            }                             // case 20

            // Satellite Product

            case 30:  {
                parameterCategory = raf.ReadByte();
                //System.out.println( "PDS parameterCategory=" + parameterCategory );

                parameterNumber = raf.ReadByte();
                //System.out.println( "PDS parameterNumber=" + parameterNumber );

                typeGenProcess = raf.ReadByte();
                //System.out.println( "PDS typeGenProcess=" + typeGenProcess );

                backGenProcess = raf.ReadByte();
                //System.out.println( "PDS backGenProcess=" + backGenProcess );

                nb = raf.ReadByte();
                //System.out.println( "PDS nb =" + nb );
                for (int j = 0; j < nb; j++)
                {
                    SupportClass.Skip(raf, 10);
                }
                break;
            }                             // case 30

            // CCITTIA5 character string

            case 254:  {
                parameterCategory = raf.ReadByte();
                //System.out.println( "PDS parameterCategory=" +
                //parameterCategory );

                parameterNumber = raf.ReadByte();
                //System.out.println( "PDS parameterNumber=" + parameterNumber );

                //numberOfChars = GribNumbers.int4( raf );
                //System.out.println( "PDS numberOfChars=" +
                //numberOfChars );
                break;
            }                             // case 254


            default:
                break;
            }             // end switch
        }
Пример #7
0
        /// <summary> Constructs a Grib1BinaryDataSection object from a raf.
        /// A bit map is defined.
        ///
        /// </summary>
        /// <param name="raf">raf with BDS content
        /// </param>
        /// <param name="decimalscale">the exponent of the decimal scale
        /// </param>
        /// <param name="bms">bit map section of GRIB record
        ///
        /// </param>
        /// <throws>  NotSupportedException  if stream contains no valid GRIB file </throws>
        //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'"
        public Grib1BinaryDataSection(System.IO.FileStream raf, int decimalscale, Grib1BitMapSection bms)
        {
            // octets 1-3 (section length)
            length = (int)GribNumbers.uint3(raf);
            //System.out.println( "BDS length = " + length );

            // octet 4, 1st half (packing flag)
            int unusedbits = raf.ReadByte();

            // TODO Check this!!!
            if ((unusedbits & 192) != 0)
            {
                throw new NotSupportedException("BDS: (octet 4, 1st half) not grid point data and simple packing ");
            }

            // octet 4, 2nd half (number of unused bits at end of this section)
            unusedbits = unusedbits & 15;
            //System.out.println( "BDS unusedbits = " + unusedbits );

            // octets 5-6 (binary scale factor)
            int binscale = GribNumbers.int2(raf);

            // octets 7-10 (reference point = minimum value)
            float refvalue = GribNumbers.float4(raf);

            // octet 11 (number of bits per value)
            int numbits = raf.ReadByte();

            //System.out.println( "BDS numbits = " + numbits );
            if (numbits == 0)
            {
                isConstant = true;
            }
            //System.out.println( "BDS isConstant = " + isConstant );

            // *** read values *******************************************************

            //UPGRADE_WARNING: Data types in Visual C# might be different.  Verify the accuracy of narrowing conversions. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1042'"
            float ref_Renamed = (float)(System.Math.Pow(10.0, -decimalscale) * refvalue);
            //UPGRADE_WARNING: Data types in Visual C# might be different.  Verify the accuracy of narrowing conversions. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1042'"
            float scale = (float)(System.Math.Pow(10.0, -decimalscale) * System.Math.Pow(2.0, binscale));

            if (bms != null)
            {
                bool[] bitmap = bms.Bitmap;

                values = new float[bitmap.Length];
                for (int i = 0; i < bitmap.Length; i++)
                {
                    if (bitmap[i])
                    {
                        if (!isConstant)
                        {
                            values[i] = ref_Renamed + scale * bits2UInt(numbits, raf);
                        }
                        else
                        {
                            // rdg - added this to handle a constant valued parameter
                            values[i] = ref_Renamed;
                        }
                    }
                    else
                    {
                        values[i] = Grib1BinaryDataSection.UNDEFINED;
                    }
                }
            }
            else
            {
                // bms is null
                if (!isConstant)
                {
                    //System.out.println( "BDS values.size = " +
                    //(((length - 11) * 8 - unusedbits) /  numbits));
                    values = new float[((length - 11) * 8 - unusedbits) / numbits];

                    for (int i = 0; i < values.Length; i++)
                    {
                        values[i] = ref_Renamed + scale * bits2UInt(numbits, raf);
                    }
                }
                else
                {
                    // constant valued - same min and max
                    int x = 0, y = 0;
                    raf.Seek(raf.Position - 53, System.IO.SeekOrigin.Begin);                     // return to start of GDS
                    length = (int)GribNumbers.uint3(raf);
                    if (length == 42)
                    {
                        // Lambert/Mercator offset
                        SupportClass.Skip(raf, 3);
                        x = GribNumbers.int2(raf);
                        y = GribNumbers.int2(raf);
                    }
                    else
                    {
                        SupportClass.Skip(raf, 7);
                        length = (int)GribNumbers.uint3(raf);
                        if (length == 32)
                        {
                            // Polar sterographic
                            SupportClass.Skip(raf, 3);
                            x = GribNumbers.int2(raf);
                            y = GribNumbers.int2(raf);
                        }
                        else
                        {
                            x = y = 1;
                            System.Console.Out.WriteLine("BDS constant value, can't determine array size");
                        }
                    }
                    values = new float[x * y];
                    for (int i = 0; i < values.Length; i++)
                    {
                        values[i] = ref_Renamed;
                    }
                }
            }
        }         // end Grib1BinaryDataSection
Пример #8
0
        /// <summary> Constructs a <tt>Grib1GridDefinitionSection</tt> object from a raf.
        ///
        /// </summary>
        /// <param name="raf">RandomAccessFile with GDS content
        ///
        /// </param>
        /// <throws>  NoValidGribException  if raf contains no valid GRIB info </throws>
        //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'"
        public Grib1GridDefinitionSection(System.IO.FileStream raf)
        {
            int reserved;             // used to read empty space

            // octets 1-3 (Length of GDS)
            length = (int)GribNumbers.uint3(raf);
            if (length == 0)
            {
                // there's a extra byte between PDS and GDS
                SupportClass.Skip(raf, -2);
                length = (int)GribNumbers.uint3(raf);
            }
            //System.out.println( "GDS length = " + length );

            // TODO Fix Checksum stuff
            // get byte array for this gds, then reset raf to same position
            // calculate checksum for this gds via the byte array

            /*
             * long mark = raf.Position;
             * sbyte[] dst = new sbyte[length - 3];
             * SupportClass.ReadInput(raf, dst, 0, dst.Length);
             * raf.Seek(mark, System.IO.SeekOrigin.Begin);
             * //UPGRADE_ISSUE: Class 'java.util.zip.CRC32' was not converted. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1000_javautilzipCRC32'"
             * //UPGRADE_ISSUE: Constructor 'java.util.zip.CRC32.CRC32' was not converted. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1000_javautilzipCRC32'"
             * CRC32 cs = new CRC32();
             * //UPGRADE_ISSUE: Method 'java.util.zip.CRC32.update' was not converted. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1000_javautilzipCRC32'"
             * cs.update(dst);
             * //UPGRADE_ISSUE: Method 'java.util.zip.CRC32.getValue' was not converted. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1000_javautilzipCRC32'"
             * checksum = System.Convert.ToString(cs.getValue());
             * //System.out.println( "GDS checksum =" + checksum );
             */

            // octets 4 NV
            int NV = raf.ReadByte();

            //System.out.println( "GDS NV = " + NV );

            // octet 5 PL the location (octet number) of the list of numbers of points in each row
            P_VorL = raf.ReadByte();
            //System.out.println( "GDS PL = " + P_VorL );

            // octet 6 (grid type)
            type = raf.ReadByte();
            //System.out.println( "GDS grid type = " + type );
            name = getName(type);

            if (type != 50)
            {
                // values same up to resolution

                // octets 7-8 (Nx - number of points along x-axis)
                nx = GribNumbers.int2(raf);
                nx = (nx == -1)?1:nx;

                // octets 9-10 (Ny - number of points along y-axis)
                ny = GribNumbers.int2(raf);
                ny = (ny == -1)?1:ny;

                // octets 11-13 (La1 - latitude of first grid point)
                lat1 = GribNumbers.int3(raf) / 1000.0;

                // octets 14-16 (Lo1 - longitude of first grid point)
                lon1 = GribNumbers.int3(raf) / 1000.0;

                // octet 17 (resolution and component flags).  See Table 7
                resolution = raf.ReadByte();
            }

            switch (type)
            {
            //  Latitude/Longitude  grids ,  Arakawa semi-staggered e-grid rotated
            //  Arakawa filled e-grid rotated
            case 0:
            case 4:
            case 40:
            case 201:
            case 202:

                // octets 18-20 (La2 - latitude of last grid point)
                lat2 = GribNumbers.int3(raf) / 1000.0;

                // octets 21-23 (Lo2 - longitude of last grid point)
                lon2 = GribNumbers.int3(raf) / 1000.0;

                // octets 24-25 (Dx - Longitudinal Direction Increment )
                dx = GribNumbers.int2(raf) / 1000.0;

                // octets 26-27 (Dy - Latitudinal Direction Increment )
                //               Np - parallels between a pole and the equator
                if (type == 4)
                {
                    np = GribNumbers.int2(raf);
                }
                else
                {
                    dy = GribNumbers.int2(raf) / 1000.0;
                }

                // octet 28 (Scanning mode)  See Table 8
                scan = raf.ReadByte();

                // octet 29-32 reserved
                reserved = GribNumbers.int4(raf);

                if (length > 32)
                {
                    // getP_VorL(raf);
                    // Vertical coordinates (NV) and thinned grids (PL) not supported - skip this
                    SupportClass.Skip(raf, length - 32);
                }

                break;                         // end Latitude/Longitude grids


            case 1:                      //  Mercator grids

                // octets 18-20 (La2 - latitude of last grid point)
                lat2 = GribNumbers.int3(raf) / 1000.0;

                // octets 21-23 (Lo2 - longitude of last grid point)
                lon2 = GribNumbers.int3(raf) / 1000.0;

                // octets 24-26 (Latin - latitude where cylinder intersects the earth
                latin1 = GribNumbers.int3(raf) / 1000.0;

                // octet 27 reserved
                reserved = raf.ReadByte();

                // octet 28 (Scanning mode)  See Table 8
                scan = raf.ReadByte();

                // octets 29-31 (Dx - Longitudinal Direction Increment )
                dx = GribNumbers.int3(raf);

                // octets 32-34 (Dx - Longitudinal Direction Increment )
                dy = GribNumbers.int3(raf);

                // octet 35-42 reserved
                reserved = GribNumbers.int4(raf);
                reserved = GribNumbers.int4(raf);

                if (length > 42)
                {
                    // getP_VorL(raf);
                    // Vertical coordinates (NV) and thinned grids (PL) not supported - skip this
                    SupportClass.Skip(raf, length - 42);
                }

                break;                         // end Mercator grids


            case 3:                      // Lambert Conformal

                // octets 18-20 (Lov - Orientation of the grid - east lon parallel to y axis)
                lov = GribNumbers.int3(raf) / 1000.0;

                // octets 21-23 (Dx - the X-direction grid length) See Note 2 of Table D
                dx = GribNumbers.int3(raf);

                // octets 24-26 (Dy - the Y-direction grid length) See Note 2 of Table D
                dy = GribNumbers.int3(raf);

                // octets 27 (Projection Center flag) See Note 5 of Table D
                proj_center = raf.ReadByte();

                // octet 28 (Scanning mode)  See Table 8
                scan = raf.ReadByte();

                // octets 29-31 (Latin1 - first lat where secant cone cuts spherical earth
                latin1 = GribNumbers.int3(raf) / 1000.0;

                // octets 32-34 (Latin2 - second lat where secant cone cuts spherical earth)
                latin2 = GribNumbers.int3(raf) / 1000.0;

                // octets 35-37 (lat of southern pole)
                latsp = GribNumbers.int3(raf) / 1000.0;

                // octets 38-40 (lon of southern pole)
                lonsp = GribNumbers.int3(raf) / 1000.0;

                // octets 41-42
                reserved = GribNumbers.int2(raf);

                if (length > 42)
                {
                    // getP_VorL(raf);
                    // Vertical coordinates (NV) and thinned grids (PL) not supported - skip this
                    SupportClass.Skip(raf, length - 42);
                }

                break;                         // end Lambert Conformal


            case 5:                      //  Polar Stereographic grids

                // octets 18-20 (Lov - Orientation of the grid - east lon parallel to y axis)
                lov = GribNumbers.int3(raf) / 1000.0;

                // octets 21-23 (Dx - Longitudinal Direction Increment )
                dx = GribNumbers.int3(raf);

                // octets 24-26(Dy - Latitudinal Direction Increment )
                dy = GribNumbers.int3(raf);

                // octets 27 (Projection Center flag) See Note 5 of Table D
                proj_center = raf.ReadByte();

                // octet 28 (Scanning mode)  See Table 8
                scan = raf.ReadByte();

                // octet 29-32 reserved
                reserved = GribNumbers.int4(raf);

                if (length > 32)
                {
                    // getP_VorL(raf);
                    // Vertical coordinates (NV) and thinned grids (PL) not supported - skip this
                    SupportClass.Skip(raf, length - 32);
                }

                break;                         // end Polar Stereographic grids

            //rotated grid
            case 10:
                //18-20
                this.lat2 = GribNumbers.int3(raf) / 1000.0;
                //21-23
                this.lon2 = GribNumbers.int3(raf) / 1000.0;

                //24-25
                var tmpDx = GribNumbers.int2(raf);
                this.dx = tmpDx == -1 ? -9999.0 : tmpDx / 1000.0;

                //26-27
                var tmpDy = GribNumbers.int2(raf);
                this.dy = tmpDy == -1 ? -9999.0 : tmpDy / 1000.0;

                //28
                scan = raf.ReadByte();
                //29-32
                GribNumbers.int4(raf);

                //33-35
                this.latsp = GribNumbers.int3(raf) / 1000.0;
                //36-38
                this.lonsp = GribNumbers.int3(raf) / 1000.0;

                //39-42 - angle
                GribNumbers.float4(raf);

                if (length > 42)
                {
                    SupportClass.Skip(raf, length - 42);
                }

                break;


            default:
                System.Console.Out.WriteLine("Unknown Grid Type : " + type);
                break;
            }             // end switch grid_type

            //System.out.println( "scan=" + scan );
            if ((scan & 63) != 0)
            {
                throw new NoValidGribException("GDS: This scanning mode (" + scan + ") is not supported.");
            }
        }         // end Grib1GridDefinitionSection( raf )