Greetings,
Just FYI - I imported an image into it’s own location with this command:
r.in.gdal input=m_4511664.jp2 output=m_4511664 location=m_4511664
In the new location:
g.region rast= m_4511664; v.in.region output=bounds
When I tried re-projecting the from the created location, I got this error:
WARNING: pj_transform() failed: failed to load datum shift file
ERROR: Unable to re-project vector map <bounds@PERMANENT> from
<m_4511664>
The cause of the error was this line in the PROJ_INFO file in the location created by r.in.gdal (first command shown in this email):
nadgrids: @null
After removing that line, reprojection works.
Cheers,
Jeshua Lacock
Founder/Engineer
<3DTOPO.com>
GlassPrinted.com
On Thu, Jul 20, 2017 at 3:05 AM, Jeshua Lacock <jeshua@3dtopo.com> wrote:
Greetings,
Just FYI - I imported an image into it’s own location with this command:
r.in.gdal input=m_4511664.jp2 output=m_4511664 location=m_4511664
In the new location:
g.region rast= m_4511664; v.in.region output=bounds
When I tried re-projecting the from the created location, I got this error:
WARNING: pj_transform() failed: failed to load datum shift file
ERROR: Unable to re-project vector map <bounds@PERMANENT> from
<m_4511664>
The cause of the error was this line in the PROJ_INFO file in the location created by r.in.gdal (first command shown in this email):
nadgrids: @null
After removing that line, reprojection works.
Probably it is this ticket?
https://trac.osgeo.org/grass/ticket/2229
Best,
Markus
On 20/07/17 12:11, Markus Neteler wrote:
On Thu, Jul 20, 2017 at 3:05 AM, Jeshua Lacock <jeshua@3dtopo.com> wrote:
Greetings,
Just FYI - I imported an image into it’s own location with this command:
r.in.gdal input=m_4511664.jp2 output=m_4511664 location=m_4511664
In the new location:
g.region rast= m_4511664; v.in.region output=bounds
When I tried re-projecting the from the created location, I got this error:
WARNING: pj_transform() failed: failed to load datum shift file
ERROR: Unable to re-project vector map <bounds@PERMANENT> from
<m_4511664>
The cause of the error was this line in the PROJ_INFO file in the location created by r.in.gdal (first command shown in this email):
nadgrids: @null
After removing that line, reprojection works.
Probably it is this ticket?
https://trac.osgeo.org/grass/ticket/2229
The error is also described in this ticket, but the issue seems to be the difficulty handling the 'nadgrids: @null' hack GDAL uses (IIUC).
See also [1].
Maybe we can create a workaround for this specific issue in v./r.proj ?
Moritz
[1] https://lists.osgeo.org/pipermail/grass-user/2017-June/076709.html
On Jul 20, 2017, at 4:11 AM, Markus Neteler <neteler@osgeo.org> wrote:
Probably it is this ticket?
https://trac.osgeo.org/grass/ticket/2229
Hi Markus,
Yes; it looks like the same issue.
When I googled for the issue, I came across an older different ticket for the same "failed to load datum shift file” error but I didn’t see that one.
Cheers,
Jeshua Lacock
Founder/Engineer
<3DTOPO.com>
GlassPrinted.com