[GRASS5] Re: [Gdal-dev] ERDAS/IMG driver error for NULL values

[Otto Dassau]

On Wed, 30 Jun 2004 15:02:55 +0200
Markus Neteler <neteler at itc.it> wrote:

> On Wed, Jun 30, 2004 at 11:14:45AM +0100, Paul Kelly wrote:
> > Hello Otto
> > 
> > Otto Dassau wrote:
> > > 
> > > 
> > > Also there is no coordinate system specified.
> 
> Otto, are you using the original GRASS driver or the new driver
> of Radim Blazek (not in CVS yet)?

I compiled GDAL with the patches from Radim for new raster support also in QGIS. I didn't use "old" libgrass.  The patches were:

grassdataset.cpp3
configure.in.patch2


> > This should be:
> > > 
> > > > name: Transverse Mercator
> > > > datum: potsdam
> > > > towgs84: 590.5,69.5,411.6,-0.796,-0.052,-3.601,8.30
> > > > proj: tmerc
> > > > ellps: bessel
> > > > a: 6377397.1550000003
> > > > es: 0.0066743722
> > > > f: 299.1528128000
> > > > lat_0: 0.0000000000
> > > > lon_0: 9.0000000000
> > > > k_0: 1.0000000000
> > > > x_0: 3500000.0000000000
> > > > y_0: 0.0000000000
> > 
> > This looks like a deficiency in the r.out.gdal script.
> 
> I don't think so. Because the r.out.gdal script is simply a wrapper
> around gdal_translate.
> 
> In fact gdalinfo doesn't print the projection:
> 
> gdalinfo ~/grassdata/spearfish57/PERMANENT/cellhd/elevation.dem | head
> Driver: GRASS/GRASS Database Rasters
> Size is 633, 466
> Coordinate System is `'
> Origin = (590010.000000,4928000.000000)
> Pixel Size = (30.00000000,-30.00000000)
> Corner Coordinates:
> Upper Left  (  590010.000, 4928000.000)
> Lower Left  (  590010.000, 4914020.000)
> Upper Right (  609000.000, 4928000.000)
> Lower Right (  609000.000, 4914020.000)
> [...]
> 
> I am using the new driver which contains the relevant code (I'm no
> expert). Somewhere this info get's lost. Or should the driver make
> use of the new GRASS 5.7 projection functions (such as GPJ_grass_to_wkt
> in grass57/lib/proj/convert.c )?
> 
> Once solved in GDAL also r.out.gdal will write out the projection.
> 
> Markus
> 
> _______________________________________________
> grass5 mailing list
> grass5 at grass.itc.it
> http://grass.itc.it/mailman/listinfo/grass5
>