问题标签 [rasterize]

我正在尝试平滑 R 中不规则多边形的边缘，即将其尖角变成圆形边缘。我正在尝试使用 来执行此操作smoothr::smooth，但是此函数对包中的对象进行操作，sf或者sp我所拥有的只是一组坐标。不知何故，将 mydata.frame变成一个SpatialPolygonsDataFrame对象（来自 package 的对象类sp）的结果是一个矩形，其极限是原始多边形的极限。有谁知道如何smoothr::smooth在保持原始多边形形状的同时将我的一组坐标转换为兼容的类的对象？这是我所做的，部分遵循此页面上的说明：

红色是我的原始多边形data.frame dd，绿色是对象sp，蓝色是平滑版本sp，sps。函数smooth完成了这项工作，问题出在转换dd为sp兼容对象的某个地方。我怀疑问题是由raster()但我不知道为什么或如何解决它。提前谢谢了。

我最近一直在使用 fasterize 包将 sf 多边形转换为栅格：

https://cran.rstudio.com/web/packages/fasterize/fasterize.pdf

当我处理大文件时，我最好直接写入磁盘而不是内存。因此，例如，而不是这样做：

更快（多边形文件，光栅模板，字段='值'）

我会这样做：

更快（多边形文件，光栅模板，字段='值'，文件='output.tif'）

这似乎不可能。有没有人有关于如何做到这一点的建议？

谢谢。

如何构建一个循环来逐步光栅化几个简单的特征？

我正在寻找的内容可以通过以下方式创建：

但是，我确信有一种方法可以通过循环接收相同的结果。结果可能是 Rasterstack。

我正在使用gdal.RasterizeLayer()GeoTiff 模板将 shapefile 转换为 GeoTiff，同时通过ATTRIBUTE. 我要输出的是一个 .tif，其中的 burn 值对应于给定属性的值。我发现它gdal.RasterizeLayer()正在燃烧为与我的属性字段中的值不对应的奇怪值。这是我目前拥有的：

发生的情况是输出 .tif 为每个属性正确分配了不同的烧录值，但该值与属性值不对应。例如，输入属性值 FCode=46006 变成燃烧值 182（不清楚为什么！）。我尝试添加和删除该'COMPRESS=DEFLATE'选项，并添加和删除gdal.RasterizeLayer(). 无影响输出燃烧值。

您可以在此处查看输入 shapefile 和属性值：输入 .shp

输出的值不正确，在这里：输出栅格

我正在寻找的是一种在一定范围内逐步光栅化（或加速）几何图形并将光栅化几何图形的各个部分连接到整个光栅对象的方法。

我猜，裁剪每个图块内的多边形（或多边形的一部分）效果很好。我的实际目标是逐块处理数据块。看起来下面的循环执行此操作，并将多边形的部分连接到整个多边形。好吧，输出（cr）与输入（p）相同......我认为这一步可能是必要的，以便在下一步中光栅化生成的（部分）多边形。我没想到多边形会被连接起来。我试图根据中间结果构建一个循环，因此添加了以下几行：

但是，我很难栅格化（中间）结果。

我想计算每个栅格单元内线串的总长度。我有一个繁琐的解决方法，但我认为有一种更有效的方法可以直接使用该rasterize函数执行此操作。下面的 MWE 显示了我想要什么以及我迄今为止的工作解决方案。

我在使用 Raster* 对象时遇到了问题。我想从栅格的顶行和底行中提取坐标，然后使用它们创建一个新栅格，其值仅在这些单元格中。这可能是一个不好的解释，所以我将提供一个示例：

这完美地工作。但是，我想在更大的范围内进行。因此，我将范围和ncol参数替换为nrow与我正在使用的真实世界栅格等效的值。见下文：

除了，现在只改变形状和分辨率后，rasterize函数只返回第一行！我很困惑...

我认为 肯定有问题rasterize，因为r_NS_nodes包含顶行和底行的坐标：

请，如果有人可以向我解释这一点，或提供解决方案，我将不胜感激。

-亚历克斯。

我在 R 中很新，所以我希望你们能帮助我完成这一步，因为我很卡住。

我有一张表，上面有监测点（坐标）和每个监测点的水质指数。我正在尝试使用 rasterize 函数从该表创建一个栅格：

rasterize(x, raster, field, fun=mean, backgroud=0, update=TRUE)

x =带坐标的数据框

raster=使用来自参考栅格的范围、nrows、ncols 等创建。

field=使用水质指数数据创建的向量。每个值都属于 x 表中的一个点。

fun=mean=我需要同一个单元格中所有点的平均值。

background=0以完成没有从字段派生的值的单元格。

update=TRUE因此单元格中的值会使用向量值进行更新

运行此程序后，我得到了一张带有正确扩展但没有数据（全为空白）和下一个摘要的地图（有时我得到的是 NA 或 2 而不是 5，不知道为什么）：

我认为问题在于功能，但我不确定。我尝试按如下方式更改功能，但结果相同。

非常感谢您提前。干杯。

在尝试栅格化伦敦病房 (london_ward) 的空间多边形并用这个新栅格掩盖每个人口栅格（pop_2011 和 pop_2001）时，我遇到了栅格化函数的问题。下面是我被告知要使用的代码：

但它继续给我以下错误代码： h(simpleError(msg, call)) 中的错误：在为函数“addAttrToGeom”选择方法时评估参数“x”时出错：初始化列表中没有参数

我正在尝试获取在伦敦的空间范围内裁剪的地图，如图 1 所示，但由于此错误而无法这样做。不完全确定我哪里出错了，非常感谢任何帮助。谢谢你。

我想在 QGIS 3.16 Windows 10 中获取我的街道多边形的中心线。我按照在 stackoverflow 和其他地方找到的步骤进行操作。

1. 在 QGIS 中使用来自 GDAL 的栅格对多边形进行栅格化。（请看结果） 光栅化多边形的结果。数据保存在 INT32 中，并且具有“无数据”单元格或 205684。

2. 在下一步中，我想使用 r.thin 创建栅格数据的中心线。我没有得到任何结果，协议中的错误如下：

来自协议的错误

我也尝试定义输出文件。但我得到同样的错误。有没有人有任何经验可能是什么问题？

太感谢了