问题标签 [sf]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
r - 如何获得带有轮廓的相交线的坐标 - R
我有一个从互联网上传的形状。基于其计算出的质心,我想从中绘制一条 50 度线并找到他与轮廓相交的坐标。知道我该怎么做吗?
谢谢。
脚本:
r - 无缝拟合两个 sf 多边形
问题
假设我们有两个无缝边界的 shapefile。只是,他们没有。有没有一种方法可以迫使它们彼此无间隙地粘在一起?
具体案例
我有两个 shapefile:一个用于欧洲地区 - REG
,另一个用于邻国 - NEI
。两个 shapefile 都取自Eurostat 存储库,应该可以很好地组合在一起;但有一些小的差距。另外,我需要简化多边形,然后差距变得非常明显。
我能想到的最好的
我尝试了几种方法,但没有成功。实现我看到的预期结果的唯一方法需要以下步骤:
- 用我的 shapefile 之间的边界创建一条线 sf;
- 从这条线创建一个足够大的缓冲区多边形以覆盖所有间隙;
- 将此缓冲区加入并解散到邻居的 shapefile 中 -
NEI
; NEI
用REG
shapefile剪掉扩展的部分。
显然,这是一个相当笨拙的解决方法。
有更好的方法吗?
此要点中可重现的示例
一个最小的例子
r - 通过R中的sf将经度和纬度序列转换为多边形
我有五个经度和纬度,形成这样的形状。
我怎样才能使用sf
这样的包轻松地将它们转换为 sf 多边形数据框?
r - 如何使用 R sf 包从底层网格中裁剪多边形?
我想知道它是否比下面使用 rstats sf 的示例更简单,以计算底层网格中多边形 shapefile 的(分数)面积覆盖率。到目前为止,我还没有找到可以为我完成这项工作的功能。这是我的例子:
现在开始计算,我首先通过从边界框('bpol')中减去'mpol'来反转多边形('ipol')并计算新的覆盖区域。
如果一个网格单元被覆盖两次,然后从网格中断处裁剪倒多边形,它会变得更加复杂
错误消息:“CPL_geos_op2(op, st_geometry(x), st_geometry(y)) 中的错误:评估错误:TopologyException:输入几何 1 无效:在 9.75 50.5 处或附近的嵌套壳在 9.75 50.5。
现在我的问题是,是否已经有一个功能可以完成这项工作?我错过了文档中的某些内容吗?当网格单元被裁剪多次时,如何解决错误?
谢谢和亲切的问候约尔格
r - R中的错误:没有适用于“st_centroid”的方法应用于“NULL”类的对象
编辑:我后来意识到错误出在我的完整数据集中,而不是代码中。我很抱歉。这对其他人可能有用也可能没有用。
我有一个小标题,有一个名为“几何”的列,它是多边形(我使用的是 sf 空间包)。我正在尝试向作为每个多边形质心的小标题添加一列。但是,我收到一个错误:“mutate_impl(.data, dots) 中的错误:评估错误:没有适用于 'st_centroid' 的方法应用于“NULL”类的对象。” 我看过一些关于类似主题的其他帖子,但不是很相似,以至于我无法弄清楚。如果我错过了什么,我道歉!
我不熟悉使用管道和 tidyverse,所以那里可能出了点问题。任何帮助,将不胜感激!
这是我的代码示例:
center1 = function(df){
centroid = st_centroid(df)
}
table2 = x %>%
mutate(center.1 =map(geometry, center1))
这是我的数据的一个非常小的样本(只有一行)——x
代码中的变量(对不起长度!):我加载了 tidyverse、tibble、dplyr、purrr 和 sf。