[GRASS-user] r.proj ERROR

GRASS GRASS user
1

Hi,
I want to re-project DEM data from a lambert-location into ETRS89/utm zone
33N
I did that many times before but now I get the following message:
                                              
r.proj location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
input=dhm_10m_ogd_austria output=dhm_10m_ogd_dachstein_NEU method=bilinear
resolution=10
ERROR: Unable to open element file <> for <WIND@dachstein>
(Sat Jan 11 14:24:16 2020) Befehl ausgeführt (0 Sek)

The WIND file is ok. When I open it there are the correct region settings.

Whats wrong?
thx!
Nik

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

2

Hi,

On Sat, Jan 11, 2020 at 2:28 PM nik <n.hafner@carto.at> wrote:

Hi,
I want to re-project DEM data from a lambert-location into ETRS89/utm zone
33N
I did that many times before but now I get the following message:

r.proj location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
input=dhm_10m_ogd_austria output=dhm_10m_ogd_dachstein_NEU method=bilinear
resolution=10
ERROR: Unable to open element file <> for <WIND@dachstein>
(Sat Jan 11 14:24:16 2020) Befehl ausgeführt (0 Sek)

The WIND file is ok. When I open it there are the correct region settings.

Whats wrong?

We need to know GRASS GIS version and operating system...

Best
Markus

3

nik wrote

Hi,
I want to re-project DEM data from a lambert-location into ETRS89/utm zone
33N
I did that many times before but now I get the following message:
                                              
r.proj location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
input=dhm_10m_ogd_austria output=dhm_10m_ogd_dachstein_NEU method=bilinear
resolution=10
ERROR: Unable to open element file <> for &lt;WIND@dachstein&gt;
(Sat Jan 11 14:24:16 2020) Befehl ausgeführt (0 Sek)

The WIND file is ok. When I open it there are the correct region settings.

Whats wrong?
thx!
Nik

which GRASS version and operating system are you using?

testing here with:

System Info
GRASS version: 7.8.2
Code revision: 3900fb114
Build date: 2019-12-11
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 tried here with a raster of the Dachstein region :wink:

g.proj -p
-PROJ_INFO-------------------------------------------------
name : MGI / Austria Lambert
datum : hermannskogel
ellps : bessel
proj : lcc
lat_0 : 47.5
lon_0 : 13.3333333333333
lat_1 : 49
lat_2 : 46
x_0 : 400000
y_0 : 400000
no_defs : defined
towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232
-PROJ_EPSG-------------------------------------------------
epsg : 31287
-PROJ_UNITS------------------------------------------------
unit : meter
units : meters
meters : 1

----------------------------

g.region -p
projection: 99 (MGI / Austria Lambert)
zone: 0
datum: hermannskogel
ellipsoid: bessel
north: 413800
south: 391100
west: 407600
east: 441400
nsres: 100
ewres: 100
rows: 227
cols: 338
cells: 76726

#################################################################

g.proj -p
-PROJ_INFO-------------------------------------------------
name : ETRS89 / UTM zone 33N
datum : etrs89
ellps : grs80
proj : utm
zone : 33
no_defs : defined
towgs84 : 0.000,0.000,0.000
-PROJ_EPSG-------------------------------------------------
epsg : 25833
-PROJ_UNITS------------------------------------------------
unit : meter
units : meters
meters : 1

----------------------------

r.proj -g location=dachstein mapset=data input=raster2proj
WARNING: Found 2 possible transformations
************************
Operation 1:
Description: axis order change (2D) + Inverse of Austria Lambert + MGI to
ETRS89 (1) + UTM zone 33N

Area of use: Austria

Accuracy within area of use: 1.5 m

PROJ string:
+proj=pipeline +step +inv +proj=lcc +lat_0=47.5 +lon_0=13.3333333333333
+lat_1=49 +lat_2=46 +x_0=400000 +y_0=400000 +ellps=bessel +step +proj=push
+v_3 +step +proj=cart +ellps=bessel +step +proj=helmert +x=577.326 +y=90.129
+z=463.919 +rx=5.137 +ry=1.474 +rz=5.297 +s=2.4232
+convention=position_vector +step +inv +proj=cart +ellps=GRS80 +step
+proj=pop +v_3 +step +proj=utm +zone=33 +ellps=GRS80
************************
Operation 2:
Description: axis order change (2D) + Inverse of Austria Lambert + MGI to
ETRS89 (5) + UTM zone 33N

