Luciano e Tito,<div>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.</div><div><br></div><div>Eu usei um pouco diferente e deu certo:</div>

<div><br></div><div><div>coord <- coord1 <- coord2 <- matrix(NA,480,216)</div><div><br></div><div>for(i in 1:480){</div><div>    for(j in 1:216){</div><div>        if (lat[i,j] < latmax</div><div>        & lat[i,j] > latmin)</div>

<div>          { if(lon[i,j] < lonmax</div><div>            & lon[i,j] > lonmin)</div><div>                    coord[i,j] <- 1</div><div>                    coord1[i,j] <- lat[i,j]</div><div>                    coord2[i,j] <- lon[i,j]</div>

<div>    }</div><div>   }</div><div>  }</div></div><div><br></div><div>Agora a minha resposta são matrizes de NA e 1 quando a lat/lon está contida na área do polígono.</div><div><br></div><div>Obrigada pelas dicas,</div>
<div>
<br></div><div>Abc,</div><div>Heloise</div><div><br></div><div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">---</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
Oc. Heloise Pavanato<br>Pós-graduação em Oceanogra Biológica<br>Instituto de Oceanografia<br>Universidade Federal do Rio Grande</div></div>