在Python 3.2 basemap中从形状文件绘制地图时出错
我正在使用Python 3.2 basemap从形状文件绘制地图 但是,底部轴上的经度值仅部分显示。 此外,缺少所有纬度值 这是我的python代码在Python 3.2 basemap中从形状文件绘制地图时出错,python,python-3.x,shapefile,matplotlib-basemap,Python,Python 3.x,Shapefile,Matplotlib Basemap,我正在使用Python 3.2 basemap从形状文件绘制地图 但是,底部轴上的经度值仅部分显示。 此外,缺少所有纬度值 这是我的python代码 import shapefile as sf import sys import numpy as np import matplotlib.pylab as plt from mpl_toolkits.basemap import Basemap map = Basemap(projection='stere', lon_0=-106.4,
import shapefile as sf
import sys
import numpy as np
import matplotlib.pylab as plt
from mpl_toolkits.basemap import Basemap
map = Basemap(projection='stere', lon_0=-106.4, lat_0= 31.9, lat_ts = 31.9, \
llcrnrlat=31.7, urcrnrlat= 31.85, \
llcrnrlon=-106.5 , urcrnrlon= -106.1, \
rsphere=6371200., resolution='l', area_thresh=1000)
plt.figure(num=None, figsize=(10, 8), dpi=80, facecolor='w', edgecolor='k')
parallels = np.arange(31.7, 31.85, 0.25)
map.drawparallels(parallels, labels=[0, 0, 0, 1] , fontsize=10, labelstyle='+/-', dashes=[2, 2])
meridians = np.arange (-106.5, -106.1, 0.25)
map.drawmeridians(meridians, labels=[0, 0, 0, 1], fontsize=10, labelstyle='+/-' , dashes=[2, 2])
无论我如何更改标签,纬度/经度图例值仍然丢失 我想这个问题的答案已经无关紧要了,但可能对其他人有用。您需要更改代码中的行
parallels = np.arange(31.7, 31.85, 0.25)
map.drawparallels(parallels, labels=[0, 0, 0, 1],...)
在
或更短
map.drawparallels(np.arange(31.7, 31.85, 0.025), labels=[1, 1, 0, 0],...)
您的纬度值丢失,因为它们未包含在区间31.7-31.85中(31.7+0.25=31.95,31.95>31.85)
map.drawparallels(np.arange(31.7, 31.85, 0.025), labels=[1, 1, 0, 0],...)