[GRASS-user] Datum not recognized by Grass
Helmut Kudrnovsky
hellik at web.de
Tue Dec 3 11:51:59 PST 2019
Markus Neteler wrote
> Hi,
>
> On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei <
> zoltans at .co
> > wrote:
>>
>> Hi,
>> I'm using EPSG 2932 and QGIS etc all are OK with it.
>>
>> I have installed everything using OSGeo4W64, so how come Grass does not
>> use the same projections database (as QGis etc)?
>> (and please can I have a pointer as to how to introduce EPSG 2932 to
>> Grass 7.8.1)
>>
>> Thanks and regards,
>> Zoltan
>>
>> ******** Running against:
>> D:\GDBroad_tiles\shp_road_poly_few\23103770_road_poly.shp
>> WARNING: Datum
> <Qatar_National_Datum_1995>
> not recognised by GRASS and
>> no parameters found
>
> Is it possible that you have a PROJ software version mixup?
>
> I tried on my Linux box:
>
> grass78 -c epsg:2932 ~/grassdata/test_2932
> Starting GRASS GIS...
> Creating new GRASS GIS location
> <test_2932>
> ...
> Cleaning up temporary files...
>
> __________ ___ __________ _______________
> / ____/ __ \/ | / ___/ ___/ / ____/ _/ ___/
> / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \
> / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ /
> \____/_/ |_/_/ |_/____/____/ \____/___//____/
>
> Welcome to GRASS GIS 7.8.2dev (ad4836c73)
> ...
>
> GRASS 7.8.2dev (test_2932):~ > g.proj -w
> PROJCS["QND95 / Qatar National Grid",
> GEOGCS["QND95",
> DATUM["Qatar_National_Datum_1995",
> SPHEROID["International 1924",6378388,297,
> AUTHORITY["EPSG","7022"]],
>
> TOWGS84[-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065],
> AUTHORITY["EPSG","6614"]],
> PRIMEM["Greenwich",0,
> AUTHORITY["EPSG","8901"]],
> UNIT["degree",0.0174532925199433,
> AUTHORITY["EPSG","9122"]],
> AUTHORITY["EPSG","4614"]],
> PROJECTION["Transverse_Mercator"],
> PARAMETER["latitude_of_origin",24.45],
> PARAMETER["central_meridian",51.21666666666667],
> PARAMETER["scale_factor",0.99999],
> PARAMETER["false_easting",200000],
> PARAMETER["false_northing",300000],
> UNIT["metre",1,
> AUTHORITY["EPSG","9001"]],
> AXIS["Easting",EAST],
> AXIS["Northing",NORTH],
> AUTHORITY["EPSG","2932"]]
when importing a shapefile with the same epsg projection, I get:
v.in.ogr --verbose input=D:\wd\test_epsg2932.shp
Using OGR driver 'ESRI Shapefile/ESRI Shapefile'
WARNING: Datum <Qatar_National_Datum_1995> von GRASS nicht erkannt und keine
Parameter gefunden.
Die Projektionsinformationen des Eingabedatensatzes und der aktuellen
Location scheinen übereinzustimmen.
D:\wd>ogrinfo --version
GDAL 3.0.2, released 2019/10/28
D:\wd>ogrinfo test_epsg2932.shp -al -so
INFO: Open of `test_epsg2932.shp'
using driver `ESRI Shapefile' successful.
Layer name: test_epsg2932
Metadata:
DBF_DATE_LAST_UPDATE=2019-12-03
Geometry: Point
Feature Count: 1
Extent: (823205.438352, 25769.568359) - (823205.438352, 25769.568359)
Layer SRS WKT:
PROJCS["QND95 / Qatar National Grid",
GEOGCS["QND95",
DATUM["Qatar_National_Datum_1995",
SPHEROID["International 1924",6378388,297,
AUTHORITY["EPSG","7022"]],
AUTHORITY["EPSG","6614"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4614"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",24.45],
PARAMETER["central_meridian",51.2166666666667],
PARAMETER["scale_factor",0.99999],
PARAMETER["false_easting",200000],
PARAMETER["false_northing",300000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","2932"]]
Data axis to CRS axis mapping: 1,2
id: Integer64 (10.0)
test: Integer (9.0)
D:\wd>testepsg test_epsg2932.prj
Validate Succeeds.
WKT[test_epsg2932.prj] =
PROJCS["QND95 / Qatar National Grid",
GEOGCS["QND95",
DATUM["Qatar_National_Datum_1995",
SPHEROID["International 1924",6378388,297,
AUTHORITY["EPSG","7022"]],
AUTHORITY["EPSG","6614"]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",24.45],
PARAMETER["central_meridian",51.2166666666667],
PARAMETER["scale_factor",0.99999],
PARAMETER["false_easting",200000],
PARAMETER["false_northing",300000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH]]
Simplified WKT[test_epsg2932.prj] =
PROJCS["QND95 / Qatar National Grid",
GEOGCS["QND95",
DATUM["Qatar_National_Datum_1995",
SPHEROID["International 1924",6378388,297]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",24.45],
PARAMETER["central_meridian",51.2166666666667],
PARAMETER["scale_factor",0.99999],
PARAMETER["false_easting",200000],
PARAMETER["false_northing",300000],
UNIT["metre",1]]
Old Style WKT[test_epsg2932.prj] = PROJCS["QND95 / Qatar National
Grid",GEOGCS["QND95",DATUM["Qatar_National_Datum_1995",SPHEROID["International
1924",6378388,297]],PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",24.45],PARAMETER["central_meridian",51.2166666666667],PARAMETER["scale_factor",0.99999],PARAMETER["false_easting",200000],PARAMETER["false_northing",300000],UNIT["metre",1]]
ESRI'ified WKT[test_epsg2932.prj] =
PROJCS["QND_1995_Qatar_National_Grid",
GEOGCS["GCS_QND_1995",
DATUM["D_QND_1995",
SPHEROID["International_1924",6378388.0,297.0]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["False_Easting",200000.0],
PARAMETER["False_Northing",300000.0],
PARAMETER["Central_Meridian",51.2166666666667],
PARAMETER["Scale_Factor",0.99999],
PARAMETER["Latitude_Of_Origin",24.45],
UNIT["Meter",1.0]]
PROJ.4 rendering of [test_epsg2932.prj] = +proj=tmerc +lat_0=24.45
+lon_0=51.2166666666667 +k=0.99999 +x_0=200000 +y_0=300000 +ellps=intl
+units=m +no_defs
maybe worth to open a ticket
-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html
More information about the grass-user
mailing list