back to first page[..][[BR]]

== Récupération de la surface en eau (couche surfacique de Bdcarthage) sur le RHT ==

'''Calcul de la surface en eau de la couche surfacique'''[[BR]]
Plan d'eau, bassin, réservoir  + Ensemble de petits plans d'eau +Ecoulement d'eau = 3431.711333+15.37315377+173.866688[[BR]]
Soit 3620.95 km² (dans PGA : 3637km²)[[BR]]
{{{
#!sql
select * from spatial_ref_sys where srid=3035;  ---unit en mètres
alter table bd_carthage2011.hydrographie_surfacique add column surface numeric;
update bd_carthage2011.hydrographie_surfacique set surface=sub.surface from 
(select id_bdcarth, area2d(the_geom) as surface from bd_carthage2011.hydrographie_surfacique) as sub
where sub.id_bdcarth=hydrographie_surfacique.id_bdcarth;

COMMENT ON COLUMN bd_carthage2011.hydrographie_surfacique.surface IS 'Surface en m² calculée avec requête spatiale area2d';

select sum(area2d(the_geom)), type from bd_carthage2011.hydrographie_surfacique group by type;
}}}

'''Création d'un buffer 100m pour récupérer les tronçons inclus dans la couche surfacique bdcarthage'''
{{{
#!sql
drop table if exists bd_carthage2011.hydrographie_surfacique_100;
create table bd_carthage2011.hydrographie_surfacique_100 as 
	select id_bdcarth, nature, type, toponyme, surface, st_buffer(the_geom,100) as the_geom from bd_carthage2011.hydrographie_surfacique;
CREATE INDEX indexhydro_100 ON bd_carthage2011.hydrographie_surfacique_100
  USING GIST ( the_geom GIST_GEOMETRY_OPS );
}}}

'''Jointure couche surfacique bdcarthage - RHT'''
{{{
#!sql
drop table if exists rht.rhtvs2_hydrographie_surfacique;
create table rht.rhtvs2_hydrographie_surfacique as select id_drain, id_bdcarth, surface, type, nature from rht.rhtvs2 r inner join bd_carthage2011.hydrographie_surfacique_100 bd on st_contains(bd.the_geom, 
r.the_geom) order by id_bdcarth;

alter table rht.rhtvs2_hydrographie_surfacique add column a_conserver boolean default TRUE;
update rht.rhtvs2_hydrographie_surfacique set a_conserver=FALSE where type in ('Sable humide','Zone rocheuse','Pleine mer','Traitement des eaux','Bassin portuaire fluvial','Bassin portuaire',
'Gravier, galet','Inconnu','Sable, gravier');

alter table rht.rhtvs2_hydrographie_surfacique rename column "type" to type;

select id_drain, id_bdcarth from rht.rhtvs2 r inner join bd_carthage2011.hydrographie_surfacique_100 bd on st_contains(bd.the_geom, 
r.the_geom) order by id_bdcarth;

select id_drain, id_bdcarth from rht.rhtvs2 r inner join bd_carthage2011.hydrographie_surfacique_100 bd on st_within(r.the_geom, 
bd.the_geom) order by id_bdcarth;
}}}

{{{
select * from bd_carthage2011.hydrographie_surfacique where type='Nappe d''eau' limit 4;
}}}

La couche surfacique étant relative à plusieurs tronçons, la densité moyenne pour chaque tronçon sera utilisée.[[BR]]

||Type||Nature||a conserver||
||Cours d'eau||Eau douce permanente||OUI||
||Plan d'eau, bassin, réservoir||Eau douce permanente||OUI||
||Ensemble de petits plans d'eau||Eau douce permanente||OUI||
||Bassin portuaire fluvial||Eau douce permanente||??||
||Traitement des eaux||Eau douce permanente||NON||
||Zone recouverte d'eau||Eau douce non permanente||??||
||Sable, gravier||Eau douce non permanente||NON||
||Sable humide||Eau Salée non permanente||NON||
||Marais salant||Eau Salée non permanente||??||
||Vase||Eau Salée non permanente||OUI||
||Gravier, galet||Eau Salée non permanente||NON||
||Zone rocheuse||Eau Salée non permanente||NON||
||Rocher, sable||Eau Salée non permanente||NON||
||Nappe d'eau||Eau salée permanente||NON||
||Ecoulement d'eau||Eau salée permanente||??||
||Bassin portuaire||Eau salée permanente||||
||Pleine mer||Eau salée permanente||NON||
||Inconnu||Névé, glacier||NON||