title: Reprojection of .tiff file using gdal
body: <p>I need to re-project a .tiff file from EPSG:4326 - WGS 84 - Geographic to EPSG: 32638 (UTM Zone 38N). When I run the code below the output file is a re-projected .tiff half the size of the original one.</p>
  
  <pre><code># 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")
  
  </code></pre>