[GRASS-user] g.proj fails on winGRASS trunk with AIG/Arc/Info Binary Grid input

[ creating a new thread because this is a different issue ]

On Mon, Jun 11, 2018 at 7:47 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

strange, works for me on Linux.

can you provide the winGRASS 7.5 output of
g.proj georef=D:\temp\rich\Bare_Earth\bh45122c3\hdr.adf -p
and
g.proj georef=D:\temp\rich\Bare_Earth\bh45122c3 -p


System Info
GRASS version: 7.4.0
GRASS SVN revision: r72154
Build date: 2018-01-27
Build platform: x86_64-w64-mingw32
GDAL: 2.2.3
PROJ.4: 4.9.3
GEOS: 3.5.0
SQLite: 3.17.0
Python: 2.7.5
wxPython: 2.8.12.1
Platform: Windows-8-6.2.9200 (OSGeo4W)

g.proj -p georef=D:\temp\rich\bareearth\Bare_Earth\bh45122c3
-PROJ_INFO-------------------------------------------------
name : unnamed
ellps : grs80
proj : lcc
lat_1 : 43
lat_2 : 45.5
lat_0 : 41.75
lon_0 : -120.5
x_0 : 399999.9999999999
y_0 : 0
no_defs : defined
-PROJ_UNITS------------------------------------------------
unit : user-defined
units : user-defineds
meters : 0.3048
Trying to open with OGR…
Trying to open with GDAL…
…succeeded.
WARNING: Datum <Not_specified_based_on_GRS_1980_ellipsoid> not recognised by
GRASS and no parameters found


g.proj -p georef=D:\temp\rich\bareearth\Bare_Earth\bh45122c3\hdr.adf
Trying to open with OGR…
Trying to open with GDAL…
…succeeded.
WARNING: Datum <Not_specified_based_on_GRS_1980_ellipsoid> not recognised by
GRASS and no parameters found
-PROJ_INFO-------------------------------------------------
name : unnamed
ellps : grs80
proj : lcc
lat_1 : 43
lat_2 : 45.5
lat_0 : 41.75
lon_0 : -120.5
x_0 : 399999.9999999999
y_0 : 0
no_defs : defined
-PROJ_UNITS------------------------------------------------
unit : user-defined
units : user-defineds
meters : 0.3048

GRASS version: 7.5.svn
GRASS SVN revision: r72772
Build date: 2018-06-06
Build platform: x86_64-w64-mingw32
GDAL: 2.2.3
PROJ.4: 4.9.3
GEOS: 3.5.0
SQLite: 3.17.0
Python: 2.7.5
wxPython: 2.8.12.1
Platform: Windows-8-6.2.9200 (OSGeo4W)


g.proj -p georef=D:\temp\rich\bareearth\Bare_Earth\bh45122c3
(Mon Jun 11 19:39:41 2018) Command finished (2 sec)


g.proj -p georef=D:\temp\rich\bareearth\Bare_Earth\bh45122c3\hdr.adf
(Mon Jun 11 19:40:49 2018) Command finished (4 sec)

No information: can you please try this on the true CLI?


now DEBUG=3 in 7.5.svn


g.proj -p georef=D:\temp\rich\bareearth\Bare_Earth\bh45122c3\hdr.adf
D1/3: G_set_program_name(): g.proj
D1/3: Trying to open
<D:\temp\rich\bareearth\Bare_Earth\bh45122c3\hdr.adf> with
OGR…
D1/3: Trying to open with GDAL…
D1/3: …succeeded.
D3/3: DatumNameMassage: Raw string found
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: DatumNameMassage: Search for datum equivalences of
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: GPJ_osr_to_grass: pszDatumNameConst:
<Not_specified_based_on_GRS_1980_ellipsoid>
(Mon Jun 11 19:42:53 2018) Command finished (4 sec)


g.proj -p georef=D:\temp\rich\bareearth\Bare_Earth\bh45122c3
D1/3: G_set_program_name(): g.proj
D1/3: Trying to open
<D:\temp\rich\bareearth\Bare_Earth\bh45122c3> with OGR…
D1/3: Trying to open with GDAL…
D1/3: …succeeded.
D3/3: DatumNameMassage: Raw string found
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: DatumNameMassage: Search for datum equivalences of
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: GPJ_osr_to_grass: pszDatumNameConst:
<Not_specified_based_on_GRS_1980_ellipsoid>


not much information.

Indeed not. It seems like the process dies somewhere in GPJ_osr_to_grass(), but only in winGRASS. This is probably worth a new ticket.

Markus M


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

No information: can you please try this on the true CLI?

C:\>g.proj -p georef=D:\temp\rich\Bare_Earth\bh45122c3
D1/3: G_set_program_name(): g.proj
D1/3: Trying to open <D:\temp\rich\Bare_Earth\bh45122c3> with OGR...
D1/3: Trying to open with GDAL...
D1/3: ...succeeded.
D3/3: DatumNameMassage: Raw string found
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: DatumNameMassage: Search for datum equivalences of
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: GPJ_osr_to_grass: pszDatumNameConst:
<Not_specified_based_on_GRS_1980_ellipsoid>
WARNING: Datum <Not_specified_based_on_GRS_1980_ellipsoid> not recognised
         by GRASS and no parameters found

