python basemap_在python中使用matplotlib basemap生成globe

用一些基础代码制作地图。

它们在我看来真的很好(我是basemap的初学者),但每次运行它们时都会遇到相同的错误:

GEOS_错误:b'IllegalArgumentException:RobustPreliminate遇到非有限数'

这是我最感兴趣的网站代码:

from mpl_toolkits.basemap import Basemap

import matplotlib.pyplot as plt

import numpy as np

fig = plt.figure(figsize=(7,6))

# set perspective angle

lat_viewing_angle = 45.4642

lon_viewing_angle = 9.1900

# define color maps for water and land

ocean_map = (plt.get_cmap('ocean'))(210)

cmap = plt.get_cmap('gist_earth')

# call the basemap and use orthographic projection at viewing angle

m1 = Basemap(projection='ortho',

lat_0=lat_viewing_angle,lon_0=lon_viewing_angle,resolution=None)

# define map coordinates from full-scale globe

map_coords_xy = [m1.llcrnrx,m1.llcrnry,m1.urcrnrx,m1.urcrnry]

map_coords_geo = [m1.llcrnrlat,m1.llcrnrlon,m1.urcrnrlat,m1.urcrnrlon]

#zoom proportion and re-plot map

zoom_prop = 7.0 # use 1.0 for full-scale map

m = Basemap(projection='ortho',resolution='l',

lat_0=lat_viewing_angle,lon_0=lon_viewing_angle,llcrnrx=-map_coords_xy[2]/zoom_prop,

llcrnry=-map_coords_xy[3]/zoom_prop,urcrnrx=map_coords_xy[2]/zoom_prop,

urcrnry=map_coords_xy[3]/zoom_prop)

# coastlines, map boundary, fill continents/water, fill ocean, draw countries

m.drawmapboundary(fill_color=ocean_map)

m.fillcontinents(color=cmap(200),lake_color=ocean_map)

m.drawcoastlines()

m.drawcountries()

# latitude/longitude line vectors

lat_line_range = [-90,90]

lat_lines = 8

lat_line_count = (lat_line_range[1]-lat_line_range[0])/lat_lines

merid_range = [-180,180]

merid_lines = 8

merid_count = (merid_range[1]-merid_range[0])/merid_lines

m.drawparallels(np.arange(lat_line_range[0],lat_line_range[1],lat_line_count))

m.drawmeridians(np.arange(merid_range[0],merid_range[1],merid_count))

# scatter to indicate lat/lon point

x,y = m(lon_viewing_angle,lat_viewing_angle)

m.scatter(x,y,marker='o',color='#DDDDDD',s=3000,zorder=10,alpha=0.7,\

edgecolor='#000000')

m.scatter(x,y,marker='o',color='#000000',s=100,zorder=10,alpha=0.7,\

edgecolor='#000000')

plt.annotate('Milan, Italy', xy=(x, y), xycoords='data',

xytext=(-110, -10), textcoords='offset points',

color='k',fontsize=12,bbox=dict(facecolor='w', alpha=0.5),

arrowprops=dict(arrowstyle="fancy", color='k'),

zorder=20)

# save figure at 150 dpi and show it

plt.savefig('ortho_zoom_example.png',dpi=150,transparent=True)

plt.show()

错误甚至连哪一行的问题都没有说出来,所以我忍不住。。。