Table des matières
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
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
