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@yahoo.com.br
                   alexandre.santos@cas.ifmt.edu.br
            ======================================================================

            _______________________________________________
            R-br mailing list
            R-br@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@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@yahoo.com.br
        alexandre.santos@cas.ifmt.edu.br ======================================================================




_______________________________________________
R-br mailing list
R-br@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@yahoo.com.br 
        alexandre.santos@cas.ifmt.edu.br                    
======================================================================