Bernese单点定位数据准备及处理

Posted 2020-08-14 王小会

原创作者 blog ：http://yifeiyao.blog.163.com/blog/static/2058932752012669731170/
1.准备所需用的数据文件，如下： 原始观测.O文件，from ftp://cddis.gsfc.nasa.gov/pub/gps/data/daily/wwww
精密星历文件，下载后为.SP3格式，更改格式为IGSwwwwd.PRE fromftp://igscb.jpl.nasa.gov/igscb/product/wwww
地球自转参数信息文件，下载后为.ERP格式，更改格式为IGSwwww7.IEP fromftp://igscb.jpl.nasa.gov/igscb/product/wwww
卫星精密钟差IGSwwwwd.CLK from ftp://igscb.jpl.nasa.gov/igscb/product/wwww
参考框架相关的参考基准IGS_08_R.CRD from ftp://ftp.unibe.ch/aiub/BSWUSER/STA
速度场IGS_08_R.VEL from ftp://ftp.unibe.ch/aiub/BSWUSER/STA 
电离层CODwwwwd.ION from ftp://ftp.unibe.ch/aiub/CODE/yyyy
(for QIF ambiguity resolution)
码偏差P1P2yymm.DCB from ftp://ftp.unibe.ch/aiub/CODE/yyyy
(for ionosphere determination)
码偏差P1C1yymm.DCB from ftp://ftp.unibe.ch/aiub/CODE/yyyy
(for clock estimation)
注：精密星历、精密钟差、地球自转参数信息也均可在ftp://cddis.gsfc.nasa.gov/pub/gps/products/wwww
以上所有的文件名均需大写，包括后缀名。
海潮可不考虑。
其中：wwww ... GPS week
d ... day of GPS week
yyyy ... four-digit year
yy ... two-digit year
mm ... month of year
2.编辑一个项目名：‘Campaign > Edit list of campaigns‘
3.选择一个可用的项目：‘Campaign > Select active campaign‘
4.新建一个工程，‘Campaign > Create new campaign‘会在项目目录下(GPSDATA/项目名)自动生成所需的文件夹。
5.更改所对应的各种日期：‘Configure > set session/compute date，包括GPS周、年积日。
6.编辑对话列表，‘Campaign > Edit session table‘，会在GPSDATA/项目名/STA下自动生成对话列表文件。
7.将第一步下载的文件按照下列路径复制到对应文件夹中，
ORB/IGSwwwwd.PRE
ORB/IGSwwww7.IEP
ORB/P1P2yymm.DCB
ORB/P1C1yymm.DCB
OUT/IGSwwwwd.CLK
STA/IGS_00_R.CRD
STA/IGS_00_R.VEL
ATM/CODwwwwd.ION
海潮改正影响不大，可以不考虑。
8. 将观测文件复制到ORX下，Copy your RINEX observation files into the ‘ORX‘
9. 更新输入文件（非常重要），‘Configure >Update input files‘，更新所有的文件。
10.‘RINEX > RINEX utilities > Extract station information‘从观测文件中提取测站相关信息在STA文件夹中生成.STA文件，包括：测站名、接收机类型、天线类型、天线相位改正等等，可以在软件中添加测站信息，也可以在文件中直接添加。
11.Create a PLD-file ‘Campaign > Station files > Tectonic plate assignment‘，在STA文件夹中生成板块.PLD文件。
12.Create a FIX-file ‘Campaign > Station files > Station selection lists‘，在STA文件夹中生成基准.fix文件。
13.Create a ABB-file ‘Campaign > Station files > Abbreviation table‘，在STA文件夹中生成站点缩写.ABB文件。
注：PLD、FIX、ABB文件也可以直接参考EXAMPLE格式直接生成或者拷贝，这样更简单。
14.‘BPE > Edit PCF program input files‘，更改PPP.PCF(解算先验坐标)。更改Phase center variations的PHAS_IGS.REL为PHAS_COD.I08（PHAS_COD.I05或者PHAS_COD.I01）；Satellite information的SATELLITE.为SATELLITE.I08(SATELLITE.I05或者SATELLITE.I01)。其它参数自行设置。
15.‘BPE >Start BPE process‘进行BPE批处理。注意：变量PLD、ABB、FIX、BLQ与编辑的文件名PLD、ABB、FIX、BLQ相同，否则运行出错。
16.如果运行出错，查找路径为GPSDATA/项目/BPE/.PRT文件下，查找最后几行，寻找错误。一般错误有：接收机类型，天线相位类型（更新这些对应文件即可）、无法找到某个文件（软件对应变量与自行编译的文件名不一样，更改文件名即可）、找不到某个变量（变量没进行更改或者变量版本不对）；错误分析详见文献《bernesegps4.2版本数据处理软件的使用与常见错误分析》。