wiki:CookBook join ROE_CCM

Version 14 (modified by cedric, 15 years ago) (diff)

--

back to first page ..
back to CookBook Eda

Join ROE - CCM

Puttin geobs into the CCM

Initially, the database was located in another table, below we just copy the table shema roe into the ccm21_eda database

REM pour Cédric
CD C:\Documents and Settings\cedric\Mes documents\Migrateur\eda\Barrages
C:\"Program Files"\PostgreSQL\8.4\bin\pg_dump -U postgres -n geobs ROE> roe_schema.sql
C:\"Program Files"\PostgreSQL\8.4\bin\psql -d ccm21_eda -h localhost -p 5433 -U postgres -f roe_schema.sql
REM====================================
REM Pour Céline
CD D:\Celine Jouanin\ROE
C:\"Program Files"\PostgreSQL\8.4\bin\pg_dump -U postgres  -n geobs ROE> roe_schema.sql
C:\"Program Files"\PostgreSQL\8.4\bin\psql -d ccm21_eda -h localhost -p 5432 -U postgres -f roe_schema.sql

We have to reproject our data

Select only the dams within France

Extract the data within the french boundaries

this is ticket #45

With ArcGis

Use the french limits (see French boundaries) in RGF93
In !QuantumGis, connect the ROE database and add the obstacle_referentiel layer
Clic right on the layer "Sauvegarder comme shapefile" with Fichiers de type : Shapefiles (*.shp), Codage=UTF-8 > ok and choose "Système de coordonnées de référence" = RGF93/Lambert93 (because the french boundary layer is in this system).
In !ArcGis, add the two layers :

  • ROE_EDA.shape
  • GEOFLADept_FR_Corse_AV_L93\DEPARTEMENT.shp

Use Analysis Tools > Extract > Clip
Input Features : ROE_EDA
Clip Features : DEPARTEMENT

If you want you can choose a "XY Tolerance", because some points in the ROE are in France but they are at the limit. (maybe 300meters?).

With QGis

first change projection to ERTRS 89 write ogr2ogr -f formatdst_datasource_name src_datasource_name -a_srs <actual projection> -t_srs <reprojected data> Launch in OSGeo4W command shell

CD C:\Documents and Settings\cedric\Mes documents\Migrateur\eda\couches_SIG\france
C:\OSGeo4W\bin\ogr2ogr DEPARTEMENT_3035.shp DEPARTEMENT.shp -s_srs "EPSG:2154" -t_srs "EPSG:3035"

The shape file DEPARTEMENT_3035.shp is created with right projection Connect the ROE database and add this newly projected obstacle_referential layer
Process Add the file GEOFLADept_FR_Corse_AV_L93\DEPARTEMENT.shp (layer>add a vector layer)
Tools > Geoprocessing tools > Clip

Build a view allowing to know the score of the different dams

this is ticket #46

Joining CCM and Dams

this is ticket #47