Rの各ポリゴンの面積を計算する

多くのポリゴンを持つシェープファイルを使用しています。

「area_sqkm」という名前のフィールドをもう1つ追加して、シェープファイルの各ポリゴンの面積を計算するにはどうすればよいですか？

shapefile r area

— HLe
ソース

どのライブラリをインストールし、何を試しましたか？

— GISKid

rgeos :: gArea（x、byid = TRUE）

— mdsumner 2016年

gis.stackexchangeへようこそ！このサイトで良い質問をすると、ある程度の調査結果が表示されることが期待されます。つまり、これまでに試したことと、該当する場合はコードを記述します。詳細については、FAQを確認してください。

— アンダーダーク

できるよ

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

crsが経度/緯度、または距離単位としてメートルを使用していると仮定します

— ロバート・ハイマンス
ソース

あるraster::area()よりも正確、非投影データにrgeos::gArea()投影されたデータと？

— Richard DiSalvo

かもね。投影が歪むからです。ポリゴンのノード間距離が大きい場合、投影によって領域がさらに歪む可能性があります（geosphereビネットを参照）。

— Robert Hijmans