问题标签 [rgeo-shapefile]

我有数百个没有坐标参考系统的 shapefile。我的目标是在 WorldClim 栅格图层上叠加空间多边形。我以前使用过这种方法，没有任何问题。但是，这次我的 shapefile 中的坐标对我来说很奇怪。多边形内边界框和坐标的每个坐标由 8 位数字组成，没有逗号或点来分隔小数点。这是其中一种形状的边界框：

这与 WorldClim 栅格图层的边界框明显不同。

当我尝试使用 plot 命令将 shapefile 覆盖在栅格图层上时，没有任何反应。

我明白这是一个投影问题。然后我尝试使用 CRS 函数在 shapefile 中分配 WorldClim 栅格图层的相同坐标系。但是，结果保持不变（即 shapefile 不在光栅上）。在序列中，我尝试使用 rgdal 包中的 spTransform 函数重新投影 shapefile 坐标。但是，由于 shapefile 没有任何参考系统，因此该功能不起作用，我不知道如何重新投影 shapefile 以匹配栅格图层。几天来我一直在研究如何处理这个问题，我相信没有参考系统是问题的关键。但是，我无法克服这个问题，我想知道是否有人可以帮助如何处理这种情况。

我正在尝试扩展/放大我从人口普查数据集下载的切罗基民族的 shapefile：https : //www2.census.gov/geo/tiger/TIGER2019/AIANNH/ 使用 rgeos 包。我目前拥有的代码是：

生成的对象丢失了原始 SPDF 的数据框，并且几乎只是一个与原始形状完全不同的圆圈。有什么我想念的吗？

我根据我之前的帖子解决方案使用光栅包擦除功能来裁剪和溶解重叠多边形 -使用 R 中的差异和联合溶解重叠多边形

对于某些多边形，擦除功能出现以下错误：

RGEOSBinTopoFunc 中的错误（spgeom1，spgeom2，byid，id，drop_lower_td，unaryUnion_if_byid_false，：TopologyException：输入几何 1 无效：在点 1.1197332302192855 或附近的自相交点 1.1197332302192855 47.203098020153668 在 1.1197332302319268 09825.

以下代码检查坐标对是否属于特定多边形。我需要使用循环或其他东西检查更多坐标数据。

现在当它调用代码时，它会检查每个值，但如果在国外只有一个坐标，则列表中所有点的状态都是'false'。

我正在尝试为空间数据课程安装软件包，但是我尝试安装到 RI 的每个软件包都会收到相同的错误消息

我在 3.4 的文档中的 R 下创建了一个文件夹。

我也收到此错误消息；

我不明白在安装 R 包之前是否必须创建一个新目录，或者我是否有阻止它安装包的设置？

我之前已经安装了软件包/有 R 2 年了，从来没有遇到过这个问题....

我有一个空间点数据框-> spatial_points

和一个多边形 - > spatial_poly

我可以使用多边形内的所有点进行子集化

但是如果我想对多边形外的所有点进行子集化，我就不能使用

之前有人问过这个问题，答案是gDifference()从rgeos包中使用。美好的。

我的问题是，为什么 [ ] 适用于内部选择，而不是相反？我真的不明白错误信息

h(simpleError(msg, call)) 中的错误：在为函数“[”选择方法时评估参数“i”时出错：一元运算符的参数无效

只是好奇。谢谢。

编辑

这是从带有多边形的子集空间点借来的示例

在这个简单的示例中，可以使用gDifference()，它在此示例中有效。但是，我的 SpatialPointDataframe 非常大，使用 gDifference 会使 R 崩溃。

I would like to create evenly spaced polylines going North to South with 50 mile spacing between each line and 10 miles long. Not sure if this is possible using sf package. In the example below, I would like to have the lines filling the counties across the state of Washington.

I've updated to include an image of what I'd like to do below.

我有我的研究区域“吉尔吉特巴尔蒂斯坦”的 shapefile。我想在地图上绘制地区名称（shapefile）

使用此代码，我可以绘制研究区域地图。如何在这张地图上显示地区名称？

我需要能够在不丢失点的情况下缩小纬度/经度数据的多边形；更重要的是，我需要在正确的方向上有效地“平滑”这些点。通常，gBuffer工作正常，但不能保证点的数量和它们的相对间距。最终，我需要保留每个点的属性，而样条曲线、平滑和其他“良好的效率”gBuffer以及多边形的增长/收缩不允许我以对 1 对 1 映射的足够置信度来保留这些属性.

例子：

我们立即看到点数的变化。我认识到，大多数时候，这是 的理想特性gBuffer，因此也许rgeos不是这种转换的最佳工具。

这个图像对我想要的整体形状有效果，但是增加了点数，这意味着我不能再依赖与原始点的一对一关系。

通常，多边形不是对称的，并且许多多边形都有这样的内切，其中许多在特定方向上“拉”点的方法会出现偏差或方向错误。

我找不到任何选项，gBuffer也找不到其他功能rgeos来保留点的数量和基本空间关系。我不需要“完美”的缩小，如果这改变了事情，但它不应该有很大的偏差。

我发现readWKT封装rgeos可能是解决此问题的最佳方法。但是，我在此过程中遇到了几个问题。

我似乎无法解决的问题是我需要接收多个多面体，因此在将空间数据添加回我的主数据帧时the_geom_sp，我可以在命令中将其作为向量调用。但是，在我有一个预先存在的多面体向量之前SpatialPolygonsDataFrame，我无法制作the_geom_sp成 a 。所以我将它格式化为一个列表 - 但是当我在命令中调用它时，它无法识别列表中的任何多面体对象。运行循环时，我还收到一条错误消息，表明我需要被归类为 a才能无损地接收来自循环的条目。我也受到 的约束，因为它一次只能应用于一个观察。SpatialPolygonsDataFrameSpatialPolygonsDataFrame" In `[<-`(`*tmp*`, n, value = the_geom_sp_entry) : implicit list embedding of S4 objects is deprecatedthe_sp_geomSpatialPolygonsDataFramereadWKT

有谁知道我可以如何解决这个问题？我已经继续并为项目的这一部分切换到 Python，但希望你们能提供任何帮助。