[R-br] Rotação de matrizes com coordenadas geográficas em UTM

[ASANTOS]

Bom dia pessoal,

           Estou cá de volta com o problema de rotacionar matrizes para 
rotacionar um grid criado com expand.grid(), o exemplo numérico que 
havia postado na lista esta OK, porém agora quando emprego coordenadas 
geográficas em UTM não estou conseguindo plotar, pois os valores no eixo 
x estão estranhos e o eixo y esta negativo, alguém saberia me dar alguma 
sugestão, segue CRM abaixo, obrigado

#
p0<-c(524237,7978856)### Primeiro ponto (Coordenada geográfica conhecida 
em UTM)
points(p0[1],p0[2], col="red")
#
#
c0 <- seq(p0[1],p0[1]+400,by=100)### Criando o grid
d0 <- seq(p0[2],p0[2]+500,by=100)
d1 <- expand.grid(x=c0, y=d0)
plot(d1[,1],d1[,2])
#
d2<-as.matrix(d1)
ang<- 60 #### Angulo a ser rotacionado
Q<-rbind(c(cos(ang),sin(ang)),c(-sin(ang), cos(ang))) ## Matriz de rotação
drt<-d2 %*% Q  # Multiplicação de matrizes
points(drt[,1],drt[,2], col="red")
#

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
Coordenador do curso Técnico em Florestas
Vice Coordenador do curso de Engenharia Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
======================================================================