问题标签 [maptools]

我发现了两种在 R 中打开 shapefile 的基本方法——使用rgdal和maptools：

两种情况下的数据结构似乎完全相同（SpatialLinesDataFrame 类，sp 包）。但是，虽然rgdal正确读取投影，maptools但没有（您可能必须手动分配 CRS）：

那么，我为什么要使用maptools打开形状文件呢？我可能只会犯手动分配 CRS 的错误！

我可以得出结论，这两种方法是等效的，但是使用rgdal打开 shapefile 总是更安全的方法吗？

可以writeSpatialShape { maptools }以某种方式强制使用 CRS 编写 shapefile 吗？

似乎没有这个选项。

PS：我知道可以使用writeOGR { rgdal }.

伙计们，我正在 Mac 中使用 R 处理一些数据，首先让我向您展示我的代码。

然后，作为函数的结果plot，我看到了这个：

现在，让我解释一下我的问题......这只是一个测试，它只有2分，但我面临的任务有1000-1200分左右，我也必须做十次左右，多亏了我无法使用我的计算机进行此模拟，因为它会花费大量时间，而不是我计划租用以下服务之一：Amazon EC2（Amazon Elastic Compute Cloud）。

在那里，我必须通过 Linux/UNIX 外壳（所以：Ubuntu，我猜。），现在我的问题与我的代码中的最后两行有关，我如何导出Kenv到文件或其他东西以便导入它，然后plot它在我的电脑。

我希望你能理解我，如果你有任何问题，请告诉我。

我刚刚使用名为indiemapper的映射工具创建了世界地图的墨卡托投影，并使用该工具将该地图导出为 SVG，现在导出的 SVG 的路径数据是这样的负点，

如果我使用上述路径数据来创建地图，那么地图就会消失，有什么理由将负点转换为正点并使该地图适合视图？

我对在 R 中使用 Shapefile 比较陌生。我想为印度创建国家和地区级别数据的可视化，我在 R 中使用 maptools。对于 shapefile，我从gadm.org获得了数据，其中包括 state/印度的区/街道级数据。

但是，我试图在 R 中加载文件，但它没有发生。这是我的代码：

这里 IND_adm2.shp 是印度的区级形状文件。

R 无法读取文件，给我这个错误：

我不太确定问题出在哪里。我有 R 3.0.1，将文件放在适当的目录中，但有些地方不对劲。我也是在 R 中制作地图的新手，所以有点烦人。有想法该怎么解决这个吗？

假设我正在使用地图工具在世界地图上绘制国家，如果我要绘制一个国家，有没有办法以不同的颜色绘制与这个国家接壤的国家？我正在使用wrld_simplmaptools 附带的 shapefile，所以说我绘制中国：

有没有办法让我把所有与中国接壤的国家都画出来。我希望能够为许多不同的国家做到这一点，所以理想情况下我想要一个通用的解决方案，而不仅仅是针对中国的解决方案。

我有大伦敦地区的形状文件。我使用包中的readShapePoly函数将maptools其作为SpatialPolygonDataFrame.

我想绘制这些多边形。我已经通过使用plotR 中的基本函数完成了。输出如下图所示：

现在，我正在尝试使用绘制相同的形状文件，ggplot2但它对我不起作用。我在图中得到了一些奇怪的线条，如下所示：

我使用的代码是：

其中 projString 是描述用于输入形状文件的投影的投影字符串。

与我之前的问题类似，假设我正在使用地图工具在世界地图上绘制国家，如果我要绘制一个国家，有没有办法找到这个国家的中心点并从该点绘制径向距离？我正在使用wrld_simpl附带的 shapefile maptools，所以如果我绘制德国：

我想找到德国的中心并从这一点绘制一个圆圈，以显示任何位于 100 公里距离内的区域。同样，我希望能够为许多不同的国家/地区执行此操作，因此理想情况下，我想要一个通用的解决方案，而不是仅针对德国的解决方案。

我有意大利各省相对生产力的数据。我根据生产力水平建立了5个组。最后，我想使用意大利各省的地图来表示这种情况，因为我不能使用颜色，因为地图的打印将是灰色的。我使用了以下内容：

其中，为简单起见，五个组是 a,b,c,d,e

结果在这里http://s9.postimg.org/m5co5g3nz/Rplot.png

我想做两件我无法弄清楚的事情。

首先我需要添加图例，我尝试过

我无法找到一种方法来插入组 (c) 的代码来代替 ??。我试过了

但他们不工作

其次，为了让事情看起来更清晰，我可能想用与 (c) 组相同的样式填充 (e) 组省份，但线条走向相反方向甚至垂直。即使在这种情况下，我也无法处理它。

谢谢

我正在从 Arcmap 导入一个多边形 shapefile，该文件已经有一个投影集并附带了它的所有文件（sbn、sbx、prj 等）。

但是，在我使用该readShapePoly功能后，当我进行摘要时，投影信息看起来像是空白的。投影是否已包含或未被识别？

我知道有一个proj4string属性，但不清楚如何使用它以及它已经附加到 .shp 文件的 prj 文件。如果有另一个功能，我可以考虑另一个做得更好的功能。不确定rgdalwithreadOGR是否符合我的要求。

编辑跟进：我尝试readOGR感谢您的回复。我正在使用此代码 test<-readOGR(dsn=getwd(), layer="grid") 并且 shapefile 在这里 speedy.sh/Ry8rU/grid.zip 并且它仍然不读取投影。