Zilefac,
I tried to run your "reproducible" code, but you didn't provide an object
"xy".
Jean
On Tue, Nov 18, 2014 at 10:20 AM, Zilefac Elvis <zilefacel...@yahoo.com>
wrote:
> Hi,
> I am trying to convert lat/lon to UTM but my results are extremely flawed.
> Is it because of wrong UTM zone or my lon points are negative?
>
> I read a shapefile from ESRI using:
> library(rgdal)
> stations <- readOGR(".","stations")
> print(proj4string(stations))
> ## +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
> How to convert to UTM the following lat/lon points? Reproducible example
> at the end.
>
>
> library(sp)
>
> coordinates(xy) <- c("X", "Y")
> proj4string(xy) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84
> +no_defs") ## for example
> res <- spTransform(xy, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84
> +no_defs"))
> res
> as(res, "SpatialPoints")
> x<-as(res, "SpatialPoints")
> xx<-as.data.frame(x)
>
> #OR
>
> project(as.matrix(xy[,c("X","Y")]), "+proj=longlat +ellps=WGS84
> +datum=WGS84 +no_defs")
> Everything goes wrong with the values.
>
> Please help.
>
>
>
>
>
> Zilefac.
>
> ______________________________________________
> R-help@r-project.org 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.
>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org 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.
