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

[Paulo Justiniano]

Alexandre

Então o correto não é que voce nao consegue plotar, e sium não consegue 
sobrepor a um outro mapa.

de fato a rotação vai criar novas coordenadas que nao tem a mesma projeção 
para ser superposta no mapa original

No caso entao sugiro voce criar um grid retangular e depois deslocar as 
coordanadas X ou Y de acordo com algum angulo desejado
por uma simples operacao trigonometrica.
a grade pode depois ser "recortada" para retar apenas os pontos na área.

note que gráficos como image() persp() contour() nao poderao ser feitos 
com esta grade.
Outra possibilidade é rotacionar a area de base

Ao final disto talvez velha a pena pensar se voce de fato precisa fazer 
tal rotacao...



On Tue, 2 Apr 2013, ASANTOS wrote:

> Obrigado Prof. Paulo,
> 
>          Mas estar fora de escala é justamente meu problema, quando tento plotar estes pontos sobre um shapefile, quando uso o grid
> regular (plot(d1[,1],d1[,2]) ) os pontos aparecem na minha imagem, porém quando uso o grid rotacionado (points(drt[,1],drt[,2],
> col="red")), não aparece, pois não tenho coordenadas UTM negativas e um Mod() não resolve o problema,
> 
> Obrigado,
> 
> 
> Alexandre
> 
> 
> Em 02/04/2013 06:58, Paulo Justiniano escreveu:
>       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.
> 
> 
> 
> _______________________________________________
> 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.
> 
> 
> -- 
> ======================================================================
> 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 
> ======================================================================
> 
>