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.man...@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 _______________________________________________ Qgis-developer mailing list Qgis-developer@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-developer