On Thu, 9 Jun 2005, Ernesto Jardim wrote:
> Hi,
>
> I'm trying to use transform my mercator locations into utm but I'm
doing
> something wrong because only x is transformed ... see
Two points: I think you are getting the transform you specify, from
+proj=merc to +proj=utm with a different x origin - I suspect you meant
+proj=latlong for the input data. If the coordinates below are somewhere
off Cabo de Sines, they are latlong.
Second, spproj is not on CRAN, so either the R-spatial-dev mailing list on
sourceforge (where spproj and friends are living while they grow up) or
the R-sig-geo list might have been a good choice. Please consider
continuing this on R-sig-geo.
Best wishes,
Roger
>
> > xy.sp
> SpatialPoints:
> loni lati
> [1,] -8.85 38.16
> [2,] -9.19 37.99
> [3,] -9.11 37.97
> [4,] -9.06 38.15
> [5,] -9.03 37.87
> [6,] -9.14 37.81
> [7,] -9.09 37.70
> [8,] -8.95 37.45
> [9,] -9.17 37.37
> [10,] -9.00 37.34
> Coordinate Reference System (CRS) arguments: +proj=merc +datum=WGS84
> > transform(xy.sp, CRS("+proj=utm +zone=29"))
> SpatialPoints:
> loni lati
> [1,] 1505638 38.62347
> [2,] 1505638 38.45141
> [3,] 1505638 38.43116
> [4,] 1505638 38.61335
> [5,] 1505638 38.32995
> [6,] 1505638 38.26922
> [7,] 1505638 38.15788
> [8,] 1505638 37.90485
> [9,] 1505638 37.82388
> [10,] 1505638 37.79351
> Coordinate Reference System (CRS) arguments: +proj=utm +zone=29
> +ellps=WGS84
>
> What am I doing wrong ?
>
> Thanks
>
> EJ
>
> ______________________________________________
> R-help at 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
>
--
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand at nhh.no