# Python代码实现：RINEX导航电文文件读取与解析 ## 代码概述 我们将编写一个Python程序，用于读取并解析名为 `brdc1290.24n` 的RINEX格式导航电文文件。代码将完成以下任务： 1. 读取文件内容 2. 解析头部信息（如电离层参数、UTC校正等） 3. 提取每颗卫星的广播星历数据 4. 将结果结构化存储为字典或DataFrame 5. 支持后续保存为CSV或其他格式 --- ## 代码实现 ```python import re from datetime import datetime, timedelta import pandas as pd def read_rinex_nav_file(file_path): """ 读取并解析 RINEX 导航电文文件 (.24n) 参数: file_path (str): 文件路径 返回: dict: 包含头部信息和卫星星历数据的字典 """ with open(file_path, 'r', encoding='utf-8') as f: lines = f.readlines() # 初始化变量 header = {} ephemeris_data = [] current_satellite = None epoch_pattern = re.compile( r'^(\d{2}) (\d{2}) (\d{2}) (\d{2}) (\d{2}) (\d{2}\.\d+)' r'\s+(-?\d+\.\d+D[+-]\d+)\s+(-?\d+\.\d+D[+-]\d+)\s+(-?\d+\.\d+D[+-]\d+)', re.IGNORECASE ) in_header = True line_iter = iter(lines) for line in line_iter: stripped_line = line.rstrip() # 检查是否仍在头部区域 if "END OF HEADER" in stripped_line: in_header = False continue if in_header: # 提取关键头部字段 if "ION ALPHA" in stripped_line: header['ion_alpha'] = list(map(lambda x: float(x.replace('D', 'E')), stripped_line[:56].split())) elif "ION BETA" in stripped_line: header['ion_beta'] = list(map(lambda x: float(x.replace('D', 'E')), stripped_line[:56].split())) elif "DELTA-UTC" in stripped_line: parts = stripped_line[:56].split() header['delta_utc_A0'] = float(parts[0].replace('D', 'E')) header['delta_utc_A1'] = float(parts[1].replace('D', 'E')) header['delta_utc_T'] = int(parts[2]) header['delta_utc_W'] = int(parts[3]) elif "LEAP SECONDS" in stripped_line: header['leap_seconds'] = int(stripped_line[:6]) else: # 开始解析星历数据 line = line.rstrip() if len(line) < 8: continue # 识别新卫星（前两个字符为整数） try: prn = int(line[0:2]) # 卫星PRN号 year = int(line[2:4]) month = int(line[4:6]) day = int(line[6:8]) hour = int(line[8:10]) minute = int(line[10:12]) second_float = float(line[12:20]) second = int(second_float) microsecond = int((second_float - second) * 1e6) # 构造完整时间对象（2000年后） base_year = 2000 if year < 80 else 1900 epoch_time = datetime(base_year + year, month, day, hour, minute, second, microsecond) # 读取第一行其余部分（a0, a1, a2） data_part = line[22:56] + ' ' + line[56:60] # 钟差相关参数 values = [data_part[i:i+19] for i in range(0, len(data_part), 19)] values = [float(v.strip().replace('D', 'E')) for v in values if v.strip()] current_satellite = { 'prn': prn, 'epoch': epoch_time, 'clock_bias': values[0], # 卫星钟偏差 $a_0$ 'clock_drift': values[1], # 卫星钟速 $a_1$ 'clock_drift_rate': values[2], # 卫星钟漂移 $a_2$ } # 读取接下来的三行（共四行一组） line1 = next(line_iter).rstrip() line2 = next(line_iter).rstrip() line3 = next(line_iter).rstrip() # 解析第二行：IODE, C_rs, Δn, M₀ parse_line_into_sat(line1, current_satellite, [ 'iode_m0', 'crs', 'delta_n', 'm0' ]) # 第三行：C_uc, e, C_us, sqrt_a parse_line_into_sat(line2, current_satellite, [ 'cuc', 'eccentricity', 'cus', 'sqrt_a' ]) # 第四行：toe, C_ic, Ω₀, C_is parse_line_into_sat(line3, current_satellite, [ 'toe', 'cic', 'omega0', 'cis' ]) # 继续读第五行（可选）：C_rc, ω, Ω̇, IDOT try: line4 = next(line_iter).rstrip() parse_line_into_sat(line4, current_satellite, [ 'crc', 'omega', 'omegadot', 'idot' ]) except StopIteration: pass ephemeris_data.append(current_satellite) except ValueError: # 忽略无法解析的行（如空行或注释） continue return { 'header': header, 'ephemeris': ephemeris_data } def parse_line_into_sat(line, sat_dict, keys): """ 将一行文本解析为多个浮点数并存入字典 """ values = [line[i*19:(i+1)*19].strip() for i in range(4)] values = [float(v.replace('D', 'E')) for v in values] for k, v in zip(keys, values): sat_dict[k] = v # === 主程序执行 === if __name__ == "__main__": # 设置文件路径 file_path = "brdc1290.24n" # 请确保此文件位于当前目录下 print("正在读取和解析 RINEX 导航文件...") result = read_rinex_nav_file(file_path) print(f"✅ 成功解析头部信息") print(f"📊 共提取 {len(result['ephemeris'])} 条卫星星历记录") # 转换为 DataFrame df = pd.DataFrame(result['ephemeris']) print("\n📈 前5条数据预览：") print(df.head()) # 存储结果到 CSV output_csv = "parsed_rinex_nav.csv" df.to_csv(output_csv, index=False) print(f"\n💾 数据已保存至 '{output_csv}'") # 可选：查看头部信息 print("\n📄 头部信息摘要：") for k, v in result['header'].items(): print(f" {k}: {v}") ``` --- ## 使用说明 1. **准备环境**： ```bash pip install pandas ``` 2. **放置文件**：将 `brdc1290.24n` 文件与脚本放在同一目录下。 3. **运行脚本**：执行后会输出解析结果，并生成 `parsed_rinex_nav.csv` 文件。 4. **输出字段说明**： - `prn`: 卫星编号 - `epoch`: 观测时刻 - 各类轨道参数（如半长轴平方根、偏心率、近地点角距等） --- ## 知识点 ### 1. RINEX 导航电文格式 标准文本格式，每颗卫星占四行，包含钟差、开普勒参数及摄动修正项。 ### 2. 科学记数法转换（D→E） Fortran风格的指数表示法（如 `0.123D+04`）需转为 Python 可识别的 `E` 形式。 ### 3. 字符串固定宽度切片解析 按每字段19列宽度截取，适用于无分隔符的规则文本数据解析。