Présentation du

Système d'Information outils:grass

outils:grass:cas_d_utilisation

Wed, 18 May 2011 14:38:17 +0200

Grille statistique à partir du MNT

Création de la grille (100x100m) pour le Département de l'Aude

v.mkgrid map=grille_lr_100m grid=2600,3000 --o position=coor coor=588600,6137100 box=100,100

Comment réaliser une vue 3D de la RNR

Données disponibles

topo

points xyz, formats dxf, lambert 3 sud
MNT IGN BD TOPO lambert 2 etendu

fonds de référence

SCAN 25 (tif) en lambert 93
orthophoto IGN (ecw) en lambert 93
orthophoto API (ecw) resolution de 20cm, lambert 3 sud

données vecteur

BD PArcellaire Lambert 93

Ressources en ligne

Procédure pour import dans un secteur en Lambert 2 étendu

r.in.xyz -s in=D:/FORMATIONS/GRASS/donnees_RNRGG/API_MNT/points_l2c.txt fs=; out=test
# x: 764157.410000 769484.330000
# y: 1882750.910000 1885657.260000
# z:   28.421000  204.329000
# Range:     min         max
 
g.region n=1885657.260000 s=1882750.910000 w=764157.410000 e=769484.330000 res=20 -p
g.region n=n+10 s=s-10 w=w-10 e=e+10 -p
wc -l D:/FORMATIONS/GRASS/donnees_RNRGG/API_MNT/points_l2c.txt
 
r.in.xyz input=D:/FORMATIONS/GRASS/donnees_RNRGG/API_MNT/points_l2c.txt output=mnt_api method=mean type=FCELL fs=; x=1 y=2 z=3 zscale=1.0 percent=100 --overwrite

Probleme de définition du drapé : recréation d'un mnt + fin et d'un scan + fin

r.resamp.rst input=mnt_20cm@gardon ew_res=1 ns_res=1 elev=mnt_1 overlap=3 zmult=1.0 tension=40.

Procédure pour import dans un secteur en Lambert 93

r.in.xyz -s in=D:/points_l93.txt fs=; out=test
# x: 810732.480000 816053.130000
# y: 6315171.660000 6318057.890000
# z:   28.421000  204.329000
# Range:     min         max
 
g.region n=6318057.890000 s=6315171.660000 w=810732.480000 e=816053.130000 res=20 -p
g.region n=n+10 s=s-10 w=w-10 e=e+10 -p
wc -l D:/points_l93.txt
 
r.in.xyz input=D:/points_l93.txt output=mnt_api method=mean type=FCELL fs=; x=1 y=2 z=3 zscale=1.0 percent=100 --overwrite

Probleme de définition du drapé : recréation d'un mnt + fin et d'un scan + fin

r.resamp.rst input=mnt_20cm@gardon ew_res=1 ns_res=1 elev=mnt_1 overlap=3 zmult=1.0 tension=40.

Réalisation d'une carte ombrée

r.shaded.relief map=mnt_250m@region_lr shadedmap=ombre_lr_275_exag_4 altitude=60 azimuth=275 zmult=4 scale=1 units=none --overwrite
r.shaded.relief map=mnt_250m@region_lr shadedmap=ombre_lr_315_exag_4 altitude=60 azimuth=315 zmult=4 scale=1 units=none --overwrite
r.shaded.relief map=mnt_250m@region_lr shadedmap=ombre_lr_355_exag_4 altitude=60 azimuth=355 zmult=4 scale=1 units=none --overwrite

r.out.png input=ombre_lr_275_exag_4@region_lr output=ombre_lr_275_exag_4.png 
r.out.png input=ombre_lr_315_exag_4@region_lr output=ombre_lr_315_exag_4.png 
r.out.png input=ombre_lr_355_exag_4@region_lr output=ombre_lr_355_exag_4.png

Changer les couleurs du MNT en modifiant les classes

r.reclass input=mnt_250m@region_lr output=mnt_0_200_m 'rules=C:/Documents and Settings/Mathieu/Bureau/classe.txt' –overwrite””

contenu du ficher classe.txt

