On 03/01/2014 01:53 AM, Ciara O'Hara wrote:
Hi everyone,
I am using the script below to create a map plot of Ireland.
library(ncdf)library(maps)library(mapdata)library(maptools)library(RColorBrewer)shapefile<-readShapeLines("C:\\Users\\sophysics\\Desktop\\IRL_adm\\IRL_adm1.shp")file<-open.ncdf("C:\\Users\\sophysics\\Desktop\\TempData\\average\\average_1961.nc")data<-get.var.ncdf(file,"var61")lat<-get.var.ncdf(file,"lat")lon<-get.var.ncdf(file,"lon")map("worldHires","Ireland",col="purple",fill=FALSE,box=NULL)image(lon,lat,data,col="purple",xaxt="n",yaxt="n",ann=FALSE,box=NULL,xlab=NULL,ylab=NULL,asp=-1)plot(shapefile,add=TRUE,border=FALSE)
I want to get that box around to fit snugly around the coastline, but I can't figure out how to do this. The
box("inner") command doesn't seem to be doing anything for me, and box("outer"),
box("figure") create boxes outside. I've tried various combinations of the par(mar) and
par(omi/oma) commands, but these refer to the boxes outside the inner plot.
Is it possible to make the map fit snugly inside a box?
Hi Ciara,
As I don't have the map data that you used, I can only suggest that you
try fiddling with the xlim and ylim arguments to get as little space
around Ireland as you can. Add axes to your original plot to see roughly
where these limits would be.
Jim
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
