Modélisation 3D


Concept général

Principe de base

Dans le cadre du processus de diagnostic du danger issu d'une zone de permafrost alpin il est nécessaire de connaître le volume de matériau mobilisable ainsi que le flux déterminant la vitesse de recharge. Ce flux est approché par la vitesse de surface et par les inclinomètres. Le volume est lui déterminé à partir d'une modélisation du fond stable. Cette modélisation des coordonnées 3D du fond stable ne peut pas s'appuyer sur une mesure précise. Il s'agit de combiner au mieux l'ensemble des informations, cartographie de terrain, géophysique, forage, mesure des déplacements, disponibles et d'en tirer les données de bases nécessaire à l'interpolation pondérée de la surface stable. La préparation de ces données d'origine diverse demande une grande rigueur. Une fois celle-ci déterminée, l'utilisation des logiciels SIG permet de déterminer les volumes, ou de faire les coupes nécessaires au calcul des volumes.

Figure 6.1 Exemple d'intégration des données en 3D

Implémentation

  • But : Gérer le processus de validation, correction, mise en forme, des données nécessaires au calcul de la topographie du bedrock.
  • Quand : en deux temps, avant la campagne de forage et à la fin de l'étude.
  • Résultat : un flow-chart.
  • Prérequis : cartographie de terrain, premiers profiles de géophysique, toutes données utiles.
  • Utile pour : planifier les forages, calculer les volumes.

Exploitation des données

  • Définition des points de mesure de profondeur du fond stable/bedrock sur la base des forages, de la géophysiques, des affleurements (définissant leur niveau de précision). Ces valeurs sont confrontées entre-elles et mise en cohérence dans un référentiel unique.
  • Les données sont ensuite validées et stockées dans une base de données flexible selon le flowchart ci-dessous. Le but étant de pouvoir faire des aller et retour en affinant les données à chaque itération.

Figure 6.2 Flux des données et de leur stockage dans l'optique de fournir un modèle tridimensionnel.

Forces et faiblesses

Forces Faiblesses
Gestion cohérente des données. Nécessite d'avoir une grande quantité de donnée.
Peu nécessiter l'emploi de logiciels onéreux.

Alternatives et perspectives

Sans une grande quantité de données, la procédure indiquée ici n'a pas de grand sens. D'autre méthode comme le SLBL (Jaboyedoff et al., 2004, Derron et al. 2005) peuvent être alors utilisée.




Préparation des données

Principe de base

Lors de toute modélisation il est important de bien contraindre les limites du modèle. Dans le cas d'un modélisation de la topographie du bedrock, il est nécessaire que le résultat final du modèle corresponde aux affleurements et falaise entourant la zone d'intérêt (AOI). De même, dans la zone d'intérêt central les données issues de différentes méthodes peuvent être entachées d'erreurs méthodologiques (cf. cf. Figure 4.3 configuration défavorable du sous-sol).

Figure 6.3 Confrontations des données de différentes origines et mise en cohérence 3D

Implémentation

  • But : Voir les données dans un contexte tridimensionnel met en évidence les zones cohérentes/incohérentes (par ex. croisement entre des coupes). Cette inspection visuelle permet de fournir des données de la meilleure qualité possible sur l'ensemble de l'AOI.
  • Quand : avant la campagne de forage.
  • Résultat : une première topographie du bedrock indiquant les lacunes de donnée.
  • Prérequis : cartographie de terrain, 1ère données géophysique.
  • Utile pour : définir les données manquantes et le modèle 3D final

Spécifications techniques

  • Pour réaliser cette préparation un logiciel permettant la mise en cohérence tridimensionnelle des données issues de source diverse est nécessaire. Après le test de plusieurs logiciels de 3D (Carlson, GeoMensura, GSI3D) nous avons opté pour le logiciel AdHoc qui était le seul à permettre une préparation des données avant la procédure de modélisation proprement dite.
  • Connaître les limites des différentes méthodes utilisées

