Hi,
I wanted to plot a coastline of Scandinavia with basemap in intermediate resolution without lakes.
The method i found was to draw the coastlines with the plot command using the coastline polygons from the projection.
(see example below, the lakes are filled with a slightly darker gray for illustration).
Please reply, if a more condensed version (e.g. proj.drawcoastlines(draw_lakes=False ?) is existing.
from pylab import *
from mpl_toolkits.basemap import Basemap
proj = Basemap(projection='lcc',
resolution='i',
llcrnrlon=9.0,
llcrnrlat=53.0,
urcrnrlon=35.0,
urcrnrlat=65.3,
lat_0=55.0,
lon_0=14.0)
figure()
for i,cp in enumerate(proj.coastpolygons):
if proj.coastpolygontypes[i]<2:
proj.plot(cp[0],cp[1],'k-')
proj.fillcontinents([0.8,0.8,0.8],lake_color=[0.75,0.75,0.75])
proj.drawmeridians(arange(12,30,4),linewidth=0.3,labels=[0,0,0,1])
proj.drawparallels(arange(54,67,2),linewidth=0.3,labels=[1,0,0,0])
show()
Richard
(i am using basemap v0.99.5)