Le portail francophone de la géomatique

Bonjour à tous,

Dans un module écrit en PYTHON, je ne parviens pas à accéder à la géometrie de couches vecteur :

SQL_req = 'SELECT id_parcel, id_melba FROM parcels, network WHERE st_intersects(parcels.geometry, network.geometry)'
grass.write_command("db.select", input='-', stdin=SQL_req, stdout = DEVNULL, stderr = DEVNULL)


me retourne :

DBMI-SQLite driver error:
Error in sqlite3_prepare():
SELECT id_parcel, id_melba FROM parcels, network WHERE st_intersects(parcels.geometry, network.geometry)
no such column: parcels.geometry

Les shapefiles ont été importés avec v.in.ogr et d'autres requêtes sont exécutées correctement. C'est quand je cherche à faire des requêtes spatiales que ça ne marche plus.

Cette requête passée dans QGIS fonctionne bien.

Avez-vous une idée?

Merci,
Christine