[gdal-dev] Results of gdalwarp from UTM32 to EPSG:4326 are not aligned?

[Pham Huu Bang]

Hello,

I've a Sentinel 2 UTM32 scene (JPEG2000 format).

What I want to do is:

+ Warp this scene from UTM32 to EPSG:4326 by gdalwarp
gdalwarp T32UQD_20200421T102021_TCI_60m.jp2 -t_srs EPSG:4326 full_warp.tif

+ Clip a small subset of this file by gdal_translate
gdal_translate -projwin 798326.68347 5795956.80269 799609.577326
5794397.95642 T32UQD_20200421T102021_TCI_60m.jp2 crop.tif

Then warp this small subset also from UTM32 to EPSG:4326 by gdalwarp:

gdalwarp crop.tif -t_srs EPSG:4326 small_warp.tif -dstnodata 0

However, when I overlay the small_warp.tif on the full_warp.tif by
QGIS 2.14.9, I see the small file is not aligned on the big one.

Illustration: https://i.stack.imgur.com/U5Vql.png

What can be the reason for it?

More details given on this link:
https://gis.stackexchange.com/questions/361410/results-of-gdalwarp-from-utm32-to-epsg4326-are-not-aligned

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