Hi R helpers,
I tried to convert a list of LatLong coordinates (in DD format) into UTM
zone 11 NAD 27. I first tried this from PBSmapping:
library(PBSmapping)
LatLong<-cbind(c(56.85359, 56.85478),c(-118.4109, -118.4035))
colnames(LatLong)<-c("X","Y")
attr(LatLong, "projection") <- "UTM"
attr(LatLong, "zone") <- 11
UTM<-convUL(LatLong)
#and that's what I get
> UTM
X Y
1 -120.9799 -1.068699
2 -120.9799 -1.068632
Now, the UTM values are supposed to be around:
X Y
1 414040 6301764.2
2 414493 6301888.39
# I don't know what is wrong, So, I then tried with rgdal
library(rgdal)
UTM<-project(LatLong, "+proj=utm +zone=11+ellps=NAD27")
# and that's what i get:
projected point not finite
projected point not finite
# Errors might come from the fact that I do have DD format in the LatLong
file. Maybe I should be able to specify that it is zone 11 north. What do
you think?
Colin
[[alternative HTML version deleted]]
______________________________________________
[email protected] 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.
