[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