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: hodgesse@uhd.edu


	[[alternative HTML version deleted]]

