问题标签 [ogr]

我正在尝试开发（使用 C++ - MSVS 12.0）一个函数，该函数发现哪些像素（来自光栅图像）的中心位于多边形内（之前使用 shapefile 填充）。我正在使用从源代码构建并使用选项 INCLUDE_OGR_FRMTS=YES 的 GDAL 1.11.0（刚刚安装，使用 devinstall）。我可以毫无问题地使用 GDAL 和大多数 OGR 功能。但是，当我使用以下代码时：

if (polygon->Contains(tmpPoint))

我收到错误消息：错误 6：未启用 GEOS 支持

有人知道如何解决这个问题吗？

我在用着：

#include "ogrsf_frmts.h"

我的函数被声明：

void FindPixels(GDALDataset *image, OGRLayer *poLayer, OGRPolygon *polygon)

我的部分代码是：

OGRPoint *tmpPoint = NULL

OGRSpatialReference *spatialReference = NULL;

spatialReference = polygon->getSpatialReference();

tmpPoint = new OGRPoint();

tmpPoint->assignSpatialReference(spatialReference);

循环开始：

tmpPoint->setX(imgTLX + (j * imgRes) + imgResHalf);

tmpPoint->setY(imgTLY - (i * imgRes) - imgResHalf);

if (polygon->Contains(tmpPoint))

提前致谢！

MB

我试图查看 GDAL 中的 OGR api，并在此处遵循教程：http ://www.gdal.org/ogr_apitut.html 。我不断得到对几乎所有内容的未定义引用。我已经在我的 linux 发行版上从源代码安装了最新版本 (1.11)：http: //trac.osgeo.org/gdal/wiki/DownloadSource。查看头文件，特别是在 gcore/gdal.h 中，我看到例如 GDALOpenEx 不存在，也不是在此之前的版本中。但是，它位于 subversion 存储库中。那么，什么给了？

我有所有这些包含数据的文件，我想将所有数据转换成一个 geojson 文件。

我正在尝试这样做：

但我认为它没有做它应该做的事情。我确实得到了一个 pga.geojson 文件，但它只包含一个矩形：

我认为该文件还应该包含一些似乎没有出现的行。为了让它使用所有文件中的数据，我还应该做什么？

我正在尝试使用 Python 中的 OGR 将 csv 文件转换为 ESRI Shapefile。问题显然是我的标题条目对于默认字段大小来说太长了：

我想遍历所有标题元素（我的 csv 文件中的列数未知）并在 shapefile 中创建新字段 - “即时”，可以这么说：

工作正常，但标题条目被切断（这是有问题的，因为我需要稍后解决它们）。有没有办法在同一步骤中定义标题条目大小？我试过了SetField(32)，但我担心这只会改变表格条目的大小。任何人都可以帮忙吗？提前致谢！

我有一个想要转换为 KML 文件的 SQL Server 表。我有点让它工作，但不完全。这是我的命令：

我得到了一个有效的 kml 文件，但它缺少 POINT。以下是所有<Placemark>标签的外观：

我应该在那里有这样的东西：

当我在 Google 地球中查看 KML 时，没有要显示的点。

“geog”列是地理数据类型（SRID：4269）。

我也尝试过使用-f LIBKML，但遇到了同样的基本问题（没意义）。我看到这个驱动程序有很多附加选项，但我不知道如何调整对ogr2ogr.

一旦它工作，我想从预定的 Powershell 脚本中调用它。

想法？

我正在尝试使用 ogr2ogr 将 GML 文件导入 SQL Server Spatial。我成功地导入了几何特征，但没有几何列的我很少。我怎样才能将它们全部导入？

编辑：

我重新安装了 GDAL，安装了最新版本，工作正常，但再次无法编写非空间特征。

不断报错：

错误 1：创建层时出错：[Microsoft][ODBC SQL Server Driver][SQL Server]“NULLCONSTRAINT”附近的语法不正确。

我正在使用 D3 创建美国地图，我必须实施三级缩放，从国家地图到州，从州到县，从县到邮编边界。

我已经完成了前两个级别，我使用了 AngularJS 指令。

GIST URL：D3 二级缩放

现在我需要实施进入县的第三级并显示所有 zip 级区域边界。我已经对 zip level geometry 进行了很多搜索，我在 stackoverflow 中发现了一个有用的问题，但是这个数据只关联了州名，我还需要县名，因为我只需要为被点击的特定县绘制边界，任何线索如何正确将是有帮助的。我发现US Atlas project，虽然希望它会帮助我，但无法正确理解它。

有关如何将此 shapefile剪辑到适当的边界框的任何建议？特别是，我只想查看这些界限：[(-105.53 , 39.9) (-104.93,40.27)]。ogr2ogr的站点有一个使用-clipdst看起来很有希望的选项的示例。我使用Anaconda来管理我的 python 安装。我安装了GDAL，这样conda update gdal会产生：

跑步which ogr2ogr

下载并解压缩 Natural Earth 子单元文件后，我运行以下命令：

不幸的是，我收到此错误：

我查看了这个错误，这里的建议是启用geos（我使用安装的一个包conda install geos）。所以我有geos版本 3.3.3。但是，我不知道如何启用ogr2ogr访问geos库。

任何帮助是极大的赞赏。我不打算使用ogr2ogr，但我没有看到很多替代品。

更新：2015-02-26

重新安装 anaconda 后，我安装了 gdal 和 geos。这是目前的情况：

我对开源 GIS 软件很陌生。我正在尝试在我的 MAC 终端中使用 GDAL 和 ogr 剪辑和投影数据层。我可以在剪辑时获取要剪辑的数据，但是当我尝试将两者（剪辑和项目）结合起来时，出现语法错误。我的语法是为命令行编写的，我想我可能需要为终端调整它。任何帮助将不胜感激！如果有人愿意提供帮助，我会在必要时尝试提供更多信息。

这是我的代码行：

当我运行代码时，我收到以下错误消息：

伙计们，

刚刚开始使用 OGR 和 Python 完成各种地理空间任务。我在 OSGEO4w 之外工作，并在我的机器上安装了带有 Python 绑定的 GDAL 以及 Python v. 2.7.8。

也就是说，我可以从命令行界面运行 python 并导入 gdal，但在使用 IDLE 环境时无法导入模块。它返回一个错误，告诉我该模块不存在。鉴于它在 cmd 提示符下工作，我的安装必须是正确的，所以有什么关系？

这对我来说是全新的，我将不胜感激您能提供的任何帮助或建议；我读过的所有以前的问题都涉及仅安装（我已经成功完成）并让 Python 从命令提示符中识别模块，它确实如此。

干杯，

麦克风