[GRASSLIST:1889] Re: r.in.gdal

[Frank Warmerdam]

Markus Neteler wrote:
> > markus,
> > g.setproj doesn`t work in  xy-locations.
> > -> XY-location cannot be projected
> > my tifs should be in GK-location.
> > is r.in.gdal not able to create GK-locations ?
> > regards
> > steph
> 
> Hi Stephan,
> 
> o.k. I didn't think of that. So another try:
> Generate a .prj file for your data and generate a new
> location with r.in.gdal which *should* accept this .prj file.
> Here a sample for Gauss-Krueger (hope it is correct):
> 
> PROJECTION TRANSVERSE
> UNITS METERS
> DATUM POTSDAM
> PARAMETERS
> 1.000
> 9 0 0
> 0 0 0
> 3500000
> 0.0

Guys,

While this might solve this problem, it seems that a better approach is
to setup the output location in advance, and then run r.in.gdal with the
-o flag to override it's coordinate system checks.  Then r.in.gdal should
preserve the original coordinates from the .tfw while the user has to
setup the projection in advance with g.setproj since the files don't contain
any projection information. 

Best regards,

---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent