[R-sig-Geo] Different result from sptransform() and ptransform()

Agustin Lobo alobolistas at gmail.com
Sat Jan 28 19:30:57 CET 2012


Hi!
I'm reprojecting and get different results from sptransform() than
from ptransform():
(objec cents in https://sites.google.com/site/filestemp2/home/cents.rda)

> str(cents,max.level=2)
Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
  ..@ data       :'data.frame':	234 obs. of  4 variables:
  ..@ coords.nrs : num(0)
  ..@ coords     : num [1:234, 1:2] 3603500 3604500 3605500 3606500 3607500 ...
  .. ..- attr(*, "dimnames")=List of 2
  ..@ bbox       : num [1:2, 1:2] 3603500 2155500 3620500 2167500
  .. ..- attr(*, "dimnames")=List of 2
  ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots

> proj4string(cents)
[1] " +init=epsg:3035 +proj=laea +lat_0=52 +lon_0=10 +x_0=4321000
+y_0=3210000 +ellps=GRS80 +units=m +no_defs"
> centslonlat <- spTransform(cents,CRS("+init=epsg:4326"))
> coordinates(centslonlat)[1:3,]
            X        Y
[1,] 1.315311 42.23956
[2,] 1.327326 42.24060
[3,] 1.339342 42.24163

> a <- ptransform(data = as.data.frame(coordinates(grid1km)), src.proj = CRSargs(CRS("+init=epsg:3035")), dst.proj = CRSargs(CRS("+init=epsg:4326")))
> coordinates(a)[1:3,]
              x         y z
[1,] 0.02295650 0.7372194 0
[2,] 0.02316621 0.7372375 0
[3,] 0.02337593 0.7372556 0

The correct results are those from sptransform(), but I would like to
know what I'm doing wrong with ptansform()

Thanks

Agus



More information about the R-sig-Geo mailing list