如何计算面积的多边形从一个大shapefile

Question

总结:

我试图在r中计算大量多边形的面积。我读过几篇关于如何做的文章(例1

扩展说明：

我实际上是在计算澳大利亚维多利亚州的房产面积。多边形表示这些属性。我从Spatial Datamart下载了所有维多利亚州的VicMaps简化模型1和2。然而，考虑到形状文件的大小，我不得不将搜索范围缩小到一个地方政府区域（LGA），并计算多边形区域（仅用于测试）。形状文件为15.5MB。

library(raster)
x <- shapefile("D:/Downloads/SDM616230/ll_gda94/shape/lga_polygon/ballarat/VMPROP/PROPERTY_PRIMARY_APPROVED.shp")
crs(x)
x$area_sqkm <- area(x) / 1000000

这是有效的，但它不是我的问题的实际解决方案，因为维多利亚州有许多地方政府机构，我计划最终在昆士兰州和新南威尔士州遵循同样的程序。

然而，试图加载一个更大的shapefile不起作用，并导致以下错误代码“错误:内存耗尽(达到极限？)".

我尝试过使用readShapePoly、readogr、st_read和read_sf将大形状文件转换为R，但它们都不起作用。我觉得文件太大了。我尝试在read_sf中使用select查询，以减小正在读取的文件的大小，但这也不起作用。我在网上读到，我应该设法将shapefile分割成我需要的数据，以缩小大小，但我不知道如何做到这一点。

希望能帮到你。

Answer 1

显然，文件太大，无法容纳单个盒子。我认为选择是

1)将文件分割成较小的文件，逐个处理。看见

https://gis.stackexchange.com/questions/195508/split-a-shapefile-into-smaller-files-on-linux-command-line

2)用一些dbms或者数据仓库来做，他们自动做这样的批处理。