odd behaviour of llgridlines() regarding axes labels

1) is ok with xpd clipping turned off

op <- par(xpd = NA)
plot(a,axes=FALSE)
llgridlines(a, plotLines=TRUE, plotLabels=TRUE)
par(op)

2)  same
op <- par(xpd = NA)
plot(a,axes=FALSE,mar=c(4, 3, 2, 3)+0.1)
llgridlines(a, plotLines=TRUE, plotLabels=TRUE)
par(op)

3)  Consider using graticule package instead? You can control the extents
independently of the meridian/parallel lines (and optionally provided a
target projection. )

You have to give it the longs/lats you want, there's no automatic detection
from an object or plot - but for me that's the price of control over what I
get.

install.packages("graticule")
library(graticule)
plot(a,axes=TRUE)
library(graticule)
g <- graticule(seq(-90, -82, by = 2), seq(10, 16, by = 2), xlim =
par("usr")[1:2], ylim = par("usr")[3:4])
plot(g, add = TRUE, lty = 2)

Cheers, Mike.

I see an odd behaviour of llgridlines() regarding axes labels:

require(rgdal)
require(raster)
a <- extent(c(-90,-81.5625,8.407168,16.63619))
a <- as(a, 'SpatialPolygons')
projection(a) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84
+towgs84=0,0,0"

No axis labels on x axis:
plot(a,axes=FALSE)
llgridlines(a, plotLines=TRUE, plotLabels=TRUE)

#same with:
plot(a,axes=FALSE,mar=c(4, 3, 2, 3)+0.1)
llgridlines(a, plotLines=TRUE, plotLabels=TRUE)

#and 2 vertical axes with:
plot(a,axes=TRUE)
llgridlines(a, plotLines=TRUE, plotLabels=TRUE)

Any hint?

Thanks
Agus

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Dr. Michael Sumner
Software and Database Engineer
Australian Antarctic Division
203 Channel Highway
Kingston Tasmania 7050 Australia

	[[alternative HTML version deleted]]