Caspian sea disappears on world map plot

[Kristoffer Magnusson]

Hello,

I'm using this shape file from Natural Earth data to plot a world map with ggplot2, but for some reason the caspian sea is missing on the plot. It's there if I use plot(), but not in my ggplot.

Here's my code:

library(maptools)

library(ggplot2)

world.sp <- readShapePoly("50m-land/50m_land", proj4string=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"))

world.sp <- fortify(world.sp)

ggplot(world.sp, aes(long, lat, group=group)) + geom_polygon()

Anyone know how to solve this problem? 

Best,

Kristoffer

[Winston Chang]

This looks like a problem with holes in polygons - the great lakes in North America don't show up either. Could you file an issue about this on github?

https://github.com/hadley/ggplot2/issues

Thanks!

-Winston

--
You received this message because you are subscribed to the ggplot2 mailing list.
Please provide a reproducible example: https://github.com/hadley/devtools/wiki/Reproducibility
 
To post: email ggp...@googlegroups.com
To unsubscribe: email ggplot2+u...@googlegroups.com
More options: http://groups.google.com/group/ggplot2

[Mark Connolly]

You might take a look at https://github.com/hadley/ggplot2/wiki/plotting-polygon-shapefiles and see if it helps.  Look down at the Problems heading.

[Kristoffer Magnusson]

Hi,

I figured out that my problem was caused by how ggplot handles holes in polygons, as Winston mentioned. I wasn't aware that ggplot plots holes as a separate polygon above the land polygon (in this case) and not as a true hole, I can get the caspian sea to show up if I set "fill = hole". However, It's not an ideal solutions as I want all water area to be transparent. 

Thanks!

Kristoffer