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).//Carolina/
/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,/
/ /
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br