[GRASS-user] r.proj ERROR

[Markus Metz]

On Mon, Jan 13, 2020 at 1:00 PM Helmut Kudrnovsky <hellik at 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 at 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20200113/e661b86a/attachment.html>