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
_______________________________________________
R-br mailing list
[email protected]
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.
