
Hi R-geo users, 
I am trying to convert a data set in UTM (from Spain, zone=30) to the corresponding 
data set in longitude-latitude. Using convUL, it gives me strange results. Here I send my code, just in case somebody can tell me what i am missing. I suspect it has something to do with + proj or something that it is needed. 
Thank you so much in advance,
Berta

library(PBSmapping)
datainUTM<-data.frame(cbind(c(559994, 559673),c(4781905,4781583), c(0, 0), c(1,2)))
names(datainUTM)<-c("X", "Y", "PID", "POS")
attr(datainUTM, "projection") <-"UTM"
attr(datainUTM, "zone") <-30
datosLL<-convUL(datainUTM)

	[[alternative HTML version deleted]]


