public static NexradScan ParseNexrad(MemoryStream stream) { var nexrad = new NexradScan(); try { var headerData = new byte[30]; byte[] descriptorData; var bytes = new byte[stream.Length]; descriptorData = new byte[stream.Length - 30]; var descriptor = new Descriptor(); var symbology = new Symbology(); // Header stream.Position = 0; stream.Read(headerData, 0, 30); var headerString = System.Text.Encoding.UTF8.GetString(headerData); headerString = headerString.Replace("\n", string.Empty).Replace("\r", " ").Replace(" ", " "); var headerStrings = headerString.Split(' '); var header = new Header { Code = headerStrings[0], Site = headerStrings[1], Date = headerStrings[2], FullProductCode = headerStrings[3] }; header.BodyCode = ReadHalfWord(stream); header.BodyDate = ReadHalfWord(stream); header.BodyTime = ReadWord(stream); header.Length = ReadWord(stream); header.SourceId = ReadHalfWord(stream); header.DestinationId = ReadHalfWord(stream); header.BlockCount = ReadHalfWord(stream); // Descriptor descriptor.Divider = ReadHalfWord(stream); descriptor.Latitude = ReadWord(stream); descriptor.Longitude = ReadWord(stream); descriptor.Height = ReadHalfWord(stream); descriptor.Code = ReadHalfWord(stream); descriptor.Mode = ReadHalfWord(stream); descriptor.VolumeCoveragePattern = ReadHalfWord(stream); descriptor.SequenceNumber = ReadHalfWord(stream); descriptor.ScanNumber = ReadHalfWord(stream); descriptor.ScanDate = ReadHalfWord(stream); descriptor.ScanTime = ReadWord(stream); descriptor.GenerationDate = ReadHalfWord(stream); descriptor.GenerationTime = ReadWord(stream); descriptor.ProductSpecific1 = ReadHalfWord(stream); descriptor.ProductSpecific2 = ReadHalfWord(stream); descriptor.ElevationNumber = ReadHalfWord(stream); descriptor.ProductSpecific3 = ReadHalfWord(stream); descriptor.Threshold1 = ReadHalfWord(stream); descriptor.Threshold2 = ReadHalfWord(stream); descriptor.Threshold3 = ReadHalfWord(stream); descriptor.Threshold4 = ReadHalfWord(stream); descriptor.Threshold5 = ReadHalfWord(stream); descriptor.Threshold6 = ReadHalfWord(stream); descriptor.Threshold7 = ReadHalfWord(stream); descriptor.Threshold8 = ReadHalfWord(stream); descriptor.Threshold9 = ReadHalfWord(stream); descriptor.Threshold10 = ReadHalfWord(stream); descriptor.Threshold11 = ReadHalfWord(stream); descriptor.Threshold12 = ReadHalfWord(stream); descriptor.Threshold13 = ReadHalfWord(stream); descriptor.Threshold14 = ReadHalfWord(stream); descriptor.Threshold15 = ReadHalfWord(stream); descriptor.Threshold16 = ReadHalfWord(stream); descriptor.ProductSpecific4 = ReadHalfWord(stream); descriptor.ProductSpecific5 = ReadHalfWord(stream); descriptor.ProductSpecific6 = ReadHalfWord(stream); descriptor.ProductSpecific7 = ReadHalfWord(stream); descriptor.ProductSpecific8 = ReadHalfWord(stream); descriptor.ProductSpecific9 = ReadHalfWord(stream); descriptor.ProductSpecific10 = ReadHalfWord(stream); descriptor.Version = ReadByte(stream); descriptor.SpotBlank = ReadByte(stream); descriptor.SymbologyOffset = ReadWord(stream); descriptor.GraphicOffset = ReadWord(stream); descriptor.TabularOffset = ReadWord(stream); if (descriptor.SymbologyOffset != 0) { symbology.Divider = ReadHalfWord(stream); symbology.BlockId = ReadHalfWord(stream); symbology.Length = ReadWord(stream); symbology.Layers = ReadHalfWord(stream); symbology.LayerDivider = ReadHalfWord(stream); symbology.DataLength = ReadWord(stream); symbology.PacketCode = ReadHalfWord(stream); symbology.IndexFirstRangeBin = ReadHalfWord(stream); symbology.RangeBinCount = ReadHalfWord(stream); // dataHeader[1] symbology.SweepCenterI = ReadHalfWord(stream); // dataHeader[2] symbology.SweepCenterJ = ReadHalfWord(stream); // dataHeader[3] symbology.ScaleFactor = ReadHalfWord(stream); // dataHeader[4] symbology.RadialCount = ReadHalfWord(stream); // dataHeader[5] symbology.RadialData = new List<Radial>(); var radialCount = symbology.RadialCount; for (var i = 0; i < radialCount; i++) { var radial = new Radial(); radial.ColorValues = new List<Int32>(); radial.RleCount = ReadHalfWord(stream); radial.StartAngle = ReadHalfWord(stream) / 10D; radial.AngleDelta = ReadHalfWord(stream) / 10D; for (var j = 0; j < radial.RleCount * 2; j++) { var values = DecodeRLES(stream, descriptor); radial.ColorValues.AddRange(values); } symbology.RadialData.Add(radial); } } nexrad = new NexradScan { Header = header, Descriptor = descriptor, Symbology = symbology }; } catch (Exception ex) { Console.WriteLine("Generic error while parsing Nexrad data", ex.ToString()); } return nexrad; }
private static List<Int32> DecodeRLES(Stream ms, Descriptor descriptor) { var data = ReadByte(ms); // read a byte, MSB is the RUN value, LSB is the COLOR CODE var value = data & 15; // get MSB var length = data >> 4; // get LSB var values = new List<Int32>(); if (descriptor.Mode == 1 && descriptor.Code >= 16 && descriptor.Code <= 21) { if (value >= 8) value = value - 8; else value = 0; } for (var i = 0; i < length; i++) { values.Add(value); } return values; }