Changes between Version 10 and Version 11 of CookBook join ROE_CCM

Timestamp:

Mar 15, 2010 10:19:39 PM (15 years ago)

Author:

cedric

Comment:

--

CookBook join ROE_CCM

@@ -19 +19 @@
 == 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[[BR]]
+In __!QuantumGis__, connect the ROE database and add the obstacle_referentiel layer[[BR]]
+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).[[BR]]
+In __!ArcGis__, add the two layers :[[BR]]
+- ROE_EDA.shape[[BR]]
+- GEOFLADept_FR_Corse_AV_L93\DEPARTEMENT.shp[[BR]]
+Use Analysis Tools > Extract > Clip[[BR]]
+Input Features : ROE_EDA[[BR]]
+Clip Features : DEPARTEMENT[[BR]]
+
+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?).[[BR]]
+
+==== With QGis ====
+You have to reproject the data 
+{{{
+C:\OSGeo4W\bin\ogr2ogr -- help
+}}}
+gets :
+
+{{{
+Usage: ogr2ogr [--help-general] [-skipfailures] [-append] [-update] [-gt n]
+               [-select field_list] [-where restricted_where]
+               [-sql <sql statement>]
+               [-spat xmin ymin xmax ymax] [-preserve_fid] [-fid FID]
+               [-a_srs srs_def] [-t_srs srs_def] [-s_srs srs_def]
+               [-f format_name] [-overwrite] [[-dsco NAME=VALUE] ...]
+               dst_datasource_name src_datasource_name
+               [-lco NAME=VALUE] [-nln name] [-nlt type] [layer [layer ...]]
+
+ -f format_name: output file format name, possible values are:
+     -f "ESRI Shapefile"
+     -f "MapInfo File"
+     -f "TIGER"
+     -f "S57"
+     -f "DGN"
+     -f "Memory"
+     -f "BNA"
+     -f "CSV"
+     -f "GML"
+     -f "GPX"
+     -f "KML"
+     -f "GeoJSON"
+     -f "Interlis 1"
+     -f "Interlis 2"
+     -f "GMT"
+     -f "SQLite"
+     -f "ODBC"
+     -f "PostgreSQL"
+     -f "MySQL"
+ -append: Append to existing layer instead of creating new if it exists
+ -overwrite: delete the output layer and recreate it empty
+ -update: Open existing output datasource in update mode
+ -select field_list: Comma-delimited list of fields from input layer to
+                     copy to the new layer (defaults to all)
+ -where restricted_where: Attribute query (like SQL WHERE)
+ -sql statement: Execute given SQL statement and save result.
+ -skipfailures: skip features or layers that fail to convert
+ -gt n: group n features per transaction (default 200)
+ -spat xmin ymin xmax ymax: spatial query extents
+ -dsco NAME=VALUE: Dataset creation option (format specific)
+ -lco  NAME=VALUE: Layer creation option (format specific)
+ -nln name: Assign an alternate name to the new layer
+ -nlt type: Force a geometry type for new layer.  One of NONE, GEOMETRY,
+      POINT, LINESTRING, POLYGON, GEOMETRYCOLLECTION, MULTIPOINT, MULTILINE,
+      MULTIPOLYGON, or MULTILINESTRING.  Add "25D" for 3D layers.
+      Default is type of source layer.
+ -a_srs srs_def: Assign an output SRS
+ -t_srs srs_def: Reproject/transform to this SRS on output
+ -s_srs srs_def: Override source SRS
+
+ Srs_def can be a full WKT definition (hard to escape properly),
+ or a well known definition (ie. EPSG:4326) or a file with a WKT
+ definition.
+}}}
+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 [[BR]]
+Process 
+Add the file GEOFLADept_FR_Corse_AV_L93\DEPARTEMENT.shp (layer>add a vector layer) [[BR]]
+Tools > Geoprocessing tools > Clip
 == Build a view allowing to know the score of the different dams ==
 this is ticket #46