0 thru 5 = 1
6 thru 10 = 2
11 thru 15 = 3
16 thru 20 = 4
21 thru 25 = 5
26 thru 30 = 6
31 thru 35 = 7
36 thru 40 = 8
41 thru 45 = 9
46 thru 50 = 10
51 thru 55 = 11
56 thru 60 = 12
61 thru 65 = 13
66 thru 70 = 14
71 thru 75 = 15
76 thru 80 = 16
81 thru 85 = 17
86 thru 90 = 18
91 thru 95 = 19
96 thru 100 = 20
101 thru 105 = 21
106 thru 110 = 22
111 thru 115 = 23
116 thru 120 = 24
121 thru 125 = 25
126 thru 130 = 26
131 thru 135 = 27
136 thru 140 = 28
141 thru 145 = 29
146 thru 150 = 30
151 thru 155 = 31
156 thru 160 = 32
161 thru 165 = 33
166 thru 170 = 34
171 thru 175 = 35
176 thru 180 = 36
181 thru 185 = 37
186 thru 190 = 38
191 thru 195 = 39
196 thru 200 = 40
* = NULL

Créer des courbes de niveaux

Création des courbes

r.contour -q -n input=mnt_250m@region_lr output=courbes_de_niveau minlevel=0 maxlevel=200 step=25 cut=0 --overwrite

Export vers un shapefile

v.out.ogr input=courbes_de_niveau@region_lr type=line,boundary dsn=courbes_0_200 layer=1 format=ESRI_Shapefile lco=\"SHPT=ARCZ\"

Fusionner des ratsers

Fusion des MNT départementaux dans un mnt régional

r.patch input=mnt_11@PERMANENT,mnt_30@PERMANENT,mnt_34@PERMANENT,mnt_48@PERMANENT,mnt_66@PERMANENT output=mnt_lr --overwrite

Extraire une partie d'un raster

gdal_translate -projwin 814586 6352530 843101 6323630 DEPT30.ASC chusclan_laudun.ASC

diagnostic

Diagnostics topo

http://grass.osgeo.org/gdp/html_grass64/r.slope.aspect.html
r.rescale

outils:grass:formation

Tue, 25 May 2010 22:58:59 +0200

Formation GRASS d'avril 2009

Attentes

Prise en main du logiciel
Mise en place d'un SIG sous Grass
Utilisation du logiciel au retour de formation
- Opérations traditionnelles
  - Requêtes spatiales
  - digitalisation
  - verification de topologies
  - reprojection de données
  - import/export de données
  - Lien avec PostgreSQL

Cas concrets

Vendres

cartes d'inondation du site en fonction d'une échelle limni lue
source unique à considérée : eaux de pluie

Gardon

Profils topo d'itinéraires
http://casoilresource.lawr.ucdavis.edu/drupal/node/375
Coupe transversale perpendiculaire aux gorges (pour carto des habitats par exempe; ou stratification des cavités)
Bloc Diagramme
Drapage de l'ortho sur le MNT
Simulation de vol ou de randonnée
Calculs de pente, d'exposition

r.IN.xyz -s IN=D:/FORMATIONS/GRASS/donnees_RNRGG/API_MNT/points_l2c.txt fs=; out=test
# x: 764157.410000 769484.330000
# y: 1882750.910000 1885657.260000
# z:   28.421000  204.329000
# Range:     min         max
 
g.region n=1885657.260000 s=1882750.910000 w=764157.410000 e=769484.330000 res=20 -p
g.region n=n+10 s=s-10 w=w-10 e=e+10 -p
wc -l D:/FORMATIONS/GRASS/donnees_RNRGG/API_MNT/points_l2c.txt
 
r.IN.xyz input=D:/FORMATIONS/GRASS/donnees_RNRGG/API_MNT/points_l2c.txt output=mnt_20cm method=mean type=FCELL fs=; x=1 y=2 z=3 zscale=1.0 percent=100 --overwrite

Problème de définition du drapé : ressource en ligne recréation d'un mnt + fin et d'un scan + fin

r.resamp.rst input=mnt_20cm@gardon ew_res=1 ns_res=1 elev=mnt_1 overlap=3 zmult=1.0 tension=40.

outils:grass:requetes_spatiales

Tue, 25 May 2010 22:56:29 +0200

Sélection des villes de Midi-Pyrénées qu'ont met dans une nouvelle couche

seul l'opérateur overlap est disponible actuellement

v.SELECT ainput=villes@PERMANENT atype=point,line,boundary,centroid,area alayer=1 binput=departements_mp@mb btype=point,line,boundary,centroid,area blayer=1 output=villes_mp operator=overlap