Exploitation des données

  • Les données sont intégrées au modèle 3D défini par le MNT de la surface.
  • Les corrections et adaptations, basées sur des hypothèses raisonnables sont apportées.
  • Dans le cas où on suspecte des artefacts, mais qu'on n'arrive pas à les corriger sur la base d'hypothèses raisonnables, on augmente l'incertitude sur la valeur mesurée.

Figure 6.4 La divergence entre les profils de sismique réfraction peut être due à la configuration du sous-sol

Figure 6.5 Le prolongement, raisonable, de la topographie des affleurements permet d'ajouter des données utiles pour maitriser les effets de bord

Figure 6.6 Mise en cohérence des informations tirées de éa géophysique et des forages

Forces et faiblesses

Forces Faiblesses
Permet d'avoir une bonne densité de points dans tous les secteurs de l'AOI.
Permet de gérer les différents points de façon simples à l'extérieur du logiciel d'interpolation
Une partie des informations sont basées sur des hypothèses.
Système expert.
Nécessite une bonne maîtrise des logiciels CAD/SIG.

Interpolation

Principe de base

Les multiples données à notre disposition peuvent toutes être exprimée comme un triplets de coordonnées, une erreur de mesure inhérente à la méthode et un code de qualité. Afin de pouvoir spatialiser ces données ponctuelles nous avons utilisé des techniques d'interpolation géostatistiques. Ces méthodes utilisent la relation spatiale existant entre les données afin de reproduire le phénomène grâce à un modèle. Ce modèle d'interpolation, qui respecte les valeurs mesurées, est ensuite appliqué aux points d'une grille à maille régulière (Matheron, 1963; Burrough, 2001).

Cette technique d'interpolation respecte les valeurs aux points de données. Basiquement elle est similaire à celle utilisée pour spatialiser les données de déplacements. Ce calcul est in fine, le but à atteindre dans l'étude du danger.

Implémentation

  • But : définir l'altitude du fond stable.
  • Quand : à la fin de l'étude.
  • Résultat : topographie du fond stable.
  • Prérequis : préparation des données.
  • Utile pour : le calcul des volumes disponibles

Spécifications techniques

  • Pour modéliser les 28 ha du “Glacier” Bonnard, 1120 points de mesures de la profondeur du fond rocheux stable ont été nécessaires. Ces points doivent être réparti de façon la plus homogène possible sur la surface à étudier (ce qui n'est pas toujours possibles d'un point de vue logistique). Ces 1120 points sont de précision et de qualité différentes (sismique réfraction et réflexion, forage, affleurement et interprétation cartographique).

Figure 6.7 Données utilisées pour l'interpolation du fond stable


  • Un MNT de haute précision est nécessaire afin d'obtenir des courbes de niveau aux limites de l'AOI. Ces données sont nécessaires pour contraindre au mieux la surface finale.

Exploitation des données

Les éléments de géostatistiques sont identiques à ceux présentés dans le chapitre sur la spatialisation des déplacements. Une fois déterminée, la surface interpolée du bedrock a été confrontée aux données des forages. Elle peut aussi être analysée comme une surface topographique traditionnelle dans un SIG pour en dériver des indices comme la pente, la courbure, etc.

Figure 6.8 Mise en évidence des talwegs sur la surface du bedrock avec en vert les géophones utilisés pour la sismique réfraction

Forces et faiblesses

Forces Faiblesses
Permet une interpolation robuste.
Défini le degré de certitude de façon différenciée sur l'ensemble de l'AOI
traitement itératif nécessitant un opératuer expérimenté.
Nécessite une grande puissance de calcul.

Difficultés et écueils

Lors de la modélisation de la topographie dans une zone pentue (ce qui est le cas lors de l'analyse d'un mouvement de masse) un biais systématique est introduit. De façon évidente les points situés à l'aval se trouvent à une altitude plus faible! Ce biais impacte la structure spatiale des données qui sera introduite dans le modèle de variogramme.

