问题标签 [rgdal]

我下载了一个数字高程数据集。该文件为“.tif”格式。我已经加载了 R。但是当我尝试将坐标参考系统从经度/纬度转换为本地 UTM 时，就会出现问题。代码和错误如下所示。

我怎么解决这个问题？

我是 R 中的 GIS 新手，我正在尝试将 lat/long 或 UTM 坐标添加到 shapefile。我从这里下载了芝加哥市的边界 (City_Boundaries.shp)：http ://www.cityofchicago.org/city/en/depts/doit/supp_info/gis_data.html

我加载了 maptools 和 rgeos 库：

我将数据带入 R 并尝试为 16T 区添加 UTM 代码：

但是，结果数据对我来说没有意义。当我查看“city”中的“coords”槽时，坐标值是，例如，X=1092925 和 Y=1944820。我使用了一个外部工具（http://home.hiwaay.net/~taylorc/toolbox/geography/geoutm.html）来查找 GPS 坐标。结果是：lat=17.511，long=-81.42。这是牙买加和洪都拉斯海岸之间的某个地方。在我看来，shapefile 坐标存在于它们自己的宇宙中；正如 shapefile 的名称所暗示的那样，坐标提供了对城市形状的准确描述，但这些坐标不会自动映射到地球上。

我的最终目标是确定是否在芝加哥发生了一些地理标记事件。我很乐意将按纬度/经度列出的这些点转换为 UTM，如下所示：

如果这更有意义，我也愿意尝试以原始格式处理事件数据。最终，我需要帮助将芝加哥城市边界变成一个多边形，我可以对照我的事件数据进行检查。像这样的东西：

如何将我的 shapefile 和事件数据放入一个通用（且准确）的参考框架中，以便我可以检查点是否在城市中？您可以提供的任何帮助将不胜感激。

我希望能够从 R 中的等高线创建一个高程图。我对使用形状文件非常陌生

目前我已经从这里下载了数据， 这些数据提供了适用于整个英国的 .shp 文件。

它还提供了等高线，总结了英国的拓扑结构。

对于高程图，我想要一个data.frame或data.table多个均匀间隔的点（彼此相距 100m）来产生一个数据输出，给出一个 x、y 和 z 值。其中 x 和 y 代表纬度和经度（或东向和北向），z 代表高度（以米为单位）。

我认为可能有一些工具会自动为您执行插值，但不确定它如何处理地理空间数据。

这是我的基本开始...

我有一个我下载的光栅，无论出于何种原因，在读入它时都没有“proj4string”。

文档说它是地理参考的，我想知道我是否可以创建一个自定义的“proj4string”来定义数据。我也可以访问 ARCGIS，但是我对它不是很熟悉，也不知道如何重置投影的中心。

提供的地理参考信息是：

我很难使用我的 R 代码。我正在尝试根据我已经拥有的数据框创建一个新的数据框，其中每个重复值单独乘以 1000 并按顺序加 1。例如，我的数据框中的值范围从 3869014 到 4524673，每个数字有多个值（最多 100 个）。例如：[3869014，3869014，3869014，3869014，3869014，3869014，3869014，3869014，3869014，3869015，3869015，3869015，3869015 What I want is: [3869014001, 3869014002, 3869014003, 3869014004, 3869014005, 3869014006, 3869014007, 3869014008, 3869015001, 3869015002, 3869015003, 3869015004, 3869016001, 3869016002, 3869016003, 3869016004, etc...]

我尝试了以下代码，但它将每个数字乘以 1000 并添加一个而不考虑重复项。它也只添加一个，而不是添加一个计数（例如：1、2、3、4 等......）。所以输出是 [3869014001, 3869014001, 3869014001, 3869014001, 等等......这不是我想要的。我对在 R 数据帧中循环有点陌生。谢谢您的帮助。

我已经安装了 geos、proj-4.8.0 和 gdal-1.9.2，但是，在安装 R 包 rgeos 之后，rgdal 无法安装。

这是我的主目录：

. . . 和我运行的错误install.packages("rgdal", type="source")：

** 构建包索引 ** 安装小插图 ** 测试是否可以加载已安装的包 dyn.load(file, DLLpath = DLLpath, ...) 中的错误：无法加载共享对象 '/mainstorage/michaeld/R/x86_64- redhat-linux-gnu-library/3.1/rgdal/libs/rgdal.so'：libgdal.so.1：无法打开共享对象文件：没有这样的文件或目录错误：加载失败执行停止错误：加载失败*删除'/ mainstorage/michaeld/R/x86_64-redhat-linux-gnu-library/3.1/rgdal'</p>

下载的源包位于'/tmp/Rtmpt51iFG/downloaded_pa​​ckages' 警告消息：在 install.packages("rgdal", type = "source") 中：安装包 'rgdal' 的退出状态非零

我推断这个错误与依赖库有关。. . 但这就是我收集的全部内容。

我应该提到，按照这里sudo yum install libgdal1 libgdal1-dev libgeos libgeos-dev的建议运行会产生以下结果：

没有可用的软件包 libgdal1。没有可用的软件包 libgdal1-dev。软件包 geos-3.4.2-1.3.amzn1.x86_64 已安装且最新版本 没有可用的软件包 libgeos-dev。

所以，也许更好的问题是我如何安装这些dev包？

感谢您的任何指导。

这可能是一个简单的问题，但是，我想查看地图上的每个输出区域，从这里下载了寄宿生。选择 2011 OA 并获得完整的分辨率，领域的范围。

我已经对其进行了压缩...并使用以下命令读取了 shp 文件

如何绘制或查看每个输出区域的坐标？

我想使用提到的两个包“rgdal”和“raster”来裁剪 GeoTiff 光栅文件。一切正常，除了生成的输出 tif 的质量非常差并且是灰度而不是彩色。原始数据是来自瑞士联邦地形办公室的高质量栅格地图，示例文件可在此处下载。

这是我的代码：

为了重现此示例，请下载示例数据并将其解压缩到文件夹“c:/files/”。奇怪的是，使用样本数据，裁剪图像的质量还可以，但仍然是灰度的。

我使用“数据类型”，“格式”选项进行了尝试，但没有得到任何结果。有人可以指出解决方案吗？我应该提供输入数据的更多信息吗？

编辑： Josh 的示例与示例数据2配合得非常好。不幸的是，我拥有的数据似乎更旧并且有些不同。如果您阅读以下 GDALinfo，您能告诉我我选择了什么选项：

I'm trying to read in a shapefile using the rgdal library, and am having no luck.

When I try to import using the following syntax:

I'm met with the error

How can I go about diagnosing and fixing the problem? Many thanks.

与R: Crop GeoTiff Raster using packages "rgdal" and "raster"中提出的问题类似，我正在尝试使用包“rgdal”和“raster”从瑞士联邦地形局裁剪地图，同时保留原始颜色表. 对于单个带状 *.tif 文件，裁剪后的图像会丢失颜色表信息，因此无法正确显示（生成的图像几乎是黑色的）。

输入文件可在此处下载，应提取到文件夹“C:/files”中。这是代码：

前面提到的帖子中提出的解决方案仅适用于 3 波段 *.tif，但不适用于 1 波段 *.tif（例如示例文件）。

一个可行的解决方案是将包含颜色表的单带 rasterLayer 转换为 3 带 RGB rasterStack（如前面提到的帖子中的评论中所述），这显然可以保存颜色表。

但是，我不知道如何在保留颜色表的同时将单波段 *.tif 转换为 3 波段 RGB rasterStack。有谁知道如何进行这种转换，或者有没有人有更好的想法来解决这个问题？