Hello again.
I'm having a little bit of trouble with spTransform (probably I'm doing it
wrong), but here is what I'm doing:
Original data frame with UTM locations near Phuket, Thailand:
> test1.df
Loc east north
1 a 748.168 602.861
2 b 754.302 615.609
> proj1 <- CRS("+proj=utm +zone=47N")
> coordinates(test1.df) <- c("east","north")
> proj4string(test1.df) <- proj1
Set up a spatial data frame
> test1.df
coordinates Loc
1 (748.168, 602.861) a
2 (754.302, 615.609) b
> spTransform(test1.df, CRS("+proj=longlat +datum=WGS84"))
coordinates Loc
1 (94.518, 0.00543748) a
2 (94.518, 0.00555246) b
>
Now the answers should be about 8 degrees north and 98 degrees East.
Any suggestions would be most appreciated.
Sincerely,
Erin
Erin M. Hodgess, PhD
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
mailto: hodge...@uhd.edu
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
