Missing Islands?

[StephGagne]

Hello,

First of all, thank you for the great module! I find it incredibly useful.

I surveyed previous issues and haven't found one corresponding to my problem. Here is my issue:
Using Basemap, my colleagues and I noticed that some islands tend to be missing from the map. Given the amount of detail in the mainland shore nearby, it doesn't appear to be a resolution problem. They seem to be missing when using drawcoastlines and fillcontinents but they do appear (although without coastlines) when using drawlsmask. The first missing island we noticed is the (entire) Canadian province of Prince Edward Island. Then, we also noticed the absence of Anticosti Island, and others. They are fairly big (see figures).

Is that intended? Let me know if I can be of any help investigating this.

Thanks!
Steph

[jswhit]

Basemap uses the GSHHS coastline data (http://www.soest.hawaii.edu/pwessel/gshhg/). I would suggest lowering the area_thresh keyword. From the docstring:

area_thresh coastline or lake with an area smaller than area_thresh in km^2 will not be plotted. Default 10000,1000,100,10,1 for resolution c, l, i, h, f.

My guess is that those islands are being filtered out (as are lakes of similar size). Try area_thresh=0.

[StephGagne]

That solved it! Thank you!
Sorry for not noticing that myself.