[R] filled.contour ON TOP of a base map
Ray Brownrigg
Ray.Brownrigg at ecs.vuw.ac.nz
Thu Apr 15 23:45:20 CEST 2010
On Thu, 15 Apr 2010, Víctor Homar Santaner wrote:
> Anyone?
>
The answer appears to be "no".
Going back to your original question though, you could just not fill the map. I.e.
filled.contour(x,y,z,levels=c(2,4,6,8,10),plot.axes=map("worldHires",add=T))
Otherwise, you may have to do things the 'hard way': set up two regions, draw a map in the
first, (here you can then use polygon to colour the ocean green, then overlay the map
again), overlay your filled contour (perhaps again using polygon), then set up the legend
in the other region.
Hope this helps,
Ray Brownrigg
On Wed, 14 Apr 2010, Víctor Homar Santaner originally wrote:
> Dear R helpers,
>
> Any suggestion on how to make the 2D gaussian to show up ON TOP of the map,
> instead of below the continents?
>
> Example code:
>
> library(mapdata)
>
> x<-seq(-15,15,.2);y<-seq(30,55,.2)
>
> expone<-function(x,y){10*exp(-sqrt(x*x+(y-40)*(y-40))/2)}
> z<-outer(x,y,expone)
>
> filled.contour(x,y,z,levels=c(2,4,6,8,10),plot.axes=map("worldHires",add=T,
>col="lightgrey",fill=T))
>
>
> Also, any idea on how to fill the ocean with a background color (e.g.
> green)?
>
>
> Best regards and thank you,
>
> Víctor.
More information about the R-help
mailing list