[R-br] Rotação de matrizes com coordenadas geográficas em UTM
Paulo Justiniano
paulojus em leg.ufpr.br
Terça Abril 2 07:58:52 BRT 2013
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 em yahoo.com.br
> alexandre.santos em cas.ifmt.edu.br
> ======================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
Mais detalhes sobre a lista de discussão R-br