<div>Kin, boa tarde!</div><div><br></div><div>O erro listado é por atribuir ao parâmetro pch o objeto 'membros' que não existe ou não está acessível (veja seu arquivo de dados e declarações).</div><div>Acredito que a ideia seria diferenciar os pontos (tipos) de acordo com o valor em membros...</div>
<div><br></div><div>Acredito que o exemplo que segue faça algo parecido com o que você deseja, embora usando somente a biblioteca sp(). </div><div><br></div><div>Inclui os pontos que você forneceu.</div><div><br></div><div>
###########</div><div><br></div><div><div>require(sp)</div><div>br.df <- read.table("<a href="http://dl.dropbox.com/u/117618178/brlite.dat">http://dl.dropbox.com/u/117618178/brlite.dat</a>", sep=";", head=T)</div>
<div>class(br.df); head(br.df); str(br.df)</div><div><br></div><div>br.lin <- Lines(list(Line(br.df)), ID = "Brasil")</div><div><a href="http://br.sl">br.sl</a> <- SpatialLines(list(br.lin))</div><div>proj4string(<a href="http://br.sl">br.sl</a>) <- CRS('+proj=longlat +datum=WGS84')</div>
<div>summary(<a href="http://br.sl">br.sl</a>); class(<a href="http://br.sl">br.sl</a>)</div><div><br></div><div>est.x <- c(-53.35,-55.60,-51.80,-51.05,-50.13,-52.41,-52.10,-52.51)</div><div>est.y <- c(-33.51,-30.83,-29.45,-29.68,-29.05,-31.78,-32.03,-30.53)</div>
<div><a href="http://est.id">est.id</a> <- c(1:8)</div><div>est.df <- data.frame(est.x,est.y,<a href="http://est.id">est.id</a>)</div><div><br></div><div>est.sp <- est.df</div><div>coordinates(est.sp) <- c('est.x','est.y')</div>
<div>proj4string(<a href="http://br.sl">br.sl</a>) <- CRS('+proj=longlat +datum=WGS84')</div><div>est.sp</div><div><br></div><div>#png(filename="c:/temp/mapa.png",width=600,height=800)</div><div>plot(<a href="http://br.sl">br.sl</a>, axes=T)</div>
<div>plot(est.sp, pch=21, bg=3, add=T)</div><div>#dev.off()</div></div><div><br></div><div>###########<br></div><div><br></div>-- <br><font face="arial, helvetica, sans-serif"><br>================================================<br>
Éder Comunello<br><br>Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ)<br>Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]<br><br>Researcher at Embrapa Western Region Agriculture</font><div><font face="arial, helvetica, sans-serif">Dourados, MS, Brazil [</font>22 16.5'S, 54 49.0'W<span style="font-family:arial,helvetica,sans-serif">]</span></div>
<div><font face="arial, helvetica, sans-serif">================================================<br>UTC-03:00</font></div><div class="gmail_extra"><br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<br>
Date: Tue, 30 Oct 2012 21:09:26 -0200<br>
From: kin Huts <<a href="https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=kim-met@hotmail.com" target="_blank">kim-met@hotmail.com</a>><br>
To: <<a href="https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: [R-br] em plot.xy(xy.coords(x, y), type = type, ...) :<br>
objeto 'membros' não encontrado<br>
Message-ID: <SNT141-W57AFF1FD679C8C33FDDE18F1620@phx.gbl><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
<br>
<br>
Olá,<br>
Gostaria de uma ajuda para entender o erro encontrado e como corrigir esse script,Estou tentando gerar um mapa para o RS com a localização das estações que estou utilizando.<br>
<br>
<br>
require(MASS)require(maps)<br>
<br>
# GERAÇÃO DO MAPA COM AS ESTAÇÕES<br>
# carregando divisão política brasileiramapaBrasil <- read.table("BASE_BR.txt",sep="")<br>
# carregando informações de localização das estações meteorológicaslatlonRS <- read.csv2("coordenadas.txt",h=F,sep="")<br>
# gerando mapa em PNGpng(filename="mapa_grupos.png",width=600,height=800)<br>
# criando grafico em branco p/ o estado do RSmap(xlim=c(-80,-35),ylim=c(-45,10), type="n")map.axes()<br>
# plotando mapalines(mapaBrasil[,2], mapaBrasil[,1])<br>
# plotando estaçõespoints( latlonRS[,2], latlonRS[,1], pch=membros)<br>
# fechando PNGdev.off()<br>
<br>
dados que estou usando são: V1 V21 -33.51 -53.352 -30.83 -55.63 -29.45 -51.84 -29.68 -51.055 -29.05 -50.136 -31.78 -52.417 -32.03 -52.18 -30.53 -52.51<br>
e o erro que aparece é:<br>
<br>
Carregando pacotes exigidos: MASSCarregando pacotes exigidos: mapsErro em plot.xy(xy.coords(x, y), type = type, ...) : objeto 'membros' não encontrado<br>
<br></blockquote></div>
</div>