Simple Basemap Example

You can do a lot with Python basemaps. Here is an example of a few things you can do...

	# Brian Blaylock
	# January 5, 2017 It's snowing outside
	"""
	Simple Basemap examples
	Example gallery: http://matplotlib.org/basemap/users/examples.html
	Documentation: https://basemaptutorial.readthedocs.io/en/latest/
	"""
	from mpl_toolkits.basemap import Basemap
	import matplotlib.pyplot as plt
	from matplotlib.collections import PatchCollection
	from matplotlib.patches import Polygon
	import numpy as np
	# Make the figure
	fig = plt.figure()
	ax = fig.add_subplot(111)
	# Easiest way to make a basemap is to use the cylidrical projection and
	# define the bottom left lat/lon and top right lat/lon corners
	# Map of Utah
	bot_left_lat =36.5
	bot_left_lon =-114.5
	top_right_lat =42.5
	top_right_lon = -108.5
	# create the map object, m
	m = Basemap(resolution='i', projection='cyl', \
	llcrnrlon=bot_left_lon, llcrnrlat=bot_left_lat, \
	urcrnrlon=top_right_lon, urcrnrlat=top_right_lat)
	# Note: You can define the resolution of the map you just created. Higher
	# resolutions take longer to create.
	# 'c' - crude
	# 'l' - low
	# 'i' - intermediate
	# 'h' - high
	# 'f' - full
	# Draw some map elements on the map
	m.drawcoastlines()
	m.drawstates()
	m.drawcountries()
	m.drawrivers(color='blue')
	# Drawing ArcGIS Basemap (only works with cylc projections??)
	# Examples of what each map looks like can be found here:
	# http://kbkb-wx-python.blogspot.com/2016/04/python-basemap-background-image-from.html
	maps = ['ESRI_Imagery_World_2D', # 0
	'ESRI_StreetMap_World_2D', # 1
	'NatGeo_World_Map', # 2
	'NGS_Topo_US_2D', # 3
	'Ocean_Basemap', # 4
	'USA_Topo_Maps', # 5
	'World_Imagery', # 6
	'World_Physical_Map', # 7
	'World_Shaded_Relief', # 8
	'World_Street_Map', # 9
	'World_Terrain_Base', # 10
	'World_Topo_Map' # 11
	]
	print "drawing image from arcGIS server...",
	m.arcgisimage(service=maps[8], xpixels=1000, verbose=False)
	print "...finished"
	# Plot a scatter point at WBB on the map object
	lon = -111.85
	lat = 40.77
	m.scatter(lon,lat,c='r',s=150)
	# Plot some wind barbs
	lons = np.arange(-115,-100,.5)
	lats = np.arange(33,48,.5)
	u = np.arange(-5,10,.5)
	v = np.arange(5,20,.5)
	m.barbs(lons, lats, u, v, color='fuchsia')
	# Plot line between two points
	# (can also use greatcircle function to be more accurate)
	x = [-110, -112]
	y = [40, 42]
	m.plot(x, y, color='navy', lw=5)
	# Fill two polygon shapes
	patches = []
	homeplate = np.array([[-114,38],[-113,37],[-112,38],[-112,40],[-114,40]])
	patches.append(Polygon(homeplate))
	triangle = np.array([[-111,38],[-110,37],[-110,42]])
	patches.append(Polygon(triangle))
	ax.add_collection(PatchCollection(patches, facecolor='lightgreen', edgecolor='k', linewidths=1.5))
	# Plot shapefiles: see here: http://basemaptutorial.readthedocs.io/en/latest/shapefile.html
	# Plot contours
	# m.contour(lons2D, lats2D, values2D) # contour lines
	# m.contourf(lons2D, lats2D, values2D) # contour color filled, can specify a cmap
	# Plot gridded data
	# m.pcolormesh(lons2D, lats2D, values2D) # can specify a cmap
	# Add plot title and other plot elements the normal way
	plt.title('Map of Utah Basemap Example')
	plt.xlabel('this is the x label')
	plt.show()

view raw simple_basemap_example.py hosted with ❤ by GitHub