Message-ID: <D0924B67.112BF6%macqueen1@llnl.gov>
Date: 2014-11-19T21:28:52Z
From: MacQueen, Don
Subject: LAT/LON to UTM in R
In-Reply-To: <939705898.2064766.1416327605067.JavaMail.yahoo@jws106125.mail.bf1.yahoo.com>
By the way, this question is more appropriate for R-sig-geo.
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 11/18/14, 8:20 AM, "Zilefac Elvis" <zilefacelvis at 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 at 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.