问题标签 [rgdal]

我想知道如何根据特定颜色的矩阵编写 GTiff？这可能是当颜色不代表定量变量，而是定性空间现象（例如人口统计集群）时。例如，如果我有一个 1000x1000 的颜色值矩阵（字符串，如 '#FF0000'），我可以轻松地绘制这些值，grid.raster(mtrx, interpolate=FALSE)但这很难与其他 shapefile 重叠。我还可以将它们转换为 SpatialPixelsDataFrame，字段名称为“x”、“y”、“col”。但我很难过如何将这些对象中的任何一个写入 GTiff。例如：

产量

非常感谢任何指点。

例子：

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

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

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

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

我有以下数据集，其坐标自 1983 年起在州平面坐标系中。所有坐标都位于长岛区 (3104)。

我需要有纬度/经度的坐标。格式或 wgs84。

有人能告诉我我是怎么做到的吗？谢谢你的帮助。

我正在尝试做与此问题中提出的相同的事情，即 Cartogram + choropleth map in R，但从 SpatialPolygonsDataFrame 开始并希望最终得到相同类型的对象。

我可以将对象保存为 shapefile，使用scapetoad，重新打开它并转换回来，但我宁愿将它全部包含在 R 中，以便该过程完全可重现，以便我可以自动编码数十种变体。

我已经在 github 上分叉了 Rcartogram 代码，并在此处添加了我的努力。

本质上，这个演示所做的是在地图上创建一个 SpatialGrid，查找网格每个点的人口密度，并将其转换为cartogram()工作所需格式的密度矩阵。到目前为止，一切都很好。

但是，如何根据 的输出对原始地图点进行插值cartogram()？

这里有两个问题。首先是将地图和网格放入相同的单位以允许插值。第二个是访问每个多边形的每个点，对其进行插值，并使它们保持正确的顺序。

网格采用网格单位，地图采用投影单位（在示例 longlat 的情况下）。必须将网格投影到 longlat，或者将地图投影到网格单元中。我的想法是制作一个假的 CRS 并将其与 in 中的spTransform()函数一起使用package(rgdal)，因为它可以轻松处理对象中的每个点。

访问每个点都很困难，因为它们在 SpPDF 对象中是几层：我认为对象>多边形>多边形>线>坐标。任何想法如何在保持整体地图结构完整的同时访问这些？

我找到了一个相当简单的例子来说明如何做到这一点，但我无法让它为我工作。我对 R 很陌生

但这是示例数字。我有数千个坐标必须转换，但我不知道如何将它们从我的表中获取到这个脚本中

我的数据集有 3 列，ID、X 和 Y。如何使用这个等式转换它们？我已经坚持了好几个星期了

我有一个带有纬度和经度的数据集，我想使用 EPSG 2790（http://spatialreference.org/ref/epsg/2790/）或 ESRI 102672（http://spatialreference.org/ref/esri/102672/）。

以前肯定有人问过这个问题；我的代码基于此处的答案（在 rgdal R 包中的 spTransform 中的“检测到非有限变换”和http://r-sig-geo.2731867.n2.nabble.com/Converting-State-Plane-Coordinates-td5457204。 html ).

但由于某种原因，我无法让它工作：

给出：

在生成一些空间分析时，我一直在使用 IPython Notebook 将我的 Python 和 R 代码集成到一个工作流中。我以前能够将必要的基础设施导入两个命名空间来执行此任务，但我认为由于后续更新而发生了一些变化。特别是，我不再能够在%%R单元格中调用rgdal库：

这很奇怪的原因是因为我可以在基础 R 和 RStudio 中调用库（输出相同）：

RStudio 和 Notebook 中的版本信息是一样的……：

...就像 R 的位置一样：

请注意，我试图从笔记本中安装rgdal ...：

...并且它与已在 RStudio 中成功加载的版本相同：

有了这个，我对这里发生的事情不知所措。虽然我可以在 RStudio 中继续，但它会分散我尝试创建的产品。任何关于此事的建议将不胜感激。

附录：看来这可能不仅仅是 %rmagic 的问题。在尝试在 Python 中执行不同的操作时，我在尝试同时导入osgeo和gdal时收到以下错误：

我真的只是再次安装了GDAL，所以我有最新的版本。想法？

我使用 readOGR 成功导入了一个多边形 shapefile，但是 .prj 文件中的信息没有被导入。

测试类：SpatialPolygonsDataFrame nfeatures：19407 范围：35551.4、1585917、6318047、9408727（xmin、xmax、ymin、ymax）坐标。参考。: 不适用

摘要（测试）类 SpatialPolygonsDataFrame 的对象坐标：最小最大 x 35551.4 1585917 y 6318047.3 9408727 预计：NA proj4string：[NA]

投影应该是 utm21n 复杂区域

我需要使用 R 包在谷歌地图中显示我的形状文件（.shp）。我正在使用 maptools 包来读取形状文件。使用函数 ==> readShapePoly() 然后我将使用 plot() 来显示。现在它在 R 静态窗口中绘制 shapefile 内容，我无法放大/缩小，而且我需要在谷歌中显示动态地图。是否有任何包可用于在谷歌地图上绘制我的形状文件的多边形？

我想在 R 中的光栅文件上绘制一个 shapefile，但我不能让它们完美重叠：光栅似乎逆时针旋转了几度。投影有问题吗？

请考虑以下 MWE