sf：：st\U交叉口：虚拟随机错误动作

Question

我有一个相对简单的任务要在R中完成：我有两个多边形层，一个“面片”层包含39个焦点多边形，另一个“土地覆盖”层包含一个焦点土地覆盖类型的多多边形。我需要剪裁这两个图层，以便在焦点多边形中仅保留焦点土地覆盖类型的范围。有时，这适用于sf：：st\U交叉口，有时适用于sf：：st\U差异和“负”土地覆盖层（包含所有非焦点土地覆盖类型的范围），有时两种方法都不起作用。起初，我认为这些不同的行为取决于地形的复杂性，但事实似乎并非如此。

我得到的错误是

Error in CPL_geos_op2(op, st_geometry(x), st_geometry(y)) :    
Evaluation error: TopologyException: Input geom 1 is invalid: Ring Self-intersection at or near point 
4372482.6526834015 5297568.4303682083 at 4372482.6526834015 5297568.4303682083.

因此，我使用sf:: st_is_simple（'patch'）检查了土地覆盖多边形和每个焦点多边形，在所有情况下都产生TRUE。

考虑以下三种情况：

我无法制作一个可复制的示例，因此我希望可以从示例图像中找出这里发生了什么。我在那里看不到任何模式，所以可能只有对st\U交叉口和st\U差异有深入了解的人才能找到解决方案。。。

Answer 1

您描述的错误不是随机的；环形自相交表示几何体无效。您应该能够通过sf：：st\u is\u valid（）进行测试。

当使用源自ESRI产品领域的空间对象时，已知会发生此错误，ESRI产品领域使用的有效性标准与OGC领域略有不同。

要解决此问题，您有几种选择：

• 过滤出有问题的几何体（通过基于sf：：st\u is\u valid（）的结果对空间对象进行子集设置）=只保留有效的几何体
• 尝试通过sf：：st\U make\U valid（）纠正几何图形-注意，这可能会导致几何图形改变，并且可能需要安装{lwgeom}软件包
• 通过sf：：st\U缓冲区（您的\uObject，0）在无效的空间对象上应用零宽度缓冲区的“魔法尘”。这种破解将强制创建新的几何体，可能会克服原始几何体中的错误

有关更多信息，请考虑sf包留档：https://r-spatial.github.io/sf/reference/valid.html