?数据格式:
?代码:
# -*- coding: utf-8 -*-
import pandas as pd
import matplotlib.pyplot as plt
import matplotlib.font_manager as fm
myfont = fm.FontProperties(fname='C:/Windows/Fonts/simsun.ttc')
plt.rcParams['font.sans-serif'] = ['simsun']
plt.rcParams['axes.unicode_minus'] = False
# 这个会直接默认读取到这个Excel的第一个表单
excel = pd.read_excel('GNSS对比.xlsx', sheet_name='Sheet3')
data = excel.iloc[:, :].values # 读取所有行的SXYZ以及DXYZ列的值,这里需要嵌套列表
name = data[1:, 0]
X1 = data[1:, 2]
Y1 = data[1:, 3]
X2 = data[1:, 5]
Y2 = data[1:, 6]
Delta_X = (X2 - X1) * 1000
Delta_Y = (Y2 - Y1) * 1000
fig, ax = plt.subplots()
fig.set_size_inches(9, 6)
DX = Delta_X.tolist()
DY = Delta_Y.tolist()
#变化向量图,这里需要注意DX,DY的参数类型是列表
q = ax.quiver(Y1, X1, DY, DX, width=0.0012, scale=1 / 0.001)
#添加点名
for i, txt in enumerate(name):
ax.annotate(txt, (Y1[i], X1[i]), size=8)
ax.set_title('GNSS第三次复测点位变化图')
plt.ylabel('南北方向')
plt.xlabel('东西方向')
#添加起始点位
ax.scatter(Y1, X1, c='r', s=5, label='第二次复测点位')
plt.legend()
plt.savefig(fname='1111', dpi=300)
print(Delta_X)
点位效果图:?
?
?
|