[GRASS-user] strange ETRS89 settings and v.out.ogr

[Nikos Alexandris]

On Fri, 2008-08-15 at 22:00 +0100, Paul Kelly wrote:
> On Fri, 15 Aug 2008, Tim Michelsen wrote:
> 
> > Hello,
> > I defined a test location in ETRS89-LAEA using the EPSG number 3035 in the GRASS
> > startup wizard.
> >
> > Here is what g.proj -p > test_3035_proj.txt prints out:
> > -PROJ_INFO-------------------------------------------------
> > name       : Lambert Azimuthal Equal Area
> > proj       : laea
> > datum      : etrs89
> > ellps      : grs80
> > lat_0      : 52
> > lon_0      : 10
> > x_0        : 4321000
> > y_0        : 3210000
> > no_defs    : defined
> > -PROJ_UNITS------------------------------------------------
> > unit       : metre
> > units      : metres
> > meters     : 1
> >
> > Now I exported a data set to a shapefile using v.out.ogr. Here is the content of
> > the associated *.prj file:
> >
> > Layer SRS WKT:
> > PROJCS["Lambert Azimuthal Equal Area",
> >    GEOGCS["grs80",
> >        DATUM["ETRF_1989",
> 
> This datum name is very unusual and I'm not sure where it could have come 
> from. What version of GRASS? What is the exact commandline you used?

Paul, FWIW I confirm that I get the same output when I create a location
with epsg=3035 on latest grass-6.4.svn
> 
> Can you try the following commands and post output?
> g.proj -w
> g.proj -we
> g.proj -w epsg=3035
> 
> >            SPHEROID["Geodetic_Reference_System_1980",6378137,298.257222101]],
> >        PRIMEM["Greenwich",0],
> >        UNIT["Degree",0.017453292519943295]],
> >    PROJECTION["Lambert_Azimuthal_Equal_Area"],
> >    PARAMETER["latitude_of_center",52],
> >    PARAMETER["longitude_of_center",10],
> >    PARAMETER["false_easting",4321000],
> >    PARAMETER["false_northing",3210000],
> >    UNIT["Meter",1]]
> > cat: Real (11.0)
> > a_cat: Real (11.0)
> > a_CODE_00: String (80.0)
> > b_cat: Real (11.0)
> > b_cat_: Real (24.15)
> > b_na2: String (80.0)
And the same after exporting with v.out.ogr a point file created with
v.random.

That's all from my side (no access to ESRI stuff).

> >
> > If I load this into ArcGIS it does complain about wrong projection information.
> 
> How exactly does it complain?
> 
> > Here is the project file when reprojecting the data set to 3035 with ogr2ogr:
> > Layer SRS WKT:
> > PROJCS["ETRS89 / ETRS-LAEA",
> >    GEOGCS["ETRS89",
> >        DATUM["European_Terrestrial_Reference_System_1989",
> 
> This is the datum name I would expect to see from GRASS too.
> 
> >            SPHEROID["GRS_1980",6378137,298.257222101]],
> >        PRIMEM["Greenwich",0],
> >        UNIT["Degree",0.017453292519943295]],
> >    PROJECTION["Lambert_Azimuthal_Equal_Area"],
> >    PARAMETER["latitude_of_center",52],
> >    PARAMETER["longitude_of_center",10],
> >    PARAMETER["false_easting",4321000],
> >    PARAMETER["false_northing",3210000],
> >    UNIT["Meter",1]]
> > cat: Real (11.0)
> > a_cat: Real (11.0)
> > a_CODE_00: String (80.0)
> > b_cat: Real (11.0)
> > b_cat_: Real (24.15)
> > b_na2: String (80.0)
> >
> >
> > Is the GRASS projection file wrong or is it not "standard"? From the content
> > they both seem to tell the same...
> 
> The datum name looks strange. Note that you can create a .prj file in 
> "ESRI format" with
> g.proj -wef > xxx.prj
> to replace the .prj file created by v.out.ogr, if it's any use.
> 
> Paul
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user