/// <summary> /// GNSS导航文件的读取 /// </summary> /// <param name="navFilePath"></param> /// <returns></returns> public MixedNavFile ReadGnssNavFlie() { MixedNavFile f = new MixedNavFile(); f.SetHeader(ParamNavFileReader.ReadHeader(RinexFileName)); //测试版本 if (f.Header.Version == 0) { using (StreamReader sr = new StreamReader(RinexFileName, UnicodeEncoding.UTF8)) { RinexUtil.SkipHeader(sr); } } using (StreamReader sr = new StreamReader(RinexFileName, UnicodeEncoding.UTF8)) { RinexUtil.SkipHeader(sr); if (f.Header.Version < 3.0) { while (sr.Peek() != -1) { EphemerisParam record = ParamNavFileReader.ReadRecordV2(sr, f.Header); f.Add(record); } } else if (f.Header.Version >= 3.0 && f.Header.Version < 4.0) { while (sr.Peek() != -1) { string line = sr.ReadLine(); if (string.IsNullOrWhiteSpace(line)) { continue; } SatClockBias recordHeader = ParamNavFileReader.ParseFirstLineV3(line, f.Header); //轨道参数 var satType = recordHeader.Prn.SatelliteType; if (EphemerisUtil.IsEphemerisParam(satType)) { var record = new EphemerisParam(recordHeader); ParamNavFileReader.ReadRecordBodyV3(sr, record); f.Add(record); } else { var record2 = new GlonassNavRecord(recordHeader); GlonassNaviFileReader.ParseFirstLineV3(line, record2); GlonassNaviFileReader.ReadRecordBodyV3(sr, record2); f.Add(record2); } } } } return(f); }
/// <summary> /// RINEX 3.0 记录读取。 /// </summary> /// <param name="sr">数据流</param> /// <param name="header">头部信息</param> /// <returns></returns> public static GlonassNavRecord ReadRecordV3(StreamReader sr, NavFileHeader header) { GlonassNavRecord record = new GlonassNavRecord(); string line = sr.ReadLine(); GlonassNaviFileReader.ParseFirstLineV3(line, record); return(ReadRecordBodyV3(sr, record)); }
/// <summary> /// 读取GLONASS导航文件 /// </summary> /// <param name="navFilePath">路径</param> /// <returns></returns> public static GlonassNavFile ReadFileV2V3(string navFilePath) { GlonassNavFile f = new GlonassNavFile(); f.Header = ParamNavFileReader.ReadHeader(navFilePath); //测试版本 if (f.Header.Version == 0) { using (StreamReader sr = new StreamReader(navFilePath, UnicodeEncoding.UTF8)) { RinexUtil.SkipHeader(sr); try { GlonassNaviFileReader.ReadRecordV2(sr, f.Header); GlonassNaviFileReader.ReadRecordV2(sr, f.Header); // f.Header.Version = 2.0; } catch { f.Header.Version = 1.0; } } } using (StreamReader sr = new StreamReader(navFilePath, UnicodeEncoding.UTF8)) { RinexUtil.SkipHeader(sr); while (sr.Peek() != -1) { GlonassNavRecord record = null; if (f.Header.Version < 3) { record = GlonassNaviFileReader.ReadRecordV2(sr, f.Header); } else if (f.Header.Version >= 3) { record = GlonassNaviFileReader.ReadRecordV3(sr, f.Header); } else { throw new ArgumentException(" 导航文件版本 " + f.Header.Version); } f.Add(record); } } return(f); }
/// <summary> /// 读取记录 /// </summary> /// <param name="sr">数据流读取器</param> /// <param name="header">头部信息</param> /// <returns></returns> public static GlonassNavRecord ReadRecordV2(StreamReader sr, NavFileHeader header) { GlonassNavRecord record = new GlonassNavRecord(); string line = sr.ReadLine(); GlonassNaviFileReader.ParseFirstLineV2(line, header, record); line = sr.ReadLine(); line = line.Replace('D', 'E'); double x = double.Parse(line.Substring(3, 19)); double vx = double.Parse(line.Substring(22, 19)); double vvx = double.Parse(line.Substring(41, 19)); record.Health = double.Parse(line.Substring(60, 19)); line = sr.ReadLine(); line = line.Replace('D', 'E'); double y = double.Parse(line.Substring(3, 19)); double vy = double.Parse(line.Substring(22, 19)); double vvy = double.Parse(line.Substring(41, 19)); record.FrequencyNumber = double.Parse(line.Substring(60, 19)); line = sr.ReadLine(); line = line.Replace('D', 'E'); double z = double.Parse(line.Substring(3, 19)); double vz = double.Parse(line.Substring(22, 19)); double vvz = double.Parse(line.Substring(41, 19)); record.FrequencyNumber = double.Parse(line.Substring(60, 19)); record.AgeOfOper = double.Parse(line.Substring(60, 19)); record.XYZ = new XYZ(x, y, z) * 1000.0; // km -> m record.XyzVelocity = new XYZ(vx, vy, vz) * 1000.0; // km -> m record.XyzAcceleration = new XYZ(vvx, vvy, vvz) * 1000.0; // km -> m return(record); }