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]]

