[GRASS-user] importing raster from a linear projection into current location

[Sebastian P. Luque]

On Fri, 1 Mar 2013 22:42:56 +0100,
Markus Neteler <neteler at osgeo.org> wrote:

> On Fri, Mar 1, 2013 at 9:53 PM, Seb <spluque at gmail.com> 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