[R] map axis on projected shapefiles

[penguins]

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.