[R] map axis on projected shapefiles

penguins
Wed, 08 Aug 2012 09:10:31 -0000

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.