<div dir="ltr"><div><div><div><div><div><div>Group,<br></div>The following command works using GRASS command session:<br><br></div>v.extract input=const output=const_clean type=area where="(value=255) and (area_sqf>1000)"<br>Extracting features...<br>Building topology for vector map <const_cl@PERMANENT>...<br>Registering primitives...<br>7091 primitives registered<br>54239 vertices registered<br>Building areas...<br>2484 areas built<br>1143 isles built<br>Attaching islands...<br>Attaching centroids...<br>Number of nodes: 5334<br>Number of primitives: 7091<br>Number of points: 0<br>Number of lines: 0<br>Number of boundaries: 6675<br>Number of centroids: 416<br>Number of areas: 2484<br>Number of isles: 1143<br>Writing attributes...<br><br><br><br></div>But does not work in grass/python script:<br>grass.run_command('v.extract', input='const', output='const_clean', type='area',where='"(value=255) and (area_sqf>1000)"' ) <br><br>....<br>Extracting features...<br> 100%<br>Building topology for vector map <const_clean@PERMANENT>...<br>Registering primitives...<br>0 primitives registered<br>0 vertices registered<br>Building areas...<br> 100%<br>0 areas built<br>0 isles built<br>Attaching islands...<br>Attaching centroids...<br>Number of nodes: 0<br>Number of primitives: 0<br>Number of points: 0<br>Number of lines: 0<br>Number of boundaries: 0<br>Number of centroids: 0<br>Number of areas: 0<br>Number of isles: 0<br><br><br></div>The GRASS script contains the following imports<br>#!/usr/bin/env python<br>import os<br>import sys<br>import subprocess<br>import shutil<br>import binascii<br>import tempfile<br>import shutil <br><br>import grass.script as grass<br>import grass.script.setup as gsetup<br><br></div>Your help is appreciated.<br><br></div>Alassane<br><div><div><br><div><br></div></div></div></div>