Question

我有一系列二进制栅格图层（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

任何人都知道如何准确计算存在区域，考虑地球曲率？

谢谢！

Answer 1

我不知道您的文件发生了什么（为什么会收到警告＃2）。但这是一个解决方法

r <- raster(nrow=13800, ncol=9600, xmn=-20, xmx=60, ymn=-40, ymx=75)
# equivalent to r <- raster(jan)
x = area(r)
x
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. : +proj=longlat +datum=WGS84 
data source : c:\temp\R_raster_Robert\2015-01-26_213612_1208_85354.grd 
names       : layer 
values      : 0.2227891, 0.8605576  (min, max)

现在，每个单元格的面积为km2。通过将这些值与带有存在/不存在值的Raster对象相乘，然后使用cellStats（，'sum'），您可以获得存在的总面积。

计算二进制未投影栅格的占用面积

1 个答案: