写在前面
使用basemap、cartopy等包在python中绘制地图添加比例尺有投影限制,不太方便。
本博文使用最基础的matplotlib画线实现比例尺的添加,为大家提供一种可选思路。
思路:
提供了三种可选风格,根据自己需要进行选择。
原理为根据1°=111千米的原理计算需要绘制的长度,因此本代码主要适用于低纬度地区的wgs84投影的绘图。高纬度地区投影变形太大,使用这种方法会不准确。
大家也可以根据这种思路自己修改代码。
详细代码:
#-----------函数:添加比例尺--------------
def add_scalebar(lon0,lat0,length):
# style 3
plt.hlines(y=lat0, xmin = lon0, xmax = lon0+length/111, colors="black", ls="-", lw=1, label='%d km' % (length))
plt.vlines(x = lon0, ymin = lat0-0.45, ymax = lat0+0.45, colors="black", ls="-", lw=1)
plt.vlines(x = lon0+length/2/111, ymin = lat0-0.45, ymax = lat0+0.45, colors="black", ls="-", lw=1)
plt.vlines(x = lon0+length/111, ymin = lat0-0.45, ymax = lat0+0.45, colors="black", ls="-", lw=1)
plt.text(lon0+length/111,lat0+0.6,'%d' % (length),horizontalalignment = 'center')
plt.text(lon0+length/2/111,lat0+0.6,'%d' % (length/2),horizontalalignment = 'center')
plt.text(lon0,lat0+0.6,'0',horizontalalignment = 'center')
plt.text(lon0+length/111/2*3,lat0+0.6,'km',horizontalalignment = 'center')
# # style 1
# plt.hlines(y=lat0, xmin = lon0, xmax = lon0+length/111, colors="black", ls="-", lw=2, label='%d km' % (length))
# plt.vlines(x = lon0, ymin = lat0-0.45, ymax = lat0+0.45, colors="black", ls="-", lw=2)
# plt.vlines(x = lon0+length/111, ymin = lat0-0.45, ymax = lat0+0.45, colors="black", ls="-", lw=2)
# plt.text(lon0+length/2/111,lat0+0.45,'500 km',horizontalalignment = 'center')
# style 2
# plt.hlines(y=lat0, xmin = lon0, xmax = lon0+length/111, colors="black", ls="-", lw=1, label='%d km' % (length))
# plt.vlines(x = lon0, ymin = lat0-0.45, ymax = lat0+0.45, colors="black", ls="-", lw=1)
# plt.vlines(x = lon0+length/111, ymin = lat0-0.45, ymax = lat0+0.45, colors="black", ls="-", lw=1)
# plt.text(lon0+length/111,lat0+0.7,'%d km' % (length),horizontalalignment = 'center')
# plt.text(lon0,lat0+0.7,'0',horizontalalignment = 'center')
代码调用
#-----------添加比例尺------------
add_scalebar(91,-12.5,1200) # 左下角经度,左下角纬度,比例尺长度
效果:
风格1:
风格2:
风格3: