====== 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 =====
* http://georezo.net/forum/viewtopic.php?pid=169094#p169094
* http://www.portailsig.org/content/grass-gis-region-resolution-resolution-3d-drapage-ou-l-apparente-quadrature-du-cercle
* [[http://lambert93.ign.fr|La gestion du lambert 93 et des reprojections]]
* http://e-geo.fr/mastersig/wiki/doku.php/fonctions:visualisation:3d:grass
===== 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 ======
[[http://grass.itc.it/grass64/manuals/html64_user/r.reclass.html|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
[[grass:cas_d_utilisation:diagnostic]]
====== Diagnostics topo ======
* http://grass.osgeo.org/gdp/html_grass64/r.slope.aspect.html
* r.rescale