TransWikia.com

Executing Python script for GRASS from PHP

Geographic Information Systems Asked by Dawid on January 5, 2021

I wrote the following Python script:

#!/usr/bin/env python
import os
import sys

gisbase = os.environ['GISBASE'] = "/usr/lib/grass64"
gisdbase=os.path.join("PATH")
location=LOCATION NAME
mapset=MAPSETNAME
sys.path.append(os.path.join(os.environ['GISBASE'], "etc", "python"))

import grass.script as grass
import grass.script.setup as gsetup
import grass.script as grass
 
gsetup.init(gisbase, gisdbase, location, mapset)
print "one"

grass.run_command('v.in.ascii','-o', input='/opt/lampp/htdocs/project/basic_statistics/sites_grass.csv', output='sites',fs=",",skip='1', x='9', y='10')
print "two"

grass.run_command('v.db.addcol', map='sites', columns="height INTEGER")
grass.run_command('g.region', rast="dem21")
grass.run_command('v.what.rast',vector='sites', raster='dem21', column='height')
print "three"

grass.run_command('v.out.ascii','-o', input="sites",output="sites_geography_stats.csv", columns='height')
print "four"

It runs without any problems when I simply execute it from UNIX terminal (/.FILENAME.py). However, my goal is to execute it from a PHP code:

<?php
...
echo system("./FILENAME.py");
...
?>

When I run the PHP code the following output is printed to the website:
one two three four

Sadly, the sites_geography_stats.csv file is not created. It looks that the python script is executed from php but its GRASS-related part is ignored.
Does anybody have any ideas how to fix it?

One Answer

Edited:

I have tried the script simply on command line as python script and found the errors:

grass.run_command('v.in.ascii','-o', ...

should be

grass.run_command('v.in.ascii', overwrite = True, ...

Generally flags are indicated like this:

flags = 'r'

but in case of the overwrite (--o) this does not apply due to the doubled dash.

Likewise

grass.run_command('v.out.ascii','-o',

should be

grass.run_command('v.out.ascii', overwrite = True,

Eventually you should not use integer but double precision to avoid "WARNING: Raster type is float and column type is integer, some data lost!!" in the v.what.rast call.

With these changes the script works properly with GRASS 6.

BTW: See here for a sample Python script.

Answered by markusN on January 5, 2021

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP