Calcular a área para cada polígono em R

Estou trabalhando com um shapefile que possui muitos polígonos.

Como adiciono mais um campo chamado "area_sqkm" e calculo a área para cada polígono no shapefile?

shapefile r area

— HLe
fonte

Quais bibliotecas você instalou e o que você tentou?

— GISKid

rgeos :: Garea (x, byid = TRUE)

— mdsumner

Bem-vindo ao gis.stackexchange! Observe que uma boa pergunta neste site deve mostrar algum grau de pesquisa de sua parte, ou seja, o que você tentou e, se aplicável, o código até o momento. Para mais informações, você pode conferir nossas perguntas frequentes .

— Underdark

Você pode fazer

library(raster)
x <- shapefile('file.shp')
crs(x)
x$area_sqkm <- area(x) / 1000000

Supondo que seus crs sejam longitude / latitude ou com metro como unidade de distância

— Robert Hijmans
fonte

os raster::area()dados não rgeos::gArea()projetados são mais precisos do que os dados projetados?

— Richard DiSalvo

Pode ser. Porque as projeções distorcem. Se a distância entre nós do polígono for alta, a projeção poderá distorcer ainda mais a área (consulte a geospherevinheta).

— Robert Hijmans