[R-br] ajuda loop

[Tito Conte]

Heloíse,

experimente fazer assim:

lat <- read.table(file="lat.txt")
lon <- read.table(file="lon.txt")

library(R.matlab)   # para ler .mat
sx <- readMat('SX.mat',maxLength=NULL,
fixNames=TRUE,verbose=F,sparseMatrixClass='SparseM')

# limites de longitude
lonSX <- sx$SX[3,1,1]
lonmin <- min(as.numeric(lonSX$Lon[1:130]))
lonmax <- max(as.numeric(lonSX$Lon[1:130]))

# limites de latitude
latSX <- sx$SX[4,1,1]
latmin <- min(as.numeric(latSX$Lat[1:130]))
latmax <- max(as.numeric(latSX$Lat[1:130]))

indice=which(lat<latmax) #indices das latitudes menores do que a máxima
indice=which(lat[indice]>latmin) #indices das latitudes maiores do que a
minima
indice=which(lon[indice]<lonmax) #indices das longitudes menores do que a
máxima
indice=which(lon[indice]>lonmin) #indices das longitudes maiores do que a
minima

no final vc terá os indices das coordenadas que encerram a área de interesse





Atenciosamente,
Tito Conte

Instituto Oceanográfico
Universidade de São Paulo
LabPesq - Laboratório Oceanografia Pesqueira
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120821/0f359a6e/attachment.html>