编辑经过更多研究但仍然没有解决方案，我正在添加大量编辑以及指向 .shp 文件的链接。

形状文件包含在这里

我有一个包含 9 个多边形的 SpatialPolygonsDataFrame，每个多边形还包含多个嵌套多边形 - “洞”。数据摘要在这里。

数据样本structure如下。

如下例所示（九个之一），父多边形有多个孔。

这是我第一次涉足sp()、rgdal()、rgeos()和其他空间包，我发现了许多关于使用运算符提取区域等的有用帖子，但仍然存在问题。虽然这篇文章提供了接近的解决方案，但我似乎无法调整代码以满足此处描述的需求。

我想获得一个 SpatialPolygonsDataFrame，它只包含来自每组子列表（即 data@polygon）的父（最大）多边形。似乎我应该能够仅提取父多边形或“溶解”这些孔。

最终结果将是 9 个多边形，每个都是 9 个列表的父级，我可以将其导出为 ESRI shapefile。

任何建议，将不胜感激。

I'm sorry this is going to be a for loop 101 question. I am struggling to write a simple for loop to generate a table of distances between cities based upon longitude-latitude data

locations returns the following table:

My goal in this particular part of the task is to produce a table of the distances (in kilometres) between each of the cities in the nature of a correlation matrix, with the diagonal being 0 (ie all cities are zero distance from themselves).

To do this I am using the sp package, which requires a matrix of long-lat values, so I can remove the text as follows:

The tool spDistsN1 allows me to get this information by comparing the distance all cities in the matrix are from one individual city. Clearly, I can use the following expression to obtain the distances of all cities from Sheffield (city or row# 1):

This correctly gives:

However, to achieve my desired correlation matrix style output, I want to achieve this for each of the cities, so I tried to write a for loop:

This gives me the correct values for NY:

So I presume I have overwritten one city by another throughout the loop. I appreciate the help in showing me where I am going wrong. Many thanks.

我想使用spplot+ sp.lines（格）而不是plot+ segments。你知道实现这一点的简单方法吗，例如R：计算两个点层之间的最短距离

形象化这个

可能可以转换为类似的东西

但是有没有一种简单的方法可以转换segments为SpatialLinesielist("sp.lines", Lines(...))

我正在尝试为空间多边形数据框中的多边形获取伪重心。今天我偶然发现了coordinates实际上为SpatialPolygonsDataFrame返回某些内容的函数。

不幸的是，我在SpatialPolygonsDataFramecoordinates的值的帮助中一无所获。有人能告诉我这些坐标是什么吗？

我在 SpatialPolygonsDataFrame 中有一个多边形列表，需要将其中一个设置为另一个孔。

我发现set_Polypath如何在新创建的多边形上定义一个洞，但如何在现有多边形上设置“洞”标志？

编辑：

正如 Edzer Pebesma 在评论中所建议的那样，向多边形添加孔的推荐方法不是修改槽而是重建多边形，如相关问题所示。

原始问题

在SpatialPolygons-class我尝试将多边形修改为另一个多边形的孔的帮助下，但是当“孔多边形”显示为另一个多边形的边框时，其内部与其余部分一样着色。

我究竟做错了什么？

在这里使用defPunched和defHole定义：

lm在 stats 中有很好的功能可以对对象进行残差和回归删除诊断glm，例如：

stats::plot.lm和stats::influence.measures

例子：

我想对spdep::spautolm （类型“SAR”）对象进行这些（或类似的）诊断，但既不plot.lm也不influence.measures接受 spdep::spautolm/sarlm对象。

你能推荐一个解决方法吗？

还是这些诊断根本不适合 SAR 类型回归？

似乎这个问题已经以不同的形式被问过几次，但我找不到正确的解决方案。我有一个带有多个多边形的 SpatialPoint 对象，并且想使用插槽“ID”对一个多边形进行子集化和绘制。

使用这个问题的例子：

我可以提取 SpatialPolygons 对象的 ID

但是我怎样才能使用这些信息来子集和绘制一个多边形呢？很高兴得到任何帮助，在此先感谢！

我正在尝试使用 GADM 数据并使用 R 绘制印度各州/地区的详细地图。我使用了以下代码

生成这张地图

我现在正在尝试输入几个标记（或点加文本）来显示地区总部，“Purulia”，位于 lon = 86.36521 lat = 23.33208 但不知何故我无法sp.layout正确使用语法。对于一组已知其长、纬度的城镇，我将需要它。如果有人可以帮助我解决这个问题，将不胜感激。

嗨，我是一个真正的自学 R 新手：我正在关注 GWR 教程， 但已经在第一页我的代码失败了

创建一个调色板并在 spplot

错误：“colors = c(�”中的意外输入

和 spplot 一样

错误：“spplot（Nocars，�”中的意外输入

你能帮忙吗，我确定它很傻。