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

[ASANTOS]

Muito obrigado Hélio,

           Muito perspicaz a solução encontrada, vou ate apagar o autocad,

Abraço,

Alexandre


-- 
======================================================================
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
======================================================================



Em 08/04/2013 18:25, Hélio Gallo Rocha escreveu:
> p0<-c(0,0)
> plot(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)
> plot(d2[,1],d2[,2])
> ang<- (pi/3) #### Angulo a ser rotacionado
> Q<-rbind(c(cos(ang),sin(ang)),c(-sin(ang), cos(ang)))  ## Matriz de 
> rotação
> Q
> drt<-d2 %*% Q  # Multiplicação de matrizes
> plot(drt[,1],drt[,2], col="red")
>  x1=  drt[,1]+524237
>  x2=  drt[,2]+ 7978856
>  d3=cbind(x1,x2)
>  d3=as.matrix(d3)
>  d3
>  plot(d3, col="red")