Мне нужно перепроецировать файл .tiff из EPSG: 4326 - WGS 84 - Geographi c в EPSG: 32638 (UTM Zone 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")