Python调用GFZRNX和RTKLIB实现GNSS数据批处理

导语: 小时解数据是全球定位系统（GPS）和其他卫星导航系统的输出结果之一。它提供了每个小时的接收机位置和时钟偏移量等关键信息。然而，当需要处理大量小时解数据时，手动操作变得费时费力。为了提高处理效率，我们可以使用Python编程语言编写代码来自动化处理这些数据。本文将介绍一个处理小时解数据的Python代码示例，帮助您快速处理大规模的数据集。

介绍: 小时解数据处理是一个常见的任务，特别是在地理测量、导航和地理空间分析等领域。它涉及对每个小时解文件进行处理，并生成定位结果。为了简化这一过程，我们将使用Python编程语言编写代码，以自动化处理小时解数据。

代码示例中的主要步骤包括：

1.设置文件和路径：根据实际情况设置小时解文件、RTKLIB可执行文件、配置文件和观测数据文件夹的路径。

2.匹配广播星历文件：根据小时解数据的起始时间，匹配对应的广播星历文件，以提供卫星轨道信息。

3.循环处理每个小时解文件：对于每个小时解文件，设置相关参数并运行RTKLIB的定位命令。

生成定位结果：根据定位命令的设置，运行RTKLIB生成定位结果文件。

这段代码示例是一个基本框架，您可以根据自己的需求进行适当的修改和扩展。例如，您可以添加数据预处理、结果分析和可视化等步骤来进一步完善处理流程。

代码示例中使用的是RTKLIB，它是一个开源的实时定位系统库，提供了一系列处理GPS和其他导航卫星系统数据的功能。通过调用RTKLIB的命令行工具，我们可以在Python中实现自动化的小时解数据处理。

首先是我们代码的文件声明信息

@Description :   Library for processing hourly solution data",
@Author      :   ZhaiWei",
@Version     :   1.0",
@Contact     :   navsense_support@163.com",
@Time        :   2023/06/23 17：21：23

 第一步，我们导入所需的模块：

import subprocess
import os
from os.path import isfile
import datetime

ymd2doy()函数用于将年月日转换为年积日。

def ymd2doy(year, month, day):
    # 将年月日转换为年积日
    dn = datetime.datetime(year, month, day, 0, 0, 0)
    return int(dn.strftime("%j"))

 process_rtklib()函数用于处理RTKLIB的相关操作。首先，它匹配广播星历文件。

def process_rtklib(rovpath, refpath, binpath_rtklib, cfgfile_rtklib, folder1):
    # 匹配广播星历文件
    stime1 = "20210101"  # 根据实际情况设置开始时间
    doy = ymd2doy(int(stime1[0:4]), int(stime1[4:6]), int(stime1[6:]))
    navfiles = [nav for nav in os.listdir(folder1) if nav.startswith("BRDC00GOP_R_2023" + str(doy).zfill(3))]
    navfile = os.path.join(folder1, navfiles[0])

设置一些变量和参数。包括设置文件路径、设置频率数量、接收机类型和RTKLIB的模式参数。其中：

1. ts、te、tint、refCoor分别被初始化为空字符串，用于存储时间戳、观测时间范围、采样间隔和参考坐标。这些变量将在后续的代码中被使用和赋值。

2. numfreq被设置为3，表示使用的频率数量。其中，1表示L1频率，2表示L1+L2频率，3表示L1+L2+L5频率。这个值将用于后续的RTKLIB处理中。

3. resType被设置为2，表示接收机类型。其中，1表示RCV_SMARTPHONE，2表示RCV_LOWCOST，3表示RCV_GEODETIC。这个值也将用于后续的RTKLIB处理中。

4. ambMode被初始化为空字符串，用于存储RTKLIB的模式参数。在这里，空字符串表示未设置模式。可以通过赋值-i来设置瞬时模式（instantaneous），或者赋值-h来设置固定并保持模式（fix and hold）。这里将ambMode设置为-h，表示使用固定并保持模式。

