[R] sp: proj4string has no impact

From: Albrecht Kauffmann <alkauffm_at_rz.uni-potsdam.de>
Date: Wed 17 Jan 2007 - 16:34:07 GMT


Hi all,

I'm faced with a problem applying the sp package: The projection argument in

readShapePoly(Shapefile,proj4string="CRS class argument")

e.g.: CRS("+proj=aea +lat_1=46 +lat_2=73 +lat_0=60 +lon_0=84 +x_0=0 +y_0=0 +ellps=clrk66 +units=m +no_defs")

doesn't have any impact on the plotted object. I also tested the simple example:

xy = cbind(x = 2 * runif(100) - 1, y = 2 * runif(100) - 1) plot(SpatialPoints(xy, proj4string =
CRS("+proj=longlat")),xlim=c(-1,1),ylim=c(-1,1))

looks exactly like

plot(SpatialPoints(xy, proj4string =CRS("+proj=stere +lon_0=98 +over"))

or

plot(SpatialPoints(xy))

without any projection.

What I'm doing wrong?

I use the latest versions of sp and rgdal.

With many thanks for any hint,

Albrecht Kauffmann
Potsdam University



R-help@stat.math.ethz.ch 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. Received on Thu Jan 18 03:39:01 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 17 Jan 2007 - 19:31:35 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.