Il est utile d'enlever cette tendance avant la modélisation du variogramme. Après plusieurs tests une surface obtenues à partir d'un MNT avec une maille de 100 m (BK trend topo 100) a été utilisé et soustrait à la valeur d'altitude des points retenus pour l'analyse. Le modèle géostatistique est ensuite développé à partir de la structure spatiale des résidus obtenus de cette soustraction.

Alternatives et perspectives

Si la densité de points est très importante, par exemple dans le cas de glacier rocheux s'avançant sur des terrains connus au préalable, des algorithmes d'interpolation plus simples et plus rapides peuvent être utilisé.

Une fois la partie inférieure du “sandwich” connu, il est intéressant d'utiliser les informations obtenues sur la profondeur (géophysique et forage) pour modéliser le remplissage du “sandwich”. Des méthodes de géostatistique en 3D peuvent alors être utilisées, comme les plurigaussiennes.

Procédure et commentaires

La procédure étant relativement longue et fastidieuse, on pourra utiliser le fichier journal ci-dessous. Pour cela, en changeant le nom des fichiers sources, utilisez la fonction batch –> Edit journal file –> Journal file –> reload.

toit_bedrock_ISATIS.ijnl
%VAR selection = sel data
#
# Lissage du MNT :
 ******* Bulletin Name *******          =B= Grid Operator
 ***** Bulletin Version ******          =N= 1000
 Directory Name                         =A= Base data
 File Name                              =A= mnt2m-tot
 1 bit Variable: v1                     =A= None
 1 bit Variable: v2                     =A= None
 1 bit Variable: v3                     =A= None
 1 bit Variable: v4                     =A= None
 1 bit Variable: v5                     =A= None
 Real Variable: w1                      =A= alt-topo
 Real Variable: w2                      =A= alt_topo liss100m
 Real Variable: w3                      =A= None
 Real Variable: w4                      =A= None
 Real Variable: w5                      =A= None
 Structural Element Definition          =A= s1  = B[25 , 25 , 1]
 Morphological Transformation           =M= w2=moving_average(s1,w1)
#
# Calcul des residus et initialisation du modele de variogramme :
 ******* Bulletin Name *******          =B= Global Trend Modeling
 ***** Bulletin Version ******          =N= 320
 Maximum number of External Drifts      =N= 1
 Data Directory Name                    =A= Base data
 Data File Name                         =A= data_BK-2012
 Data Selection Name                    =A= None
 Variable                               =A= toit BK
 Weight                                 =A= None
 External Drift #1                      =A= alt_topo liss100m
 Trend at Data                          =A= toit BK trend topo liss100m
 Residuals at Data                      =A= toit BK residus trend topo liss100m
 Grid Directory Name                    =A= Base data
 Grid File Name                         =A= mnt2m-tot
 Grid Selection Name                    =A= None
 Trend                                  =A= toit BK trend topo liss100m
 External Drift on Target #1            =A= alt_topo liss100m
 Trend Option                           =A= 1 f1
 Model Pset Name                        =A= toit BK NonStat topo_liss100m
 Save in Parameter File                 =A= Y
 Use Grid File                          =A= Y
#
# Calcul du variogramme experimental :
 ******* Bulletin Name *******          =B= Logical Operations on Selections
 ***** Bulletin Version ******          =N= 320
 Directory Name                         =A= Base data
 File Name                              =A= data_BK-2012
 First Selection Name                   =A= sel data
 Selection Operator                     =A= .AND.
 Second Selection Name                  =A= Sans code 5 et 6
 NEW Selection Name                     =A= sel data sans code 5 et 6
