保姆级教程：手把手教你从CDDIS和武汉大学IGS数据中心下载RTKLIB所需数据（附文件命名规则详解）-平芜编程栈

高精度定位数据获取实战指南：从CDDIS到IGS的完整流程解析

刚接触RTKLIB时，最让人头疼的莫过于面对各种数据源网站和复杂的文件名。记得我第一次打开CDDIS官网时，满屏的字母数字组合文件名让我完全摸不着头脑——什么"ABPO202000101D0.01O"、"gfz20864.sp3"这些看似随机的字符串，到底哪个才是我需要的？本文将带你一步步破解这些"密码"，让你从数据获取小白变成游刃有余的实践者。

1. 数据获取前的准备工作

在开始下载数据前，我们需要先了解RTKLIB解算所需的核心数据类型。不同于简单的GPS定位，高精度定位需要多种数据协同工作：

观测值数据：接收机实际记录的卫星信号数据，通常以.o或.rnx为扩展名
导航电文：包含卫星轨道和时钟信息，常见格式有.n、.g和.rnx
精密星历：包括.sp3(轨道)和.clk(钟差)文件，用于提升定位精度
辅助改正文件：如海潮负荷(.blq)、天线相位中心(.atx)等

建议新手在电脑上创建如下目录结构：

RTKLIB_Data/ ├── Observation/ ├── Navigation/ ├── SP3/ ├── CLK/ └── Correction/

提示：所有数据文件的时间系统必须保持一致，常见的三种时间表示法需要特别注意：年积日(DOY)、GPS周和周内日、标准日历日期。

2. CDDIS数据下载详解

NASA的CDDIS(卫星数据信息系统)是最全面的GNSS数据源之一，但它的文件命名规则也是最让初学者困惑的。让我们以2023年5月1日(年积日121)的数据为例，解析典型文件的下载路径。

2.1 导航电文获取

访问CDDIS的每日数据目录：

https://cddis.nasa.gov/archive/gnss/data/daily/2023/121/

你会看到类似这样的文件：

BRDC00IGS_R_20231210000_01D_GN.rnx BRDM00DLR_R_20231210000_01D_MN.rnx

文件命名规则解析：

文件名部分	示例	含义
测站/类型	BRDC00IGS	广播星历(IGS综合)
系统标识	R	RINEX 3.x格式
日期	2023121	2023年第121天
时长	0000_01D	1天数据
内容类型	_GN	GPS导航电文

多系统导航电文通常更实用，推荐下载_MN结尾的文件

2.2 精密星历下载

精密星历位于产品目录，需要先确定GPS周数。2023年5月1日是GPS周2242的第1天，因此访问：

https://cddis.nasa.gov/archive/gnss/products/2242/

关键文件类型：

gfs22421.sp3：GFZ分析中心的精密轨道
igs22421.clk：IGS综合的精密钟差
igs22421.erp：地球定向参数

# 使用wget下载示例 wget --user=[注册邮箱] --password=[密码] https://cddis.nasa.gov/archive/gnss/products/2242/igs22421.sp3

注意：CDDIS从2020年起需要注册才能下载，建议使用与学术机构关联的邮箱注册。

3. 武汉大学IGS中心实战指南

对于国内用户，武汉大学的IGS数据中心(http://www.igs.gnsswhu.cn/)提供了更稳定的访问体验。其界面更加友好，支持按测站或日期检索。

3.1 观测值数据获取

进入"数据下载"→"观测数据"
选择日期(如2023-05-01)
输入测站名(如WUHN)
选择数据采样率(通常选30秒)
点击"检索"后下载.rnx文件

常见测站命名规律：

前4字符通常代表地点(如WUHN=武汉)
后3字符可能表示接收机编号

3.2 精密产品下载

在"产品下载"区域可以获取：

精密轨道(SP3)
精密钟差(CLK)
电离层格网(ION)
差分码偏差(DCB)

武汉中心特有的"快速产品"通常在3小时内发布，适合实时应用。

4. 辅助文件获取与处理

高精度定位还需要多种改正文件，这些通常不随日期变化，可以一次性下载备用。

4.1 天线相位中心改正

从IGS官网获取最新版天线文件：

https://files.igs.org/pub/station/general/igs14.atx

在RTKLIB中配置路径：

Options → Files → Antenna/Receiver Parameters File

4.2 海潮负荷改正

推荐使用FES2004模型：

# Python下载示例 import urllib.request url = "http://ftp.aiub.unibe.ch/BSWUSER52/STA/FES2004.BLQ" urllib.request.urlretrieve(url, "FES2004.BLQ")

在RTKPOST中配置：

Options → Files → OTL BLQ File

5. 文件转换与质量控制

下载的文件有时需要格式转换才能使用。常见转换场景包括：

CRX转RINEX：使用crx2rnx工具

crx2rnx -f CAS10000A.crx > CAS10000A.rnx

Hatanaka压缩：减小观测文件体积

rnx2crx -f WUHN00CHN_R_20231210000_01D_30S_MO.rnx

数据质量检查要点：

文件大小异常小可能是下载不完整
用文本编辑器检查文件头信息
确认时间范围覆盖所需时段
检查多系统数据的卫星系统标识

6. 自动化脚本实践

对于经常需要下载数据的用户，可以编写简单的自动化脚本。以下是Python示例：

import ftplib import datetime def download_cddis(day, save_path): date = datetime.datetime.strptime(day, "%Y-%m-%d") doy = date.strftime("%j") gps_week = (date - datetime.datetime(1980,1,6)).days // 7 with ftplib.FTP("gdc.cddis.eosdis.nasa.gov") as ftp: ftp.login("your_email", "password") # 下载导航电文 ftp.cwd(f"/gnss/data/daily/{date.year}/{doy}/") ftp.retrbinary(f"RETR BRDM00DLR_R_{date.year}{doy}0000_01D_MN.rnx", open(f"{save_path}/nav.rnx", 'wb').write) # 下载精密星历 ftp.cwd(f"/gnss/products/{gps_week}/") ftp.retrbinary(f"RETR igs{gps_week}{date.weekday()}.sp3", open(f"{save_path}/orbit.sp3", 'wb').write)