Fix several broken examples

Author: molinav

examples/geos_demo.py

@@ -20,11 +20,11 @@ def get_input(prompt):
 # plot land-sea mask.
 # land red, oceans blue.
 # lakes=True means plot inland lakes with ocean color.
+m.drawmapboundary()
 m.drawlsmask(land_color='red',ocean_color='blue',lakes=True)
 # draw parallels and meridians.
 m.drawparallels(np.arange(-90.,120.,30.))
 m.drawmeridians(np.arange(0.,420.,60.))
-m.drawmapboundary()
 plt.title('Geostationary Map Centered on Lon=%s' % (lon_0))
 
 # map with continents drawn and filled.

examples/geos_demo_2.py

@@ -39,7 +39,7 @@
 m = Basemap(projection='geos', lon_0=lon_0, satellite_height=satellite_height,
             resolution='l', llcrnrlon=ll_lon, llcrnrlat=ll_lat, urcrnrlon=ur_lon, urcrnrlat=ur_lat)
 # add data
-m.imshow(data, cmap=plt.cm.gray, interpolation='nearest')
+m.imshow(data[::-1], cmap=plt.cm.gray, interpolation='nearest')
 plt.clim(0, 255)
 # draw coastlines.
 m.drawcoastlines(linewidth=0.5, color=overlay_color)

examples/ortho_demo.py

@@ -22,13 +22,13 @@ def get_input(prompt):
 # land coral, oceans aqua.
 # lakes=True means plot inland lakes with ocean color.
 # resolution = 5 (default) means use 5 min dataset (can use 2.5)
+m.drawmapboundary()
 m.drawcoastlines()
 m.drawlsmask(land_color='coral',ocean_color='aqua', lakes=True,\
         resolution=resolution,grid=grid)
 # draw parallels and meridians.
 m.drawparallels(np.arange(-90.,120.,30.))
 m.drawmeridians(np.arange(0.,420.,60.))
-m.drawmapboundary()
 plt.title('Orthographic Map Centered on Lon=%s, Lat=%s' % (lon_0,lat_0))
 
 # map with continents drawn and filled (continent filling fails for

examples/plotcities.py

@@ -1,20 +1,22 @@
 from __future__ import (absolute_import, division, print_function)
 
-from matplotlib.mlab import prctile_rank
+import numpy as np
 import matplotlib.pyplot as plt
 from mpl_toolkits.basemap import Basemap as Basemap
+prctile_rank = lambda x, p: np.searchsorted(np.percentile(x, np.atleast_1d(p)), x)
 
 # cities colored by population rank.
 
 m = Basemap()
 shp_info = m.readshapefile('cities','cities')
-x, y = zip(*m.cities)
-pop = []
-for item in m.cities_info:
+x, y, pop = [], [], []
+for item, (x_i, y_i) in zip(m.cities_info, m.cities):
     population = item['POPULATION']
-    if population < 0: continue # population missing
-    pop.append(population)
-popranks = prctile_rank(pop,100)
+    if population >= 0:
+        pop.append(population)
+        x.append(x_i)
+        y.append(y_i)
+popranks = prctile_rank(pop,np.linspace(0, 100, 101))
 colors = []
 for rank in popranks:
     colors.append(plt.cm.jet(float(rank)/100.))

examples/plotozone.py

@@ -66,7 +66,7 @@
 
 toplot = np.ma.masked_values(o3[0, 0], 0.) * 1000.
 bounds = np.percentile(toplot.compressed().ravel(), np.linspace(5, 95, 9).tolist())
-ptch = m.pcolor(X, Y, toplot, cmap = WhGrYlBu, norm = plt.matplotlib.colors.BoundaryNorm(bounds, 20), vmin = bounds[0], vmax = bounds[-1])
+ptch = m.pcolor(X, Y, toplot, cmap = WhGrYlBu, norm = plt.matplotlib.colors.BoundaryNorm(bounds, 20))
 
 # Add a colorbar using proportional spacing, but
 # colors based on 10 distinct bins

examples/warpimage.py

@@ -81,7 +81,7 @@
 # define orthographic projection centered on Europe.
 m = Basemap(projection='ortho',lat_0=40,lon_0=40,resolution='l')
 # plot a gray-scale image specified from a URL.
-im = m.warpimage("http://earthobservatory.nasa.gov/Features/BlueMarble/Images/gebco_bathy.5400x2700.jpg")
+im = m.warpimage("https://eoimages.gsfc.nasa.gov/images/imagerecords/73000/73963/gebco_08_rev_bath_3600x1800_color.jpg")
 # draw coastlines.
 m.drawcoastlines(linewidth=0.5,color='0.5')
 # draw lat/lon grid lines every 30 degrees.