[Gdal-dev] Re: problem with gdalwarp (resulting images are not georeferenced properly)

[rumicuna4-2007 at yahoo.com]

Some people reported the links to the images were not working (because of a %20 that was added at the end), here are the proper links to the images:
http://ecovagabond.com/gis/Denmark.InitialCrop.tif
http://ecovagabond.com/gis/Denmark.Warped.Gdal.tif
Thanks,
Diego

rumicuna4-2007 at yahoo.com wrote: Hello all,
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.

The command I am using is:

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

I have placed the image I am trying to warp and the resulting file here:

http://ecovagabond.com/gis/Denmark.InitialCrop.tif 
http://ecovagabond.com/gis/Denmark.Warped.Gdal.tif 

In case  this is helpful I am also including below the results of running gdalinfo on both files.
I am using a freshly installed version of FWTools1.3.0 on a WindowsXP machine.
The initial file was created by running gdal_translate on an image of the world.
Any help would be greatly appreciated.
 
Thank you,
-Diego


gdalinfo.exe Denmark.InitialCrop.tif

Driver: GTiff/GeoTIFF
Size is 988, 433
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.2572235630016,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin =  (7.520833333333343,57.937500000000000)
Pixel Size = (0.008333333333333,-0.008333333333333)
Metadata:
  AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  (   7.5208333,  57.9375000) (  7d31'15.00"E, 57d56'15.00"N)
Lower Left  (   7.5208333,  54.3291667) (  7d31'15.00"E, 54d19'45.00"N)
Upper Right (  15.7541667,  57.9375000) ( 15d45'15.00"E, 57d56'15.00"N)
Lower Right (  15.7541667,  54.3291667) ( 15d45'15.00"E, 54d19'45.00"N)
Center      (  11.6375000,  56.1333333) ( 11d38'15.00"E, 56d 8'0.00"N)
Band 1 Block=988x8 Type=Byte, ColorInterp=Red
Band 2 Block=988x8 Type=Byte, ColorInterp=Green
Band 3 Block=988x8 Type=Byte, ColorInterp=Blue

gdalinfo.exe Denmark.Warped.Gdal.tif

Driver: GTiff/GeoTIFF
Size is 888, 679
Coordinate System  is:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235630016,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Orthographic"],
    PARAMETER["latitude_of_origin",56.154163],
    PARAMETER["central_meridian",11.61861],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
     UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (-266952.302792707110000,186552.524760038390000)
Pixel Size = (603.704336812790870,-603.704336812790870)
Metadata:
  AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  ( -266952.303,  186552.525) (  7d 7'11.01"E, 57d55'29.90"N)
Lower Left  ( -266952.303, -223362.720) (  7d31'40.14"E, 54d15'29.98"N)
Upper Right (  269137.148,  186552.525) ( 16d 9'15.23"E, 57d55'25.29"N)
Lower Right (  269137.148, -223362.720) ( 15d44'34.17"E, 54d15'25.78"N)
Center      (    1092.423,  -18405.098) ( 11d38'10.15"E, 56d 9'59.87"N)
Band 1 Block=888x3 Type=Byte, ColorInterp=Red
Band 2 Block=888x3 Type=Byte, ColorInterp=Green
Band 3 Block=888x3 Type=Byte, ColorInterp=Blue






-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20070504/02199df2/attachment.html