[gdal-dev] gdal_grid inverts y axis

Christopher Dedels cdedels at gmail.com
Sun Apr 17 16:10:57 EDT 2011


I am using gdal_grid to create a grid from a point shapefile.  gdalinfo
reports inverted y-axis coordinates for the output raster.  The command I am
running is:

gdal_grid -ot Float32 -outsize 800 600 -l sites -zfield elev sites.shp
/tmp/sites.tif

I am using gdal 1.8.0.  the output of gdalinfo for sites.tif is:

Driver: GTiff/GeoTIFF
Files: /tmp/sites.tif
Size is 800, 600
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin = (-123.860031127930000,38.011638641357401)
Pixel Size = (0.033339290618897,0.026156539916992)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (-123.8600311,  38.0116386) (123d51'36.11"W, 38d 0'41.90"N)
Lower Left  (-123.8600311,  53.7055626) (123d51'36.11"W, 53d42'20.03"N)
Upper Right ( -97.1885986,  38.0116386) ( 97d11'18.96"W, 38d 0'41.90"N)
Lower Right ( -97.1885986,  53.7055626) ( 97d11'18.96"W, 53d42'20.03"N)
Center      (-110.5243149,  45.8586006) (110d31'27.53"W, 45d51'30.96"N)
Band 1 Block=800x2 Type=Float32, ColorInterp=Gray

Is there any way to invert the y-axis back?  Thanks for the help.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20110417/b07d414b/attachment.html


More information about the gdal-dev mailing list