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_surface in ('Pleine mer','Zone rocheuse','Rocher, sable','Inconnu','Sable, gravier','Sable humide','Gravier, galet','Traitement des eaux');

alter table rht.rhtvs2_hydrographie_surfacique rename column "type" to type;
select id_drain, sum(surface)as surface , count(id_drain) from rht.rhtvs2_hydrographie_surfacique where a_conserver=TRUE group by id_drain order by count(id_drain);


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 --> on ne prend pas donc pas l'estuaire de la Gironde ?||
||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||