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

[ASANTOS]

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
======================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130402/5dc8c1ca/attachment.html>