or
or
By clicking below, you agree to our terms of service.
New to HackMD? Sign up
Syntax | Example | Reference | |
---|---|---|---|
# Header | Header | 基本排版 | |
- Unordered List |
|
||
1. Ordered List |
|
||
- [ ] Todo List |
|
||
> Blockquote | Blockquote |
||
**Bold font** | Bold font | ||
*Italics font* | Italics font | ||
~~Strikethrough~~ | |||
19^th^ | 19th | ||
H~2~O | H2O | ||
++Inserted text++ | Inserted text | ||
==Marked text== | Marked text | ||
[link text](https:// "title") | Link | ||
 | Image | ||
`Code` | Code |
在筆記中貼入程式碼 | |
```javascript var i = 0; ``` |
|
||
:smile: | ![]() |
Emoji list | |
{%youtube youtube_id %} | Externals | ||
$L^aT_eX$ | LaTeX | ||
:::info This is a alert area. ::: |
This is a alert area. |
On a scale of 0-10, how likely is it that you would recommend HackMD to your friends, family or business associates?
Please give us some advice and help us improve HackMD.
Do you want to remove this version name and description?
Syncing
xxxxxxxxxx
Notes / Mémo pour les enseignants
Liens utiles
Diagramme de Flux
Rusle
Version TD
Modifications à faire pour l'application sur la save:
Version Basin Versant de la Save
Détails
Formule plus générale pour facteur R :
numpy.where(A<=850,(0.04830*(A**1.610)),(587.8-(1.249*A)+(0.004105*(A**2))))
Formule plus générale pour facteur C :
Données
Note certaines données peuvent en fait être dans les gpkg → ne pas faire attention aux formats qui distinguent ici raster / vecteur
Modèle builder
Fichier à utiliser
Versions :
Inrae :
Modele_erosion_INRAe_v3_saga7.3.0.model3
Modele_erosion_INRAe_v3_saga7.3.0_adaptee.model3
Modele_erosion_INRAe_v3_saga7.8.2.model3
Rusle :
Modele_RUSLE_INRAe_v3_saga7.3.0.model3
Modele_RUSLE_INRAe_v3_saga7.3.0_adaptee.model3
Modele_RUSLE_INRAe_v3_saga7.8.2.model3
La version adaptée à la Save n'a pas été faite pour la 7.8.2 mais il n'y a que deux traitements à changer a priori.
Compatibilité entre les versions de SAGA
Les algos
fill sink
etcatchment area
sont dispos dans la version7.3.0
mais pas dans la version7.8.2
(celle du serveur). J'ai trouvé leurs nouveaux noms :fill sink
fill sink (Planchon/Darboux, 2001)
catchment area
Accumulation de flux (haut vers le bas)
Des nouveaux paramètres sont à renseigner (ils étaient renseignés par défaut avant, de manière transparente).
Dans la
7.8.2
, pourAccumulation de flux (haut vers le bas)
, il faut changer les valeurs par défaut des paramètres.Thresholded Linear Flow
:true
Target accumalation
, je ne suis pas sûr de la signification de ce paramètre, j'ai mis le mnt "filled".On obtient alors bien les mêmes résultats qu'avec la version
7.3.0
Avec le docker
Attention, après plusieurs tests :
catchment area
ne fonctionne pas, en ayant testé les options :method
Thresholdd Linear Flow
catchment area (flow tracing)
fonctionne, bien que la méthode D8 qu'on utilisait avant ne soit pas disponible. En tous cas, les résultats en sortie ne sont pas les mêmes dépendemment du MNT corrigé en entrée qu'on lui donne (sortie de planchon/darboux ou wang et liu), celui de planchon et darboux parait visuellement plus adapté. -> TODO, à creuser.Rusle
V2 vers V3
Test exécutation sur le serveur
Environ 15s.
Détails implémentation
Pour le découpage/ré-échantillonnage de la couche OCS :
mode
) –> remarque : la sortie ré-échantillonnée à moins d'effet poivre et sel que la sortie par défaut;-tr 25.0 25.0
)Emprise de référence :
Jusqu'à présent, l'emprise des couches de sortie était écrite en dur dans les traitements du modeleur. Elle est maintenant accessible depuis "Extraire l'emprise d'une couche".
Configuration des noms des champs à utiliser dans les calculatrices de champs comme des entrées du modèles
Champs vectoriel
;Les utiliser dans des formules :
Valeur précalculée
(1 sur la figure ci-dessous)'
pour définir une chaîne de caractère et des opérateurs de concaténation||
. Enfin, attention à bien entourer les noms des champs (contenu dans la variable@myfield
dans l'exemple ci-dessus) de double guillemet"
.Vous pouvez constater que vous avez un modèle où les champs d'une couche sont entrés de votre traitement.
Propriétés du sol pour les nuls
Rusle
Facteur K
Rappel :
\(K= 2,8 \times 10{-7} M^{1,14} (12 – MO) + 4,3 \times 10^{-3} (s – 2) + 3,3 \times 10^{-3} (p – 3)\)
Avec :
L'équation actuelle est :
M = (% limons + % sable) (100 - % argile).
Or, ce n'est pas le % sable mais le % sable très fin. Puisque la base de données DONESOL ne fournit pas cette information, le % sable sera mis à 0 et l'équation de M deviendra :
M = (% limons) (100 - % argile)
Comment les obtenir ?
UCS et UTS
Quelles données ?
Légende :
Énoncé de la semaine de modélisation
Questions auxquelles on veut répondre
Question générale : quel est le risque erosif sur le BV de la Save ?
Objectifs
Au total :
Liste des taches
Spatialiser le facteur R à partir des données de précipitation
Nouvelle nomenclature
nomenclature.txt
-> "Code OSO ou Code_Group RGP + 900 (voir annexe des spécifications IGN du RPG pour connaître le Code_Group) ;Specs_IGN_RPG_2-0.pdf
pour les codes RPG ;nomenclature_oso.pdf
pour les codes oso ;nomenclature_CLC.pdf
;Nouvelles source de donnnées pour les propriétés des sols
Pré-traitement
découpage, ré-échantillonnage à 25m et reprojection peuvent se faire en renseignant bien dans l'outil reprojection (gdal-warp) :
Produire des indicateurs
Atlas
(ressources dispos dans les cours de 1A et une aide en ligne)
Oral
Guide erreurs et questions courantes
Nouveau :
SAGA : à l'issue des deux traitements, le format du fichier en sortie est sdat, souvent pour le calcul de L, le format de sortie est alors par défaut un .sdat. Pour beaucoup, si on laisse se format, la projection de la couche L en sortie est mauvaise -> forcer le format de sortie à geotiff pour régler le pb.
Sortie de A : attention à la semiologie, iels se focalisent sur les valeurs extrêmes qui sont abhérentes.
Jointure de K et rasterisation : tant que la jointure n'est pas permanente, l'opération de rastérisation va planter. Attention au type de données
Info sur R et sa spatialisation : rajouter l'étape de calcul de P.
Découpage / Réchantillonage :
Calcul de C : les infos de la manière dont les deux nomenclatures cohabitent ne sont pas trouvées facilement.
Comment débugger A si celui ci présente des erreurs ?
On fait on devrait leur faire masquer tout ce qui ne sont pas des surfaces agricoles, car en réalité, quand on fait une moyenne/médianne, on cherche à avoir un valeur moyenne/médianne des parcelles agricoles, et pas inclure les surfaces boisées ou les surfarces urbaines (par ex. )
Ancien
Notion de sémiologie : classes ordinales et nominales –> trop loin. Expliquer à l'oral a minima.
reclass avec Grass : attention peut y avoir des problèmes de chemin :
Calculatrice Raster GDAL
Si l'erreur suivante apparaît c'est qu'il y a une erreur dans la syntaxe de la formule, il ne faut pas faire copier coller, et faire attention aux espaces :
*
**
et non^
;.
et pas,
Ici parce que la couche était ouverte dans QGIS et le résultat n'a pas pu être écrit. Il suffit de retirer la couche de QGIS.
INRAE
reclass
:Rusle
Détecter les types de champs
. Sinon les pourcentages seront détecté comme des chaîne de caractères (i.e. du texte) ;