Area of use: Austria

Accuracy within area of use: 0.15 m

PROJ string:
+proj=pipeline +step +inv +proj=lcc +lat_0=47.5 +lon_0=13.3333333333333
+lat_1=49 +lat_2=46 +x_0=400000 +y_0=400000 +ellps=bessel +step
+proj=hgridshift +grids=AT_GIS_GRID.gsb +step +proj=utm +zone=33
+ellps=GRS80
************************
See also output of:
projinfo -o PROJ -s "EPSG:31287" -t "EPSG:25833"
Please provide the appropriate PROJ string with the pipeline option
************************
Selected pipeline:
+proj=pipeline +step +inv +proj=lcc +lat_0=47.5 +lon_0=13.3333333333333
+lat_1=49 +lat_2=46 +x_0=400000 +y_0=400000 +ellps=bessel +step +proj=push
+v_3 +step +proj=cart +ellps=bessel +step +proj=helmert +x=577.326 +y=90.129
+z=463.919 +rx=5.137 +ry=1.474 +rz=5.297 +s=2.4232
+convention=position_vector +step +inv +proj=cart +ellps=GRS80 +step
+proj=pop +v_3 +step +proj=utm +zone=33 +ellps=GRS80
************************
Input map <raster2proj@data> in location <dachstein>:
n=5275655.1124836 s=5252234.92874731 w=381830.9328192 e=416111.44820006
rows=227 cols=338

----------------------------

g.region n=5275655.1124836 s=5252234.92874731 w=381830.9328192
e=416111.44820006 rows=227 cols=338

g.region -p
projection: 1 (UTM)
zone: 33
datum: etrs89
ellipsoid: grs80
north: 5275655.1124836
south: 5252234.92874731
west: 381830.9328192
east: 416111.44820006
nsres: 103.17261558
ewres: 101.42164314
rows: 227
cols: 338
cells: 76726

----------------------------

r.proj location=dachstein mapset=data input=raster2proj
output=rasterreprojected pipeline=+proj=pipeline +step +inv +proj=lcc
+lat_0=47.5 +lon_0=13.3333333333333 +lat_1=49 +lat_2=46 +x_0=400000
+y_0=400000 +ellps=bessel +step +proj=push +v_3 +step +proj=cart
+ellps=bessel +step +proj=helmert +x=577.326 +y=90.129 +z=463.919 +rx=5.137
+ry=1.474 +rz=5.297 +s=2.4232 +convention=position_vector +step +inv
+proj=cart +ellps=GRS80 +step +proj=pop +v_3 +step +proj=utm +zone=33
+ellps=GRS80

Input:
Cols: 338 (original: 338)
Rows: 227 (original: 227)
North: 413800.000000 (original: 413800.000000)
South: 391100.000000 (original: 391100.000000)
West: 407600.000000 (original: 407600.000000)
East: 441400.000000 (original: 441400.000000)
EW-res: 100.000000
NS-res: 100.000000

Output:
Cols: 338 (original: 338)
Rows: 227 (original: 227)
North: 5275655.112484 (original: 5275655.112484)
South: 5252234.928747 (original: 5252234.928747)
West: 381830.932819 (original: 381830.932819)
East: 416111.448200 (original: 416111.448200)
EW-res: 101.421643
NS-res: 103.172616

Allocating memory and reading input raster map...
Projecting...
r.proj complete.

it's using the new proj library version 6

C:\>projinfo -o PROJ -s "EPSG:31287" -t "EPSG:25833"
Candidate operations found: 2
Note: using '--spatial-test intersects' would bring more results (3)
-------------------------------------
Operation n┬░1:

unknown id, Inverse of Austria Lambert + MGI to ETRS89 (1) + UTM zone 33N,
1.5 m, Austria

