On Fri, 1 Mar 2013 22:42:56 +0100, Markus Neteler <[email protected]> wrote:
> On Fri, Mar 1, 2013 at 9:53 PM, Seb <[email protected]> wrote: >> Hi, >> I am trying to import a series of rasters (HDF) in a linear >> projection, i.e. without any projection information, where x and y >> are simple indexes; gdalinfo shows: > Could you say a bit more about the HDF files? The xy output is > probably right but GCPs may be in the file itself (so gdalwarp is > needed or the like). This file was produced with WIM/WAM (wimsoft.com), and the full output of gdalinfo is: $ gdalinfo chl_48_200.hdf Driver: HDF4Image/HDF4 Dataset Files: chl_48_200.hdf Size is 1572, 1114 Coordinate System is `' Metadata: Base=10 Intercept=-2 Map Projection=WIM Linear Scaling=Logarithmic Slope=0.015 WIM Color Range=0, 255 WIM Color Stretch=48, 200 WIM Linear Coeffs=45, 0.012723, -40, -0.008977 Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( 0.0, 1114.0) Upper Right ( 1572.0, 0.0) Lower Right ( 1572.0, 1114.0) Center ( 786.0, 557.0) Band 1 Block=1572x636 Type=Byte, ColorInterp=Gray The relevant information is in the "WIM Linear Coeffs" line, where 45, -40 are the coordinates of the upper right corner of the map, and 0.012723, -0.008977 are the lon, lat resolution. Thanks, -- Seb _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
