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 > [email protected] > http://lists.osgeo.org/mailman/listinfo/grass-user _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
