gdal_rasterize can miss points on the extent edges
|Reported by:||EliL||Owned by:||warmerdam|
|Severity:||normal||Keywords:||rasterize, edge, half-pixel-shift|
Running gdal_rasterize -tr 100 100 -a Elevation -l town town.shp townraster.tif converts 12 points to 9 pixels (with values) in the output raster. It should convert it to 11 pixels (with values) in the output raster (1 point has a null value and should not get converted).
town.shp is a point layer with 12 points, 11 of those points have values in -a Elevation (1 is null), and in the output tif, 9 pixels have values. The far south and far east points do not have corresponding pixels with values.
Viewing the results in qgis, the points and pixels do not always line up correctly, perhaps it is the same cause for the extent edges to be missing pixels (i.e. those pixels would be in the next row or column of pixels).
Results should be reproducible with town.shp, attached.
Change History (8)
comment:4 Changed 5 years ago by
|Priority:||low → normal|
|Severity:||minor → normal|