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:[email protected]
[email protected]
======================================================================
_______________________________________________
R-br mailing list
[email protected]
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
[email protected]
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:[email protected]
[email protected]
======================================================================
_______________________________________________
R-br mailing list
[email protected]
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.
