Olá, <div>Seria uma função para recortar uma quadricula dentro de um intervalo de coordenadas?</div><div>Se for, talvez isso que segue abaixo possa ajudar:</div><div><br></div><div><div>#---------------------------------------------------------------------------</div>
<div><br></div><div>coords<-matrix(sample(0:100,200,r=T),100,2);colnames(coords)<-c("Lat","Lon")</div><div><br></div><div>coordinates<-c(0,50,50,100)#Sequencia de valores: </div><div>#1-Latitude Minima</div>
<div>#2- Latitude Maxima</div><div>#3- Longitude Minima</div><div>#4- Longitude Maxima</div><div><br></div><div>coords.1<-coords[c(which((coords[,"Lat"])>=(coordinates[1]))),]#.lat.min</div><div>coords.2<-coords.1[c(which((coords.1[,"Lat"])<=(coordinates[2]))),]#.lat.max</div>
<div>coords.3<-coords.2[c(which((coords.2[,"Lon"])>=(coordinates[3]))),]#.long.min</div><div>coords.final<-coords.3[c(which((coords.3[,"Lon"])<=(coordinates[4]))),]#.long.max</div></div><div>
#Aqui no final saem somente os pontos que estão no intervalo determinado acima</div><div><br></div><div>#---------------------------------------------------------------------------
</div><div><br></div><div>Espero ter ajudado.</div><div>Abraços</div><div>Luciano</div><div><br></div><div><br><br><div class="gmail_quote">2012/8/21 Heloíse Pavanato <span dir="ltr"><<a href="mailto:helopavanato@gmail.com" target="_blank">helopavanato@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Boa tarde a todos,<div><br></div><div>Estou com um problema em um loop.</div><div>Tenho duas matrizes de 480 por 216 com latitudes e longitudes e quero selecionar dentre estas somente as que estão dentro de determinado limite (latLim e lonLim), que são pares.</div>
<div><br></div><div>Estou usando:</div><div><br></div><div><div>lat <- read.table(file="lat.txt")</div><div>lon <- read.table(file="lon.txt")</div></div><div><br></div><div>library(R.matlab) # para ler .mat</div>
<div>sx <- readMat('SX.mat',maxLength=NULL, fixNames=TRUE,verbose=F,sparseMatrixClass='SparseM')</div><div><br></div><div># limite de longitude</div><div>lonSX <- sx$SX[3,1,1]
</div><div>lonlim <- as.numeric(lonSX$Lon[1:130])</div><div><br></div><div># limite de latitude</div><div>latSX <- sx$SX[4,1,1]
</div><div>latlim <- as.numeric(latSX$Lat[1:130])</div><div><br></div><div>####</div><div>latA <- lonA <- matrix(NA,480,216)</div><div><br></div><div>for(i in 1:(dim(lat)[1])){</div>
<div> for(j in 1:(dim(lon)[2])){</div><div> if (c(lat[i,j],lon[i,j]) > min(latLim,lonLim) # seleciona a lat/lon só para a área do polígono</div><div> & c(lat[i,j],lon[i,j]) < max(latLim,lonLim))</div>
<div> latA[i,j] <- lat[i,j]</div><div> lonA[i,j] <- lon[i,j]</div><div> }</div><div> } </div><div><br></div><div>Mas a sintaxe não está correta pois eu não estou sabendo dizer que cada par latLim,lonLim é um delimitador de lat e lon.</div>
<div><br></div><div>Se alguém puder me dar uma mãozinha, por favor.</div><div><br></div><div>Os arquivos estão em anexo.</div><div><br></div><div>Obrigada,</div><div>Heloise</div><div><br></div><div>---</div><div><span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Oc. Heloise Pavanato</span><br style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
<span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Pós-graduação em Oceanogra Biológica</span><br style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
<span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Instituto de Oceanografia</span><br style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
<span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Universidade Federal do Rio Grande</span>
</div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
<blockquote style="margin:0 0 0 40px;border:none;padding:0px"></blockquote><span style="color:rgb(51,51,51);font-family:'Arial Narrow',sans-serif">Luciano
F. Sgarbi</span><br><blockquote style="margin:0 0 0 40px;border:none;padding:0px"></blockquote><font color="#333333" face="'Arial Narrow', sans-serif">Mestrando em Ecologia e Evolução - UFG</font><br><blockquote style="margin:0 0 0 40px;border:none;padding:0px">
</blockquote><span style="color:rgb(51,51,51);font-family:'Arial Narrow',sans-serif">Laboratório
de Ecologia de Insetos</span><br><blockquote style="margin:0 0 0 40px;border:none;padding:0px"></blockquote><div><span style="color:rgb(51,51,51);font-family:'Arial Narrow',sans-serif">Cel.
(62)8174-2262 Lab. (62)3521-1732</span><div>
<p></p><p></p></div>
<div style="display:inline"></div></div><br>
</div>