recuperer les z d'un batiment sur un ascii grid dans postgis

première étape charger un mns dans postgis:

 ./raster2pgsql -I -C -r -t 50x50 -d -s 21781 /home/xao/postgis/MNS_zone_test.asc mns >/tmp/mns.sql

puis on insère dans notre base de données: psql -U postgres -f /tmp/mns.sql gis01 enfin on peut créer une table temporaire avec une séquence

CREATE SEQUENCE seq_bati_isolines_id MINVALUE 1 START WITH 1 INCREMENT BY 1 CACHE 10

CREATE TABLE bati_isolines as
 (SELECT nextval('seq_bati_isolines_id') as id,
  no_eca,
  (gv).geom,
  (gv).val as z
 FROM (SELECT b.no_eca , ST_Intersection(geom,rast) gv
  FROM bati_pol b, mns
  WHERE no_eca like '11330'
  AND st_intersects(geom,rast)) as t
)

une fois ce petit détail reglé on peut allez voir nos isolignes dans QGIS avec le plugin "DBManager" avec la requête SELECT id,geom,z FROM bati_isolines