Luciano e Tito,
Obrigada pela ajuda, mas infelizmente não deu certo pois eu tenho duas matrizes diferentes, e a sintaxe que vocês me passaram serviria para vetores.
Eu usei um pouco diferente e deu certo:
coord <- coord1 <- coord2 <- matrix(NA,480,216)
for(i in 1:480){
for(j in 1:216){
if (lat[i,j] < latmax
& lat[i,j] > latmin)
{ if(lon[i,j] < lonmax
& lon[i,j] > lonmin)
coord[i,j] <- 1
coord1[i,j] <- lat[i,j]
coord2[i,j] <- lon[i,j]
}
}
}
Agora a minha resposta são matrizes de NA e 1 quando a lat/lon está contida na área do polígono.
Obrigada pelas dicas,
Abc,
Heloise
---
Oc. Heloise Pavanato
Pós-graduação em Oceanogra Biológica
Instituto de Oceanografia
Universidade Federal do Rio Grande