Commit 512b2a65 authored by Benoit Urruty's avatar Benoit Urruty

LHS

parent 8913fcd4
This diff is collapsed.
......@@ -219,17 +219,17 @@
},
{
"cell_type": "code",
"execution_count": 2,
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"sf = gpd.read_file(\"/home/urrutyb/Documents/PhD_TiPACCS/topo_antarctica/InSAR_GL_Antarctica_v02/InSAR_GL_Antarctica_v02.shp\")\n",
"sf = gpd.read_file(\"/home/urrutyb/Documents/PhD_TiPACCS/data_antarctica/InSAR_GL_Antarctica_v02/InSAR_GL_Antarctica_v02.shp\")\n",
"sf=sf.to_crs({'init':'epsg:3031'})"
]
},
{
"cell_type": "code",
"execution_count": 12,
"execution_count": 4,
"metadata": {
"scrolled": false
},
......@@ -762,7 +762,7 @@
},
{
"cell_type": "code",
"execution_count": 7,
"execution_count": 5,
"metadata": {
"scrolled": true
},
This diff is collapsed.
This diff is collapsed.
......@@ -32,7 +32,7 @@ to approach as best as possible the observed surface velocity.
Viscosity exponent :
$$
Viscosity\ exponent = \frac{1}{n}\ avec\ n=
Viscosity\ exponent = \frac{1}{n}
$$
......@@ -79,13 +79,18 @@ $$
The value of the velocity $u$ and $v$ are exported to computed in the cost function.
### The cost function
$$
J_{tot}=J_v+\lambda_{1}J_{div}+\lambda_{2}J_{\alpha}+\lambda_{3}J_{eta}
$$
The cost function is the comparison of the compute surface velocity
and the observed one. We compute here the integer of all the difference
$$
J = \sum_1^{Nobs}{\frac{1}{2}(U_{SSA} - U_{Obs})^2}
J_v = \sum_1^{Nobs}{\frac{1}{2}(U_{SSA} - U_{Obs})^2}
$$
......@@ -93,7 +98,15 @@ $$
and the derivative of the cost function.
$$
velocityb = \sum_1^{Nobs}{(U_{SSA} - U_{Obs})}
velocityb = \sum_1^{Nobs}{(U_{SSA} - U_{Obs})}
$$
### The  penalization of the divergence of the flux
$$
J_{div}=\frac{1}{2}\int (\dot{a}-div(\overline uH))^2 \\
(\frac{\partial H}{\partial t}-\frac{\partial H}{\partial t}^{obs})^2=(\dot{a}-div(\overline uH)-\frac{\partial H}{\partial t}^{obs})^2
$$
......@@ -115,19 +128,12 @@ The output from this solver are the nodal derivative of the friction parameter a
### The regularization terms
The regularization can be done one the first derivative. We gave a weight $\lambda$ to this the value computed. This weight have to be choice according to the order of the optimized parameter.
The regularization can be done one the first derivative of the nodal variable $V$. We gave a weight $\lambda_{reg}$ to this the value computed. This weight have to be choice according to the order of the optimized parameter.
$$
J_{reg} = \int_{\Omega} 0.5 (|dV/dx|)^2 d\Omega
\\with\ V \ the\ nodal\ variable
$$
#### regularization of Beta
#### regularization of Eta
### The optimization M1QN3
......
This diff is collapsed.
c
# -*- mode: org -*-
#+STARTUP: overview indent inlineimages logdrawer
#+TITLE: Journal
......@@ -2816,3 +2814,129 @@ données immportante pour l'initialisation
Entered on [2020-06-04 jeu. 15:44]
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/journal.org::*reunion][reunion]]
*** 2020-06-12 vendredi
**** Inylchek
Personnes présentes:
- thomas boileau
- johan berthet
- bruno fromento
- eric larose
- Luc Moreau
-
comprendre la difleunce
bcp de publi depuis 2018 par une coopération postdam et vienne
trouver les drain qui peuvent conduire ce plux d'eau vers l'aval
geopolitique :
deux postes frontières passé. Pas de soucis pour acceder a la zone
zone tres frequente en ete par les touristes
Cooporération en chine? bien vu ou pas?
Notion de risque pour les populations? manque d'eau en chine?
Documenter l'evenements vidage du lac glaciaire car n'a pas encore ete
fait (Eric)
hydrologie sous glaciaire comment ca se propage (Eric)
Si coté speleo jamais realiser c'est bien d'aller le faire (Luc)
problématique sur l'exploration des moulins :
chute de pierre dans le moulins
fonte du glacier, petit debit
deux ou trois ans d'instrumentations...
Eric: possible de capter poser nodes mais il est d'accord qu'il faut faire attention au donnée car il faudra
du monde pour traiter les données. Et donc y faire attention.
Johan berthet: definir la stratégie. Glof pertinent. utiliser des
entrees méthodologique non explorés. Jusqu'a present les glof on
etait etudier par les glacio et amener des hydrogeologues donne une
entrée original
Luc: l'idéee principale c'est le glof sont suivi et les risques
liés. dans quel mesure la vidange ce passe et a quelle niveau. vitesse
de l'evenement
lit de la rivière très large
Glof peut arriver entre mai et octobre
Budget:
2 appel a projet :
- fondation rolex. details et temps qu'on peut consacrer à cela. Une
idée du projet et des actions qu'on peut mener autour du "glof"
(vidange du lac). Modélisation?
180 000€ dont une part dédié à de l'instrumentation
- National geographic qui ont décalé leur appel a projet a cause du
covid. volet cinéma... faire image avec scientifique (on dirait du
déjà vu ahah)
30 000€
- MAIF aussi. Pas d'idée
Necessaire d'alimenter en documentation ces deux financeurs
Je prend en charge la synthèse des idées
Date:
une partie été dernière quinzaine d'aout pour voir l'activité hydro et
une partie en novembre pour rentrer dans les moulins
Entered on [2020-06-12 ven. 19:04]
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/journal.org::*presentation][presentation]]
*** 2020-06-18 jeudi
**** Réunion 18/06
***** interpolation python flowline
erreur lié à la method cubic
***** erreur ssa/Uobs
***** espace occigen et python
***** 1200 ITER 1200 simu
regarder simulation de manière detaillé
puis faire 100 simu
tirage aléatoire entre gamme de paramètre sur un facteur 10
latine hypercube sampling
aspect assimilation de donnée :
terme d'ébauche : quel est la probabilité que je donne a l'etat
initial
matrice de covariance d'ebauche B
coef de friction
eta : hilmar regularisation
prevoir un transient de 10 ans pour voir la difference sur la
reponse. Sensibilité
critère d'arret : quand fonction cout arrive au niveau de l'erreur j'arrete
Entered on [2020-06-18 jeu. 09:22]
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/liste.org::*linear%20solver%20tester%20GCR][linear solver tester GCR]]
*** 2020-06-22 lundi
**** biblio
***** An extended L-curve method for choosing a regularization parameter in electrical resistance tomography
utilisation de method basé sur la L-curve
plot de la l curve ||x²||=f(rms)
- method sur la courbure L -curve
- method sur le derivé seconde de la L-curve
***** Model terrain correction using variational adjoint method with Tikhonov-total variation regularization
J=J_1+ \mu_1 J_2 + \mu_2 J_3
***** Aster 2019 Thikonov regularization
Entered on [2020-06-22 lun. 16:09]
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/journal.org]]
* Tasks
** TODO CR Kirgiz
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/liste.org::*Tasks][Tasks]]
** DONE CR pyrénées
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/liste.org::*Tasks][Tasks]]
** TODO create bed with noise
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/journal.org::*r%C3%A9union%20PhD][réunion PhD]]
......@@ -8,11 +14,8 @@
location voiture]]
** TODO finir workflow
** TODO finir ecrit sur inversions
** TODO biblio regularization tikhonov
** TODO régularization pénalisation dérivé première ou apriori
** TODO regarder travail de johannes sur le butressing
** TODO linear solver tester GCR
** DONE faire test interpolation de mon code python
metho cubic resultat etrange
lineaire mieux
** DONE biblio regularization tikhonov
** DONE régularization pénalisation dérivé première ou apriori
** DONE regarder travail de johannes sur le butressing
** DONE linear solver tester GCR
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment