[原创]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处理星历文件格式