1

我无法从大型形状文件（916.2 Mb，4618197 个元素 - 从这里提取点数据：https ://earthdata.nasa.gov/data/near-real-time-data/firms/active-fire-data ) 在 R 中。我正在使用readShapeSpatialinmaptools来读取形状文件，这需要一段时间但最终可以工作： 

worldmap <- readShapeSpatial("shp_file_name")

然后我有一个data.frame我想要提取数据的坐标。然而，R 真的在为此苦苦挣扎，要么失去连接，要么死机，即使只有一组坐标！

pt <-data.frame(lat=-64,long=-13.5)
pt<-SpatialPoints(pt)
e<-over(pt,worldmap)

任何人都可以建议我以更有效的方式执行此操作吗？

还是我需要在更强大的东西上运行这个脚本（目前使用的是 2.3 GHz 处理器的 mac mini）？

非常感谢！

4

1 回答 1

1

“点数据”是指经度和纬度坐标吗？如果是这种情况，您可以通过以下方式获取 shapefile 的基础数据：

worldmap@data

您可以像查看任何其他数据框一样查看它，例如：

View(worldmap@data)

您还可以以与通常相同的方式访问此数据框中的列，除非您不需要@data，例如：

worldmap$LATITUDE

最后，建议使用readOGRfrom rgdalpackage 而不是maptools::readShapeSpatial前者读取 CRS/projection 信息。

于 2015-05-11T13:04:57.293 回答