RTKPLOT导出卫星数据并进行Python天空图绘制

RTKPLOT导出卫星数、DOP值、高度角等信息

依次点击File->Open Obs File,File->Open Nav File,分别导入观测文件及导航文件。点击右上角的设置,在Satellite System中选择卫星系统,Elevation Mask中设置卫星截止高度角。

依次点击File->Save # of Sats/DOP 可对可见卫星数/DOP值进行保存。依次点击File->Save AZ/EL/SNR/MP可对卫星方位角/卫星高度角/信噪比/L1多路径进行保存,保存数据如下,可以此进行重绘制。

【Python】根据导出的AZ/EL/SNR/MP数据进行天空图绘制

单系统代码:

import math
import matplotlib.pyplot as pltf = open('cpt-base.txt', 'r') #Add the derived file pathID=[]
AZ=[]
EL=[]ln = f.readline()
while ln:ln = f.readline()if not ln:breakstr = ln.split()id = str[2]az = float(str[3])/180*math.piel = 90 - float(str[4])ID.append(id)AZ.append(az)EL.append(el)f.close()
ax = plt.subplot(111, projection='polar')
ax.set_theta_direction(-1)
ax.set_theta_zero_location('N')
ax.set_rticks([0,15,30,45,60,75,90])satnum = 32
if  (ID[0][0] == 'G'):satnum = 32
elif(ID[0][0] == 'R'):satnum = 27
elif(ID[0][0] == 'E'):satnum = 36
elif(ID[0][0] == 'C'):satnum = 51
SATAZ=[]
SATEL=[]
SATID=[]
for i in range(1,satnum):for j in range(0,len(ID)):if int(ID[j][1:3])==i:az=AZ[j]SATAZ.append(az)el=EL[j]SATEL.append(el)id=ID[j]SATID.append(id)if len(SATAZ)>0:c = ax.scatter(SATAZ,SATEL,s=1, marker=".",alpha=0.75)ax.text(SATAZ[0],SATEL[0],SATID[0])SATAZ = []SATEL = []SATID = []ax.yaxis.set_label_position('right')
ax.tick_params('y', labelleft=False)
plt.savefig('skyplot.jpg', dpi=360)
plt.show()

单系统结果:

多系统代码:

import math
import matplotlib.pyplot as pltf = open('AZ_EL_SNR_MP.txt', 'r')ID=[]
PRN=[]
AZ=[]
EL=[]ln = f.readline()
while ln:ln = f.readline()if not ln:breakstr = ln.split()id = str[2]prn = int(str[2][1:3])if (str[2][0] == 'R'):prn = prn+ 32elif (str[2][0] == 'E'):prn = prn+ 59elif (str[2][0] == 'C'):prn = prn+ 95az = float(str[3])/180*math.piel = 90 - float(str[4])ID.append(id)PRN.append(prn)AZ.append(az)EL.append(el)
f.close()
ax = plt.subplot(111, projection='polar')
ax.set_theta_direction(-1)
ax.set_theta_zero_location('N')
ax.set_rticks([0,15,30,45,60,75,90])SATAZ=[]
SATEL=[]
SATID=[]
for i in range(1,146):for j in range(0,len(ID)):if int(PRN[j])==i:az=AZ[j]SATAZ.append(az)el=EL[j]SATEL.append(el)id=ID[j]SATID.append(id)if len(SATAZ)>0:if i<=32:c = ax.scatter(SATAZ, SATEL,s=1, color='b', marker=".", alpha=0.75)elif i<=59:c = ax.scatter(SATAZ, SATEL,s=1, color='r', marker=".", alpha=0.75)elif i<=95:c = ax.scatter(SATAZ, SATEL,s=1, color='g', marker=".", alpha=0.75)else:c = ax.scatter(SATAZ, SATEL,s=1, color='y', marker=".", alpha=0.75)middle=int(len(SATAZ)/2)ax.text(SATAZ[middle],SATEL[middle],SATID[0],fontsize=8)SATAZ = []SATEL = []SATID = []ax.yaxis.set_label_position('right')
ax.tick_params('y', labelleft=False)
plt.savefig('skyplot.jpg', dpi=360)
plt.show()

多系统结果:

