[GRASS-user] r.proj ERROR

[Helmut Kudrnovsky]

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 <dhm_10m_ogd_austria at 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!

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