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
./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
Commentaires