gdalwarp on 3.4.3 behaving differently to 3.4.2

[qgib]

Author Name: Ross Edwards (Ross Edwards)
Original Redmine Issue: 20912
Affected QGIS version: 3.4.3
Redmine category:gdal_tools
Assignee: Ross Edwards

---

I found an entirely repeatable issue with gdalwarp commands that worked perfectly on 3.4.2, but on 3.4.3 they appear to generate slightly off height and width.

on 3.4.2 I run the following 2 commands separately in Python console -

gdalwarp -t_srs EPSG:32631 -wo SOURCE_EXTRA=1000 -tr 5.0 5.0 -r cubic -of GTiff "E:/Users/rosso/Arma/Terrain Heightmaps/Heightmaps Opentopo/output_srtm.asc" D:/Arma/QGIS/crm4.tif

gdalwarp -of GTiff -cutline E:\Users\rosso\Arma\QGIS\UTM31N\Test23-11-18\shape.shp -crop_to_cutline "D:/Arma/QGIS/crm4.tif" D:/Arma/QGIS/cut6.tif

which generates a perfect 8192 x 8192 tif
https://ibb.co/6PmvQ8x

However running exactly the same on 3.4.3

gdalwarp -t_srs EPSG:32631 -wo SOURCE_EXTRA=1000 -tr 5.0 5.0 -r cubic -of GTiff "E:/Users/rosso/Arma/Terrain Heightmaps/Heightmaps Opentopo/output_srtm.asc" D:/Arma/QGIS/crm5.tif

gdalwarp -of GTiff -cutline E:\Users\rosso\Arma\QGIS\UTM31N\Test23-11-18\shape.shp -crop_to_cutline "D:/Arma/QGIS/crm5.tif" D:/Arma/QGIS/cut7.tif

I get a tif with 8191 x 8191
https://ibb.co/0GcptBT

[qgib]

Author Name: Ross Edwards (Ross Edwards)

---

Just to confirm the shapefile feature used to cut the tif is a perfect square 40960 x 40960 so at 5m resolution I expect 8192 output

[qgib]

Author Name: Nyall Dawson (@nyalldawson)

---

If the difference is encountered when using the exact same gdal command then the change must be in gdal itself -- you'll need to file this on gdal's tracker.

---

• resolution was changed from to up/downstream
• status_id was changed from Open to Closed