问题标签 [rasterize]

我在将多边形转换为 R 中的栅格时遇到了很大的麻烦。我想做的是：我有 574 个物种的 shapefile（即多边形）。也就是说，在属性表中它有 574 行（即 FID 在 0 到 573 之间）。可以在此处找到数据子集：https ://drive.google.com/file/d/1AdTChjerCXopE1-PZIAPp5ZXISdq45i8/view?usp=sharing

我想将其转换为栅格。在输出栅格中，我看到最小值和最大值分别为 1 和 574。我怀疑的是：它将单元格中的字段 ID 作为不应该的像素值获取。单元格值应来自覆盖多边形。任何帮助将不胜感激。下面是示例代码：

输出栅格的属性如下：

我有一个由两个字段/属性组成的 shapefile，一个是整数，另一个是字符串。

gdal.RasterizeLayer()我可以通过首先使用函数将 shapefile 刻录为 .tiff 图像作为第一个波段来将整数提取到 Python 数组中。然后，我使用my_raster.GetRasterBand(1).ReadAsArray()将整数作为数组读取。

但是，我想从其他字段/属性中提取字符串值。我做了完全相同的事情，但我已经更改了gdal.RasterizeLayer()规范中的属性名称。但是，调用GetRasterBand(1).ReadAsArray()只会给我零。

有谁知道是否可以从栅格中读取字符串？

顺便说一句：我使用与此处完全相同的代码。

从

纯 Python 版本——gdal.RasterizeLayer

我需要将 PDF 转换为不同的 ICC 颜色配置文件。通过不同的搜索和测试，我找到了一种方法：

• 首先，我将我的 PDF 转换为 PS 文件：

  /li>

• 然后我使用以下内容将 PS 转换回 PDF（这是为了生成有效的 PDF/X-3 文件）：

  /li>

ICC 配置文件已嵌入，一切正常。唯一的问题是整个最终 PDF 都被光栅化了。在这里，我失去了起始文件中的所有路径和其他矢量元素质量。我需要保持它们是矢量的，因为这个 PDF 将有一个特定的应用程序。

硬件加速是 Direct2D 支持的一项功能。这是我的问题。

据我所知，硬件加速受 GPU 型号、驱动程序版本等限制。有人知道这件事的细节吗？也就是说，如何判断一台电脑是否支持d2d硬件加速。

下图是在 Chrome 浏览器中捕获的。

我希望为每个 ISO3（空间多边形数据帧 wrld_simpl 中的列名）提取栅格 r 的像素值的总和。我正在考虑使用函数 rasterize，然后是 zonal，但是当光栅化 wrld_simpl 时，我丢失了定义 ISO3 的字符串（例如 AUS、USA...）。非常感谢您的建议！理想情况下，我的最终输出将是一个数据帧，其中每个 ISO3 都与一个值相关联（相应的像素值总和）

我已经通过以下示例使用我自己的线样带数据创建密度表面模型。

输出是使用 ggplot2 的阴影/渐变密度图

但我不知道如何将其转换为栅格，以便我可以提取特定 GPS 点的猎物密度。据我了解，ggplots与空间工具不兼容，所以我需要转换它吗？但是怎么做？任何帮助将非常感激！

我必须对 30 个大 shapefile（>800 万个单元格）的每 62 列进行栅格化，从而产生 1860 个栅格。

我能够在 R 中以较小的形状（约 8000 个特征）轻松快速地做到这一点。然而，当只是尝试在 R 中加载这个大形状时，我花了 2 个小时，16G RAM 和 47% 的 64G 交换。调用 rasterize 函数时，由于内存问题，无法运行。但是我发现 QGIS 可以非常快地做到这一点，但是对于每列一次，这将花费我的生命来运行它。我尝试在 QGIS 中使用带有循环的 Python 控制台来遍历列，但没有成功。

我首先收到以下错误：

然后，似乎 gdalogr:rasterize 无法处理FIELD = n,。我对么？

请你帮我处理这段代码好吗？谢谢！

作为一名编码爱好者，我目前正在创建自己的“地理空间”工具，就像我希望它的工作方式一样。然而，一开始我已经面临一个问题。我的工具应该使用 GeoPandas 来提取信息，然后使用 OGR/GDAL 进行数据编辑，因为我希望它能够快速工作。我喜欢分析很多大数据！

带有问题的代码片段应该光栅化单个 GeoPandas 多边形。我尝试使用这条路径来做到这一点。- 使用 geopandas WKT-polygon 从多边形中提取 - 使用 WKT-polygon 创建 OGR 特征 - 使用 GDAL 将其栅格化。

我面临的问题是我只检索一个由 0 组成的栅格，而不是 0 和 1 ......

代码如下所示：

非常感谢能给我提示正确方向的人:-)。

干杯!

我正在使用 GDAL.rasterize 栅格化一个简单的点形状文件。shapefile 点仅由 Xco-ord、Yco-ord 和一个整数数据值组成。关于输出文件的所有内容都很好，除了仅包含 No_data 的单行，该行在栅格向下约 3/4 处神秘地插入，因此其下方的所有后续行随后出现向南 100m 未对齐。但是该异常行的数据存在于 shapefile 中？

我尝试过创建其他格式而不是 TIFF，例如 EHDr，但结果都一样

所以，认为这是与记忆有关的，我试着减少他的程度

如果我将范围缩小到仅在插入的空行的行下方栅格化，则结果输出仍然是偏移的

如果我将范围缩小到仅在插入的空数据行上方光栅化，则结果输出对于该部分来说很好，就像我做整个范围时一样

我在减少列数但保留行的范围内采取了一些措施，同样的事情发生了

所以我认为它不再与记忆有关

输出栅格是 100mx100m 网格中的简单 6256 列乘 12361 行 tiff

范围是 45080,670080,4355,1240255 CRS 是 EPSG：27700

这是我使用的 Gdal .rasterise 开关

我需要发生的是一个没有异常行 3/4 的栅格，它抵消了所有后续行，

是什么导致插入该异常行？

现在我可以通过转换为 .asc 并编辑异常行来手动纠正这个问题，但我宁愿找到程序原因

非常感谢所有帮助和考虑。这是问题的图片 绿色是创建的栅格，蓝色与原始数据点交叉，可以清楚地看到已插入的 no_data 行以及随后的向下移位enter image description here

我有一个由 3 个多边形组成的 spatialPolygonsDataFrame。第三个多边形与第一个多边形具有相同的形状，但在第二个多边形所在的位置有一个孔。

我使用另一个问题的答案（如何在 SpatialPolygonsDataFrame 中向多边形添加一个洞？）。

当我mask()用来“裁剪” rasterr1时，会创建孔，尽管三角形多边形具有值并且确实不是真正的孔。但是它被带有孔的第三个多边形“覆盖”：

当我更改多边形的顺序时，不会创建任何孔：

函数 rasterize 以相同的方式受到影响：

在我的真实数据中，我有没有“填充”多边形的洞，而那些我想保留为洞的洞。

如何修复正在创建空洞的多边形的 spatialPolygonsDataFrame？

如何在不重新排序但“转换”孔创建多边形的情况下解决此问题？