
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