No seu ultimo comando voce adicionou pontos ao grafico anterior com points() entretanto os pontos estao fora da escala do gráfico e nao aparecem nba umtima linha simplesmente troque points() por plot() On Tue, 2 Apr 2013, ASANTOS wrote:
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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.