Olá,
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 .matsx <- readMat('SX.mat',maxLength=NULL, fixNames=TRUE,verbose=F,sparseMatrixClass='SparseM')# limite de longitudelonSX <- sx$SX[3,1,1]lonlim <- as.numeric(lonSX$Lon[1:130])# limite de latitudelatSX <- 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
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Laboratório de Ecologia de Insetos