Python绘制RTKLIB导出的信噪比和多路径数值

将RTKLIB计算出来的SNR和MP值导出，运行程序，即可呈现自己绘制的SNRhe MP曲线图

import matplotlib.pyplot as plt

linedata = []
i = 1
q = 1
PRN = []
EL = []
SNR = []
MP = []
GPRN = []
GEL = []
GSNR = []
GMP = []
path = input('请输入BDS的SNR数据地址：') #
with open(path,'r+') as f1:
    linedata = f1.readlines()

# path1 = input('请输入GPS的SNR数据地址：') 
# with open(path1,'r+') as f2:
#     linedata1 = f2.readlines()

while i< len(linedata):
    cdata = linedata[i]
    cprn = cdata[25:28]
    cEL = float(cdata[42:46])
    cSNR = float(cdata[51:57])
    cMP = float(cdata[60:67])
    PRN.append(cprn)
    if cEL>=15:
        EL.append(cEL)
        MP.append(cMP)
        SNR.append(cSNR)

    i = i+1

# while q< len(linedata1):
#     cdata = linedata1[q]
#     cprn = cdata[25:28]
#     cEL = float(cdata[42:46])
#     cSNR = float(cdata[51:57])
#     cMP = float(cdata[60:67])
#     GPRN.append(cprn)
#     GEL.append(cEL)
#     GMP.append(cMP)
#     GSNR.append(cSNR)
#     q = q+1


plt.rcParams['font.sans-serif'] = ['KaiTi']
plt.rcParams['axes.unicode_minus'] = False
fig, ax = plt.subplots(figsize=(20, 10))
# plt.subplot(2,1,1)
font1 = {
'size'   : 30
}
A = ax.scatter(EL,MP, s=10, c='b', marker='o',label ='All')
# B = ax.scatter(GEL,GMP, s=10, c='r', marker='o',label ='GPS')
plt.xlabel('Elevation(°)',fontdict={'size': 30})
plt.ylabel('Multipath(m)',fontdict={'size': 30})
plt.yticks(size=25, color='k')
plt.xticks(size=25, color='k')
plt.xticks(range(5,95,5))
plt.text(83,1.3,'RMS=0.7681',fontdict={'size': 30,})
plt.grid(linestyle="-.",linewidth=0.15,color='k',axis='both')
plt.legend(handles=[A],prop=font1)
plt.show()

import matplotlib.pyplot as plt

linedata = []
i = 1
q = 1
PRN = []
EL = []
SNR = []
MP = []
GPRN = []
GEL = []
GSNR = []
GMP = []
path = input('请输入BDS的SNR数据地址：') #
with open(path,'r+') as f1:
    linedata = f1.readlines()

# path1 = input('请输入GPS的SNR数据地址：') #
# with open(path1,'r+') as f2:
#     linedata1 = f2.readlines()

while i< len(linedata):
    cdata = linedata[i]
    cprn = cdata[25:28]
    cEL = float(cdata[42:46])
    cSNR = float(cdata[51:57])
    cMP = float(cdata[60:67])
    PRN.append(cprn)
    if cEL >= 15:
        EL.append(cEL)
        MP.append(cMP)
        SNR.append(cSNR)
    i = i+1
#
# while q< len(linedata1):
#     cdata = linedata1[q]
#     cprn = cdata[25:28]
#     cEL = float(cdata[42:46])
#     cSNR = float(cdata[51:57])
#     cMP = float(cdata[60:67])
#     GPRN.append(cprn)
#     GEL.append(cEL)
#     GMP.append(cMP)
#     GSNR.append(cSNR)
#     q = q+1


plt.rcParams['font.sans-serif'] = ['KaiTi']
plt.rcParams['axes.unicode_minus'] = False
fig, ax = plt.subplots(figsize=(20, 10))
# plt.subplot(2,1,1)
font1 = {
'size'   : 25,
}
A = ax.scatter(EL,SNR, s=10, c='b', marker='o')
B = ax.scatter(GEL,GSNR, s=10, c='r', marker='o')
plt.xlabel('Elevation(°)',fontdict={'size': 30,})
plt.ylabel('SNR(dBHz)',fontdict={'size': 30,})
plt.yticks(size=25, color='k')
plt.xticks(size=25, color='k')
plt.xticks(range(5,95,5))
plt.grid(linestyle="-.",linewidth=0.15,color='k',axis='both')
# plt.legend(handles=[A],prop=font1)
# plt.subplot(2,1,2)
# plt.plot(EL,MP,'r',label = '多路径')
# plt.xticks(range(0,90,5))
# plt.grid(linestyle="-.",linewidth=0.25,color='k',axis='both')
plt.show()