问题标签 [gdal]

我正在尝试运行一个简单的脚本来手动创建一些几何图形，然后对所述几何图形执行缓冲区。但是，每当我尝试执行 Buffer 操作时，python 都会完全崩溃。

我正在运行 Windows 7 64 位。

我正在尝试使用 GDAL 和 Python 实现图像分类的最小距离算法。在计算样本区域的平均像素值并将它们存储到数组列表（“sample_array”）中后，我将图像读入一个名为“values”的数组中。使用以下代码，我遍历该数组：

对于大图像，此过程需要很长时间。这就是为什么我想问是否有人知道更快的方法。我对python中不同变量的访问速度了解不多。或者也许有人知道我可以使用的库。在此先感谢，马里奥

我正在尝试在我的 Windows 7 工作站上安装 GDAL 和 Python。

• Python 版本：2.6.6
• 格达尔 1.8
• （Visual Visual C++ Studio 2010 速成版）

我按照http://ucsb.pbworks.com/w/page/1936549/Using-GDAL-with-Python上的说明将 PYTHONPATH 和 GDAL_DATA 添加到我的系统环境变量中。此外，我通过添加 GDAL 二进制文件夹修改了我的 PATH 变量。

我从以下网址下载的 GDAL Python 绑定：http: //vbkto.dyndns.org/sdk/PackageList.aspx ?file=release-1600-gdal-1-8-mapserver-5-6.zip

“ from osgeo import gdal ” 导致以下错误消息：

在此先感谢，

我正在尝试使用 Delphi 和 GDAL18 创建等值线。为此，我使用以下代码：

GDALContourGenerate 命令返回“不支持的几何类型” - 错误。

我通过以下方式包含了 gdal18.dll：

我也尝试过其他几何类型（例如 wkbLineString25D），但这并没有帮助。如果您有任何建议，我会很高兴。非常感谢，马里奥

[编辑]我发现当我用“nil”重播“layer”（在 GDALContourGenerate 中）时会发生同样的错误。所以也许问题出在其他地方。[/edit]

我在一个脚本中遇到错误，该脚本迭代目录中的一些文本文件，进行一些数据操作并使用 R 中 rgdal 包中的 writeOGR 写入 ESRI shapefile。在一定数量的文件之后，我不断收到错误被写入（似乎是大约 43-44 个单独的文件，或大约 14 个完整的 shapefile，因为每个 shapefile 创建了许多文件）。

错误通常类似于

此外，脚本崩溃后，我收到一个错误，例如

当我尝试读取另一个文件或做任何事情时。

此时不显示任何打开的文件。

这是一些重现问题的代码：

我在 OSX 上使用 R 2.12.2。错误发生在 64 位和 32 位版本上。

有没有人在 GDAL 工具中使用过 ogrinfo [-spat xmin ymin xmax ymax] 选项？我能够在形状文件上运行 - sql 查询并获得答案/形状，但是如果我在 spat 中使用 lat/long 值，我不会得到结果（但是我没有得到错误）。我根本找不到一个例子。

我想使用 python、gdal 和 numpy 创建一些高程/高度场栅格。我被困在 numpy 上（可能还有 python 和 gdal。）

在 numpy 中，我一直在尝试以下操作：

从 osgeo 导入 gdal 从 gdalconst 导入 *

我想我错过了一些简单的东西，期待你的建议。

谢谢，

克里斯

（稍后继续）

• terragendataset.cpp,v 1.2 *
  • 项目：Terragen(tm) TER 驱动程序
  • 目的：Terragen TER 文档的阅读器
  • 作者：Ray Gardener，Daylon Graphics Ltd. *
  • 该模块的部分来自 GDAL 驱动程序
  • 弗兰克·沃默丹，见http://www.gdal.org

我提前向 Ray Gardener 和 Frank Warmerdam 道歉。

Terragen 格式说明：

对于写入：SCAL = gridpost 距离，以米为单位 hv_px = hv_m / SCAL span_px = span_m / SCAL offset = 参见 TerragenDataset::write_header() scale = see TerragenDataset::write_header() 物理 hv = (hv_px - offset) * 65536.0/scale

我们告诉来电者：

这告诉我，在上面的 WriteArray(somearray) 之前，我必须同时设置 GeoTransform 和 SetProjection 以及 SetUnitType 才能工作（可能）

来自 GDAL API 教程：Python import osr import numpy

我很抱歉创建了一个过长的帖子和一个忏悔。如果我做对了，将所有笔记放在一个地方（长篇文章）会很好。

忏悔录：

我之前拍摄了一张照片（jpeg），将其转换为 geotiff，并将其作为瓷砖导入 PostGIS 数据库。我现在正在尝试创建用于覆盖图片的高程栅格。这可能看起来很愚蠢，但我希望向一位艺术家致敬，同时又不冒犯那些辛勤工作以创造和改进这些伟大工具的人。

