Re: [GRASS-user] Datum not recognized by Grass
Within GRASS: g.proj epsg=2932 -p -PROJ_INFO- name : QND95 / Qatar National Grid ellps : international proj : tmerc lat_0 : 24.45 lon_0 : 51.21667 k : 0.9 x_0: 20 y_0: 30 towgs84: -119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065 no_defs: defined -PROJ_EPSG- epsg : 2932 -PROJ_UNITS unit : meter units : meters meters : 1 looks fine. Outside GRASS, using projinfo from PROJ 6.2.1: projinfo -o PROJ -s EPSG:2932 -t EPSG:4326 Candidate operations found: 1 - Operation n°1: unknown id, Inverse of Qatar National Grid + QND95 to WGS 84 (1), 0 m, Qatar - onshore PROJ string: +proj=pipeline +step +inv +proj=tmerc +lat_0=24.45 +lon_0=51.21667 +k=0.9 +x_0=20 +y_0=30 +ellps=intl +step +proj=push +v_3 +step +proj=cart +ellps=intl +step +proj=helmert +x=-119.4248 +y=-303.65872 +z=-11.00061 +rx=1.164298 +ry=0.174458 +rz=1.096259 +s=3.657065 +convention=position_vector +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1 Datum transformation happens with +proj=helmert which is the same as the towgs84 parameters in GRASS. This does not matter, GRASS does not use these parameters for reprojection, it relies on PROJ to select an appropriate datum transformation depending on the source and target CRS. Markus M On Mon, Dec 9, 2019 at 10:15 PM Markus Metz wrote: > > > On Wed, Dec 4, 2019 at 12:23 PM Zoltan Szecsei > wrote: > > > > Hi Helmut, > > Thanks for your comments. > > > > I installed everything with OSGeo4W64, and QGIS get the EPSG:2932 but > > Grass not. > > > > Perhaps I have a PATH or some other setting problem? > > Perhaps let me know what search paths Grass uses for proj4, and what > > proj files and locations I should scan for. > > In this case where the EPSG code is known, there is no need to do anything > but to ignore the warning. GRASS will use the EPSG code if available and > passes it to PROJ when it comes to reprojection. > > Markus M > > > > > Regards, > > Zoltan > > > > On 2019/12/04 01:30, Helmut Kudrnovsky wrote: > > >> ignore the warning and use GRASS with PROJ6, granted that authority > name > > > (e.g. EPSG) and authority code (e.g. 2932) are known for both CRS's in > case > > > of reprojection > > > > > > this issue is already by GRASS with PROJ6, see > > > > > > > > > GRASS version: 7.8.1 > > > Code revision: c865432c9 > > > Build date: 2019-11-10 > > > Build platform: x86_64-w64-mingw32 > > > GDAL: 3.0.2 > > > PROJ: 6.2.1 <= > > > GEOS: 3.8.0 > > > SQLite: 3.29.0 > > > Python: 3.7.0 > > > wxPython: 4.0.7 > > > Platform: Windows-10-10.0.18362-SP0 (OSGeo4W) > > > > > > > > > and the output from the underlying PROJ: 6.2.1 > > > > > > > > > C:\>projinfo EPSG:2932 -o PROJ,WKT2_2018 > > > PROJ.4 string: > > > +proj=tmerc +lat_0=24.45 +lon_0=51.21667 +k=0.9 +x_0=20 > > > +y_0=30 +ellps=intl > > > > +towgs84=-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065 > > > +units=m +no_defs +type=crs > > > > > > WKT2_2018 string: > > > PROJCRS["QND95 / Qatar National Grid", > > > BASEGEOGCRS["QND95", > > > DATUM["Qatar National Datum 1995", > > > ELLIPSOID["International 1924",6378388,297, > > > LENGTHUNIT["metre",1]]], > > > PRIMEM["Greenwich",0, > > > ANGLEUNIT["degree",0.0174532925199433]], > > > ID["EPSG",4614]], > > > CONVERSION["Qatar National Grid", > > > METHOD["Transverse Mercator", > > > ID["EPSG",9807]], > > > PARAMETER["Latitude of natural origin",24.45, > > > ANGLEUNIT["degree",0.0174532925199433], > > > ID["EPSG",8801]], > > > PARAMETER["Longitude of natural origin",51.21667, > > > ANGLEUNIT["degree",0.0174532925199433], > > > ID["EPSG",8802]], > > > PARAMETER["Scale factor at natural origin",0.9, > > > SCALEUNIT["unity",1], > > > ID["EPSG",8805]], > > > PARAMETER["False easting",20, > > > LENGTHUNIT["metre",1], > > > ID["EPSG",8806]], > > > PARAMETER["False northing",30, > > > LENGTHUNIT["metre",1], > > > ID["EPSG",8807]]], > > > CS[Cartesian,2], > > > AXIS["(E)",east, > > > ORDER[1], > > > LENGTHUNIT["metre",1]], > > > AXIS["(N)",north, > > > ORDER[2], > > > LENGTHUNIT["metre",1]], > > > USAGE[ > > > SCOPE["unknown"], > > > AREA["Qatar - onshore"], > > > BBOX[24.55,50.69,26.2,51.68]], > > >
Re: [GRASS-user] Datum not recognized by Grass
On Wed, Dec 4, 2019 at 12:23 PM Zoltan Szecsei wrote: > > Hi Helmut, > Thanks for your comments. > > I installed everything with OSGeo4W64, and QGIS get the EPSG:2932 but > Grass not. > > Perhaps I have a PATH or some other setting problem? > Perhaps let me know what search paths Grass uses for proj4, and what > proj files and locations I should scan for. In this case where the EPSG code is known, there is no need to do anything but to ignore the warning. GRASS will use the EPSG code if available and passes it to PROJ when it comes to reprojection. Markus M > > Regards, > Zoltan > > On 2019/12/04 01:30, Helmut Kudrnovsky wrote: > >> ignore the warning and use GRASS with PROJ6, granted that authority name > > (e.g. EPSG) and authority code (e.g. 2932) are known for both CRS's in case > > of reprojection > > > > this issue is already by GRASS with PROJ6, see > > > > > > GRASS version: 7.8.1 > > Code revision: c865432c9 > > Build date: 2019-11-10 > > Build platform: x86_64-w64-mingw32 > > GDAL: 3.0.2 > > PROJ: 6.2.1 <= > > GEOS: 3.8.0 > > SQLite: 3.29.0 > > Python: 3.7.0 > > wxPython: 4.0.7 > > Platform: Windows-10-10.0.18362-SP0 (OSGeo4W) > > > > > > and the output from the underlying PROJ: 6.2.1 > > > > > > C:\>projinfo EPSG:2932 -o PROJ,WKT2_2018 > > PROJ.4 string: > > +proj=tmerc +lat_0=24.45 +lon_0=51.21667 +k=0.9 +x_0=20 > > +y_0=30 +ellps=intl > > +towgs84=-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065 > > +units=m +no_defs +type=crs > > > > WKT2_2018 string: > > PROJCRS["QND95 / Qatar National Grid", > > BASEGEOGCRS["QND95", > > DATUM["Qatar National Datum 1995", > > ELLIPSOID["International 1924",6378388,297, > > LENGTHUNIT["metre",1]]], > > PRIMEM["Greenwich",0, > > ANGLEUNIT["degree",0.0174532925199433]], > > ID["EPSG",4614]], > > CONVERSION["Qatar National Grid", > > METHOD["Transverse Mercator", > > ID["EPSG",9807]], > > PARAMETER["Latitude of natural origin",24.45, > > ANGLEUNIT["degree",0.0174532925199433], > > ID["EPSG",8801]], > > PARAMETER["Longitude of natural origin",51.21667, > > ANGLEUNIT["degree",0.0174532925199433], > > ID["EPSG",8802]], > > PARAMETER["Scale factor at natural origin",0.9, > > SCALEUNIT["unity",1], > > ID["EPSG",8805]], > > PARAMETER["False easting",20, > > LENGTHUNIT["metre",1], > > ID["EPSG",8806]], > > PARAMETER["False northing",30, > > LENGTHUNIT["metre",1], > > ID["EPSG",8807]]], > > CS[Cartesian,2], > > AXIS["(E)",east, > > ORDER[1], > > LENGTHUNIT["metre",1]], > > AXIS["(N)",north, > > ORDER[2], > > LENGTHUNIT["metre",1]], > > USAGE[ > > SCOPE["unknown"], > > AREA["Qatar - onshore"], > > BBOX[24.55,50.69,26.2,51.68]], > > ID["EPSG",2932]] > > > > > >> The problem is that GRASS still assumes datum transformation from X to > > WGS84, whereas with PROJ6, > >> WGS84 is no longer required as pivot datum. Datum transformations from > > datum X to datum Y can > >> sometimes (often) be done without going through WGS84. The requirement is > > to have EPSG or other > >> authority codes and to use PROJ 6. > > this warning is just from a simple v.in.ogr of a EPSG:2932-based shapefile > > into a EPSG:2932-created GRASS location. > > > > so, there may be some more adaptions needed for actions like v.in.ogr and > > co? > > > > > > > > - > > best regards > > Helmut > > -- > > Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html > > ___ > > grass-user mailing list > > grass-user@lists.osgeo.org > > https://lists.osgeo.org/mailman/listinfo/grass-user > > -- > > = > Zoltan Szecsei GPrGISc 0031 > Geograph (Pty) Ltd. > GIS and Photogrammetric Services > > P.O. Box 7, Muizenberg 7950, South Africa. > > Mobile: +27-83-6004028 (WhatsApp only) > Qatar: +974 5083 2722 www.geograph.co.za > = > > ___ > grass-user mailing list > grass-user@lists.osgeo.org > https://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Datum not recognized by Grass
Hi Helmut, Thanks for your comments. I installed everything with OSGeo4W64, and QGIS get the EPSG:2932 but Grass not. Perhaps I have a PATH or some other setting problem? Perhaps let me know what search paths Grass uses for proj4, and what proj files and locations I should scan for. Regards, Zoltan On 2019/12/04 01:30, Helmut Kudrnovsky wrote: ignore the warning and use GRASS with PROJ6, granted that authority name (e.g. EPSG) and authority code (e.g. 2932) are known for both CRS's in case of reprojection this issue is already by GRASS with PROJ6, see GRASS version: 7.8.1 Code revision: c865432c9 Build date: 2019-11-10 Build platform: x86_64-w64-mingw32 GDAL: 3.0.2 PROJ: 6.2.1 <= GEOS: 3.8.0 SQLite: 3.29.0 Python: 3.7.0 wxPython: 4.0.7 Platform: Windows-10-10.0.18362-SP0 (OSGeo4W) and the output from the underlying PROJ: 6.2.1 C:\>projinfo EPSG:2932 -o PROJ,WKT2_2018 PROJ.4 string: +proj=tmerc +lat_0=24.45 +lon_0=51.21667 +k=0.9 +x_0=20 +y_0=30 +ellps=intl +towgs84=-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065 +units=m +no_defs +type=crs WKT2_2018 string: PROJCRS["QND95 / Qatar National Grid", BASEGEOGCRS["QND95", DATUM["Qatar National Datum 1995", ELLIPSOID["International 1924",6378388,297, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4614]], CONVERSION["Qatar National Grid", METHOD["Transverse Mercator", ID["EPSG",9807]], PARAMETER["Latitude of natural origin",24.45, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8801]], PARAMETER["Longitude of natural origin",51.21667, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8802]], PARAMETER["Scale factor at natural origin",0.9, SCALEUNIT["unity",1], ID["EPSG",8805]], PARAMETER["False easting",20, LENGTHUNIT["metre",1], ID["EPSG",8806]], PARAMETER["False northing",30, LENGTHUNIT["metre",1], ID["EPSG",8807]]], CS[Cartesian,2], AXIS["(E)",east, ORDER[1], LENGTHUNIT["metre",1]], AXIS["(N)",north, ORDER[2], LENGTHUNIT["metre",1]], USAGE[ SCOPE["unknown"], AREA["Qatar - onshore"], BBOX[24.55,50.69,26.2,51.68]], ID["EPSG",2932]] The problem is that GRASS still assumes datum transformation from X to WGS84, whereas with PROJ6, WGS84 is no longer required as pivot datum. Datum transformations from datum X to datum Y can sometimes (often) be done without going through WGS84. The requirement is to have EPSG or other authority codes and to use PROJ 6. this warning is just from a simple v.in.ogr of a EPSG:2932-based shapefile into a EPSG:2932-created GRASS location. so, there may be some more adaptions needed for actions like v.in.ogr and co? - best regards Helmut -- Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user -- = Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 (WhatsApp only) Qatar: +974 5083 2722 www.geograph.co.za = ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Datum not recognized by Grass
Hi Markus, Thanks for the interest in my problem. Herewith g.version output: Kind regards, Zoltan (Wed Dec 4 10:23:18 2019) D:\sandbox\\tests\PERMANENT\.tmp/unknown\1608.0.py version=7.8.1 date=2019 revision=c865432c9 build_date=2019-11-10 build_platform=x86_64-w64-mingw32 build_off_t_size=8 libgis_revision=0 libgis_date="?" proj=6.2.1 gdal=3.0.2 geos=3.8.0 sqlite=3.29.0 (Wed Dec 4 10:23:18 2019) Command finished (0 sec) On 2019/12/03 20:06, Markus Neteler wrote: On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei 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)? Please post the output of g.version -rge Best Markus -- = Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 (WhatsApp only) Qatar: +974 5083 2722 www.geograph.co.za = ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Datum not recognized by Grass
> ignore the warning and use GRASS with PROJ6, granted that authority name (e.g. EPSG) and authority code (e.g. 2932) are known for both CRS's in case of reprojection this issue is already by GRASS with PROJ6, see GRASS version: 7.8.1 Code revision: c865432c9 Build date: 2019-11-10 Build platform: x86_64-w64-mingw32 GDAL: 3.0.2 PROJ: 6.2.1 <= GEOS: 3.8.0 SQLite: 3.29.0 Python: 3.7.0 wxPython: 4.0.7 Platform: Windows-10-10.0.18362-SP0 (OSGeo4W) and the output from the underlying PROJ: 6.2.1 C:\>projinfo EPSG:2932 -o PROJ,WKT2_2018 PROJ.4 string: +proj=tmerc +lat_0=24.45 +lon_0=51.21667 +k=0.9 +x_0=20 +y_0=30 +ellps=intl +towgs84=-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065 +units=m +no_defs +type=crs WKT2_2018 string: PROJCRS["QND95 / Qatar National Grid", BASEGEOGCRS["QND95", DATUM["Qatar National Datum 1995", ELLIPSOID["International 1924",6378388,297, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4614]], CONVERSION["Qatar National Grid", METHOD["Transverse Mercator", ID["EPSG",9807]], PARAMETER["Latitude of natural origin",24.45, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8801]], PARAMETER["Longitude of natural origin",51.21667, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8802]], PARAMETER["Scale factor at natural origin",0.9, SCALEUNIT["unity",1], ID["EPSG",8805]], PARAMETER["False easting",20, LENGTHUNIT["metre",1], ID["EPSG",8806]], PARAMETER["False northing",30, LENGTHUNIT["metre",1], ID["EPSG",8807]]], CS[Cartesian,2], AXIS["(E)",east, ORDER[1], LENGTHUNIT["metre",1]], AXIS["(N)",north, ORDER[2], LENGTHUNIT["metre",1]], USAGE[ SCOPE["unknown"], AREA["Qatar - onshore"], BBOX[24.55,50.69,26.2,51.68]], ID["EPSG",2932]] >The problem is that GRASS still assumes datum transformation from X to WGS84, whereas with PROJ6, >WGS84 is no longer required as pivot datum. Datum transformations from datum X to datum Y can >sometimes (often) be done without going through WGS84. The requirement is to have EPSG or other >authority codes and to use PROJ 6. this warning is just from a simple v.in.ogr of a EPSG:2932-based shapefile into a EPSG:2932-created GRASS location. so, there may be some more adaptions needed for actions like v.in.ogr and co? - best regards Helmut -- Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Datum not recognized by Grass
There are some solutions to the problem of datum not recognized by GRASS: - add Qatar_National_Datum_1995 to the datums known by GRASS - ignore the warning and use GRASS with PROJ6, granted that authority name (e.g. EPSG) and authority code (e.g. 2932) are known for both CRS's in case of reprojection The problem is that GRASS still assumes datum transformation from X to WGS84, whereas with PROJ6, WGS84 is no longer required as pivot datum. Datum transformations from datum X to datum Y can sometimes (often) be done without going through WGS84. The requirement is to have EPSG or other authority codes and to use PROJ 6. Markus M On Tue, Dec 3, 2019 at 8:52 PM Helmut Kudrnovsky wrote: > Markus Neteler wrote > > Hi, > > > > On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei > > > zoltans@.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 > > > > 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 > > > > ... > > 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.216667], > > PARAMETER["scale_factor",0.9], > > PARAMETER["false_easting",20], > > PARAMETER["false_northing",30], > > 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 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.21667], > PARAMETER["scale_factor",0.9], > PARAMETER["false_easting",20], > PARAMETER["false_northing",30], > 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"]],
Re: [GRASS-user] Datum not recognized by Grass
Markus Neteler wrote > Hi, > > On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei > zoltans@.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 > > 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 > > ... > 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.216667], > PARAMETER["scale_factor",0.9], > PARAMETER["false_easting",20], > PARAMETER["false_northing",30], > 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 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.21667], PARAMETER["scale_factor",0.9], PARAMETER["false_easting",20], PARAMETER["false_northing",30], 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.21667], PARAMETER["scale_factor",0.9], PARAMETER["false_easting",20], PARAMETER["false_northing",30], 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.21667],
Re: [GRASS-user] Datum not recognized by Grass
Markus Neteler wrote > Hi, > > On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei > zoltans@.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 > > 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 > > ... > 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.216667], > PARAMETER["scale_factor",0.9], > PARAMETER["false_easting",20], > PARAMETER["false_northing",30], > UNIT["metre",1, > AUTHORITY["EPSG","9001"]], > AXIS["Easting",EAST], > AXIS["Northing",NORTH], > AUTHORITY["EPSG","2932"]] C:\>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.21667], PARAMETER["scale_factor",0.9], PARAMETER["false_easting",20], PARAMETER["false_northing",30], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AXIS["Easting",EAST], AXIS["Northing",NORTH], AUTHORITY["EPSG","2932"]] - best regards Helmut -- Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Datum not recognized by Grass
Markus Neteler wrote > Hi, > > On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei > zoltans@.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 > > 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 > > ... > 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.216667], > PARAMETER["scale_factor",0.9], > PARAMETER["false_easting",20], > PARAMETER["false_northing",30], > UNIT["metre",1, > AUTHORITY["EPSG","9001"]], > AXIS["Easting",EAST], > AXIS["Northing",NORTH], > AUTHORITY["EPSG","2932"]] > > # The versions here: > > GRASS 7.8.2dev (test_2932):~ > g.version -rge > version=7.8.2dev > date=2019 > revision=ad4836c73 > build_date=2019-11-28 > build_platform=x86_64-pc-linux-gnu > build_off_t_size=8 > libgis_revision=0 > libgis_date="?" > proj=5.2.0 > gdal=2.3.2 > geos=3.7.1 > sqlite=3.30.0 > > Looks all fine. > Please check if an old PROJ 4 is installed which might confuse things. with C:\>g.version -rge version=7.8.1 date=2019 revision=c865432c9 build_date=2019-11-10 build_platform=x86_64-w64-mingw32 build_off_t_size=8 libgis_revision=0 libgis_date="?" proj=6.2.1 gdal=3.0.2 geos=3.8.0 sqlite=3.29.0 I get with the location wizzard: Starting GRASS GIS... Traceback (most recent call last): File "C:\OSGEO4~1\apps\grass\grass78\gui\wxpython\location_wizard\wizard.py", line 1607, in OnPageChanging self.parent.parent, transforms=ret) File "C:\OSGEO4~1\apps\grass\grass78\gui\wxpython\location_wizard\dialogs.py", line 617, in __init__ transforms = '---\n\n0\nDo not apply any datum transformations\n\n' + transforms TypeError: can only concatenate str (not "bytes") to str - best regards Helmut -- Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Datum not recognized by Grass
On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei 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)? Please post the output of g.version -rge Best Markus ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Datum not recognized by Grass
Hi, On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei 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 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 ... 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.216667], PARAMETER["scale_factor",0.9], PARAMETER["false_easting",20], PARAMETER["false_northing",30], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AXIS["Easting",EAST], AXIS["Northing",NORTH], AUTHORITY["EPSG","2932"]] # The versions here: GRASS 7.8.2dev (test_2932):~ > g.version -rge version=7.8.2dev date=2019 revision=ad4836c73 build_date=2019-11-28 build_platform=x86_64-pc-linux-gnu build_off_t_size=8 libgis_revision=0 libgis_date="?" proj=5.2.0 gdal=2.3.2 geos=3.7.1 sqlite=3.30.0 Looks all fine. Please check if an old PROJ 4 is installed which might confuse things. Best Markus ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Datum not recognized by Grass
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 not recognised by GRASS and no parameters found Check if OGR layer <23103770_road_poly> contains polygons... Creating attribute table for layer <23103770_road_poly>... WARNING: Name <23103770_road_poly> is not SQL compliant. Must start with a letter. Importing 28 features (OGR layer <23103770_road_poly>)... -- = Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 (WhatsApp only) Qatar: +974 5083 2722 www.geograph.co.za = ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Datum not recognized by GRASS
Hi Anders, (please keep it on the list) unfortunately you need a recompiled version. I wonder if anyone is creating weekly binary snapshots of GRASS 6.4 for Mac OS X... Otherwise - you could recompile yourself (far from impossible, I am geographer :-) ) - you could amend the metadata with the gdaltools gdal_translate -a_srs ... oldfile.tif newfile.tif (see also http://gfoss.blogspot.com/2008/06/gdal-raster-data-tips-and-tricks.html - Fixing broken projection/datum info for raster data ) cheers Markus 2010/7/2 Anders Gonçalves da Silva ander...@gmail.com: Hi Markus and Paul. Thank you for the help. I am not positive, though, how to actually go about correcting this on my own version of GRASS. ... Mac OS X ... ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Datum not recognized by GRASS
On Wed, 30 Jun 2010, Anders Gonçalves da Silva wrote: Hello all. I am trying to add a GeoTIFF to a project, but I keep getting the following message: -- GRASS 6.4.0RC6 (pontal):~/Documents/Academic/Thesis/Archive/Sampling gdalinfo pontal.TIF Driver: GTiff/GeoTIFF Files: pontal.TIF Size is 4337, 2727 Coordinate System is: PROJCS[UTM Zone 23S, GEOGCS[South_American_1969, DATUM[South_American_1969, This seems to be non-standard; the EPSG name should be South_American_Datum_1969 and then GRASS would pick it up correctly. Perhaps we need to add an exception in lib/proj/convert.c; unfortunately I'm just about to head away for a few days and don't have time right now. Paul SPHEROID[South American 1969,6378160,298.250654724]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,-45], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,1000], UNIT[unknown,1]] Origin = (280686.23560014529,7573063.00580391901) Pixel Size = (30.000,-30.000) Metadata: TIFFTAG_SOFTWARE=ERDAS IMAGINE AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=PIXEL Corner Coordinates: Upper Left ( 280686.236, 7573063.006) ( 47d 7'24.05W, 21d55'58.77S) Lower Left ( 280686.236, 7491253.006) ( 47d 8'4.36W, 22d40'17.64S) Upper Right ( 410796.236, 7573063.006) ( 45d51'49.82W, 21d56'39.98S) Lower Right ( 410796.236, 7491253.006) ( 45d52'6.24W, 22d41'0.39S) Center ( 345741.236, 7532158.006) ( 46d29'51.35W, 22d18'33.65S) Band 1 Block=4337x5 Type=Byte, ColorInterp=Red Band 2 Block=4337x5 Type=Byte, ColorInterp=Green Band 3 Block=4337x5 Type=Byte, ColorInterp=Blue -- I have tried to use gdalwarp to reset it to something that GRASS can read. While I get it to run, the projections are not lining up. In particular, it seems to be off by a few degrees to the East. Any help in fixing this would be greatly appreciated. Thank you. Anders.___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Datum not recognized by GRASS
2010/6/30 Paul Kelly paul-gr...@stjohnspoint.co.uk: On Wed, 30 Jun 2010, Anders Gonçalves da Silva wrote: Hello all. I am trying to add a GeoTIFF to a project, but I keep getting the following message: -- GRASS 6.4.0RC6 (pontal):~/Documents/Academic/Thesis/Archive/Sampling gdalinfo pontal.TIF Driver: GTiff/GeoTIFF Files: pontal.TIF Size is 4337, 2727 Coordinate System is: PROJCS[UTM Zone 23S, GEOGCS[South_American_1969, DATUM[South_American_1969, This seems to be non-standard; the EPSG name should be South_American_Datum_1969 and then GRASS would pick it up correctly. Perhaps we need to add an exception in lib/proj/convert.c; unfortunately I'm just about to head away for a few days and don't have time right now. I have added it in 6.4, 6.5, 7 in: lib/proj/convert.c --- convert.c (revision 42660) +++ convert.c (working copy) @@ -744,6 +744,8 @@ Militar_Geographische_Institut, Potsdam_Datum_83, Deutsches_Hauptdreiecksnetz, +South_American_1969, +South_American_Datum_1969, NULL }; Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Datum not recognized by GRASS
Hello all. I am trying to add a GeoTIFF to a project, but I keep getting the following message: -- GRASS 6.4.0RC6 (pontal):~/Documents/Academic/Thesis/Archive/Sampling gdalinfo pontal.TIF Driver: GTiff/GeoTIFF Files: pontal.TIF Size is 4337, 2727 Coordinate System is: PROJCS[UTM Zone 23S, GEOGCS[South_American_1969, DATUM[South_American_1969, SPHEROID[South American 1969,6378160,298.250654724]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,-45], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,1000], UNIT[unknown,1]] Origin = (280686.23560014529,7573063.00580391901) Pixel Size = (30.000,-30.000) Metadata: TIFFTAG_SOFTWARE=ERDAS IMAGINE AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=PIXEL Corner Coordinates: Upper Left ( 280686.236, 7573063.006) ( 47d 7'24.05W, 21d55'58.77S) Lower Left ( 280686.236, 7491253.006) ( 47d 8'4.36W, 22d40'17.64S) Upper Right ( 410796.236, 7573063.006) ( 45d51'49.82W, 21d56'39.98S) Lower Right ( 410796.236, 7491253.006) ( 45d52'6.24W, 22d41'0.39S) Center ( 345741.236, 7532158.006) ( 46d29'51.35W, 22d18'33.65S) Band 1 Block=4337x5 Type=Byte, ColorInterp=Red Band 2 Block=4337x5 Type=Byte, ColorInterp=Green Band 3 Block=4337x5 Type=Byte, ColorInterp=Blue -- I have tried to use gdalwarp to reset it to something that GRASS can read. While I get it to run, the projections are not lining up. In particular, it seems to be off by a few degrees to the East. Any help in fixing this would be greatly appreciated. Thank you. Anders.___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
