问题标签 [ogr]

我有一个要转换为 GeoJson 的 csv 文件。

我运行以下 ogr2ogr 命令

我得到了这个结果

为什么几何为空？我如何告诉 ogr2ogr 哪些值是纬度和经度？

我一直在一个项目中使用 GDAL/OGR，它进展顺利。我一直在将各种源格式的数据重新投影到单个 CRS（使用 C# 绑定）中，但我遇到的一个问题是相关的OGR 样式信息。

源文档中的大多数字体大小都存储在“地图地面单位”中，所以如果我直接复制 OGR 样式，文本大小将根据地图单位的差异放大或缩小。

所以我的问题是，将样式规范中的地图地面单位从地图地面单位转换为像素单位的最佳方法是什么？

非常感谢您提供的任何帮助，

我想将 .DWG 文件读入 R。DWG 在 GDAL 中支持的 OGR 格式列表中：

http://www.gdal.org/ogr/ogr_formats.html

所以理论上，我应该能够在 rgdal 中使用 readOGR() 读取 .DWG 文件。但我不知道如何定义dsn。

ogrDrivers() 命令为我提供了可用驱动程序的列表，但该列表与 GDAL 网站上的可用 OGR 格式列表不匹配。

根据 DWG 文件的 OGR 页面，应该有一个 DWG 驱动程序，但它甚至不在 rgdal 中的 OGR 驱动程序列表中。

这是否意味着 rgdal 不支持 DWG 文件，即使 GDAL 中的 OGR 支持它们？

非常感谢。

给定一个 shapefile：

Natural_earth / ne_10m_admin_0_sovereignty.zip

鉴于我们想为 D3js 数据可视化重新投影它，我们可以在不同的层次上重新投影。

1.使用以下方法获取重新投影的 shapefile ( 1 ) ogr2ogr：

或2.使用 (npm)获取重新投影的 topojson ( 2 ) topojson.js：

或3.获取重新投影的 D3js 数据/SVG ( 1 )，D3js 代码包括：

概述：

Mike Bostock > Projected Topojson告诉我们，第一种和第二种方式“消除了在渲染时投影几何图形的需要，提高了性能 [...] 因为每个点的重要性是在屏幕上的区域而不是地球表面上测量的。” 总之，最终像素质量/文件权重比更好。

另一方面，在渲染时重新投影几何体可以在最后一刻进行更灵活的投影。

了解更多？

这就是我所知道的。如果有人可以对这些方式进行更多解释，分享有关参数的帮助资源（ogr 的 EPSG 列表，d3js 投影列表），以及每种方式的各自优势/劣势，这可能是一本非常有趣的并行手册。

注意：我会给出我的答案，但我只是开始挖掘它。我想周围有更多有经验的人。

我有一系列地震的经度和纬度。我希望能够将它们分成陆地上的和海上的。

是否有 r 函数可以做到这一点？

我正在尝试获取带有美国州名和 ID 的 topojson 文件。使用http://www.tnoda.com/blog/2013-12-07供参考。

我读到有人建议我们生成自己的 topojson 文件，我一直在尝试让 ogr2ogr 在 MAC OSX v 10.7.5 上工作。

在命令行中运行 ogr2ogr 给出：

我不确定它在抱怨什么。libjpeg 的最新版本不是 13 而是 9。我认为它可能需要版本 8。但是在按照http://ijg.org/files/jpegsrc.v8a.tar.gz的安装说明进行操作后，我仍然遇到同样的错误 。

任何帮助将不胜感激。谢谢！

PS：我的最终目标是获得这样的地图http://bl.ocks.org/mbostock/4060606但我需要 json 中的 2 个字母状态名称，因为这就是我的数据中的内容。

我的组织使用批处理文件对我们的数据库执行查询，我遇到了一种情况，即我需要执行一个太大而无法一次完成的查询。所以我所做的是修改一个批处理文件来循环一个文本文件并单独查询每一行并将结果附加到输出文件中。

问题是当我回显查询变量时，它的读取完全符合预期，但是当我将它传递给 ogr 代替 sql 字符串时，它似乎是空白的。我不确定我是否做错了什么，或者我想要做的事情是不可能的。谁能澄清这里发生了什么？

我正在使用python读取shapefile，但是遇到了一些问题，这是核心代码：

如果我直接运行脚本，它会按预期工作，但是一旦我在 Web 环境中运行它就会抛出错误，我会得到如下错误：

如果我将行更改为：

有效。

这是怎么回事？

似乎replacepython中的函数让事情变得很奇怪。

更新

似乎我明白了，这是由replacepython中的函数引起的，它返回一个不同的而不是str：

然后我使用：

有效。

但我不确定为什么它在 shell 环境中工作，而不是在 web 环境中工作。我希望有人能解释一下。

OK, After a bit of fiddling, I've tweaked a script from the site hyperlink in the second comment line. The purpose of the script is to clip/mask a LARGE raster (i.e. that cannot fit into a 32-bit Python 2.7.5 application) in GTiff format with a shapefile with multiple polygons (each with a "Name" record) and save the clipped rasters into a "clip" sub-directory, where each masked grid is named after each polygon's "Name". Like the original script, it assumes that the GTiff and shapefile are in the same projection and overlap correctly, and it processes ~100 masks/sec. However, I have modified that script to 1) work with a float-valued elevation grid, 2) only load the window of the larger grid into memory that is bounded by the current polygon (i.e. to reduce the memory load), 2) exports GTiff's that have the right (i.e. not shifted) geo-location and value.

HOWEVER, I having trouble with each masked grid having a what I'll call a "right-sided shadow". That is for every ~vertical line in a polygon where the right side of that line is outside of the given polygon, the masked grid will includes one raster cell to the right of that polygon-side.

Thus, my question is, what am I doing wrong that gives the masked grid a right-shadow?

I'll try to figure out how to post an example shapefile and tif so that others can reproduce. The code below also has comment lines for integer-valued satellite imagery (e.g. in as in the original code from geospatialpython.com).

在 Java 中对 OGR 做一些事情 - 特别是从 EPSG 代码中导入空间参考，我得到这个错误：

问题是（与我看到的关于相同错误的堆栈上的另一个问题不同）我的用户环境变量（win7）中已经有一个GDAL_DATA指向C:\Program Files\GDAL\gdal-datagcs.csv 存在的位置。

我确实-Djava.library.path="C:\Program Files\GDAL\"在我的项目的运行时选项中有一行以确保链接 GDAL（即使我也有一个指向那里的环境变量 PATH），我是否需要做类似的事情或在代码中强制/设置GDAL_DATA环境变量或者这不是问题，这是别的什么吗？（不会是第一次 GDAL/OGR 绑定对我来说很奇怪）