Мне нужно перепроецировать файл .tiff из EPSG: 4326 - WGS 84 - Geographic в EPSG: 32638 (зона UTM 38N). Когда я запускаю код ниже, выходной файл представляет собой перепроецированный .tiff размером вдвое меньше исходного.
# Reprojection
from osgeo import gdal
filename = r"/home/DATA/Test.tif"
input_raster = gdal.Open(filename)
output_raster = r"/home/DATA/Test_UTM38N.tif"
gdal.Warp(output_raster,input_raster,dstSRS="+init=epsg:32638")
У меня был подобный опыт при использовании версии gdalwarp для командной строки. Указание разрешения x и y устранило проблему. То есть, в моем случае добавление «-tr 1000 1000» к команде gdalwarp привело к получению перепроецированной выходной сетки примерно того же размера, что и оригинал.
Также обратите внимание, что документ gdalwarp для опции '-tr' (https://gdal.org/programs/gdalwarp.html#cmdoption-gdalwarp-tr) говорит:
Если не указано (или не получено из -te и -ts), gdalwarp сгенерирует выходной растр с xres=yres, и это даже при использовании gdalwarp в сценариях, не связанных с повторным проецированием.
