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

outils/grass/cas_d_utilisation.txt · Dernière modification: 2011/05/18 14:38 par hassan
www.chimeric.de Creative Commons License Valid CSS Propulsé par DokuWiki Get firefox!! Changements récents - flux RSS Valid XHTML 1.0 Hébergé par Alwaysdata