RTKPLOT导出卫星数据并进行Python天空图绘制相关推荐

  1. CSDN博客第一期订阅专栏:跟“风云卫星”数据工程师学Python

    在大伙的期盼中,CSDN博客第一期订阅专栏来啦!我们邀请了许向武老师,为大家分享的主题为:跟"风云卫星数据"工程师学Python. 作者简介 许向武:擅长Python语言编程,有超 ...

  2. python 导出数据并发邮件_Python自动化导出zabbix数据并发邮件脚本

    Zabbix没有报表导出的功能,于是通过编写脚本导出zabbix数据并发邮件.效果如下: 下面是脚本,可根据自己的具体情况修改: #!/usr/bin/python` `#coding:utf-8` ...

  3. python导入excel数据-Python数据处理之导入导出excel数据

    欢迎点击上方"AntDream"关注我 .Python的一大应用就是数据分析了,而数据分析中,经常碰到需要处理Excel数据的情况.这里做一个Python处理Excel数据的总结, ...

  4. python 导出大量数据到excel_怎么在python中将大量数据导出到Excel文件

    怎么在python中将大量数据导出到Excel文件 发布时间:2021-01-16 11:09:40 来源:亿速云 阅读:76 作者:Leah 这期内容当中小编将会给大家带来有关怎么在python中将 ...

  5. python wireshark_用python编写脚本从wireshark导出的数据文件中提取数据

    上篇文章搭建了一个UDP多播程序的基础,所谓基础,就是看着它,我可以写简单的多播程序了,可以在这个基础上面开始工作了. 会多播了,多播的内容从哪里来,播出什么内容呢?呵呵,有个设备,没有通讯协议,用w ...

  6. python操作hive数据库代码_python导出hive数据表的schema实例代码

    本文研究的主要问题是python语言导出hive数据表的schema,分享了实现代码,具体如下. 为了避免运营提出无穷无尽的查询需求,我们决定将有查询价值的数据从mysql导入hive中,让他们使用H ...

  7. Python Pandas导出Hbase数据到dataframe

    Python导出Hbase数据的思路: 使用happybase连接Hbase 使用table.scan()扫数据,将得到的数据整理为dataframe格式 将从Hbase中得到的byte类型的数据转为 ...

  8. python导出百万数据到excel_使用python将大量数据导出到Excel中的小技巧分享

    使用python将大量数据导出到Excel中的小技巧分享 今天小编就为大家分享一篇使用python将大量数据导出到Excel中的小技巧心得,具有很好的参考价值,希望对大家有所帮助.一起跟随小编过来看看 ...

  9. python如何导出数据库数据库_python导出数据库数据的方法

    python导出数据库数据的方法 发布时间:2020-06-26 14:09:43 来源:亿速云 阅读:138 这篇文章运用简单易懂的例子给大家介绍python导出数据库数据的方法,代码非常详细,感兴 ...

最新文章

  1. 滴滴人脸识别申诉照片怎么拍_滴滴司机理发被停账号,平台规则到底如何遵守才能避免踩坑?...
  2. linux性能调优原创翻译系列
  3. 知识图谱最新论文清单,高阶炼丹师为你逐一解读
  4. Ubuntu16.04编译RK3399:make kernel.img error
  5. 【转】C++ Vector用法深入剖析
  6. 数组正遍历,数组倒遍历
  7. 古风手机壁纸,国潮的你不可错过!
  8. 利用VBA把PowerDesigner的comment复制到name
  9. linux查看磁盘使用情况及每个文件夹占用磁盘大小
  10. matlab矩阵中每一行数除以一个数
  11. 基于Dialog的MFC程序在启动时隐藏为托盘程序(三)
  12. 文档碎片createDocumentFragment
  13. 理正深基坑弹性计算方法_【干货】关于深基坑计算的几点思考(一个实际项目的总结)...
  14. 未来20年美国AI研究路线图
  15. 如何用数据库可视化工具:DataGrip导入数据
  16. RAID5数据应该如何恢复
  17. 电话都免费了 移动IM这么牛X 运营商知道么?
  18. Log slf4j+logback配置
  19. 个人对于po、bo、vo的理解
  20. 计算机网络中 mac地址的二进制位数,《计算机网络》阶段性测试试题(二)

热门文章

  1. 36家互联网大厂真实薪资曝光!平均月薪 25k+
  2. Bilibili 网页版下载视频 B站下载MP4格式视频(不用下载工具)
  3. 如何使用模拟量转485模块
  4. 2009年12月最新的全国铁路线路图
  5. 高中毕业会考计算机试题,高中计算机会考试题及答案2
  6. webstorm破解方法(2018.3.2 64x版本)
  7. AI测谎仪:通过微表情在法庭检测谎言,性能已优于人类
  8. 手机连无线网后没有网络连接到服务器,手机无线网络连接上但上不了网
  9. Android 截屏(Screenshot)代码流程小结
  10. 如何利用知识管理软件提高员工工作效率