#
 ******* Bulletin Name *******          =B= Exploratory Data Analysis
 ***** Bulletin Version ******          =N= 1100
 Directory name                         =A= Base data
 File name                              =A= data_BK-2012
 Selection name                         =A= sel data sans code 5 et 6
 Data Variable #1                       =A= toit BK residus trend topo liss100m
 Weight Variable                        =A= None
 Use the grid organization              =A= N
 Tabulated Representation               =A= N
 Reference Variable                     =A= toit BK residus trend topo liss100m
 Vertical View Disposal                 =A= Y
 Format for Statistics - Type           =A= Decimal
 Format for Statistics - Length         =A= 10
 Format for Statistics - Digits         =A= 2
 Use the Weight Variable                =A= N
 Default Symbol size                    =N= 0.2
 Default Character size                 =N= 0.2
 Default Line Style                     =A= Solid
 Default Line Color                     =A= foreground
 Define parameters before graphic       =A= Y
# Definition of a new Statistical Page
 Type of representation                 =A= Variogram
 List of selected variables             =A= toit BK residus trend topo liss100m
 Start Graphic Definition Parameters    =A= Y
# Graphic Parameters
 Draw axes and title                    =A= Y
 Draw active points only                =A= N
 Pick points                            =A= Y
# Specific Calculation Parameters
 Draw the experimental variogram        =A= Y
 Draw the variogram cloud               =A= Y
 Draw the experimental variance         =A= Y
 Use the same distance axis             =A= Y
 Draw the number of pairs               =A= N
 Draw the histogram of pairs            =A= N
 Same scale for all directions          =A= Y
 Draw all the direction                 =A= Y
 Calculate the variogram cloud          =A= N
 Rank of the direction                  =N= 0
 Type of calculation                    =N= 0
 Type of variography                    =A= Variogram
 Number of directions                   =N= 1
 Number of regular directions           =N= 0
 Orthogonal direction                   =A= N
 ***** Rotation Version *****           =N= 800
# Reference Plane
 Convention Name                        =A= Geologist
 Convention Code                        =A= A+X-Z
 1st Angle                              =A= 40 degrees
 2nd Angle                              =A= 0 degrees
 3rd Angle                              =A= 0 degrees
# Parameters for Direction #1
 Rank of the elementary definition      =N= 0
 Label for this Direction               =A= 
 Number of lags                         =N= 20
 Special code usage                     =A= N
 Grid skipping value                    =N= 1
 Grid index along X                     =N= 0
 Grid index along Y                     =N= 0
 Grid index along Z                     =N= 0
 ***** Rotation Version *****           =N= 800
