Tenham cuidado. ângulo não é medida de distância, é de posicão.<br><br><div class="gmail_quote">Em 27 de março de 2011 12:24, Pedro Andrade <span dir="ltr"><<a href="mailto:pedro.andrade@inpe.br">pedro.andrade@inpe.br</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Prezada Carolina,<br>
<br>
todo dado espacial está (ou pelo menos deveria estar) associado a uma projeção cartográfica. Estas projeções estão divididas basicamente em dois grupos, as que trabalham com valores em metros e as que trabalham com ângulos. Se você tem os valores em metros, basta dividir a área por 10^6. Entretanto, se os seus valores estão em ângulos, as áreas certamente estão em ângulos quadrados porque praticamente todos os algoritmos espaciais não diferenciam ângulos de metros e consideram que os dados estão em metros. Se for o caso, você tem que primeiramente converter o seu dado para uma projeção que trabalha com metros para então poder calcular as áreas.<br>

<br>
O problema então se transforma em saber qual é a projeção do seu dado. Se a pessoa que gerou o dado que você está usando for cuidadosa, além dos arquivos .shp, .shx e .dbf, existe o arquivo .prj que contem tudo o que você precisa, e pode ser aberto por qualquer bloco de notas da vida. Caso contrário, você deve procurar saber com uma pessoa que já usou o dado. Se tudo mais falhar, resta a solução POG de verificar o próprio dado para se ter uma estimativa da projeção. Normalmente ângulos usam valores após a vírgula (10.03342, 10.03351, por exemplo), enquanto que metros usam valores antes da virgula (13245.2, 13230.1).<br>

<br>
um abraço,<br>
<br>
Pedro<br>
<br>
<br>
On 26-03-2011 19:51, Carolina Pena wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
/Alguém conhece algum pacote/função no R capaz de calcular a área  em km2 de um mapa (shapefile)./<div><div></div><div class="h5"><br>
<br>
/Eu li o mapa no R utilizando o comando 'readShapePoly'; consigo acessar a "area" (slots) do objeto através dos seguintes comandos:/<br>
<br>
/area = sapply(slot(map2, "polygons"),<br>
  function(i) slot(i, "area"))//<br>
  /<br>
<br>
/No entanto, não sei como transformar essa informação em km2/<br>
<br>
/Obrigada,/<br>
<br></div></div>
/Carolina/<br>
<br>
/   /<div class="im"><br>
<br>
<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
</div></blockquote><div><div></div><div class="h5">
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Luis Iván Ortiz Valencia<br>Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ<br>Estatístico Msc.<br>Spatial Analyst Msc.<br>...............................................<br>
Curriculum Lattes<br><br><a href="http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4778724J3" target="_blank">http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4778724J3</a><br>...............................................<br>
Aquarela Cusco Hostel Perú<br><br><a href="http://www.aquarelacuscohostel.com/" target="_blank">http://www.aquarelacuscohostel.com/</a><br>...............................................<br>