[Qgis-developer] Raster providers

[Jürgen E. Fischer]

On Thu, 31. Mar 2011 at 11:12:42 +0200, Radim Blazek wrote:
> On Wed, Mar 9, 2011 at 12:24 PM, Giovanni Manghi
> <giovanni.manghi at gmail.com> wrote:
> > Regardless their CRS, GRASS rasters are always added to qgis with WGS84 CRS.
> 
> For me, QGIS sets CRS to raster's CRS but some info seems to be lost
> in createFromWkt().
> 
> GRASS gives: GEOGCS["wgs84",DATUM["WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563],TOWGS84[0.000,0.000,0.000]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]]
> 
> QgsCoordinateReferenceSystem.createFromWkt().toWkt() gives:
> GEOGCS["WGS 84",DATUM["unknown",SPHEROID["WGS84",6378137,298.257223563],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]]

> QgsCoordinateReferenceSystem.authid(): EPSG:4148
> QgsCoordinateReferenceSystem.description(): Hartebeesthoek94

> Everything for CRS QgsCoordinateReferenceSystem created from WKT given
> by GRASS above.

That shouldn't be a problem projectionwise.

It found EPSG:4148 because that has the same projection parameters as EPSG:4326
and the WKT string doesn't have any information about which EPSG is desired.
Otherwise the description and authid would change, still that wouldn't change
the projection.

Or am I missing something else that is lost?


Jürgen
 
-- 
Jürgen E. Fischer         norBIT GmbH               Tel. +49-4931-918175-20
Dipl.-Inf. (FH)           Rheinstraße 13            Fax. +49-4931-918175-50
Software Engineer         D-26506 Norden               http://www.norbit.de

-- 
norBIT Gesellschaft fuer Unternehmensberatung und Informationssysteme mbH
Rheinstrasse 13, 26506 Norden
GF: Jelto Buurman, HR: Amtsgericht Emden, HRB 5502