On Mon, Feb 19, 2018 at 4:36 AM, Jonathan Reith <reith@mundialis.de> wrote:
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')
##################
You have to use r.in.xyz with gcore.parse_command as I suggested
earlier. You are now using 2 different APIs, PyGRASS (where you use
r.in_xyz) and Python Scripting Library (where you must use the real
name of the command).
Anna
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:
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