问题标签 [maptools]

给定一个 shapefile，我如何塑造和使用数据文件，以便能够使用与 shapefile 中的形状区域对应的标识符来绘制专题图？

到目前为止，我如何采取下一步来生成主题/等值线图，例如根据 Debtors.Petition 值对每个区域进行着色？

（我也刚刚注意到一个可能的问题——大写 GOR 级别不匹配：“Yorkshire and the Humber”和“Yorkshire and The Humber”）

如何在 shapefile 对象后面绘制“光栅”对象？两者都可以自己绘制，但点不会绘制在栅格上：

我正在尝试在另一个已经使用 R 的彩色多边形的 shapefile 之上绘制多边形的 shapefile，并且我正在使用 maptools 和 rgdal。我想将底部 shapefile（具有多个多边形）剪辑到由顶部 shapefile 分隔的区域，顶部 shapefile 具有四个多边形（我认为）。我知道如何填充顶部多边形的内部，当然，通过使用 plot() 中的 col= 命令。但我想要做的是将顶部多边形的外部填充为白色，并通过顶部多边形的内部使底部彩色多边形可见。在我能够做到这一点之后，我将重新绘制底部多边形的边界。我在 Windows 7 机器上使用 R 2.13.2 和最新版本的 maptools 和 rgdal 以及它们的所有依赖项。

我有依赖于这些 shapefile 的代码：ftp: //ftpext.usgs.gov/pub/er/wi/la.crosse/McKann/bcr%20arcinfo%20files/

我可以使范围图的内部成为一种颜色

但我想做的是让外面变成白色（这样我以后可以在上面画出州和省的边界），并保持里面透明，这样我就可以看到范围内的州和省的颜色。

如果您需要更多信息，请告诉我，并提前感谢您的帮助。也许有更简单的方法...

如何使用 R maptools 库将 Dataframe 转换为 SpatialGridDataFrame？我是 Rpy2 的新手，所以这可能是一个非常基本的问题。

R代码是：坐标（dataf）=~X+Y

在 Python 中：

除了这个特殊的问题，我很高兴能得到一些关于更一般性想法的反馈：我的最终目标是使用 gstat 库对空间数据进行回归克里金法。R 脚本工作正常，但我想从 Python/Arcgis 调用我的脚本。您如何看待这项任务，这可以通过 rpy2 实现吗？

非常感谢！理查德

我想使用 R 生成一个非常基本的世界地图，其中包含一组特定的国家，这些国家用红色填充，表示它们是疟疾流行国家。

我在数据框中列出了这些国家/地区，但我很难将它们覆盖在世界地图上。

I have tried using the wrld_simpl object and also the joinCountryData2Map method in the rworldmap package.

I would comment on this answer to prevent addition of a possibly redundant question but I do not have enough reputation at the moment, apologies for this.

https://stackoverflow.com/a/9102797/1470099

I am having difficulty understanding the arguments given to the plot() command - I wondered if there was just an easy way to tell R to plot all of the country NAMEs in my list on the wrld_simpl map instead of using grepl() etc. etc.

我正在使用 maptools 在 R 中打开 shapefile，但在 Rstudio 中出现以下错误：

当我在 RGui（而不是 Rstudio）中尝试完全相同的代码时，它可以工作！我想继续使用 Rstudio，因为它非常方便，所以有人知道如何解决这个问题吗？

谢谢！

这给了我一张世界地图，并显示了地图顶部和底部的大量空白。我希望这个空间消失，因为我需要打印许多（~60）张地图。我正在使用knitr将嵌入地图的报告。例如

但我不认为这是一个knitr问题。所以：

1. 我怎样才能摆脱空白？
2. 如何确保地图从左到右填满我的页面？

我试过了，，，，但mai没有这样的运气！我猜测chunk 选项会给出我的问题 2 的结果，但由于我严重卡在问题 1 上，我发现很难说。marparout.width

非常感谢！这里有很多东西要学！

我想根据 @data 数据框中的相应属性值从 SpatialPolygonsDataFrame 对象中简单地删除一些多边形，以便我可以绘制简化/子集的 shapefile。到目前为止，我还没有找到一种方法来做到这一点。

例如，假设我想从这个世界 shapefile中删除所有面积小于 30000 的多边形。我该怎么做呢？

或者，类似地，我怎样才能删除南极洲？

如果我做这样的事情，情节不会反映任何变化。

如果我这样做，结果相同：

任何帮助表示赞赏！

我尝试将具有多行的空间对象（从 OSM 检索到的河流）转换为 KML。对于具有单线的对象，使用 kmlLine 很容易。但是，对于多行，以下方法不起作用，我尝试从文档中调整示例是徒劳的：

我正在使用鸟带（欧洲人响铃）记录来研究日出方位角和迁徙航向之间的关系。我有来自欧洲和北美的几千只鸟在 1935 年至 1990 年期间的条带（振铃）位置。我的初步分析表明，日出方位角、中午太阳高度和鸟类在秋季迁徙的方向之间存在关系.

这种关系随日期而变化，因此我最初使用美国海军天文台计算器 ( http://aa.usno.navy.mil/data/docs/AltAz.php ) 的计划变得非常耗时。我查看了包maptools和moonsunin R，但无法弄清楚如何使用矩阵而不是单个日期和位置运行它们各自的命令。

另外，我不需要日出的时间，只需要方位角。这些值也不必精确到一定程度。如果有人有任何建议，我很乐意听取您的意见。