这些值将在后续的RTKLIB命令中使用。

    ts, te, tint, refCoor = "", "", "", ""
    os.chdir(os.path.dirname(rovpath))
    numfreq = 3  # 1:L1 ,2:L1+L2 3:L1+L2+L5
    resType = 2  # RCV_SMARTPHONE=1, RCV_LOWCOST=2, RCV_GEODETIC=3
    ambMode = ""  # -i: 瞬时模式     -h: 固定并保持模式

检查浮动站文件路径是否存在，并创建文件夹。

    if isfile(rovpath):
        fname = os.path.basename(rovpath)
        fname1 = os.path.splitext(fname)[0]
        folder = os.path.dirname(rovpath)
        if not os.path.exists(folder):
            os.makedirs(folder)

 使用循环和参数设置构建RTKLIB的命令，并运行。

        for weight in range(1, 2):
            for ref in range(1, 2):
                for ambMode in ["-h"]:
                    site = fname[:4]
                    solFile = "%s_spp_ref%d_w%d_amb%s_freq%d.pos" % (fname1, ref, weight, ambMode, numfreq)
                    rtkcmd = '%s %s %s  %s %s %s -dop 0 -x 3 -y 2 -f %d -rec %d -weg %d -ref %d -k %s -o %s %s %s %s' % \
                             (binpath_rtklib, ambMode, ts, te, tint, refCoor, numfreq, resType, weight, ref, cfgfile_rtklib,
                              solFile, rovpath, refpath, navfile)
                    print(rtkcmd)

                    # 运行命令
                    os.chdir(folder)
                    subprocess.run(rtkcmd, stdout=subprocess.DEVNULL, stderr=subprocess.DEVNULL)

main()函数是程序的入口点。在这里，你需要根据实际情况设置文件和路径，然后调用process_rtklib()函数进行RTKLIB处理。

def main():
    # 设置文件和路径
    rovpath = "浮动站每小时文件路径"
    refpath = "参考站每小时文件路径"
    binpath_rtklib = "RTKLIB可执行文件路径"
    cfgfile_rtklib = "RTKLIB配置文件路径"
    folder1 = "观测文件夹路径"

    process_rtklib(rovpath, refpath, binpath_rtklib, cfgfile_rtklib, folder1)

main()

以上代码是一个示例，它假设您已经分割并获得了小时解文件。在实际应用中，您需要根据自己的数据和文件组织方式进行相应的设置和调整。下面是对代码中涉及的几个关键参数的说明：

- `rovpath`和`refpath`：分别表示测站文件和参考站文件的路径，这里是对应的小时解文件路径。您需要根据实际情况设置这两个参数，确保它们指向正确的文件路径。

- `binpath_rtklib`：RTKLIB的可执行文件路径。RTKLIB是一个开源的实时定位系统库，包含多个命令行工具，用于处理GPS和其他导航卫星系统的数据。您需要指定正确的`binpath_rtklib`路径，确保能够调用RTKLIB的相关功能。

- `cfgfile_rtklib`：RTKLIB的配置文件路径。配置文件包含了RTKLIB定位命令的参数设置，如定位模式、导航文件等。根据您的需求，您可以根据示例中的`cfgfile_rtklib`设置来创建自己的配置文件，并指定正确的路径。

总结: 处理小时解数据的Python代码示例可以帮助您自动化处理大规模的小时解数据集。通过设置文件和路径，并结合RTKLIB的功能，您可以轻松地处理每个小时解文件，并生成定位结果。这种自动化的处理方式提高了效率，节省了时间和精力，并为后续的数据分析和应用提供了便利。您可以根据自己的需求和实际情况，对代码示例进行修改和扩展，以满足特定的处理要求。

注意：在使用代码示例之前，请确保已正确安装RTKLIB，并根据实际情况进行路径和参数的设置。