PROJ string:
+proj=pipeline +step +proj=axisswap +order=2,1 +step +inv +proj=lcc
+lat_0=47.5 +lon_0=13.3333333333333 +lat_1=49 +lat_2=46 +x_0=400000
+y_0=400000 +ellps=bessel +step +proj=push +v_3 +step +proj=cart
+ellps=bessel +step +proj=helmert +x=577.326 +y=90.129 +z=463.919 +rx=5.137
+ry=1.474 +rz=5.297 +s=2.4232 +convention=position_vector +step +inv
+proj=cart +ellps=GRS80 +step +proj=pop +v_3 +step +proj=utm +zone=33
+ellps=GRS80

-------------------------------------
Operation n┬░2:

unknown id, Inverse of Austria Lambert + MGI to ETRS89 (5) + UTM zone 33N,
0.15 m, Austria, at least one grid missing

PROJ string:
+proj=pipeline +step +proj=axisswap +order=2,1 +step +inv +proj=lcc
+lat_0=47.5 +lon_0=13.3333333333333 +lat_1=49 +lat_2=46 +x_0=400000
+y_0=400000 +ellps=bessel +step +proj=hgridshift +grids=AT_GIS_GRID.gsb
+step +proj=utm +zone=33 +ellps=GRS80

Grid AT_GIS_GRID.gsb needed but not found on the system.

in the test I used the operation 1 as AT_GIS_GRID.gsb isn't installed here
in my system.

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

4

On Sat, Jan 11, 2020 at 5:55 PM Helmut Kudrnovsky <hellik@web.de> wrote:
[...]

unknown id, Inverse of Austria Lambert + MGI to ETRS89 (5) + UTM zone 33N,
0.15 m, Austria, at least one grid missing

PROJ string:
+proj=pipeline +step +proj=axisswap +order=2,1 +step +inv +proj=lcc
+lat_0=47.5 +lon_0=13.3333333333333 +lat_1=49 +lat_2=46 +x_0=400000
+y_0=400000 +ellps=bessel +step +proj=hgridshift +grids=AT_GIS_GRID.gsb
+step +proj=utm +zone=33 +ellps=GRS80

Grid AT_GIS_GRID.gsb needed but not found on the system.

It has been added [1] and released only recently [2].
It may take some time to be shipped in the various distro/Windows packages.

Markus

[1] https://github.com/OSGeo/proj-datumgrid/issues/71
[2] https://github.com/OSGeo/proj-datumgrid/releases/tag/europe-1.5

5

Hi

which GRASS version and operating system are you using?

I also use:
GRASS Version: 7.8.2
Code revision: 3900fb114
Build date: 2019-12-11
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)

I still get this message:
r.proj location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
input=dhm_10m_ogd_austria output=dhm_10m_ogd_dachstein_NEU method=bilinear
resolution=10
ERROR: Unable to open element file <> for <WIND@dachstein>

and I tested it from an other location too but its the same problem:
r.proj location=ASTER_lat-long mapset=PERMANENT
input=DHM_SRTM30m_Aut-plus-Umgebung output=dhm_30m_ASTER_dachstein
method=bilinear resolution=30
ERROR: Unable to open element file <> for <WIND@dachstein>

What can be wrong?
Nik

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

6

nik wrote

Hi

which GRASS version and operating system are you using?

I also use:
GRASS Version: 7.8.2
Code revision: 3900fb114
Build date: 2019-12-11
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)

I still get this message:
r.proj location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
input=dhm_10m_ogd_austria output=dhm_10m_ogd_dachstein_NEU method=bilinear
resolution=10
ERROR: Unable to open element file <> for &lt;WIND@dachstein&gt;

and I tested it from an other location too but its the same problem:
r.proj location=ASTER_lat-long mapset=PERMANENT
input=DHM_SRTM30m_Aut-plus-Umgebung output=dhm_30m_ASTER_dachstein
method=bilinear resolution=30
ERROR: Unable to open element file <> for &lt;WIND@dachstein&gt;

What can be wrong?
Nik

--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html
_______________________________________________
grass-user mailing list

grass-user@.osgeo

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

What does r.proj -g output? With PROJ: 6.2.1, a proj pipeline has to be
given in r.proj. it's not in your cmd

Can you report all steps with g.proj, g.region, r.proj in source and
target locations as I've done?

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

7

