title: georeferencing with gdal goes wrong
body: <p>although there are several pages dedicated to this task, it seems that i can not make that work. I have a non georeferenced tif file which I provide you the output from gdalinfo</p>
  
  <pre><code>Driver: GTiff/GeoTIFF
  Files: Arctic_r05c02.2017235.terra.250m.tif
  Size is 4096, 4096
  Coordinate System is `'
  Metadata:
    TIFFTAG_SOFTWARE=ppm2geotiff v0.0.9
  Image Structure Metadata:
    INTERLEAVE=PIXEL
  Corner Coordinates:
  Upper Left  (    0.0,    0.0)
  Lower Left  (    0.0, 4096.0)
  Upper Right ( 4096.0,    0.0)
  Lower Right ( 4096.0, 4096.0)
  Center      ( 2048.0, 2048.0)
  Band 1 Block=4096x1 Type=Byte, ColorInterp=Red
  Band 2 Block=4096x1 Type=Byte, ColorInterp=Green
  Band 3 Block=4096x1 Type=Byte, ColorInterp=Blue
  </code></pre>
  
  <p>I do have the corner coordinate extention in lat/lon. Hence I am creating Ground Control Points with thee following command:</p>
  
  <pre><code>gdal_translate -of GTiff -gcp 0 0 -180.0000 +63.1066 -gcp 0 4096 -161.5651 +68.5979 -gcp 4096 0 +161.5651 +68.5979 -gcp 4096 4096 -180.0000 +76.3728 input.tif output.tif
  </code></pre>
  
  <p>which results to:</p>
  
  <pre><code>Driver: GTiff/GeoTIFF
  Files: Arctic_r05c02.2017235.terra.250mproj.tif
  Size is 4096, 4096
  Coordinate System is `'
  GCP Projection = 
  GCP[  0]: Id=1, Info=
            (0,0) -&gt; (-180,63.1066,0)
  GCP[  1]: Id=2, Info=
            (0,4096) -&gt; (-161.5651,68.5979,0)
  GCP[  2]: Id=3, Info=
            (4096,0) -&gt; (161.5651,68.5979,0)
  GCP[  3]: Id=4, Info=
            (4096,4096) -&gt; (-180,76.3728,0)
  Metadata:
    TIFFTAG_SOFTWARE=ppm2geotiff v0.0.9
  Image Structure Metadata:
    INTERLEAVE=PIXEL
  Corner Coordinates:
  Upper Left  (    0.0,    0.0)
  Lower Left  (    0.0, 4096.0)
  Upper Right ( 4096.0,    0.0)
  Lower Right ( 4096.0, 4096.0)
  Center      ( 2048.0, 2048.0)
  Band 1 Block=4096x1 Type=Byte, ColorInterp=Red
  Band 2 Block=4096x1 Type=Byte, ColorInterp=Green
  Band 3 Block=4096x1 Type=Byte, ColorInterp=Blue
  </code></pre>
  
  <p>The next step is to resample with gdalwarp in order to create a projected file.</p>
  
  <pre><code>gdalwarp -r near -s_srs epsg:4326 -t_srs epsg:3413 -tr 250 250 input.tif output.tif
  </code></pre>
  
  <p>But the file seems distort:
  <a href="https://i.stack.imgur.com/Cp4X5.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/Cp4X5.png" alt="enter image description here"></a></p>
  
  <p>Any idea what I can further do? Is there a way to directly change the corner coordinates inside the tiff file instead of adding GCP's? Would that help?</p>
  
  <p>Some more information which i forgot. The data can be downloaded at:</p>
  
  <pre><code>https://lance-modis.eosdis.nasa.gov/imagery/subsets/?subset=Arctic_r05c02.2017235.terra.250m
  </code></pre>
  
  <p>and for the bounding box definition:</p>
  
  <pre><code>https://lance-modis.eosdis.nasa.gov/imagery/subsets/?subset=Arctic_r05c02.2017235.terra.250m.met
  </code></pre>
  
  <p>Since the box definition from NASA's side looks a bit odd, i tried the same procedure for the following tile:</p>
  
  <pre><code>https://lance-modis.eosdis.nasa.gov/imagery/subsets/?subset=Arctic_r04c03.2017235.terra.250m
  </code></pre>
  
  <p>and box definition:</p>
  
  <pre><code>https://lance-modis.eosdis.nasa.gov/imagery/subsets/?subset=Arctic_r04c03.2017235.terra.250m.met
  </code></pre>
  
  <p>Although the box looks ok the following image was produced:</p>
  
  <p><a href="https://i.stack.imgur.com/i6Jcc.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/i6Jcc.png" alt="enter image description here"></a></p>