Gdal Rasterize在栅格中插入不必要的空值行,从而使后续的行错分

时间:2019-04-24 12:16:07

标签: gis gdal rasterize

我正在使用GDAL.rasterize对点的简单shapefile进行栅格化。 shapefile点仅由Xco-ord,Yco-ord和整数数据值组成。关于输出文件的所有内容都很好,除了仅包含No_data的单行,该行是在栅格向下3/4处的过程中神秘插入的,因此下面的所有后续行随后都向南偏移100m。但是shapefile中存在该异常行的数据?

我尝试创建其他格式而不是TIFF,例如EHDr,但结果都相同

所以,我认为这与记忆有关,所以我尝试减小他的程度

如果仅将范围缩小到插入的空行的行以下的栅格,则输出仍将偏移

如果我将范围缩小为仅插入插入的空数据行上方的栅格,那么对于该部分,结果输出就很好,就像我进行整个范围操作一样

我减少了cols的数量但保留了行,只是发生了细小的事情

所以我认为它的内存不再相关

输出栅格为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

0 个答案:

没有答案