What does r.proj -g output? With PROJ: 6.2.1, a proj pipeline has to be
given in r.proj. it's not in your cmd

Can you report all steps with g.proj, g.region, r.proj in source and
target locations as I've done?

-----
best regards
Helmut
--

g.proj -p
-PROJ_INFO-------------------------------------------------
name : Lambert Conformal Conic
proj : lcc
datum : hermannskogel
ellps : bessel
lat_1 : 49
lat_2 : 46
lat_0 : 47.5
lon_0 : 13.33333333333333
x_0 : 400000
y_0 : 400000
no_defs : defined
towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1

----------------------------------------------------------------------------------------------

g.region -p
projection: 99 (Lambert Conformal Conic)
zone: 0
datum: hermannskogel
ellipsoid: bessel
north: 576922.51207363
south: 273692.51207363
west: 106549.26720377
east: 694629.26720377
nsres: 10
ewres: 10
rows: 30323
cols: 58808
cells: 1783234984

============================================================

g.proj -p
-PROJ_INFO-------------------------------------------------
name : Universal Transverse Mercator
proj : utm
datum : etrs89
ellps : grs80
zone : 33
towgs84 : 0,0,0,0,0,0,0
no_defs : defined
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1

----------------------------------------------------------------------------------------------

g.region -p
projection: 1 (UTM)
zone: 33
datum: etrs89
ellipsoid: grs80
north: 5286581.58847414
south: 5248335.67358037
west: 380870.53610518
east: 429633.13185005
nsres: 9.998932
ewres: 10.00053235
rows: 3825
cols: 4876
cells: 18650700

----------------------------------------------------------------------------------------------

I downloaded and installed http://www.bev.gv.at/portal/page?_pageid=713,2157075&_dad=portal&_schema=PORTAL

r.proj -g location=MGI2wgs84_Austria_Lambert mapset=PERMANENT input=dhm_10m_ogd_austria output=dhm_10m_ogd_dachstein_NEU method=bilinear resolution=10 pipeline=+proj=pipeline +step +inv +proj=lcc +lat_0=47.5 +lon_0=13.3333333333333 +lat_1=49 +lat_2=46 +x_0=400000 +y_0=400000 +ellps=bessel +step +proj=hgridshift +grids=AT_GIS_GRID.gsb +step +proj=utm +zone=33 +ellps=GRS80
n=1.#INF s=5132938.26037436 w=78264.72048991 e=1.#INF rows=30323 cols=58808
Eingabekarte <dhm_10m_ogd_austria@PERMANENT> in Location <MGI2wgs84_Austria_Lambert>:
(Mon Jan 13 10:19:18 2020) Befehl ausgeführt (0 Sek)

But nothing happend!

And please tell me why do I have to give a pipeline in r.proj? I thought that would defined by the locations resp. regions.

thx!

8

nik wrote

What does r.proj -g output? With PROJ: 6.2.1, a proj pipeline has to be
given in r.proj. it's not in your cmd

Can you report all steps with g.proj, g.region, r.proj in source and
target locations as I've done?

-----
best regards
Helmut
--

g.proj -p
-PROJ_INFO-------------------------------------------------
name : Lambert Conformal Conic
proj : lcc
datum : hermannskogel
ellps : bessel
lat_1 : 49
lat_2 : 46
lat_0 : 47.5
lon_0 : 13.33333333333333
x_0 : 400000
y_0 : 400000
no_defs : defined
towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1

----------------------------------------------------------------------------------------------

g.region -p
projection: 99 (Lambert Conformal Conic)
zone: 0
datum: hermannskogel
ellipsoid: bessel
north: 576922.51207363
south: 273692.51207363
west: 106549.26720377
east: 694629.26720377
nsres: 10
ewres: 10
rows: 30323
cols: 58808
cells: 1783234984

============================================================

g.proj -p
-PROJ_INFO-------------------------------------------------
name : Universal Transverse Mercator
proj : utm
datum : etrs89
ellps : grs80
zone : 33
towgs84 : 0,0,0,0,0,0,0
no_defs : defined
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1

----------------------------------------------------------------------------------------------

