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