ヨーロッパとアフリカの種の有無を示す一連のバイナリ ラスター レイヤー (ASCII ファイル) があります。このファイルは、投影されていない緯度/経度 (WGS84) データに基づいています。私の目的は、R を使用してプレゼンス エリアを計算することです (ArcGIS にアクセスできません)。
ラスター パッケージに面積を計算する機能があることは知っていますが、これが投影されていないデータに対して正確ではないことが心配です。ラスター パッケージの cellStats 関数も調べました。これを使用して、占有されているセルの数を「合計」できますが、これには同じ問題があると感じています。
jan<-raster("/filelocation/file.asc")
jan
class : RasterLayer
dimensions : 13800, 9600, 132480000 (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333 (x, y)
extent : -20, 60, -40, 75 (xmin, xmax, ymin, ymax)
coord. ref. : NA
data source : "/filelocation"
names : file.asc
values : -2147483648, 2147483647 (min, max)
area(jan)
class : RasterLayer
dimensions : 13800, 9600, 132480000 (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333 (x, y)
extent : -20, 60, -40, 75 (xmin, xmax, ymin, ymax)
coord. ref. : NA
names : layer
values : 6.944444e-05, 6.944444e-05 (min, max)
Warning messages:
1: In .local(x, ...) :
This function is only useful for Raster* objects with a longitude/latitude coordinates
2: In .rasterFromRasterFile(grdfile, band = band, objecttype, ...) :
size of values file does not match the number of cells (given the data type)
cellStats(jan,"sum")
[1] 3559779
地球の曲率を考慮して、プレゼンス エリアを正確に計算する方法を知っている人はいますか?
ありがとう！