[R-br] Malha sobre poligono

[Eder David Borges da Silva]

#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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110816/2cb35c51/attachment.html>