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@yahoo.com.br alexandre.santos@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")