R语言 提取tif地图中特定面积的数据

你可以使用 R 语言中的一些地图处理库，如 rgdal 和 raster，来提取特定面积的数据。

下面是一个简单的例子，假设你有一个名为 "map.tif" 的 tif 地图文件，你想提取其中一个名为 "my_area" 的多边形区域中的数据。

首先，你需要读取地图文件：

library(raster)
map <- raster("map.tif")

然后，你需要创建一个 SpatialPolygons 对象来表示你要提取的区域。这个对象可以手动创建，也可以从文件中读取。

这里我们手动创建一个简单的矩形区域：

library(sp)
poly <- SpatialPolygons(list(Polygons(list(Polygon(cbind(c(0, 0, 1, 1), c(0, 1, 1, 0)))), "my_area")))

接下来，你可以使用 extract() 函数来提取这个区域中的数据：

area_data <- extract(map, poly)

这将返回一个包含所有像素值的向量，其中只包括位于你指定区域内的像素。

你还可以使用其他参数来进一步控制提取过程，例如使用 mask 参数来指定一个遮罩图层，只提取两个图层中重叠的像素数据。