[原创]Bernese数据处理报告

Bernese数据处理报告 

实验目的 

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

2. 2012年年积日为094的七个中国IGS基准站(shaochanbjfskunmguaolhaztwtf)数据进行分时段处理,得到四个时段的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.sp3igs16822.clkigs16822.erp),电离层文件(COD16822.ION)。ftp://unibe.ch 下载*.DCB文件(P1P21204.DCBP1C11204.DCB)。在求差定位解算和结果分析中用到的参考坐标为IGS发布的上述跟踪站周解坐标(igs12P1682.snx)。 

非差数据批处理 

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

图像   图像 

3.建立工程和子目录 

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.设置时段的时间信息 

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

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

图像 

3.设置处理时段 

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.修改STA文件 

此时会报错,找不到天线类型: 

图像 

3.找不到天线类型 

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

图像 

3.修改天线相位文件 

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_STAINFV_PLDINFV_ABBINFV_CRDREFV_CRDMRGV_PCV参数分别改为12_09412_09412_094IGS08_RIGS08_RI08,如3.11所示(这里是因为*.pld*.crd和天线类型文件都是用的08版本),点击“Save”保存。 

XI 进行数据的批处理Menu> BPE >Start BPE process下,选择处理时段数为4 

图像 

3.12 设置处理时段数 

点击“Next”,修改Process control filePPP 

图像 

3.13 修改处理控制文件 

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

图像 

3.14 非差数据批处理完成 

求差数据批处理 

IGS跟踪站BJFS SHAO LHAZ作为坐标已知点(参考站),解算从站(roverCHAN 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 fileRNX2SNX 

图像

 4.4 修改处理控制文件   

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

 图像 图像 

4.5 差数据批处理结果文件 

在实验中遇到的问题有: 

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

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

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

 图像 

3.1 星历文件格式错误 

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

图像 

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

 
posted @ 2013-01-08 18:44  lcxu2  阅读(4881)  评论(3编辑  收藏  举报