# Direction Definition
 Convention Name                        =A= Geologist
 Convention Code                        =A= A+X-Z
 1st Angle                              =A= 90 degrees
 2nd Angle                              =A= 0 degrees
 3rd Angle                              =A= 0 degrees
 Tolerance on distance                  =N= 0.5
 Tolerance on direction                 =N= 90
 Lag value                              =N= 25 m
 Last Lag to be Refined                 =N= 
 Lag Subdivision                        =N= 
 Slicing Width                          =N= 
 Slicing Height                         =N= 
 Page Name                              =A= Page 2: Variogram (toit BK residus trend topo li)
 P-Set name                             =A= None
 Save Action                            =A= Override
 Page Setup Begin                       =A= Y
 ***** Page Setup Version *****         =N= 1000
 Page Setup Print in Batch              =A= N
 Page Setup Format Name                 =A= Computer Graphic Metafile (CGM)
 Page Setup Print to File               =A= N
 Page Setup Print File                  =A= Variogram__toit_BK_residus_trend_topo_liss100m_.cgm
 Page Setup Printer Name                =A= Default
 Page Setup Printer Keep File           =A= N
 Page Setup Printer Keep Name           =A= Variogram__toit_BK_residus_trend_topo_liss100m_
 Page Setup Save in Batch               =A= N
 Page Setup Save Directory              =A= Graphiques
 Page Setup Save File Name              =A= Variogram (toit BK residus trend topo liss100m)
 Page Setup Lock Legend Size            =A= Y
 Page Setup Draw Title Block            =A= Y
 Page Setup Page Frame                  =A= Y
 Page Setup Use Gray Scale              =A= N
 Page Setup White BG Print              =A= Y
 Page Setup Draw Title                  =A= Y
 Page Setup Title On Top                =A= Y
 Page Setup Title Frame                 =A= N
 Page Setup Title Text                  =M= Variogram (toit BK residus trend topo liss100m)
 Page Setup Title Size                  =N= 0.4
 Page Setup Draw Description            =A= Y
 Page Setup Manual Description          =A= N
 Page Setup Description Frame           =A= Y
 Page Setup Description Size            =N= 0.2
 Page Setup Logo                        =A= Y
 Page Setup Logo Bitmap                 =A= N
 Page Setup Logo Bitmap File Name       =A= None
 Page Setup Logo Date                   =A= Y
 Page Setup Logo User                   =A= Y
 Page Setup Logo User Name              =A= ERIBAR
 Page Setup Logo Study                  =A= Y
 Page Setup Logo Study Name             =A= Bonnard_2012
 Page Setup Logo Size                   =N= 0.2
 Page Setup Graphic Size Mode           =N= 1
 Page Setup Page Width                  =N= 19
 Page Setup Page Height                 =N= 27.7
 Page Setup Scale Value                 =N= 
 Page Setup Magnification Factor        =N= 1
 Page Setup End                         =A= Y
 Number of Views                        =N= 1
 View Label Begin                       =A= Y
 ***** View Label Version *****         =N= 1000
 View Label Visible                     =A= Y
 View Label Manual Text                 =A= N
 View Label Text                        =M! 
 View Label Background Color            =A= transparent
 View Label Text Color                  =A= foreground
 View Label Border Color                =A= foreground
 View Label Text Size                   =N= 0.2
 View Label Position                    =A= North West
 View Label Gravity                     =A! Center
 View Label Alignment                   =A= Center
 View Label Inside                      =A= Y
 View Label X Center                    =N! 0
 View Label Y Center                    =N! 0
 View Label Left-Right Offset           =N! 0.2
 View Label Bottom-Top Offset           =N! 0.2
 View Label End                         =A= Y
 Grid Axes Begin                        =A= Y
 ***** Grid Axes Version *****          =N= 1100
 Number of Grid Axes                    =N= 1
 Grid Axis Number 1                     =A= Y
 Grid Axis Name                         =A= All View
 Draw Grid and Subgrid                  =A= N
 Number of Horizontal Subdivisions      =N= 0
 Number of Vertical Subdivisions        =N= 0
 Grid Axes End                          =A= Y
 End of Definition                      =A= Y
 Default Minimum                        =N= 0
 Default Maximum                        =N= 0
 Default Number of classes              =N= 1
 Default Use Log Scale                  =A= N
 Default Log Shift                      =N= 0
#
# Ajustement du modele :
 ******* Bulletin Name *******          =B= Variogram Fitting
 ***** Bulletin Version ******          =N= 1000
 Parameter Set Name                     =A= toit BK residus trend topo liss100m
 Fitting option                         =N= 2
 Number of iterations                   =N= 100
 Verbose Automatic Sill Fitting option  =A= N
 Override checks for KU usage           =A= Y
 Automatic Sill Fitting                 =A= N
 Advanced Parameters                    =A= N
 Masking minimum number of pairs        =A= N
 Masking maximum distance               =A= N
# Parameters for the Fitting Window
 Fitting: Representing the variogram    =A= Y
 Fitting: Option for lag representation =A= N
 Fitting: Join experimental lags        =A= Y
 Fitting: Representing number of pairs  =A= Y
 Fitting: Representing variogram names  =A= N
 Fitting: Representing pairs histogram  =A= N
 Fitting: Histograms with same scale    =A= Y
 Fitting: Representing the model        =A= Y
 Fitting: Representing the envelop      =A= N
 Fitting: Representing the variance     =A= Y
 Fitting: Representing title and axes   =A= Y
