[R-sig-Debian] spTransform

Thierry Onkelinx thierry.onkelinx at inbo.be
Mon May 15 13:32:11 CEST 2017

Dear Abed,

This is rather a question for r-gis-geo

Try two things:

1) simplify the CRS of the transformation
pointsTransformed <- spTransform(pointsutm, CRS("+proj=longlat"))

2) check the projection of the original data

Best regards,

ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht

To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to say
what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

2017-05-12 17:51 GMT+02:00 Abed W. Ayyad <abdulhadi.ayyad at uclmail.net>:

> Salutations,
> Long story short: I need to create an SPDF to superimpose on a Google map.
> I have .shp file which has a long list of locations listed in (what I guess
> must be) UTM. I need to transform these into lon, lat. My use of
> spTransform, however, keeps giving me the wrong results. For what it's
> worth, it only seems to be messing up the latitude.
> My final line of code looks like this:
> pointsTransformed <- spTransform(pointsutm, CRS("+proj=longlat +datum=WGS84
> +zone=36S"))
> The line before looks like:
> pointsutm <- SpatialPoints(cbind(pointsmapspatial$Point_Y,
> pointsmapspatial$Point_X), proj4string = CRS("+proj=utm +zone=36S
> +ellps=GRS80")
> Can you see what I did wrong?
>         [[alternative HTML version deleted]]
> _______________________________________________
> R-SIG-Debian mailing list
> R-SIG-Debian at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-debian

	[[alternative HTML version deleted]]

More information about the R-SIG-Debian mailing list