C:\>g.proj -p georef=D:\temp\rich\Bare_Earth\bh45122c3\hdr.adf
D1/3: G_set_program_name(): g.proj
D1/3: Trying to open <D:\temp\rich\Bare_Earth\bh45122c3\hdr.adf> with OGR...
D1/3: Trying to open with GDAL...
D1/3: ...succeeded.
D3/3: DatumNameMassage: Raw string found
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: DatumNameMassage: Search for datum equivalences of
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: GPJ_osr_to_grass: pszDatumNameConst:
<Not_specified_based_on_GRS_1980_ellipsoid>
WARNING: Datum <Not_specified_based_on_GRS_1980_ellipsoid> not recognised
         by GRASS and no parameters found

Indeed not. It seems like the process dies somewhere in GPJ_osr_to_grass(),

but only in winGRASS.

it seems so.

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

On Mon, 11 Jun 2018, Helmut Kudrnovsky wrote:

Indeed not. It seems like the process dies somewhere in GPJ_osr_to_grass(),
but only in winGRASS.
it seems so.

All:

   The process dies when using Slackware-14.2 with grass-7.5svn rev. 72794,
gdal-2.3.0, and proj-5.1.0.

   I thought of trying to convert the data from AIG to OPKG using
gdal_translate, but I've no idea of what options and value to select since
I've not before worked translating spatial data formats. For example, since
these are LiDAR data and not remotely sensed satallite imagery I've no idea
how to specify bands. (I've not worked with such imagery since GRASS-4.1.)

   If you think trying Open GeoPackage is worth doing, please suggest
gdal_translate options for me to use.

Regards,

Rich

  If you think trying Open GeoPackage is worth doing, please >suggest
gdal_translate options for me to use.

regarding gdal_translate options and formats, it's a questions for
consulting the GDAL manual and/or GDAL mailing list.

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

On Mon, Jun 11, 2018 at 11:01 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

No information: can you please try this on the true CLI?

C:>g.proj -p georef=D:\temp\rich\Bare_Earth\bh45122c3
D1/3: G_set_program_name(): g.proj
D1/3: Trying to open <D:\temp\rich\Bare_Earth\bh45122c3> with OGR…
D1/3: Trying to open with GDAL…
D1/3: …succeeded.
D3/3: DatumNameMassage: Raw string found
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: DatumNameMassage: Search for datum equivalences of
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: GPJ_osr_to_grass: pszDatumNameConst:
<Not_specified_based_on_GRS_1980_ellipsoid>
WARNING: Datum <Not_specified_based_on_GRS_1980_ellipsoid> not recognised
by GRASS and no parameters found

C:>g.proj -p georef=D:\temp\rich\Bare_Earth\bh45122c3\hdr.adf
D1/3: G_set_program_name(): g.proj
D1/3: Trying to open <D:\temp\rich\Bare_Earth\bh45122c3\hdr.adf> with OGR…
D1/3: Trying to open with GDAL…
D1/3: …succeeded.
D3/3: DatumNameMassage: Raw string found
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: DatumNameMassage: Search for datum equivalences of
<Not_specified_based_on_GRS_1980_ellipsoid>
D3/3: GPJ_osr_to_grass: pszDatumNameConst:
<Not_specified_based_on_GRS_1980_ellipsoid>
WARNING: Datum <Not_specified_based_on_GRS_1980_ellipsoid> not recognised
by GRASS and no parameters found

Indeed not. It seems like the process dies somewhere in GPJ_osr_to_grass(),
but only in winGRASS.

it seems so.

Gotcha! Valgrind gave me a hint, hopefully fixed in trunk r72802.

Markus M


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

On Tue, 12 Jun 2018, Markus Metz wrote:

Gotcha! Valgrind gave me a hint, hopefully fixed in trunk r72802.

Markus M,

   Here, running 7.5svn r72804 there's no change with my data set.

Rich

On Tue, Jun 12, 2018 at 4:07 PM, Rich Shepard <rshepard@appl-ecosys.com> wrote:

On Tue, 12 Jun 2018, Markus Metz wrote:

Gotcha! Valgrind gave me a hint, hopefully fixed in trunk r72802.

Markus M,

Here, running 7.5svn r72804 there’s no change with my data set.

as of r72802, it should be possible to create a location from the data, and the projection check should work.

About valid values only along the edges, that’s a different issue, most probably there is nothing else in the input data. The data provider would know if this is correct.

Markus M

Rich


grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

On Tue, 12 Jun 2018, Markus Metz wrote:

as of r72802, it should be possible to create a location from the data, and
the projection check should work.

Markus,

   I've not had an issue with creating a location, only importing the data.

About valid values only along the edges, that's a different issue, most
probably there is nothing else in the input data. The data provider would
know if this is correct.

   The other data, highest_hits, has the same problem. I'll download the data
again.

Rich

as of r72802, it should be possible to create a location from the data, and

the projection check should work.

confirmed, it works now also in winGRASS trunk.

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

On Wed, Jun 13, 2018 at 10:25 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

as of r72802, it should be possible to create a location from the data, and
the projection check should work.

confirmed, it works now also in winGRASS trunk.

Great, thanks a lot for testing on winGRASS!

Markus M


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