[gdal-dev] Difference results from QGIS project on fly and gdalwarp for EPSG:4326 of Sentinel 2 UTM zone 32?

[Pham Huu Bang]

Hello,

I always thought QGIS and gdalwarp should return same result for projecting
a 2D Sentinel 2 image (e.g: T32UQD_20200421T102021_TCI_60m.jp2) from UTM 32
to EPSG:4326.

What I did was:

   - Open the .jp2 in QGis 2.14.9 with 'Enable on the fly CRS
   transformation'.
   -

   Warp the .jp2 to geotiff by gdalwarp:

   gdalwarp T32UQD_20200421T102021_TCI_60m.jp2 -t_srs EPSG:4326
   4.warped.tiff

But the output results from them look different significantly (on the left
hand side with QGIS it looks rotated and on the right hand side with
gdalwarp it looks aligned).

Illustration image of differences: https://i.stack.imgur.com/Idisx.png

More information given described at:
https://gis.stackexchange.com/q/360919/17733

Could someone show the opinion?

Thanks,
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20200508/f1359961/attachment.html>