[GRASS-user] command works inside GRASS but not with python script outside GRASS
alassane toure
atourej at gmail.com
Sat Jul 2 18:02:22 PDT 2016
Group,
The following command works using GRASS command session:
v.extract input=const output=const_clean type=area where="(value=255) and
(area_sqf>1000)"
Extracting features...
Building topology for vector map <const_cl at PERMANENT>...
Registering primitives...
7091 primitives registered
54239 vertices registered
Building areas...
2484 areas built
1143 isles built
Attaching islands...
Attaching centroids...
Number of nodes: 5334
Number of primitives: 7091
Number of points: 0
Number of lines: 0
Number of boundaries: 6675
Number of centroids: 416
Number of areas: 2484
Number of isles: 1143
Writing attributes...
But does not work in grass/python script:
grass.run_command('v.extract', input='const', output='const_clean',
type='area',where='"(value=255) and (area_sqf>1000)"' )
....
Extracting features...
100%
Building topology for vector map <const_clean at PERMANENT>...
Registering primitives...
0 primitives registered
0 vertices registered
Building areas...
100%
0 areas built
0 isles built
Attaching islands...
Attaching centroids...
Number of nodes: 0
Number of primitives: 0
Number of points: 0
Number of lines: 0
Number of boundaries: 0
Number of centroids: 0
Number of areas: 0
Number of isles: 0
The GRASS script contains the following imports
#!/usr/bin/env python
import os
import sys
import subprocess
import shutil
import binascii
import tempfile
import shutil
import grass.script as grass
import grass.script.setup as gsetup
Your help is appreciated.
Alassane
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20160702/1664eb96/attachment.html>
More information about the grass-user
mailing list