thanks for the answers, but this did not work.

Again, here is my code, this time not just the minimum version:

#!/usr/bin/env python
import os
from grass_session import Session
from grass.script import core as gcore

from grass.pygrass.modules.shortcuts import general as g
from grass.pygrass.modules.shortcuts import raster as r

home= "/home/jreith/grassdata/nrw"
gisdb= "/home/jreith/grassdata"
loc= "nrw"
maps= "elevation"
epsg= "EPSG:4647"
xyz= "/home/jreith/geodaten/dgm1_05314000_Bonn_EPSG4647_XYZ/tmp.xyz"

if not os.path.exists(gisdb):
    os.makedirs(gisdb)

if not os.path.exists(os.path.join(gisdb, loc)):
    with Session(gisdb=gisdb, location=loc, create_opts=epsg):
        print("Created a new location!")
else:
    print("Location already exist!")

with open(xyz, "r") as t:
    print(t.read())

with Session(gisdb=gisdb, location=loc, mapset="elevation", create_opts=""):
    print("r.in_xyz works:")
    r.in_xyz(input=xyz, output="new_file", flags="s",separator="space")
    gisenvironment = gcore.parse_command("g.gisenv", flags="s")
    print(gisenvironment)
    compregion = gcore.parse_command('r.in_xyz', input=xyz, flags='s',
output="new_file",separator='space')
##################*
*

*This is my output:*

GRASSBIN: grass74
GISBASE: /home/jreith/source/grass-7.4.svn/dist.x86_64-pc-linux-gnu
Location already exist!
32360000.00 5610000.00 170.55
32360000.00 5610001.00 170.56
32360000.00 5610002.00 170.58
32360000.00 5610003.00 170.55
32360000.00 5610004.00 170.58



r.in_xyz works:
Range:     min         max
x:    32360000    32360000
y:     5610000     5610004
z:      170.55      170.58
{u'MAPSET': u"'elevation';", u'GISDBASE': u"'/home/jreith/grassdata';",
u'LOCATION_NAME': u"'nrw';"}
Traceback (most recent call last):
  File "grass_scripts/error_script.py", line 33, in <module>
    compregion = gcore.parse_command('r.in_xyz', input=xyz, flags='s',
output="new_file",separator='space')
  File
"/home/jreith/source/grass-7.4.svn/dist.x86_64-pc-linux-gnu/etc/python/grass/script/core.py",
line 516, in parse_command
    res = read_command(*args, **kwargs)
  File
"/home/jreith/source/grass-7.4.svn/dist.x86_64-pc-linux-gnu/etc/python/grass/script/core.py",
line 471, in read_command
    process = pipe_command(*args, **kwargs)
  File
"/home/jreith/source/grass-7.4.svn/dist.x86_64-pc-linux-gnu/etc/python/grass/script/core.py",
line 444, in pipe_command
    return start_command(*args, **kwargs)
  File
"/home/jreith/source/grass-7.4.svn/dist.x86_64-pc-linux-gnu/etc/python/grass/script/core.py",
line 380, in start_command
    return Popen(args, **popts)
  File
"/home/jreith/source/grass-7.4.svn/dist.x86_64-pc-linux-gnu/etc/python/grass/script/core.py",
line 74, in __init__
    subprocess.Popen.__init__(self, args, **kwargs)
  File "/usr/lib/python2.7/subprocess.py", line 711, in __init__
    errread, errwrite)
  File "/usr/lib/python2.7/subprocess.py", line 1343, in _execute_child
    raise child_exception
OSError: [Errno 2] No such file or directory
###########################

I used a minimum example for the xyz-file. As you can see, there is no
problem with the location or database.

r.in_xyz is correct and also works itself. I also can parse the g.gisenv
to a variable and print it.

The problem is, that I cannot parse r.in_xyz to a new variable.

Maybe you have some hints for me :)

best regards

Jonathan


Am 12.02.2018 um 15:35 schrieb Pietro:
> Dear Jonathan,
>
> the error is due to grass_session that is not creating the location if
> missing.
> I don't have time in this day to fix this issue in grass_session, so
> the fastest fsolution at the momenth is to check and create what is
> needed step by step.
>
> I did not have xyz file to test so I've only execute g.gisenv and it
> works, let me know if it works also with r.inxyz:
>
> ```python
> from __future__ import print_function
> import os
> from grass_session import Session
> from grass.script import core as gcore
>
> GISDBASE = "/tmp/grassdata"
> LOCATION = "nrw"
> EPSG = "EPSG:4326"
>
>
> if not os.path.exists(GISDBASE):
>     os.makedirs(GISDBASE)
>
> if not os.path.exists(os.path.join(GISDBASE, LOCATION)):
>     with Session(gisdb=GISDBASE, location=LOCATION,
>                  create_opts=EPSG):
>         print("Created a new location!")
> else:
>     print("Location already exist!")
>
>
> with Session(gisdb=GISDBASE, location=LOCATION, mapset="elevation",
>              create_opts=""):
>     gcore.run_command("g.gisenv")
> ```
>
> Best regards
>
> Pietro

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to