Hello all,<br>I am getting incorrect results when I use gdalwarp and cannot figure out what I am doing wrong. The warping looks perfect, but the coords embedded in the resulting geotiff are shifted south so the resulting image does not show up correctly when overlayed on a map. I tried with other programs and did not have the same problem, I would rather use gdal because of the great interpolation capabilities.<br><br>The command I am using is:<br><br>gdalwarp.exe -rc -s_srs "+proj=latlong +datum=WGS84" -t_srs "+proj=ortho +datum=WGS84 +lon_0=11.618610 +lat_0=56.154163" -co INTERLEAVE=PIXEL Denmark.InitialCrop.tif Denmark.Warped.tif<br><br>I have placed the image I am trying to warp and the resulting file here:<br><br><a href="http://ecovagabond.com/gis/Denmark.InitialCrop.tif%20">http://ecovagabond.com/gis/Denmark.InitialCrop.tif </a><br><a href="http://ecovagabond.com/gis/Denmark.Warped.Gdal.tif%20">http://ecovagabond.com/gis/Denmark.Warped.Gdal.tif </a><br><br>In case
this is helpful I am also including below the results of running gdalinfo on both files.<br>I am using a freshly installed version of FWTools1.3.0 on a WindowsXP machine.<br>The initial file was created by running gdal_translate on an image of the world.<br>Any help would be greatly appreciated.<br> <br>Thank you,<br>-Diego<br><br><br><span style="font-weight: bold;">gdalinfo.exe Denmark.InitialCrop.tif</span><br><br>Driver: GTiff/GeoTIFF<br>Size is 988, 433<br>Coordinate System is:<br>GEOGCS["WGS 84",<br> DATUM["WGS_1984",<br> SPHEROID["WGS 84",6378137,298.2572235630016,<br> AUTHORITY["EPSG","7030"]],<br> AUTHORITY["EPSG","6326"]],<br> PRIMEM["Greenwich",0],<br> UNIT["degree",0.0174532925199433],<br> AUTHORITY["EPSG","4326"]]<br>Origin =
(7.520833333333343,57.937500000000000)<br>Pixel Size = (0.008333333333333,-0.008333333333333)<br>Metadata:<br> AREA_OR_POINT=Area<br>Corner Coordinates:<br>Upper Left ( 7.5208333, 57.9375000) ( 7d31'15.00"E, 57d56'15.00"N)<br>Lower Left ( 7.5208333, 54.3291667) ( 7d31'15.00"E, 54d19'45.00"N)<br>Upper Right ( 15.7541667, 57.9375000) ( 15d45'15.00"E, 57d56'15.00"N)<br>Lower Right ( 15.7541667, 54.3291667) ( 15d45'15.00"E, 54d19'45.00"N)<br>Center ( 11.6375000, 56.1333333) ( 11d38'15.00"E, 56d 8'0.00"N)<br>Band 1 Block=988x8 Type=Byte, ColorInterp=Red<br>Band 2 Block=988x8 Type=Byte, ColorInterp=Green<br>Band 3 Block=988x8 Type=Byte, ColorInterp=Blue<br><br><span style="font-weight: bold;">gdalinfo.exe Denmark.Warped.Gdal.tif</span><br><br>Driver: GTiff/GeoTIFF<br>Size is 888, 679<br>Coordinate System
is:<br>PROJCS["unnamed",<br> GEOGCS["WGS 84",<br> DATUM["WGS_1984",<br> SPHEROID["WGS 84",6378137,298.2572235630016,<br> AUTHORITY["EPSG","7030"]],<br> AUTHORITY["EPSG","6326"]],<br> PRIMEM["Greenwich",0],<br> UNIT["degree",0.0174532925199433],<br> AUTHORITY["EPSG","4326"]],<br> PROJECTION["Orthographic"],<br> PARAMETER["latitude_of_origin",56.154163],<br> PARAMETER["central_meridian",11.61861],<br> PARAMETER["false_easting",0],<br> PARAMETER["false_northing",0],<br>
UNIT["metre",1,<br> AUTHORITY["EPSG","9001"]]]<br>Origin = (-266952.302792707110000,186552.524760038390000)<br>Pixel Size = (603.704336812790870,-603.704336812790870)<br>Metadata:<br> AREA_OR_POINT=Area<br>Corner Coordinates:<br>Upper Left ( -266952.303, 186552.525) ( 7d 7'11.01"E, 57d55'29.90"N)<br>Lower Left ( -266952.303, -223362.720) ( 7d31'40.14"E, 54d15'29.98"N)<br>Upper Right ( 269137.148, 186552.525) ( 16d 9'15.23"E, 57d55'25.29"N)<br>Lower Right ( 269137.148, -223362.720) ( 15d44'34.17"E, 54d15'25.78"N)<br>Center ( 1092.423, -18405.098) ( 11d38'10.15"E, 56d 9'59.87"N)<br>Band 1 Block=888x3 Type=Byte, ColorInterp=Red<br>Band 2 Block=888x3 Type=Byte, ColorInterp=Green<br>Band 3 Block=888x3 Type=Byte, ColorInterp=Blue<br><br><BR><BR>