g.region -p
projection: 1 (UTM)
zone: 33
datum: etrs89
ellipsoid: grs80
north: 5286581.58847414
south: 5248335.67358037
west: 380870.53610518
east: 429633.13185005
nsres: 9.998932
ewres: 10.00053235
rows: 3825
cols: 4876
cells: 18650700

----------------------------------------------------------------------------------------------

I downloaded and installed
http://www.bev.gv.at/portal/page?_pageid=713,2157075&_dad=portal&_schema=PORTAL

r.proj -g location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
input=dhm_10m_ogd_austria output=dhm_10m_ogd_dachstein_NEU
method=bilinear resolution=10 pipeline=+proj=pipeline +step +inv
+proj=lcc +lat_0=47.5 +lon_0=13.3333333333333 +lat_1=49 +lat_2=46
+x_0=400000 +y_0=400000 +ellps=bessel +step +proj=hgridshift
+grids=AT_GIS_GRID.gsb +step +proj=utm +zone=33 +ellps=GRS80
n=1.#INF s=5132938.26037436 w=78264.72048991 e=1.#INF rows=30323
cols=58808
Eingabekarte &lt;dhm_10m_ogd_austria@PERMANENT&gt; in Location
<MGI2wgs84_Austria_Lambert>
:
(Mon Jan 13 10:19:18 2020) Befehl ausgeführt (0 Sek)

But nothing happend!

And please tell me why do I have to give a pipeline in r.proj? I thought
that would defined by the locations resp. regions.

thx!

for r.proj -g you don't need the pipeline, try just

r.proj -g location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
input=dhm_10m_ogd_austria

see also example https://grass.osgeo.org/grass78/manuals/r.proj.html

-------------------
# same calculation, but in a form which can be cut and pasted into a
g.region call
r.proj input=elevation location=ll_wgs84 mapset=user1 -g
-------------------

the pipeline is needed later for the reprojection process.

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

9

On Mon, Jan 13, 2020 at 1:00 PM Helmut Kudrnovsky <hellik@web.de> wrote:

nik wrote

I downloaded and installed
http://www.bev.gv.at/portal/page?_pageid=713,2157075&_dad=portal&_schema=PORTAL

r.proj -g location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
input=dhm_10m_ogd_austria output=dhm_10m_ogd_dachstein_NEU
method=bilinear resolution=10 pipeline=+proj=pipeline +step +inv
+proj=lcc +lat_0=47.5 +lon_0=13.3333333333333 +lat_1=49 +lat_2=46
+x_0=400000 +y_0=400000 +ellps=bessel +step +proj=hgridshift
+grids=AT_GIS_GRID.gsb +step +proj=utm +zone=33 +ellps=GRS80
n=1.#INF s=5132938.26037436 w=78264.72048991 e=1.#INF rows=30323
cols=58808
Eingabekarte <dhm_10m_ogd_austria@PERMANENT> in Location
<MGI2wgs84_Austria_Lambert>
:
(Mon Jan 13 10:19:18 2020) Befehl ausgeführt (0 Sek)

But nothing happend!

Nik, have you set the current region correctly?

for r.proj -g you don’t need the pipeline, try just

r.proj -g location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
input=dhm_10m_ogd_austria

see also example https://grass.osgeo.org/grass78/manuals/r.proj.html

same calculation, but in a form which can be cut and pasted into a

g.region call
r.proj input=elevation location=ll_wgs84 mapset=user1 -g

the pipeline is needed later for the reprojection process.

A pipeline is not really needed, in case of several possible operations, PROJ should automagically select the best available operation. Note that also with r.proj -g coordinate transformation takes place.

Markus M

10

Am 13.01.2020 um 14:23 schrieb Markus Metz:

Nik, have you set the current region correctly?

>
> for r.proj -g you don't need the pipeline, try just
>
> r.proj -g location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
> input=dhm_10m_ogd_austria
>
> see also example https://grass.osgeo.org/grass78/manuals/r.proj.html
>
> -------------------
> # same calculation, but in a form which can be cut and pasted into a
> g.region call
> r.proj input=elevation location=ll_wgs84 mapset=user1 -g
> -------------------
>
> the pipeline is needed later for the reprojection process.

A pipeline is not really needed, in case of several possible operations, PROJ should automagically select the best available operation. Note that also with r.proj -g coordinate transformation takes place.

