2021-03-03BERNESE软件定轨批处理详细教程

BERNESE软件定轨批处理详细教程
非专业人士，工作需要又处理一天的grace卫星轨道数据，感觉自己之前写的教程自己都有些看不懂。遂在此细写一份grace卫星数据下载、linux下数据格式转换、Bernese安装、定轨批处理所需文件、文件位置与格式要求、结果输出等说明。
①grace卫星数据下载
建议下载使用软件FILAZILLA进行下载（网页下载容易中断，造成数据不完整）
下载地址：
ftp://isdcftp.gfz-potsdam.de/grace/Level-1B/JPL/INSTRUMENT/RL02
下载得到的文件格式为二进制的.dat压缩文件，解压之后，GVN与GPS开头文件为BERNESE处理时所需文件。
GPS开头的是我们要找的GPS星载接收机观测文件；
GVN开头的是JPL实验室公布的grace卫星轨道；
GRACE数据说明：其数据的处理和归档由喷气动力实验室（JPL），德克萨斯州大学空间研究中心（UTCSR）和波茨坦地学中心（GFZ）完成。以月解的球谐系数的方式提供。其科学数据产品分为四类：Level-0、Level-1A、Level-IB、Level-2。重点介绍Level-2数据。
Level-2数据包含几种不同类型的重力场模型。其中我们可以得到重力场模型的类型有:GSM,GCM,GAA,GAB,GAC,GAD等。
GSM和GCM分别表示只由GRACE数据建立的静态重力场模型和由GRACE数据联合其他数据建立的静态重力场模型；
GAA,GAB,GAC,GAD分别表示GSM中去除的非潮汐大气、非潮汐海洋、非潮汐大气和海洋、海洋底部洋枪的重力场模型。
②数据格式转换
进行数据处理时，需要先将二进制格式文件转换为RINEX格式，具体可参照链接：
https://blog.csdn.net/mr_zhouming/article/details/82700604
注意以下几点：
I、安装linux虚拟机时，如果无法在本机和虚拟机之间拖动文件，可查看如何设置建立共享文件夹，通过mv命令进行文件的互传。
II、注意修改GRACEsystempath后的路径，改为要转换的文件所在路径，不然无法转换。
III、文件转换命令有两种，不会就看readme！！！介绍的很详细。
③bernese安装
常规安装步骤，程序readme！！！关键是安装完成后环境变量配置。
I、DEP405.EPH文件拷贝到安装目录；
II、crx2rnx.exe和gzip.exe拷贝在安装目录下一个新文件下（如Bernese_me）;
III、将Bernese_me添加到参数变量的路径中，如下图：
此电脑属性-高级设置-环境变量-（找到path：******这一行）-编辑

④定轨批处理所需文件
星载GPS观测数据和姿态数据、精密星历、精密钟差、极移和DCB改正。
注:姿态数据非必须，如果没有，需在处理过程中设置相关配置；
另注：COD、IGS下载的数据均可，只是不同下载地址下载的文件需要放在对应的文件夹下，且处理时应配置统一。
⑤bernese批处理
3. BERNESE软件进行GRACE低轨卫星定轨
① 数据准备：示例采用2016.03.23数据，除转换好的O文件外，需准备以下数据：http://www.aiub.unibe.ch/download/CODE/

② 文件放置位置
将O观测文件存放在你所建立项目中的RAW文件夹（注意文件命名与后缀），姿态数据ATT放在ORB文件夹中（如果没有，后面设置为空）中；
将EPH.Z、CLK.Z和ERP.Z这三个文件存储在GPSDATA\DATAPOOL\COD文件夹中（如果是IGS数据中心下载，放入IGS文件夹中）；
将DCB.Z下载到GPSDATA\DATAPOOL\BSW52文件夹中；
因为BERNESE软件本身文件中就带有GRACE卫星各种信息，所以，不需要更改如EXP_LEO.CRD、EXP_LEO.STA和EXP_LEO.ABB之类的文件信息。
③ 数据处理

I、Campaign-Edit new campaign，创建新的项目（两次提示直接点确定）；
II、Campaign-Create new campaign；
III、Campaign-Select new campaign，选择要处理的项目；
IV、Configure-Set session Date，时段设置；
V、BPE->Start BPE processing进入批处理界面。
④ 批处理配置修改
因ATT文件非必要文件且难以下载获取（未找到可用ATT文件），但资料表明对最后的处理精度无影响，需要进行BPE设置，将ATT输入项去掉。
点击BPE，选择Edit PCF program input files，进行PCF项修改。主要将MAUPRP/GPSEST/CODSPP中ATT输入项去掉，如下所示：

选择LEOPOD.PCF，进入配置选择界面，需要修改的有：140/200/500/511/600，以140MAUPRP.INP为例

选择140MAUPRP.INP，点击OK，进入设置界面

点击NEXT，找到ATT文件输入位置，将ATT文件输入位置清空，点击SAVE即可，其他类同。
⑤ 进行BPE批处理
设置名称等相关信息，点击run，根据报错修改（项目文件夹下BPE文件夹，每一个步骤的log文件），直到处理成功。

4. 处理结果
   输出文件整理：
   ${S}/{V_RESULT}/yyyy/OUT/{V_LEO}{V_F}yyssss.PRC
   ${S}/{V_RESULT}/yyyy/OUT/RD{V_F}yyssss.RES
   ${S}/{V_RESULT}/yyyy/OUT/KN{V_F}yyssss.RES
   ${S}/{V_RESULT}/yyyy/OUT/KN{V_F}yyssss.CLK.gz
   ${S}/{V_RESULT}/yyyy/STA/KN{V_F}yyssss.KIN
   ${S}/{V_RESULT}/yyyy/ORB/RD{V_F}yyssss.PRE
   ${S}/{V_RESULT}/yyyy/ORB/RD{V_F}yyssss.ELE
   ${S}/{V_RESULT}/yyyy/ORB/RD{V_F}yyssss.STD
   ${S}/{V_RESULT}/yyyy/ORB/KN{V_F}yyssss.PRE
   具体解算结果看软件说明书Page641开始所示。