Hi,
I have overlayed 2 projected shapefiles using the "plot" function.
When I plot them the numbers next to the x-axis ticks are printed in
scientific format (e.g. -4e+05, -3e+05,...,), and the problem is these dont
match up to the projected point locations I also wish to overlay from my
dataset. For example, according to the x-axis ticks the prime meridian (0°)
is in totally the wrong position.
I have tried "degAxis" which only changes the axis to 4e+05°W, 3e+05°W,...,
and I cant seem to get map.grid to work.
Does anyone know why this would happen or what I have done wrong?
#Code is
shelfline <- readOGR(dsn="C:\\ArcGIS\\phy_bathy_contours.shp", layer =
"phy_bathy_contours")
world <- readOGR(dsn="C:\\ArcGIS\\phy_coastline_plus_ssi_so.shp", layer =
"phy_coastline_plus_ssi_so")
plot(shelfline, axes=T,xlab=expression("Longitude"^o),
ylab=expression("Latitude"^o), yaxt="n", xaxt="n")
plot(world, col="lightslategray", add=T)
#attempts to alter axes
degAxis(1)
map.grid(col="black", labels=TRUE, pretty=TRUE, lty=2)
Many thanks
--
View this message in context:
http://r.789695.n4.nabble.com/map-axis-on-projected-shapefiles-tp4639554.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
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.
