[R-br] Calcular a área em Km2 a partir de um shapefile

[Pedro Andrade]

Prezada Carolina,

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.

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).

um abraço,

Pedro


On 26-03-2011 19:51, Carolina Pena wrote:
> /Alguém conhece algum pacote/função no R capaz de calcular a área  em km2 de um mapa (shapefile)./
>
> /Eu li o mapa no R utilizando o comando 'readShapePoly'; consigo acessar a "area" (slots) do objeto através dos seguintes comandos:/
>
> /area = sapply(slot(map2, "polygons"),
>   function(i) slot(i, "area"))//
>   /
>
> /No entanto, não sei como transformar essa informação em km2/
>
> /Obrigada,/
>
> /Carolina/
>
> /   /
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br