这位艺术家是比利时人，所以米是有序的。她在纽约曼哈顿下城工作，因此，UTM 18。现在一些合理的 SetGeoTransform。上面提到的图片是w=3649/h=2736，我得考虑一下。

另一个尝试：

显然越来越近，但不清楚 SetUTM(18,1) 是否被拾取。这是哈德逊河中的 4x4 还是 Local_CS（坐标系）？什么是无声的失败？

更多阅读

4x4 米是一个相当小的逻辑跨度。

所以，也许这是最好的。SetGeoTransform 得到正确的单位，设置比例，你就有了你的 Terragen 世界空间。

最后的想法：我不编程，但在某种程度上我可以跟随。也就是说，我首先想知道我的小 Terragen 世界空间中的数据是否以及是什么样的（以下感谢http://www.gis.usu.edu/~chrisg/python/2009/第 4 周）：

所以这是令人欣慰的。我想上面使用的 numpy c 和这个结果之间的区别在于在这个非常小的逻辑跨度上应用 Float16 的操作。

然后到'hf2'：

几乎完全令人满意，虽然看起来我在 La Concordia 秘鲁。所以我必须弄清楚如何说-like more north, like New York North。SetUTM 是否采用第三个元素来暗示北或南“多远”。似乎我昨天遇到了一个带有字母标签区域的 UTM 图表，比如赤道的 C 和纽约地区的 T 或 S。

我实际上认为 SetGeoTransform 本质上是建立左上角的北向和东向，这影响了 SetUTM 的“北/南多远”部分。转到 gdal-dev。

后来还是：

帕丁顿熊去秘鲁是因为他有一张票。我到那里是因为我说那是我想去的地方。Terragen 以它的方式工作，给了我我的像素空间。随后对 srs 的调用作用于 hf2 (SetUTM)，但东移和北移是在 Terragen 下建立的，因此 UTM 18 已设置，但位于赤道的边界框内。够好了。

gdal_translate 带我去了纽约。我在windows中，所以是命令行。和结果；

所以，我们回到了纽约。可能有更好的方法来处理这一切。我必须有一个接受 Create 的目标，因为我也从 numpy 假设/即兴创作数据集。我需要查看其他允许创建的格式。GeoTiff 中的海拔是一种可能性（我认为。）

感谢所有评论、建议和对适当阅读的温和推动。用python制作地图很有趣！

克里斯

我已经下载了 USDA 作物数据层（数据可在此处获得，请注意大文件： http: //www.nass.usda.gov/research/Cropland/SARS1a.htm）并希望将其导入 R 进行分析。我导入了其他大型 GTIFF 文件（例如 LANDSAT 数据），但在调用该函数时出现以下错误：

其中“文件”是美国农业部提供的任何 .tif 文件（我使用的是中西部地区汇编）。

我怀疑问题是 GTiff 是分类数据（土地利用类别）而不是真正的序数/光谱数据，但我不明白为什么它会在图像编辑器（例如预览）中呈现良好并在 R / GDAL 中退缩.

谢谢您的帮助。

编辑：如果这应该在 GIS 论坛上进行，请告诉我。一般来说，空间问题的编程方法似乎就在这里，不是吗？

注意：我已经从每个命令中删除了其他 5 个文件名，以使其更具可读性。我正在开发一个 Python 工具，它将为我构建一个批处理文件，但是当我尝试运行这个 .bat 时，我得到的只是第一个成功运行的命令，然后返回到控制台。我在 FWTools shell 中运行它。我已经在十六进制编辑器中打开了批处理文件并搜索了额外的换行符并返回字符，但我没有看到任何异常，所以 Python 似乎正在正确编写。

简要地：

如何参数化 .prj WKT 文件，以便我可以执行 7 点转换（wiki）。我知道false_easting and false_northing参数是如何工作的，但是如何调整比例？我不是说 scale_factor'

这就是问题描述：

我在非 GIS 环境（交通建模软件）中保存了交通网络（矢量图层）。网络由节点（点）和折线（道路链接）组成。它主要是从随机背景完成的，不管任何投影、坐标等。

我需要为网络设置适当的投影。

我可以访问 .prj 文件（如果我在说 WGS84 投影，我可以切换到任何其他投影）

所以这就是我正在尝试的：
我尝试 7 点 Helmert 转换（http://proj.maptools.org/gen_parms.html）。我在 .prj 文件中使用 towgs84 转换作为 WKT 参数，我假设旋转矩阵为零（我可以这样做吗？）并且我只计算 delta_x、delta_y 和缩放参数。

但是它不会起作用。这是我的 .prj ，TOWGS84 中的参数不影响转换：

所以我尝试使用false_norting和false_easting参数，这些工作很好，并适当地改变我的网络，但是：

它不会改变我的网络规模，只会改变位置。那么如何使用 .prj 文件重新调整我的网络呢？

感谢您的任何提示