<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">Alexandre,</span><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div>
<div style="font-family:arial,sans-serif;font-size:13px">Andei olhando o seu CRM sobre rotação.</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">
Tava achando que era a sua matriz de multiplicação, mas não é..</div><div style="font-family:arial,sans-serif;font-size:13px">fiz tres mudanças, veja se é isso</div><div style="font-family:arial,sans-serif;font-size:13px">
<br></div><div style="font-family:arial,sans-serif;font-size:13px">1. Ao invés de usar as coordenadas UTM, que não davam certo, coloquei inicial (0,0)</div><div style="font-family:arial,sans-serif;font-size:13px">2. ao invés de usar ang=60, usei ang=pi/3</div>
<div style="font-family:arial,sans-serif;font-size:13px">3. Depois somei o resultado com as coordenadas iniciais em UTM,<br></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">
<br></div><div style="font-family:arial,sans-serif;font-size:13px">é uma "gambiarra", pode ser , mas parece que deu certo</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">
<div>p0<-c(0,0)</div><div>plot(p0[1],p0[2], col="red")</div><div class="im"><div><br></div><div>c0 <- seq(p0[1],p0[1]+400,by=100)### Criando o grid</div><div>d0 <- seq(p0[2],p0[2]+500,by=100)</div><div>
d1 <- expand.grid(x=c0, y=d0)</div><div><br></div><div>plot(d1[,1],d1[,2])</div></div><div>d2<-as.matrix(d1)</div><div>plot(d2[,1],d2[,2])</div><div>ang<- (pi/3) #### Angulo a ser rotacionado</div><div class="im">
Q<-rbind(c(cos(ang),sin(ang)),c(-sin(ang), cos(ang))) ## Matriz de rotação</div><div>Q</div><div class="im">drt<-d2 %*% Q # Multiplicação de matrizes</div><div>plot(drt[,1],drt[,2], col="red")</div><div>
x1= drt[,1]+524237</div><div> x2= drt[,2]+ 7978856</div><div> d3=cbind(x1,x2)</div><div> d3=as.matrix(d3)</div><div> d3</div><div> plot(d3, col="red")</div><div><br></div><div><br></div></div><div class="gmail_extra" style="font-family:arial,sans-serif;font-size:13px">
Hélio<br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 3 de abril de 2013 20:35, ASANTOS [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4658838h46@n4.nabble.com" target="_blank">ml-node+s2285057n4658838h46@n4.nabble.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Prof. Paulo,<br>
<br>
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,<br>
<br>
Redobrados agradecimentos,<br>
<br>
Alexandre<br>
<br>
<br>
<div>Em 03/04/2013 16:15, Paulo Justiniano
escreveu:<br>
</div>
<blockquote style="border-left:2px solid #cccccc;padding:0 1em" type="cite"><div><div class="h5">Alexandre
<br>
<br>
Então o correto não é que voce nao consegue plotar, e sium não
consegue sobrepor a um outro mapa.
<br>
<br>
de fato a rotação vai criar novas coordenadas que nao tem a mesma
projeção para ser superposta no mapa original
<br>
<br>
No caso entao sugiro voce criar um grid retangular e depois
deslocar as coordanadas X ou Y de acordo com algum angulo desejado
<br>
por uma simples operacao trigonometrica.
<br>
a grade pode depois ser "recortada" para retar apenas os pontos na
área.
<br>
<br>
note que gráficos como image() persp() contour() nao poderao ser
feitos com esta grade.
<br>
Outra possibilidade é rotacionar a area de base
<br>
<br>
Ao final disto talvez velha a pena pensar se voce de fato precisa
fazer tal rotacao...
<br>
<br>
<br>
<br>
On Tue, 2 Apr 2013, ASANTOS wrote:
<br>
<br>
</div></div><blockquote style="border-left:2px solid #cccccc;padding:0 1em" type="cite"><div><div class="h5">Obrigado Prof. Paulo,
<br>
<br>
Mas estar fora de escala é justamente meu problema,
quando tento plotar estes pontos sobre um shapefile, quando uso
o grid
<br>
regular (plot(d1[,1],d1[,2]) ) os pontos aparecem na minha
imagem, porém quando uso o grid rotacionado
(points(drt[,1],drt[,2],
<br>
col="red")), não aparece, pois não tenho coordenadas UTM
negativas e um Mod() não resolve o problema,
<br>
<br>
Obrigado,
<br>
<br>
<br>
Alexandre
<br>
<br>
<br>
Em 02/04/2013 06:58, Paulo Justiniano escreveu:
<br>
No seu ultimo comando voce adicionou pontos ao grafico
anterior com points() entretanto os pontos estao fora da escala
do
<br>
gráfico e nao aparecem
<br>
<br>
nba umtima linha simplesmente troque points() por plot()
<br>
<br>
On Tue, 2 Apr 2013, ASANTOS wrote:
<br>
<br>
Bom dia pessoal,
<br>
<br>
Estou cá de volta com o problema de
rotacionar matrizes para rotacionar um grid criado com
<br>
expand.grid(), o exemplo numérico que havia postado
na lista esta OK, porém agora quando emprego coordenadas
<br>
geográficas em UTM não estou conseguindo plotar,
pois os valores no eixo x estão estranhos e o eixo y esta
<br>
negativo, alguém saberia me dar alguma sugestão,
segue CRM abaixo, obrigado
<br>
<br>
#
<br>
p0<-c(524237,7978856)### Primeiro ponto
(Coordenada geográfica conhecida em UTM)
<br>
points(p0[1],p0[2], col="red")
<br>
#
<br>
#
<br>
c0 <- seq(p0[1],p0[1]+400,by=100)### Criando o
grid
<br>
d0 <- seq(p0[2],p0[2]+500,by=100)
<br>
d1 <- expand.grid(x=c0, y=d0)
<br>
plot(d1[,1],d1[,2])
<br>
#
<br>
d2<-as.matrix(d1)
<br>
ang<- 60 #### Angulo a ser rotacionado
<br>
Q<-rbind(c(cos(ang),sin(ang)),c(-sin(ang),
cos(ang))) ## Matriz de rotação
<br>
drt<-d2 %*% Q # Multiplicação de matrizes
<br>
points(drt[,1],drt[,2], col="red")
<br>
#
<br>
<br>
--
<br>
======================================================================
<br>
Alexandre dos Santos
<br>
Proteção Florestal
<br>
Coordenador do curso Técnico em Florestas
<br>
Vice Coordenador do curso de Engenharia Florestal
<br>
IFMT - Instituto Federal de Educação, Ciência e
Tecnologia de Mato Grosso
<br>
Campus Cáceres
<br>
Caixa Postal 244
<br>
Avenida dos Ramires, s/n
<br>
Bairro: Distrito Industrial
<br>
Cáceres - MT CEP: 78.200-000
<br>
Fone: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a>
(VIVO)
<br></div></div><div><div class="h5">
<a href="http://user/SendEmail.jtp?type=node&node=4658838&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br>
<a href="http://user/SendEmail.jtp?type=node&node=4658838&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br>
======================================================================
<br>
<br>
_______________________________________________
<br>
R-br mailing list
<br>
<a href="http://user/SendEmail.jtp?type=node&node=4658838&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
<br>
Leia o guia de postagem
(<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo
reproduzível.
<br>
<br>
<br>
<br>
_______________________________________________
<br>
R-br mailing list
<br>
<a href="http://user/SendEmail.jtp?type=node&node=4658838&i=3" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
<br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e
forneça código mínimo reproduzível.
<br>
<br>
<br>
-- <br>
======================================================================
<br>
Alexandre dos Santos
<br>
Proteção Florestal
<br>
Coordenador do curso Técnico em Florestas
<br>
Vice Coordenador do curso de Engenharia Florestal IFMT -
Instituto Federal de Educação, Ciência e Tecnologia de Mato
Grosso
<br>
Campus Cáceres
<br>
Caixa Postal 244
<br>
Avenida dos Ramires, s/n
<br>
Bairro: Distrito Industrial Cáceres - MT
CEP: 78.200-000
<br>
Fone: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)
<br>
<a href="http://user/SendEmail.jtp?type=node&node=4658838&i=4" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br>
<a href="http://user/SendEmail.jtp?type=node&node=4658838&i=5" rel="nofollow" link="external" target="_blank">[hidden email]</a>
======================================================================
<br>
<br>
<br>
</div></div></blockquote><div><div class="h5">
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
R-br mailing list
<a href="http://user/SendEmail.jtp?type=node&node=4658838&i=6" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</div></div></blockquote>
<br>
<pre cols="72"><div class="im">--
======================================================================
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: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)
</div><div class="im"><a href="http://user/SendEmail.jtp?type=node&node=4658838&i=7" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<a href="http://user/SendEmail.jtp?type=node&node=4658838&i=8" rel="nofollow" link="external" target="_blank">[hidden email]</a>
======================================================================
</div></pre><div class="im">
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4658838&i=9" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr noshade size="1" color="#cccccc">
</div><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div class="im">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</div><a href="http://r-br.2285057.n4.nabble.com/R-br-Rotacao-de-matrizes-com-coordenadas-geograficas-em-UTM-tp4658805p4658838.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Rotacao-de-matrizes-com-coordenadas-geograficas-em-UTM-tp4658805p4658838.html</a>
</div><div class="HOEnZb"><div class="h5">
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4" target="_blank">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>