# Parameters for the Global Window
 Global: Representing the variogram     =A= Y
 Global: Option for lag representation  =A= N
 Global: Join experimental lags         =A= Y
 Global: Representing number of pairs   =A= N
 Global: Representing variogram names   =A= Y
 Global: Representing pairs histogram   =A= N
 Global: Histograms with same scale     =A= Y
 Global: Representing the model         =A= Y
 Global: Representing the envelop       =A= Y
 Global: Representing the variance      =A= Y
 Global: Representing title and axes    =A= Y
 Global Representation mode             =N= 0
 Number of rows in Global Window        =N= 1
 Number of cols in Global Window        =N= 1
 View rank of the Window Global         =N= 1
# Matrix of covariances at the origin
 C00[ 1, 1]                             =N= 0
 Model Normation mode                   =N= 0
 Constrained sill for variable  1       =N= 
 Constraint on nugget effect            =A= N
 Nugget[ 1, 1]                          =N= 0
# Masking Rules
 Minimum number of pairs in dir[ 1]     =N= 
 Maximum distance in dir[ 1]            =N= 
 Graphic Bounds along X-axis            =A= N
 X-bounds for dir[ 1]                   =N= 
 Graphic bounds along Y-axis            =A= N
 Ymax-Bounds for vars[ 1, 1]            =N= 
 Ymin-Bounds for vars[ 1, 1]            =N= 
 ***** Rotation Version *****           =N= 800
# Simulation Rotation
 Convention Name                        =A= Mathematician
 Convention Code                        =A= +Z+Y+X
 1st Angle                              =A= 0 degrees
 2nd Angle                              =A= 0 degrees
 3rd Angle                              =A= 0 degrees
 Simulation Orientation                 =N= 0
 Simulation mesh along X                =N= 0.01 m
 Simulation mesh along Y                =N= 0.01 m
 Simulation nodes along X               =N= 100
 Simulation nodes along Y               =N= 100
 Simulation Number of bands             =N= 2000
 ***** Rotation Version *****           =N= 800
