问题标签 [rgdal]

我有一个包含 900K 随机采样坐标点的数据集，每个点都有一个值。我想在地图上制作一个网格，并为每个单元格分配单元格内所有点的平均值。最后，我想将网格绘制为热图。

网格的分辨率可以是 100 平方公里。小区可达 10000 平方公里。

This code works fine on a Macbook, but fails on my Windows desktop with the error:

Any ideas what's causing it?

Full code:

我有一个数据框，其中包含我想与 Google 地球一起使用的各个位置的坐标。这是一个显示结构的简单示例：

我按照这个 SO 帖子和本指南使用包中的函数成功创建KML输出，但是我在调​​整属性时遇到了麻烦。这是代码：writeOGR()rgdal

结果仅适用于查看位置，但我想更改<styleUrl>属性并<name>填充属性。没有它，Google 地球会显示带有[no name]属性的位置：

这是生成的.kml文件：

我需要获取一个<name>元素来填充SimpleField name="city"内容，或者将<name>City</name>标签添加到每个<Placemark>. 我想要的是这样的最终结果（注意添加的<Style>定义、<styleUrl>属性<Placemark>和<name>添加的属性）：

这是结果的样子（类似于我的目标）：

该rgdal文档提到了一个layer_options属性，但没有什么直观地突出给我......

• layer_options = c("<name>????</name>")?
• layer_options = c("<styleUrl>#custom</styleUrl")?
• 还有什么？

上面直接传递标签的尝试似乎不会影响输出。

除了从 中创建默认输出之外，我在谷歌搜索中发现的示例并不多writeOGR()，如上所示。感谢您的任何建议。

rgdal adn rgoes在我们的 redhat linux 机器上编译软件包时出现错误。我试图做一些研究，但找不到可能的解决方案。你能帮我解决这个问题吗，因为这对我来说非常重要。

我正在尝试使用以下代码将全球地图 shapefile 从地理坐标系（纬度/经度，WGS84）转换为 UTM 投影：

不幸的是，我不断收到以下错误代码：

我已阅读有关此主题的以下两个问题，但都与小数据示例有关，而不是与具有大量数据点的 shapefile 有关： rgdal package lat/long -> UTM； rgdal R Package 中 spTransform 中的“检测到非有限变换”。我有一些 UTM 格式的多边形文件希望添加到地图中，我真的需要 UTM 投影。任何关于我可能做错了什么或纠正问题的方法的建议将不胜感激。

我已经在 python 中看到了此代码，但无法弄清楚如何在 R 中执行此操作。我有一些代码，但它对我不起作用。

我的数据有一个纬度列和一个经度列，但它们在国家网格系统中。我不知道如何发布数据集以提供帮助

我正在尝试使用 rgdal。我取 lat 和 long 并制作一个 lat, long 变量。我发现应该有帮助的代码是这个

我不知道我是否正确设置了威士忌。坐标是基本问题。当我运行第三个命令时，出现此错误

(function (classes, fdef, mtable) 中的错误：无法为签名 '"character", "CRS"' 的函数 'proj4string<-' 找到继承的方法</p>

RowID 酒厂纬度经度 1 Aberfeldy 286580 749680 2 Aberlour 326340 842570 3 AnCnoc 352960 839320

我正在处理一个包含所有省份的阿富汗 shapefile。我最终想通过区域命令绘制一些数据，每个命令都包含几个省份（这就是我所说的：http ://en.wikipedia.org/wiki/File:Afghanistan_ISAF_Sept2008.jpg ）。我对在 R 中处理 shapefile 不是很熟悉，但我在这个问题上的教程中没有找到太多的方法。是否可以重新绘制多边形或覆盖 RC 边界并以这种方式绘制等值线？谢谢。

如何在 shapefile 之后保留多边形的信息？让我试着解释一下：

我有一个包含这些数据的 shapefile：

这个数据：

我可以合并data和mapa@data绘制它plot()，但我想使用 ggplot2。但是fortify()的输出没有任何原始变量。例如：

两者（fortify 和 mapa@data）都有一个 id 变量，但值不同。所以，再一次，我的问题是：如何将mapa@data' 的信息传递给fortify()' 的输出（或另一个允许使用 ggplot2 的函数）

我已经阅读了包中使用的 shapefile readShapePoly，maptools但无法使用readOGR. 我希望有人可以帮助我阅读 shapefile readOGR。

orcounty.shp我从这里下载了文件：http: //geography.uoregon.edu/geogr/topics/maps.htm

我还下载了相关文件：orcounty.shx, orcounty.sbx, orcounty.sbn,orcounty.dbf并将所有五个文件放入文件夹中：c:/users/mark w miller/gis_in_R/shapefile_example/

以下代码读取 shapefile 并显示一些属性：

但是，当我尝试使用以下代码读取相同的 shapefile 时，我收到一个错误：

错误消息说：

我可以阅读 Natural Earth http://www.naturalearthdata.com/ shapefiles 使用：

因此，自然地球 shapefile 和俄勒冈 shapefile 之间显然存在差异orcounty.shp。

感谢您对如何阅读的任何orcounty.shp建议readOGR。我的问题类似于这里的问题：rgdal / readOGR - 无法从 .zip 读取 shapefile