[R-br] Rotação de matrizes com coordenadas geográficas em UTM
ASANTOS
alexandresantosbr em yahoo.com.br
Quarta Abril 3 20:34:30 BRT 2013
Prof. Paulo,
Muito obrigado, mas desisti de rotacionar a área da base apesar
de necessário, minhas áreas ou seja talhões são irregulares e a criação
de um expand.grid() com posterior over() para criar uma máscara de
pontos que estão contidos no interior da área não aproveitam muito bem a
área disponível. Por ora, vou me rebaixar a ter que jogar o grid no
autocad e rotacionar, mas futuramente pretendo achar uma solução para isto,
Redobrados agradecimentos,
Alexandre
Em 03/04/2013 16:15, Paulo Justiniano escreveu:
> 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
>> ======================================================================
>>
>>
>
>
> _______________________________________________
> 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/20130403/7909d1e9/attachment.html>
Mais detalhes sobre a lista de discussão R-br