sp相关内容

空间多边形列表中的R-空间多边形数据帧

我正在寻找从SpatialPolygons列表创建SpatialPolygonsDataFrame的方法? 下面是一个多边形列表的示例,应从该列表创建包含该列表的所有多边形的SpatialPolygonsDataFrame。 编辑:SpatialPolygonsDataFrame必须从SpatialPolygons列表创建!因为我的原始数据不包含SpatialPolygons作为单独的值, ..
发布时间:2022-05-29 21:32:15 其他开发

包含面的空间点的子集

我有一个SpatialPolygonsDataFrame(Spolydf)和一个SpatialPointsDataFrame(Spintdf)。这些层具有不同的范围,但重叠。 我可以使用 选择位于多边形内的点 fall.within.poly ..
发布时间:2022-05-29 21:21:42 其他开发

SpatialPolygonsDataFrame对象的坐标()值?

我正在尝试获取空间多边形数据框中多边形的伪重心。今天我偶然发现了coordinates函数,它实际上为SpatialPolygonsDataFrame返回了一些东西。 不幸的是,我在coordinates的帮助下没有找到有关SpatialPolygonsDataFrame的值。谁能告诉我这些坐标是什么? 推荐答案 它是多边形质心。找到源代码here,查找函数FindCG。计算的方 ..
发布时间:2022-05-29 18:26:51 其他开发

如何将多边形修改为更改其槽的洞(SpatialPolygons)

编辑: 正如Edzer Pebesma在评论中所建议的,向多边形添加洞的推荐方法不是修改槽,而是重建多边形,如in this related question所示。 原始问题 在SpatialPolygons-class的帮助下,我试图将一个多边形修改为另一个多边形洞,但当该“洞多边形”显示为另一个多边形的边界时,其内部的颜色与其他多边形相同。 我做错了什么? 使用d ..
发布时间:2022-05-29 18:18:41 其他开发

R:如何计算网格的每个单元格中有多少个点?

我制作了一个参考网格,单元格50x50m,基于一只项圈动物的GPS位置。我想在ArcGIS中执行相当于空间连接的操作,并计算每个单元格中的点数。 我已经使用SpatialPointsDataFrame对象创建了一个参考网格(数据框已经使用UTM坐标系投影) RESO ..
发布时间:2022-05-29 17:10:39 其他开发

如何绘制多个经纬度坐标的多边形并计算曲面面积?

我有以下数据:https://ufile.io/p9s0le6g ...它由485个经度/经度观测组成。我最终想要计算这个数字的表面积(以平方米或平方公里为单位)。我想用数据点(外部点)创建一个多边形,然后计算曲面。然而,我想不出一个好的多边形。 我的坐标数据点如下所示: 我想要一个看起来像那样(或类似)的多边形: 以下是我尝试的内容: library(sp) df ..
发布时间:2022-04-11 16:44:26 其他开发

如何使用 data.table 有效地计算坐标对之间的距离:=

我想找到最有效(最快)的方法来计算经纬度坐标对之间的距离. 提出了一个不太有效的解决方案(这里) 使用 sapply 和 spDistsN1{sp}.我相信如果在 data.table 中使用 spDistsN1{sp} 和 := 运算符,这可以更快,但我没有无法做到这一点.有什么建议吗? 这是一个可重现的示例: # 加载库库(数据表)图书馆(dplyr)图书馆(SP)图书馆(rge ..
发布时间:2022-01-13 18:40:12 其他开发

R sp:Polygon的面积单位

我用R读入了一个shape文件,用sp包分析了油田的多边形(WGS84的经纬度)及其各自的区域.不幸的是我不知道面积输出的单位.例如.面积输出平均为 0.85 个单位(最大 4.34 个),这可能不是以平方公里为单位,因为这对于油田来说太小了. 有谁知道sp包中Polygons的面积输出单位是多少?非常感谢! 解决方案 要在纬度坐标系中正确计算多边形的面积,最好使用“spTransf ..
发布时间:2021-11-17 04:20:43 其他开发

