Boa tarde a todos,

Estou com um problema em um loop.
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.

Estou usando:

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')

# limite de longitude
lonSX <- sx$SX[3,1,1]
lonlim <- as.numeric(lonSX$Lon[1:130])

# limite de latitude
latSX <- sx$SX[4,1,1]
latlim <- as.numeric(latSX$Lat[1:130])

####
latA <- lonA <- matrix(NA,480,216)

for(i in 1:(dim(lat)[1])){
    for(j in 1:(dim(lon)[2])){
        if (c(lat[i,j],lon[i,j]) > min(latLim,lonLim)      # seleciona a lat/lon só para a área do polígono
        & c(lat[i,j],lon[i,j]) < max(latLim,lonLim))
        latA[i,j] <- lat[i,j]
        lonA[i,j] <- lon[i,j]
      }
    } 

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.

Se alguém puder me dar uma mãozinha, por favor.

Os arquivos estão em anexo.

Obrigada,
Heloise

---
Oc. Heloise Pavanato
Pós-graduação em Oceanogra Biológica
Instituto de Oceanografia
Universidade Federal do Rio Grande