0

我正在使用 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 开关

gdal_rasterize -l !fileOUT!  -a OP_DATAFIELD -tr 100.0 100.0 -a_nodata -9999 -te 44780.0 4155.0 670380.0 1240255.0 -ot Int16  C:\WorkingMDT\!SHPfileIN!.shp C:\WorkingMDT\!fileOUT!.tiff

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

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

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

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

4

0 回答 0