所有点到第一个点的距离

我在街道上有一些点,我正在尝试计算所有点与街道上第一个点之间的距离 这是我的 sf 对象: library(rgeos)图书馆(rgdal)图书馆(SF)#我实际上在 5 公里的一段上有多个点df 我试过 st_distance 但距离不正确.我得到的是让每一个点都到街道的起点直到终点.所以基本上它必须从 0 到 5000m. 解决方案 如果在 st_distance() 中加 ..
发布时间:2021-07-23 19:44:19 其他开发

R中的空间线起点和终点

我正在尝试使用 sp 包来访问线串的起点和终点,类似于 ST_StartPoint 和 ST_EndPoint使用 psql 生成. 无论我如何尝试访问该行,我都会收到错误或 NULL 值: >onetrip@lines[[1]][1]onetrip@lines[[1]][1] 中的错误:“S4"类型的对象不可子集>onetrip@lines@Lines@coords错误:尝试从没有插槽的 ..
发布时间:2021-07-23 19:44:12 其他开发

更改 SF 对象的 CRS

我有一些“道路"作为 sp 对象: class : SpatialLinesDataFrame特点:17360范围:490176.4、567680.9、148639.1、212821(xmin、xmax、ymin、ymax)坐标.参考: +init=epsg:27700 +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_ ..
发布时间:2021-07-23 19:43:43 其他开发

使用 sf 和 sp 在 R 中的投影差异

我有一个从 GeoTIFF 转换为 shapefile 的网格.我想将 shapefile 转换和导出为 GeoPackage 并更改投影,以便在 GIS 中打开时使用英国国家网格作为地理坐标系.然而,这似乎只适用于 sp 而不是 sf(它似乎没有保留像数据这样的方面). 这是一个问题,因为我想导出包含多个图层的 GeoPackages,您目前只能在 sf 而不是 sp.我做错了什么吗? ..
发布时间:2021-06-22 18:30:58 其他开发

改变形状文件的投影

我正在尝试将德国形状文件的投影从 NA 更改或分配到 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0,但不知何故它不能正常工作. 可重现的示例:Shapefile 和其他文件可以在此处下载: 我尝试的是以下内容: 库(地图工具)图书馆(sp)图书馆(rgeos)图书馆(rgdal)投影.x 所以,问题 ..
发布时间:2021-06-22 18:30:29 其他开发

在城市上创建六边形网格,并与lon/lat点关联(在R中)

我已经对此进行了一段时间的研究,但是还没有找到适合我需要的解决方案,或者我可以进行充分的调整以适合我的情况: 我有多个城市的大型汽车共享数据集,在这些城市中每个位置我都有充电需求(例如,行= carID,55.63405、12.58818,充电需求).我现在想将整个城市(上面的例子是哥本哈根)的区域划分为一个六边形网格,并用ID标记每个停车位置(例如,行= carID,55.63405、12 ..
发布时间:2021-05-10 19:18:12 其他开发

对于一个数据集中的每个点,计算到第二个数据集中的最近点的距离

尝试为 SpatialPointsDataFrame 中的每个点查找到第二个 SpatialPointsDataFrame 中最接近点的距离(等同于ArcGIS for中的“最近"工具)两个 SpatialPointDataFrames ). 我可以通过使用 gDistance 计算所有成对距离并采用 min ( 在这里,您可以使用它快速为 SpatialPoints 对象 b 中的每个点 ..
发布时间:2021-05-10 19:17:43 其他开发

计算分组数据中点之间的欧式距离

在下面的数据中(包含在 dput 中),我对三个人( IndIDII )进行了重复观察(经度和经度).请注意,每个人都有不同数量的位置. >达特印度印度WintLat WintLong1 BHS_265 BHS_265-2015 47.61025 -112.72102 BHS_265 BHS_265-2016 47.59884 -112.70893 BHS_770 BHS_770-2016 42 ..
发布时间:2021-04-30 20:54:43 其他开发