[GRASS-user] r.proj ERROR
Helmut Kudrnovsky
hellik at web.de
Sat Jan 11 08:55:45 PST 2020
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 <WIND at 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
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 ;-)
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 at 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
More information about the grass-user
mailing list