thx Markus and Helmut!

Yes that is what I would expect and what I´ve done many times before but since shortly I get this message:

r.proj location=MGI2wgs84_Austria_Lambert mapset=PERMANENT
input=dhm_10m_ogd_austria output=dhm_10m_ogd_dachstein_NEU method=bilinear
resolution=10
ERROR: Unable to open element file <> for <WIND@dachstein>

I tested it from an other location too but it´s the same problem:
r.proj location=ASTER_lat-long mapset=PERMANENT
input=DHM_SRTM30m_Aut-plus-Umgebung output=dhm_30m_ASTER_dachstein
method=bilinear resolution=30
ERROR: Unable to open element file <> for <WIND@dachstein>

What does that mean?

the pojection and region settings of both input and output are als follows:

g.proj -p
-PROJ_INFO-------------------------------------------------
name : Lambert Conformal Conic
proj : lcc
datum : hermannskogel
ellps : bessel
lat_1 : 49
lat_2 : 46
lat_0 : 47.5
lon_0 : 13.33333333333333
x_0 : 400000
y_0 : 400000
no_defs : defined
towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1

----------------------------------------------------------------------------------------------

g.region -p
projection: 99 (Lambert Conformal Conic)
zone: 0
datum: hermannskogel
ellipsoid: bessel
north: 576922.51207363
south: 273692.51207363
west: 106549.26720377
east: 694629.26720377
nsres: 10
ewres: 10
rows: 30323
cols: 58808
cells: 1783234984

============================================================

g.proj -p
-PROJ_INFO-------------------------------------------------
name : Universal Transverse Mercator
proj : utm
datum : etrs89
ellps : grs80
zone : 33
towgs84 : 0,0,0,0,0,0,0
no_defs : defined
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1

----------------------------------------------------------------------------------------------

g.region -p
projection: 1 (UTM)
zone: 33
datum: etrs89
ellipsoid: grs80
north: 5286581.58847414
south: 5248335.67358037
west: 380870.53610518
east: 429633.13185005
nsres: 9.998932
ewres: 10.00053235
rows: 3825
cols: 4876
cells: 18650700

11

On Mon, Jan 13, 2020 at 3:01 PM Nikolai Hafner <n.hafner@carto.at> wrote:

I tested it from an other location too but it´s the same problem:
r.proj location=ASTER_lat-long mapset=PERMANENT
input=DHM_SRTM30m_Aut-plus-Umgebung output=dhm_30m_ASTER_dachstein
method=bilinear resolution=30
ERROR: Unable to open element file <> for WIND@dachstein

What does that mean?

This is a bug in r.proj: the mapset dachstein does not exist in the source location, therefore also no WIND file in the mapset dachstein because this mapset does not exist. The new PR #286 https://github.com/OSGeo/grass/pull/286 fixes this bug, but it exposes another bug in r.import.

Markus M

the pojection and region settings of both input and output are als follows:

g.proj -p
-PROJ_INFO-------------------------------------------------
name : Lambert Conformal Conic
proj : lcc
datum : hermannskogel
ellps : bessel
lat_1 : 49
lat_2 : 46
lat_0 : 47.5
lon_0 : 13.33333333333333
x_0 : 400000
y_0 : 400000
no_defs : defined
towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1

g.region -p
projection: 99 (Lambert Conformal Conic)
zone: 0
datum: hermannskogel
ellipsoid: bessel
north: 576922.51207363
south: 273692.51207363
west: 106549.26720377
east: 694629.26720377
nsres: 10
ewres: 10
rows: 30323
cols: 58808
cells: 1783234984

============================================================

g.proj -p
-PROJ_INFO-------------------------------------------------
name : Universal Transverse Mercator
proj : utm
datum : etrs89
ellps : grs80
zone : 33
towgs84 : 0,0,0,0,0,0,0
no_defs : defined
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1

g.region -p
projection: 1 (UTM)
zone: 33
datum: etrs89
ellipsoid: grs80
north: 5286581.58847414
south: 5248335.67358037
west: 380870.53610518
east: 429633.13185005
nsres: 9.998932
ewres: 10.00053235
rows: 3825
cols: 4876
cells: 18650700

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