Difference between revisions of "File:Parallel Mesh Generation Around the world.png"

import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt
fig=plt.figure(figsize=[15,10])

ax = plt.axes(projection=ccrs.Robinson())
ax.add_feature(cfeature.LAND,facecolor='lightgrey')
#ax.add_feature(cfeature.OCEAN)
ax.add_feature(cfeature.COASTLINE)
ax.add_feature(cfeature.BORDERS)

# make the map global rather than have it zoom in to
# the extents of any plotted data
#ax.set_global()

#ax.set_extent([-160,160,-90,90])


#https://stackoverflow.com/a/25421922
transform = ccrs.PlateCarree()._as_mpl_transform(ax)
for key, value in places.items():
    print(key)
    lon,lat,x,y = value
    # add marker
    plt.plot(lon,lat, marker='D',color='red',markersize=3,transform=ccrs.Geodetic())
    # add text
    plt.annotate(key,xy=(lon,lat), xytext=(x,y),arrowprops=dict(arrowstyle="->",color='red'),xycoords=transform)

# Save the plot by calling plt.savefig() BEFORE plt.show()
#plt.savefig('map.svg')
plt.savefig('map.png',dpi=200,bbox_inches='tight')

plt.show()