Dear rSigGeo,

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") <- "LL"
attr(LatLong, "zone") <-11
UTM<-convUL(LatLong)


#and that's what I get
> UTM
         X         Y
1 14088.69 -16092.76
2 14082.61 -16096.32

Now, the UTM values are supposed to be around:
    X    Y
1 414040     6301764.2
2 414493           6301888.39

So, I then tried with rgdal
library(rgdal)

UTM<-project(LatLong, "+proj=utm +zone=11 ellps=WGS84")

#and that's what i get:

projected point not finite
projected point not finite

# The possible errors I see are DD format in the LatLong file. Maybe I
should be able to specify that it is zone 11 north or I think NAD 27 should
also be specified somewhere. What do you think?

	[[alternative HTML version deleted]]

