Bernese数据处理报告

Bernese数据处理报告 

实验目的 

1. 学习Bernese软件的安装和使用。 

2. 对2012年年积日为094的七个中国IGS基准站（shao、chan、bjfs、kunm、guao、lhaz、twtf）数据进行分时段处理，得到四个时段的PPP定位结果；固定某三个测站坐标，求两个时段的求差定位结果。 

实验步骤 

安装软件 

①解压软件压缩包。 

②添加系统环境变量。 

③安装Activeperl软件。 

数据准备 

从ftp://cddis.gsfc.nasa.gov/pub/gps/data/daily 下载上述7个站2012-4-3（对应年积日为94）整天的原始观测值文件（Rinex 格式O文件），原始观测值文件采样间隔为30s。进行非差定位解算时，从ftp://ftp.unibe.ch/aiub/CODE/2012/下载对应的精密轨道、钟差文件和对应的地球自转参数文件（igs16822.sp3，igs16822.clk，igs16822.erp），电离层文件（COD16822.ION）。从ftp://unibe.ch 下载*.DCB文件（P1P21204.DCB，P1C11204.DCB）。在求差定位解算和结果分析中用到的参考坐标为IGS发布的上述跟踪站周解坐标（igs12P1682.snx）。 

非差数据批处理 

I 建立项目。在“Menu> Campaign> Edit list of campaigns”中定义新项目的名字“Xu”，点击“Save”保存项目。 

    

图3.1 建立工程和子目录 

II 创建项目相关的子目录。“Menu> Campaign > Create new campaign”点击“Run”。 

III 激活工程。在“Menu> Campaign > Select active campaign”选择活动工程“Xu”。 

 

图3.2 选择活动工程 

IV设置时段的时间信息。在“Menu > Configure > Set session/compute_date”下设置“年 年积日”为“2012 94”，点击“Compute”、“Set”， 然后保存。 

 图3.3 设置时段的时间信息 

V 处理数据准备。将RINEX格式的观测值文件（O文件）放入项目的RAW和ORX目录下； 轨道相关文件（包括精密轨道.sp3和钟差文件.clk）以及与它相对应的地球自转参数文件（后缀名由.erp 改成 .iep）及DCB文件放入项目ORB目录下。 

VI 设置处理时段。“Menu> Campaign > Edit session table”下分为均等四个时段。 

 

图3.4 设置处理时段 

VII 创建测站信息文件。“Menu> RINEX > RINEX utilities > Extract station information”，

选择Original RINEX observation files，并选择所有o文件作为测站信息生成的原始观测文件；

将STA命名为“12_094”。 

 

图3.5创建测站信息文件 

完成后，打开STA目录，编辑项目对应的.sta文件，将“TYPE 002: STATION INFORMATION”下的测站名复制到“TYPE 001: RENAMING OF STATIONS”下，删除“OLD STATION NAME REMARK”之后的字符。 

 

图3.6 编辑STA文件 

VIII 修改测站信息文件。“Menu> Campaign >Edit station files> Station information”，在“Old station name”一列中填入“*+原始站名”，点击“Save”保存。

  

图3.7 修改STA文件 

此时会报错，找不到天线类型： 

 

图3.8 找不到天线类型 

此时，点击“OK”、“Next”，修改天线类型文件为“PHAS_COD.I08”，再保存即可通过。 

 

图3.9 修改天线相位文件 

IX BPE数据准备。将精密钟差文件“*.clk”拷贝至 OUT目录下；电离层文件“*.ION”拷入ATM文件夹下；从“D:\bern50\BERN50\ GPS\DOC\” 目录下将 “*.CRD, *.BLQ, *.PLD *.VEL”等文件拷贝至STA目录下，重命名为“12_094.*”。 

X 编辑批处理控制文件。“Menu> BPE >Edit process control file”，选择PPP.PCF，  

 

图3.10 选择处理控制文件 

点击next至出现 

 

图3.11 修改处理参数 

将V_STAINF、V_PLDINF、V_ABBINF、V_CRDREF、V_CRDMRG、V_PCV参数分别改为12_094、12_094、12_094、IGS08_R、IGS08_R、I08，如图3.11所示（这里是因为*.pld、*.crd和天线类型文件都是用的08版本），点击“Save”保存。 

XI 进行数据的批处理。在“Menu> BPE >Start BPE process”下，选择处理时段数为4， 

 

图3.12 设置处理时段数 

点击“Next”，修改Process control file为PPP。 

 

图3.13 修改处理控制文件 

    点击“Run”运行，即可得到处理结果，处理完成界面如下： 

 

图3.14 非差数据批处理完成 

求差数据批处理 

将IGS跟踪站BJFS SHAO LHAZ作为坐标已知点（参考站），解算从站（rover）CHAN GUAO KUNM TWTF站的坐标。其中参考站坐标来自IGS公布的上述周解算坐标（1682周）文件igs12P1682.snx。解算步骤如下。 

I设置处理时段。“Menu> Campaign > Edit session table”下分为均等两个时段。 

 

图4.1 设置处理时段 

II 设置解算固定点。将 D:\bern50\BERN50\GPS\DOC 下的EXAMPLE.FIX 拷贝到SAT目录下，并改名为COORFIX.FIX，该文件为固定点配置文件，修改其内容，Station Name即要设置为参考站的站点名，设置如下。 

 

图4.2 固定点配置 

III进行数据的批处理。在“Menu> BPE >Start BPE process”下，选择处理时段数为2， 

 

 

 

图4.3 设置处理时段数

点击“Next”，修改Process control file为RNX2SNX。 

 图4.4 修改处理控制文件   

点击“Run”运行，即可得到处理结果，处理完成得到输出结果文件如下：

   

图4.5 求差数据批处理结果文件 

在实验中遇到的问题有： 

钟差文件格式错误。在进行非差数据处理时，钟差文件应使用igs格式，cod格式无法进行正常处理。 

下载的IGS测站周结果文件应该为COD*.CRD，而不是COE*.CRD，两者是不一样的周结果坐标文件。 

在求差数据处理过程中，进行BPE解算时，使用默认配置会出现以下错误：

  

图3.1 星历文件格式错误 

此时，应在“Menu> BPE >Start BPE process”下翻至第四页，将“Orbit/ERP，DCB，ION information”这一行第二列改为IGS，即igs格式的星历文件即可运行。 

 

图3.2 更改BPE处理星历文件格式