使用R中的图像进行区域估计

Question

所以我正在寻找一种仅使用地图图像估算区域面积的方法。我这样做的原因是我想计算在海平面上有一定增加时会丢失的区域，我找不到任何类型的元数据，只有地图（图像格式）。以下是此类地图的链接：

所以我想到的是使用EBimage包将此图像转换为灰度图像，然后使用像素强度作为标准来计算代表潜在威胁区域的像素数。

我的问题是可能的吗？我们怎样才能将像素强度作为标准？如果还有其他方法可以解决这个问题？

此外，如果有办法获取用于绘制我不知道的地图的元数据，请告诉我。

谢谢大家。

编辑：

感谢hrbrmstr我能够使用rgdal包读取网格数据int R.因此，为了计算我尝试使用rgeos包的区域，但来自CRESIS的数据集没有形状文件。那么我们如何定义多边形并计算面积呢？ 对不起，这个问题似乎很愚蠢。这是我第一次处理空间数据和分析

Answer 1

数据位于ESRI文件中：

library(rgdal)

grid_file_1m <- new("GDALReadOnlyDataset", "/full/path/to/inund1/w001001.adf")
grid_1m <- asSGDF_GROD(grid_file_1m, output.dim=c(1000, 1000))

plot(grid_1m, bg="black")

grid_1m_df <- as.data.frame(grid_1m)

str(grid_1m_df)
## 'data.frame': 2081 obs. of  3 variables:
##  $ band1: int  1 1 1 1 1 1 1 1 1 1 ...
##  $ x    : num  -77.2 -76.9 -76.5 -76.1 -75.8 ...
##  $ y    : num  83.1 83.1 83.1 83.1 83.1 ...