问题标签 [maptools]

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

形状文件包含在这里

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

数据样本structure如下。

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

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

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

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

任何建议，将不胜感激。

我无法从大型形状文件（916.2 Mb，4618197 个元素 - 从这里提取点数据：https ://earthdata.nasa.gov/data/near-real-time-data/firms/active-fire-data ) 在 R 中。我正在使用readShapeSpatialinmaptools来读取形状文件，这需要一段时间但最终可以工作：

然后我有一个data.frame我想要提取数据的坐标。然而，R 真的在为此苦苦挣扎，要么失去连接，要么死机，即使只有一组坐标！

任何人都可以建议我以更有效的方式执行此操作吗？

还是我需要在更强大的东西上运行这个脚本（目前使用的是 2.3 GHz 处理器的 mac mini）？

非常感谢！

我正在检查qtm helpfiles，但它没有说明如何更改图例百分比。

此行qtm(doe, fill = "sep12_PERC", fill.title = "Installs")创建此地图，但我们希望 Installs 为

如何解决这个问题？

谢谢！

这可能是一个显而易见的问题。

我有一个澳大利亚 Shapefile （可以在这里找到http://www.abs.gov.au/AUSSTATS/abs@.nsf/DetailsPage/1259.0.30.001July%202011?OpenDocument）和一些包含纬度和经度的自定义数据一些人。

如何使用纬度和经度将形状文件连接到自定义数据，以便我可以创建人口分布图并相应地为地图的各个部分着色？

谢谢你。

我设法创建了一张地图，甚至包括了一个指北针，但无法map.scale正常工作并出现这种错误：

map.scale 中的错误（x = -83，y = 12，ratio = FALSE，relwidth = 0.2，cex = 0.6）：未使用的参数（ratio = FALSE，relwidth = 0.2，cex = 0.6）

这是代码：

我计算了一个变量“final_factor”，我确信它是正确的。当我尝试在地图上绘制这些数据以显示不同的值时，我知道加利福尼亚拥有最高值之一，但它显示为最低值之一。我认为这与阿拉斯加和夏威夷没有出现在地图上有关吗？但我在向量中有两个州的名称：

打印状态向量：

编辑

我正在尝试使用加载 Shapefile，readShapeSpatial()但出现错误

“readShapeSpatial(shp) 中的错误：无法读取文件类型”。

关于这里出了什么问题的任何想法？这是我正在使用的 shapefile 的代码，它来自 2010 年巴西人口普查。

奇怪的是我可以使用readOGRfrom package加载文件rgdal

我目前正在尝试使用 R 在空间地图上绘制沿海的弯曲贸易路线。结果应如下所示： http ://spatialanalysis.co.uk/2012/03/mapped-british-shipping-1750-1800 / 或者这个

有人可以向我推荐一个可以做到这一点的包裹吗？主要问题是路线及其“曲线”。

我对 R 不是很熟练（刚开始使用 ggplot2），但我需要知道它是否可能，它有多难，以及它是否做得更好，以其他方式学习什么才能进行这种类型的可视化。

提前谢谢帕维尔

我有一个列表（从 Stata dta 文件导入），其中包含为地图定义空间多边形的条目，以提供一个想法：

我想把它转换成一个数据结构，我可以用它在 R 中创建一个地图（没有更大的难度，也不需要广泛的专业知识，以前从未做过）。我推断 sp 包的 SpacialPolygons 类型是最简单的选择。此外，从定义来看，似乎 SpacialPolygons 方法（在此包中定义，请参阅包的文档，第 79 页）似乎是从列表转换为此数据类型的正确方法。

不幸的是，该方法并不那么简单，我需要一些帮助。我的（天真）尝试产生了一个我不理解的错误，并且在谷歌搜索中没有出现任何有趣的结果：

上面的列表可以转换为 SpacialPolygon 吗？如果是这样，怎么做？如果不是，我应该选择什么格式？谢谢。

我正在分析一个复杂数据集中的时间模式，该数据集中包含几个环境变量以及来自各种动物物种的活动数据。这些数据由多个实验装置收集，每个装置的数据每分钟存储一次。该项目已经运行了几年，所以我的数据集相当大。

我的一个数据集的前几行如下所示：

因为我想将这些变量与整个季节的日出和日落等各种自然周期联系起来，所以我使用了这个包maptools来计算日出和日落时间

由于大多数动物的行为不同，这取决于它是白天还是晚上，我使用日落/日出时间来计算数据来计算一个新变量，它在晚上为 0，在白天为 1：

到目前为止，一切都很好......甚至可以maptools使用民用/航海/天文黄昏和黎明的开始而不是日出和日落。

然而，这就是我的问题开始的地方。我想对实验中的所有日子进行编号。而不是像往常一样在午夜增加日计数器，而且很容易做，我想在日落时增加日计数器（或者可能在未来的实验中另一个可移动的时间，如日出，航海黄昏和黎明，......） . 由于日落并非每天都在同一时间发生，所以对我来说，这不是一个简单的问题要解决。

我只提出了一个for-loop，这不是一种很好的做事方式。此外，考虑到我在几个设置中每分钟收集一次价值超过 6 年的数据点，我可以坐下来观察构造板块的移动，而 R 会通过这样的一大堆循环：

除了丑陋和缓慢之外，这段代码还有一个大问题：它不处理缺失值。有时，由于设备故障，数小时或数天都没有记录数据。如果在日落期间没有记录数据，则上述代码不会增加日计数器。这意味着我还需要 - 以某种方式 - 合并日期/时间代码。很容易创建自实验开始以来的天数变量：

也许可以使用这些数字，可能与Herka 的nice rle-algorithm 结合使用。

我曾经dput从一个设置中获取几个月的数据，包括几大块缺失的数据，以及新创建的变量（如本文和Herka 的回答中所述）可在此处获得。

我一直在寻找更好，更好，尤其是更快的东西，但一直无法想出一个好的技巧。我摆弄了我的数据框的子集，但得出的结论是这可能是一种愚蠢的方法。我看过maptools,lubridate和GeoLight. 我搜索了 Google、Stack Overflow 和各种书籍，比如 Hadley Wickham 的出色的Advanced R。一切都无济于事。也许我错过了一些明显的东西。我希望这里有人可以帮助我。