basemap icon indicating copy to clipboard operation
basemap copied to clipboard

Using mercator projection: I want to color only lands/countries, and all land(except Antarctica) plus all ocean gets colored

Open Code-Sage opened this issue 5 years ago • 0 comments

This code is from Basemap docs pages, with some changing to focus on my issue at hand. I want to color(in this case, magenta) all lands(Antarctica can be overlooked), and all ocean should be white like Antarctica is right now. The shapefile is from here

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from matplotlib.patches import Polygon
from matplotlib.collections import PatchCollection
from matplotlib.patches import PathPatch
import numpy as np

fig = plt.figure(figsize=(20,20))
ax = fig.add_subplot(111)

map = Basemap(llcrnrlon=-180,llcrnrlat=-85,urcrnrlon=180.,urcrnrlat=85, projection='merc')

map.readshapefile('countries_lakes/ne_10m_admin_0_countries_lakes', 'units')

patches = []

for info, shape in zip(map.units_info, map.units):
    patches.append(Polygon(np.array(shape), True))
        
ax.add_collection(PatchCollection(patches, facecolor= 'm', edgecolor='k', linewidths=1.))

plt.savefig('dem.png', bbox_inches='tight')`

dem

Code-Sage avatar Feb 09 '20 19:02 Code-Sage