在给定点的多边形内创建线
给定一个空间多边形 library(sf) df % st_as_sf(coords = c("lon", "lat"), cr
..
sf相关内容
查找多线串和多边形Sf r之间的交点
我在地图上有一组空间坐标，并有一个穿过它们的多行字符串。我需要计算出这条线在每种颜色的多边形中花费的长度。 注意事项： 我的最小表示法不是很好，我的实际数据是一个数据。Frame有一个颜色列(在真实数据中是组)，一些标识列，以及一列多边形几何图形。我真的不知道如何制作，所以我的Reprex只有3个单独的多边形。 有时两种颜色的多边形会重叠，在这种情况下，我希望将这条线与重叠形状相交的
..
有没有一种方法可以删除ggploy外部的空格？
我使用代码创建一个静态地图，并将其包含在一个闪亮的应用程序中。我有一个问题，那里的地块大小比实际的地块大得多，所以布局看起来不太对劲。有什么方法可以减少绘图外部的空白吗？ 如果可能，我希望输出仅为绘图本身。 # reprex library(terra) library(sf) library(magrittr) library(ggplot2) lux
..
使用SF聚合栅格中的值
使用SF聚合栅格中的值 我需要的是为每个栅格聚合一些度量值。 假设我们有一些数据坐标和值，并且我想创建一个热图。 首先，我使用Simple Feature框架创建了一个网格和栅格。 现在我需要获取df中的每个坐标，并检查它们是否在一个栅格中。然后为每个栅格计算任何其他聚合函数的平均值。 # Packages -------------------------------------
..
如何在R中为.shp多边形文件设置CRS
我对在R中使用地理空间数据非常陌生。我在我的研究区域中有一个由几个象限组成的Shapefile，我根据“QUIDID”对其进行了分类。我可以很好地使用gglot和sf包来显示我的数据，但是，我的数据似乎不能正确显示。在形状文件中检查任何空间数据信息时，似乎没有： > proj4string(quads) Error in (function (classes, fdef, mtable)
..
R RASTER：：EXTRACT无法创建数据。帧
我希望使用RASTER：：EXTRACT和SF对象从栅格中提取值。然而，我遇到了一些困难。我希望创建一个数据框，其中包含单元格号(以派生XY值)、提取的单元格值和ID字段(如果我从多行提取)。我通过以下几点做到了这一点： dtm_values
..
SF&；dplyr：按相同坐标分组不起作用
我有一个从S2卫星栅格(10x10米)提取的数据集，其中有12个值(ras.df.ll)，但6个在一个块(T21JYG)中，第二个在另一个(T21JYG)中。 我想计算瓷砖之间相同(x，y坐标)的平均值，但没有成功。我找不到任何方法来识别第一个块中的第一行与第二个块中的第一行的坐标相同，只是我的数据集的末尾。在我的示例中： library(sf) library(sfheaders) lib
..
用geom_sf绘制基于点的核图
有人知道如何使用ggplot2和sf绘制基于事件发生的内核密度图吗？ 例如，考虑来自sp包的meuse数据集(假设每个点都是一个事件发生) library(sf) library(ggplot2) # data data(meuse, package = "sp") # as_sf meuse_sf
..
缓冲区重叠时计数点数
我在下面包含了我的所有代码和指向示例数据的链接。 简要说明：我的缓冲区重叠；我想计算离学校一定距离内的商店数量。 我特别想知道离学校1000米以内有多少家商店，离学校2000米以内有多少家商店，因为我想比较一下差异。当然，其中一些学校的缓冲是重叠的。因此，一家商店可能距离A校1500米，但距离B校只有750米。因此，它被算作距离一所学校1000米以内，应该只被计算在B校的1000米内，
..
为什么使用st_intersects而不是st_intersects？
st_intersection与st_intersects相比非常慢。那么，为什么不使用后者而不是前者呢？下面是一个带有小玩具数据集的示例，但是执行时间上的差异是巨大的，因为我的实际集只有62,020个点，与实际的地理区域边界相交。我有24 GB的RAM，st_intersects代码需要几秒钟，而st_intersection代码需要15分钟以上(可能更多，我还没有耐心等待...)。st_int
..
R|SF：我有点，每个点周围有2个缓冲区。如果较大的缓冲区重叠(但较小的缓冲区不重叠)，如何将这些点合并为单个单元？
我目前在地图上有点(学校)。每个人在点(学校)周围都有两个缓冲区。一个是450米，一个是250米。如果点重叠，我希望将它们视为单个单元(因为否则事情会变得复杂)，但我希望它们保持其覆盖的几何图形/面积。 因此，在这里给出的示例地图上，我希望将排名前三的学校/点合并为一个单位。&我希望它们保留其覆盖的面积，但仅将R作为一个单位计算。如果我使用"；st_Union"；函数，则必须
..
通过非SF列内部连接两个SF对象
我正在尝试使用内联接或左联接来联接两个SF数据帧。这些数据帧内部都有几何图形列。我一直收到错误： CHECK_JOIN(x，y)出错： y应为data.frame；对于空间连接，请使用st_joinfalse 可重现示例如下： df1
..
R中传单中聚类标记的聚合加权线串
我正在尝试绘制位置和加权连接线串.当我放大或缩小标记的聚类时，可以很好地调整.显示的集群标签是标记的聚合 node_val. 我想对线串做类似的事情，这样 该图不显示连接单个标记的蓝线，而是显示连接标记簇的线，并且 连接标记簇的新线串的宽度根据 wgt 变量进行自定义. 我希望下面的代码能说明问题: 库(dplyr)图书馆(传单)图书馆(旧金山)set.seed(123)N
..
使用传单中的 map_click 选择多个项目，链接到闪亮应用程序 (R) 中的 selectizeInput()
我想创建一个传单地图，您可以在其中选择多个多边形，这将更新闪亮应用程序中的 selectizeInput().这包括在 selectizeInput() 中删除选定的多边形时. 我稍微更改/更新了 the此处答案中的代码(使用 sf 而不是 sp 和更多 dplyr，我可以在其中计算出基本 R 是什么). 多边形可能会通过与 input$clicked_locations 绑定的 ob
..
如何使用 st_intersects() 等 sf 方法过滤 R 简单特征集合?
SF 是 R-Spatial 包，旨在处理诸如 dyplr 和管道之类的整洁语法. 我想对一个简单的特征集合对象做一个简单的空间过滤.给定一个简单的特征集合，我想从集合中返回满足某些几何条件的所有特征.特别是，我想找到与另一个对象相交的特征. SF 提供了函数 st_intersects(x,y,...) 来执行此操作，但我无法让它与 dplyr 一起使用. 我正在使用 R 3.
..
如何将 geom_sf 生成的地图放在 ggmap 生成的栅格之上
我尝试了以下代码: 库(ggplot2)图书馆(ggmap)图书馆(SF)nc 两层没有正确重叠；我尝试使用 coord_sf() 更改投影，但没有成功. 有什么建议吗?谢谢 解决方案 我自己也在努力解决这个问题，并在 这篇文章 我想出了一个解决方案.ggmap 对象的边界框在 WGS84 (EPSG:4326) 中，但实际栅格在 EPSG:3857 中.您必须破解 ggmap
..
R sf:LinearRing 的点不形成闭合线串
我正在尝试计算一组多边形的质心. 我的数据集geodata 包含五列，其中包括一列sfc_GEOMETRY 的geometry 列，有45759 行. 当我运行 sf::st_centroid(geodata) 时，我收到以下消息 CPL_geos_op("centroid", x, numeric(0), integer(0), numeric(0), Error in CPL
..
如何从st_sample(R包sf)中获取六边形样本?
我想从一个区域创建一些样本点.这些点必须给人一种密度的印象.我希望它们不是随机的，以避免人们认为它们是“真实"的观察结果.我希望它们在整个区域呈六边形分布.如何获得这样的样本?st_sample 带有 type = "hexagonal" 并不能解决问题. 一个可重复的例子: 库(sf)nc 最后一行代码给出了这个错误信息: seq_len(nrow(xy)) 中的错误:参数必须强制
..
根据另一个 sf 对象计算 sf 对象的面积
我有一个研究区(左)，还有一条道路(中间).我想计算由两种模式(右)相交产生的面积——这将产生 5 个子区域，这些子区域的总和为研究区域对象的总面积. 两个对象都是 sf 并且具有相同的投影(见下文). 如何使用 sf 库计算子区域? >str(myarea)“sf"和“data.frame"类:1 obs.共 4 个变量:$ 编号:编号 0$ 栖息地 : num 1$ 地区 : n
..
将 sf 转换为标记的 ppp
我使用下面的代码成功地将一个 sf 对象转换为 ppp: sf_owin = maptools::as.ppp.SpatialPointsDataFrame(as_Spatial(sf_points__flat)) sf_points__flat 看起来像这样: 简单的特征集合，包含 131 个特征和 3 个字段几何类型:多点尺寸:XYbbox:xmin:-1932934 ymin:
..