# Variogram Map Rotation
 Convention Name                        =A= Mathematician
 Convention Code                        =A= +Z+Y+X
 1st Angle                              =A= 0 degrees
 2nd Angle                              =A= 0 degrees
 3rd Angle                              =A= 0 degrees
 Vmap Orientation                       =N= 0
 Vmap Number of sectors                 =N= 90
 Vmap Number of lags                    =N= 100
 Vmap Lag                               =N= 1 m
 Draw the Fitting window                =A= Y
 Page Setup Begin                       =A= Y
 ***** Page Setup Version *****         =N= 1000
 Page Setup Print in Batch              =A= N
 Page Setup Format Name                 =A= Computer Graphic Metafile (CGM)
 Page Setup Print to File               =A= N
 Page Setup Print File                  =A= Variogram_Model_-_Fitting_Window.cgm
 Page Setup Printer Name                =A= Default
 Page Setup Printer Keep File           =A= N
 Page Setup Printer Keep Name           =A= Variogram_Model_-_Fitting_Window
 Page Setup Save in Batch               =A= N
 Page Setup Save Directory              =A= Graphiques
 Page Setup Save File Name              =A= V-Model (toit BK residus trend topo liss100m)(Bas
 Page Setup Lock Legend Size            =A= Y
 Page Setup Draw Title Block            =A= Y
 Page Setup Page Frame                  =A= Y
 Page Setup Use Gray Scale              =A= N
 Page Setup White BG Print              =A= Y
 Page Setup Draw Title                  =A= Y
 Page Setup Title On Top                =A= Y
 Page Setup Title Frame                 =A= N
 Page Setup Title Text                  =M= Variogram Model - Fitting Window
 Page Setup Title Size                  =N= 0.4
 Page Setup Draw Description            =A= Y
 Page Setup Manual Description          =A= N
 Page Setup Description Frame           =A= Y
 Page Setup Description Size            =N= 0.2
 Page Setup Logo                        =A= Y
 Page Setup Logo Bitmap                 =A= N
 Page Setup Logo Bitmap File Name       =A= None
 Page Setup Logo Date                   =A= Y
 Page Setup Logo User                   =A= Y
 Page Setup Logo User Name              =A= ERIBAR
 Page Setup Logo Study                  =A= Y
 Page Setup Logo Study Name             =A= Bonnard_2012
 Page Setup Logo Size                   =N= 0.2
 Page Setup Graphic Size Mode           =N= 1
 Page Setup Page Width                  =N= 19
 Page Setup Page Height                 =N= 27.7
 Page Setup Scale Value                 =N= 
 Page Setup Magnification Factor        =N= 1
 Page Setup End                         =A= Y
 Draw the Display window                =A= N
 Draw the Simulation Window             =A= N
 Draw the Variogram Map Window          =A= N
# Automatic Variogram Fitting
 Use Automatic Fitting                  =A= Y
 Initial Combination                    =A= Spherical
 Add Nugget                             =A= N
 Auto Fit in Batch                      =A= N
# Definition of the Model
 ***** Model Version *****              =N= 903
 Modification status                    =A= Y
 Name of the Model Set                  =A= toit BK NonStat topo_liss100m
 Switch Rotation                        =A= Y
 Global Rotation                        =A= N
 ***** Rotation Version *****           =N= 800
# Global Rotation
 Convention Name                        =A= Geologist
 Convention Code                        =A= A+X-Z
 1st Angle                              =A= 90 degrees
 2nd Angle                              =A= 0 degrees
 3rd Angle                              =A= 0 degrees
 Drift type                             =A= 1 f1
 Increment status                       =A= N
 Order for increments                   =N= 0
 Increment Direction                    =A= X Direction
 Lag for increments                     =N= 0 m
 Convolution defined                    =A= N
 Truncation defined                     =A= N
 Covariance type                        =A= Spherical Model
 Third parameter                        =N= 1
 Switch Anisotropy                      =A= N
 Scale Factor                           =N= 240 m
 ***** Rotation Version *****           =N= 800
# Local Rotation
 Convention Name                        =A= Geologist
 Convention Code                        =A= A+X-Z
 1st Angle                              =A= 90 degrees
 2nd Angle                              =A= 0 degrees
 3rd Angle                              =A= 0 degrees
 Coregionalization                      =N= 190
#
# Krigeage avec derive externe :
 ******* Bulletin Name *******          =B= Standard (Co-)Kriging
 ***** Bulletin Version ******          =N= 1100
 Kriging Option                         =A= Punctual
 Number of variables                    =N= 1
 Maximum number of External Drifts      =N= 1
 Full set of output variables           =A= N
 Rank of the extracted drift            =N= 1
 Data Directory Name                    =A= Base data
 Data File Name                         =A= data_BK-2012
 Data Selection Name                    =A= sel data sans code 5 et 6
 Variable #1                            =A= toit BK
 External Drift on Data #1              =A= alt_topo liss100m
 Variance of Measurement error          =A= None
 Kriging Weights                        =A= None
 Target Directory Name                  =A= Base data
 Target File Name                       =A= mnt2m-tot
 Target Selection Name                  =A= AOI
 Estimation for Variable #1             =A= toit BK KDE topo liss100m
 St dev for Variable #1                 =A= toit BK KDE topo liss100m EcT
 External Drift on Target #1            =A= alt_topo liss100m
 Collocated Variable                    =A= None
 Local mean for Variable #1             =A= None
 Radius of the convolution Ball         =N= 0 m
 Gradient or Laplacian Unit             =A= Meter
 Measurement error option               =A= N
 Name of the Erroneous Variable         =A= toit BK
 Collocated CoKriging option            =A= N
 Collocated Variable in Data            =A= toit BK
# Definition of the Model
 ***** Model Version *****              =N= 903
 Modification status                    =A= N
 Name of the Model Set                  =A= toit BK NonStat topo_liss100m
# Definition of the Neighborhood
 ***** Neighborhood Version *****       =N= 1000
 Modification status                    =A= N
 Name of the Neighborhood Set           =A= unique
# Definition of the Model Special Options
 ***** Model Special Version *****      =N= 400
 Simple Kriging                         =A= N
 Filtering Option                       =A= N
 Factorial Kriging Option               =A= N
 Deconvolution Option                   =A= N
 Stop at first inversion problem        =A= Y
 ***** Variable Model Version *****     =N= 1100
 LGS Directory Name                     =A= None
 LGS File Name                          =A= None
 LGS Base Name                          =A= 
 LGS Global Use Rotation                =A= N
 LGS Global 2D                          =A= N
 LGS Global Rotation / Z                =A= None
 LGS Global Rotation / Y                =A= None
 LGS Global Rotation / X                =A= None
 LGS Global Use Sill                    =A= N
 LGS Global Sill                        =A= None
 LGS Model Number of Variables          =N= 1
 LGS Number of Structures               =N= 1
 LGS Structure #1 Use Rotation          =A= N
 LGS Structure #1 2D                    =A= N
 LGS Structure #1 Rotation / Z          =A= None
 LGS Structure #1 Rotation / Y          =A= None
 LGS Structure #1 Rotation / X          =A= None
 LGS Structure #1 Use Range             =A= N
 LGS Structure #1 Proportional          =A= N
 LGS Structure #1 Range Factor          =A= None
 LGS Structure #1 Range / X             =A= None
 LGS Structure #1 Range / Y             =A= None
 LGS Structure #1 Range / Z             =A= None
 LGS Structure #1 Use Sill              =A= N
 LGS Structure #1 Sill Factor           =A= None
 LGS Structure #1 Use Third             =A= N
 LGS Structure #1 Third Parameter       =A= None
 LGS Neighborhood Use Rotation          =A= N
 LGS Neighborhood 2D                    =A= Y
 LGS Neighborhood Rotation / Z          =A= None
 LGS Neighborhood Rotation / Y          =A= None
 LGS Neighborhood Rotation / X          =A= None
 LGS Neighborhood Use Radius            =A= N
 LGS Neighborhood Proportional          =A= N
 LGS Neighborhood Radius / Z            =A= None
 LGS Neighborhood Radius / Y            =A= None
 LGS Neighborhood Radius / X            =A= None
 LGS Neighborhood Factor                =A= None
#
 ******* Bulletin Name *******          =B= Calculator
 ***** Bulletin Version ******          =N= 1000
 Directory Name                         =A= Base data
 File Name                              =A= mnt2m-tot
 Selection Name                         =A= None
 Float Variable Number                  =N= 5
 Selection Variable Number              =N= 2
 Macro Variable Number                  =N= 2
 Variable Name                          =A= v1  =alt-topo
 Variable Name                          =A= v2  =toit BK KDE topo liss100m
 Transformation                         =M= v2=ifelse(v1>v2,v2,v1)
#
 ******* Bulletin Name *******          =B= Migrate Grid to Point
 ***** Bulletin Version ******          =N= 320
 Number of variables                    =N= 1
 Input Directory Name                   =A= Base data
 Input File Name                        =A= mnt2m-tot
 Input Selection Name                   =A= None
 Input Variable #1                      =A= toit BK KDE topo liss100m
 Output Directory Name                  =A= Base data
 Output File Name                       =A= data_BK-2012
 Output Selection Name                  =A= None
 Migrated Variable #1                   =A= toit BK KDE topo liss100m
 Maximum Migration Distance             =N= 
#
 ******* Bulletin Name *******          =B= ArcView Export
 ***** Bulletin Version ******          =N= 320
 Isatis Directory Name                  =A= Base data
 Isatis File Name                       =A= mnt2m-tot
 Isatis Variable Name                   =A= toit BK KDE topo liss100m
 File Type                              =A= Ascii format
 Data File Name                         =A= K:\CREALP\LT+GR\BONNARD\18-ModelisationGeol_3D\2012\SIG\ToitBK_kde100liss.asc




model-3d.txt · Dernière modification: 2014/09/09 17:20 (modification externe)
Recent changes RSS feed Creative Commons License Donate Minima Template by Wikidesign Driven by DokuWiki