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@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@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 _______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
