问题标签 [rgdal]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
r - 在 shapefile 后面绘制栅格
如何在 shapefile 对象后面绘制“光栅”对象？两者都可以自己绘制，但点不会绘制在栅格上：
r - R maptools 或 rgdal 包——如何填充多边形的外部？
我正在尝试在另一个已经使用 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 - 在 rgdal R 中使用 spTransform 重新投影空间点时出错
G'day，
我有大量的经纬度坐标，它们位于 CRS 澳大利亚大地基准 66/84（简称 AGD66）中。我想将这些坐标从 AGD66 更改为 WGS84，因为它们之间大约有 200m 的差异，而且我在 WGS84 中有其他坐标和图层。我试图通过以下方式做到这一点：
有谁知道我为什么会收到此错误或对如何将这些坐标从 AGD66 更改为 WGS84 有任何建议？提前感谢您的帮助。
干杯，
亚当
r - 如何使用 rgdal 包将 NAD 83 坐标转换为纬度和经度？
我有坐标，所有坐标都应该位于 DC，但我不知道如何将它们从 NAD 83 转换为 R 中的纬度和经度。我正在使用spTransform()rgdal 包中的函数并收到关于不合格的错误数据。
proj4strings 的其他组合产生相同的错误。我相信错误是因为坐标太大，但我不确定为什么会这样。坐标的文档如下：
值在马里兰州平面米 NAD 83 地图投影中。
我对映射和投影非常陌生，感谢您提供任何帮助。
r - rgdal / readOGR - 无法从 .zip 读取 shapefile
我正在尝试加载一个 zip 级别的 shapefile 来进行一些绘图，每个： https ://github.com/hadley/ggplot2/wiki/plotting-polygon-shapefiles http://www.nceas.ucsb.edu/scicomp/用例/ReadWriteESRIShapeFiles 等
我的代码：
我收到此错误：
我检查了驱动程序，但老实说我无法解释输出：
file.info 给出：
鉴于我尝试过的互联网搜索，看起来我不是唯一一个遇到此问题的人，但我一直无法找到答案。我不确定问题是否与 .zip 文件夹中的 shapefile 有关。由于它是一台工作计算机，我必须等待 IT 安装 WinZip，以便我可以提取 shapefile 并在 .zip 文件夹之外尝试。希望那是明天。
此外，软件包帮助声明“请注意，数据源目录中的杂散文件（例如 *.dbf）可能会导致伴随 *.shp 丢失的异常错误。” .zip 文件夹具有以下内容：a.dbf、b.prj、c.shp、d.shp.xml、e.shx。
你能提供的任何帮助都会很棒！-亚历克斯
r - 为什么gdal给了我来自R的变量的不同最大值和最小值？
这三行将打开并从 netcdf 文件中读取一个变量：
当我查看提取的数据时，我发现我有非常小的负值。然后我使用 gdalinfo 发现最大值为 400，最小值为 0（这是正常的，这就是值的样子）。任何想法为什么我得到错误的号码？
r - R spCbind 错误
我之前已经成功地将信息添加到 shapefile（请参阅我在http: //rusergroup.swansea.ac.uk/Healthmap.ashx?HL=map 上的帖子）。
但是，我只是尝试使用稍微不同的 shapefile（威尔士的新当地健康委员会）再次执行此操作，并且代码失败spCbind并出现“行名不相同错误”
我的行之前确实有不同的名称，这没有造成任何问题。我重新标记了 wds.xtra 中的列以匹配“NEW_LABEL”，但这没有帮助。
wales.lonlat 和 wds.xtra 之间的标签和标签顺序完全匹配。
（我使用的是基于 R 2.13.2 构建的 Revolution R 5.0）
r - shapefile（R）中的只读数据槽？
我有一些非常大的 shapefile。SpatialPolygonsDataFrame我可以使用rgdal函数将它们读入readOGR，但是每个文件都需要很长时间。我实际上只对插槽data.frame中显示的内容感兴趣。@data有没有办法只读取数据，跳过资源密集型多边形？
示例代码：
r - 在 SpatialPixelsDataFrame 上覆盖 shapefile 时产生的 NA
将一系列流域多边形（来自 shapefile）叠加到 SpatialPixelsDataFrame 对象时，我遇到了一个奇怪的问题。有 70 个集水区，但我得到NA了 60 个集水区的 s。
数据可在 Dropbox 上获得，以下脚本将下载、解压缩等文件并创建可重现的示例。
问题是NAs：
不应该有任何NAs 因为网格覆盖了整个英国地区并且集水区很小，所以我可能做错了什么！如果是这样，我做错了什么，我该如何解决？我想最终得到 70 个流域中每个流域的Deposition（数据）平均值。SpatialPixelsDataFrameshp
[注意，网格是虚拟数据；我无法分发实际的沉积网格，但数据和对象本质上（完全）与depo.sp.
附加信息
r - 以下错误是什么意思：TopologyException：发现LINESTRING之间的非非非交集
我正在尝试强化使用 rgdal 加载到 R 中的形状文件，但出现以下错误：
我正在使用来自 maplibrary.org 的非洲大陆的形状文件。它可以从我的保管箱中获得： https ://www.dropbox.com/s/etqdw3nky52czv4/Africa%20map.zip
我正在使用以下代码：
我得到了我之前提到的错误。我认为 R 的多边形存在一些问题 - 有没有办法解决这个问题？