[R-br] Malha sobre poligono

[Henrique Dallazuanna]

Você pode fazer algo como:

l <- replicate(scan(), locator(4))

Onde o scan() recebe seu n, depois é so tratar para armazenar da forma
que vc precisa

2011/8/16 Eder David Borges da Silva <eder em leg.ufpr.br>:
> #Caros colegas,
> #minha duvida é uma implementação para gerar uma malha da seguinte forma:
> # 1) Imaginem que tenho um poligono como esse?
> x <- c(-0.512320711, -0.503592775, -0.669423554, -0.841072957, -0.459953097,
> -0.186477777, -0.221389520, -0.000281815,
> 0.319742495,  0.569943319,  0.503029145,  0.209188642,  0.022992680,
> -0.169021906, -0.241754703, -0.305759565, -0.352308556, -0.390129611,
> -0.422132042, -0.422132042,-0.512320711)
> y <- c( 0.73217213,  0.38635022,  0.08640264, -0.43585901, -0.80638249,
> -0.44997419, -0.03357638,
> 0.12169059,-0.13238264,  0.07228747,  0.48162768,  0.81333441 , 1.10622439,
>  1.13092596,  1.14504114,  1.13445475,  1.12739716,  1.12033957,
>  1.06035006,  1.00388934, 0.73217213)
> plot(x,y,type='l')
> # 2) dentro deste poligono quero fazer quadrados irregulares como este
> px1 <- c(-0.6131326, -0.4633484, -0.3762646, -0.5034070)
> py1 <- c(-0.6037305,-0.7381249, -0.5774359, -0.3992173)
> polygon(px1,py1)
> # 3) quero fazer isso com o mouser clicando 4 vezes e guardando as
> coordenadas  e cada coordenada ganha um ID, ficando assim:
> data.frame(id=1:length(px1),px1,py1)
> # 4) este quatros pontos forman um elemento (quadrado irregular), onde o
> elemento guarda o ID os pontos clicados.
> elemento  <- c(1,1,2,3,4)#Idelemento - 4 ids dos pontos
> # 5) Esse procedimento é repetido n  vezes cobrindo toda a região com
> poligonos irregulares, como resultado final:
> # um data frame com id e x, y do ponto
> # um data frame com  id do elemento
> # restrições para o segundo elemento eu devo clicar por exemplo no ponto 2 e
> 3 e fazer mais o 5 e 6 para fechar ele, desta forma,
> # os elementos tem de ficar plotados para aproveitar os pontos já feitos.
> #OBS: para os deterministicos de plantão isso é uma malha para modelos
> numericos de mecanica de fluidos.
> #OBS1: Sugestoes de softwares livres que faz isso são bem vindo.
> # Diante do problema alguem tem sugestoes de funçoes que posso utilizar, na
> verdade é um mix de locator + identifique, pois não
> # consigui fechar todo o raciocinio de como fazer isso.
> # Obrigado pessoal
> #Att
> _______________________________________________
> 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.
>



-- 
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40" S 49° 16' 22" O