0

我有一系列二进制栅格图层（ascii 文件），显示欧洲和非洲某个物种的存在/不存在。该文件基于未投影的纬度/经度 (WGS84) 数据。我的目标是使用 R 计算存在区域（我无权访问 ArcGIS）。

我知道 raster 包有计算面积的功能，但我担心这对于未投影的数据不准确。我也看过 raster 包中的 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

有谁知道准确计算存在区域的方法，考虑地球曲率？

谢谢！

4

1 回答 1

0

我不知道您的文件发生了什么（为什么会收到警告＃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') 您可以获得存在的总面积。

于 2015-01-27T05:40:06.417 回答