matlab unistd.h,GitHub - tianruivpn1001/C-project-for-solving-tropospheric-delay-using-GGOS-troposph...

C project for solving tropospheric delay using GGOS tropospheric products

项目作者：田睿 本科毕业于北京航空航天大学，现在读研，GNSS领域入门两年。

请注意！

该项目原始路径为：E:\Thesis\RTKLIB_Workspace\WorkSpace\For_Thesis\For_Thesis

项目名称为For_Thesis(即上传的文件夹中包含的内容)

在笔者Win10系统、VS2017上成功运行，亲测可用！！！

请务必仔细阅读此文件！！！

==============程序主要功能==============

本程序的主要功能是基于GGOS对流层产品进行对流层延迟解算，还可用于rtklib开源项目的二次开发及改进。 总之，该程序的目的是方便其他使用rtklib开源项目的GNSS研究者，提高其科研效率。

==============程序基本信息==============

本程序基于rtklib2.4.3 b33版本进行二次开发，采用C语言编写，并在Win10、Visual Studio 2017平台上完成开发及调试，项目路径为： ‪E:\Thesis\RTKLIB_Workspace\WorkSpace\For_Thesis\For_Thesis 相比于原项目，主要改动点为：

1.增加了Demo_Tropo.c文件(最主要的改进)。

2.增加了测试用的main.cpp文件，作为程序入口，应注意！主程序文件必须为cpp文件，不能使用c文件。

3.在rtklib.h 文件中增加了部分函数的声明。

4.在Windows系统下，采用VS2017编译rtklib项目会遇到各种问题。对此，对原项目进行了必要的改动， 如增加了unistd.h文件、去掉main函数的rnx2rtkp.c文件等，具体编译方法详见后文。

5.增加了数据文件orography_ell.txt等，应注意！数据文件必须放在项目目录内，因为主程序中使用的是相对路径。

==============Windows下采用Visual Studio IDE的项目编译方法==============

Windows下常用的C/C++ IDE即Visual Studio，功能丰富全面，适合GNSS研究者使用。 当然，对计算机比较精通的研究者可以采用linux系统，不使用IDE，采用cmake等工具进行编译。事实上，这样做更适合rtklib的二次开发。 然而，对于刚接触GNSS领域的硕/博研究生而言，学习linux系统及cmake等编译工具，无疑会大大拉长学习周期， 耽误宝贵的科研时间，降低科研效率。对此，笔者深有体会！因此，笔者坚持在windows系统下用Visual Studio IDE开发此项目， 遇到的编译问题很多，可参考如下博客解决：

https://blog.csdn.net/wuwuku123/article/details/100030177

https://blog.csdn.net/sd28you28/article/details/82911273

https://blog.csdn.net/zhangtao_heu/article/details/79536427

https://blog.csdn.net/qq_35363018/article/details/101317869

https://blog.csdn.net/baixia3551/article/details/101085788?utm_medium=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-4.nonecase&depth_1-utm_source=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-4.nonecase https://blog.csdn.net/WalterBrien/article/details/80754215

此外，笔者在编译过程中还遇到一些没有解决的问题，在此总结发布以飨读者：

要通过Project——Manage NuGet Packages添加pthread、dirent等

一定记得切换为64位 注意预编译器添加的是WIN32而非Win32！！！

预编译器中记得添加

WIN32

_DEBUG _CONSOLE

_CRT_SECURE_NO_WARNINGS

_WINSOCK_DEPRECATED_NO_WARNINGS

ENAGLO

DLL

主程序必须是.cpp文件！！！

记得添加unistd.h文件 经上述操作后仍出现：

unresolved external symbol showmsg等错误

经测试，发现下述方法可行： 将RTKLIB-rtklib_2.4.3\app\rnx2rtkp文件夹下的rnx2rtkp.c文件放到项目的src文件夹下并加入到工程，并去掉其中的主程序(也可以把这个当做主程序)。

==============数据源==============

所用数据均可在网站https://vmf.geo.tuwien.ac.at/上下载。 应注意！VMFG_FC产品必须注册账户后才能获取，注册方法即通过英文邮件向相关管理人员提出申请，管理人员会通过邮件提供给你一个账户。 ==============开发参考==============

参考了官方的matlab程序vmf1_grid.m及vmf1_ht.m

参考了如下文献：

《GGOS对流层延迟产品精度分析及在PPP中的应用》

《不同全球对流层天顶延迟产品在中国区域的比较》

《Troposphere mapping functions for GPS and VLBI from ECMWF operational analysis data》

《Implementation and testing of the gridded Vienna Mapping Function 1 (VMF1)》

《Generation and Assessment of VMF1-Type Grids Using North-American Numerical Weather Models》

《Discussion and recommendations about the height correction for a priori zenith hydrostatic delays derived from ECMWF data》

还有一些官方网站上提供的资料(请参见官网)

==============关于程序结果的一点讨论==============

经测试，本项目运行结果与官方matlab程序有一定出入，差距在厘米级。 经检查，笔者自认为这并非是编程的问题，而是必然存在的数值计算误差，因为C程序与matlab程序的数值精度不同。 笔者测试的结果是本程序更接近IGS发布的天顶对流层延迟数据。 也可采用C与matlab联合编程应用GGOS产品，但运行较慢，且移植性不好。 如有发现程序中问题的，欢迎在GitHub库上及时发布新版本。

==============重要声明==============

鉴于笔者入门不久，水平有限，有所疏漏在所难免。 该项目仅供研究者参考，并非TU Wien发布的官方代码，可能存在笔者尚未发现的未知错误。 如各位研究者发现问题，欢迎在GitHub上发布修正版本。非常欢迎各位研究者对本项目查漏补缺，系统测试！ 在该项目尚未成熟之前，建议各位研究者审慎地考虑在科研中到底是应用本项目，还是使用C与matlab联合编程，后者虽慢，但毕竟是官方反复验证过的matlab代码。 当然，该项目注释详尽，参考价值较高，可作为学习GGOS产品应用的重要参考！ 在此开源发布，权作抛砖引玉，期冀各位研究者对其进一步改进与完善！