# # Import shape helping to spatialize/analyse further data on map # ******* Bulletin Name ******* =B= ArcView Import ***** Bulletin Version ****** =N= 1100 Arcview File Type =A= SHAPEFILE GRID File Type =A= Ascii format GRID Data File Name =A= K:\CREALP\Projets\51320-LT+GR\BONNARD\01-Deplacements\analyse 2012\data\slope2m.asc GRID Isatis Directory Name =A= grid-2m GRID Isatis File Name =A= 2006-07 GRID Isatis Variable Name =A= slope 2 SHAPEFILE Shapefile Name =A= D:\Projets\LT+GR\GIS\morphometry\Couches\interpretation_eb_2012.shp SHAPEFILE Imported Fields =A= ID SHAPEFILE Imported Fields =A= TYPE SHAPEFILE Imported Fields =A= Ice SHAPEFILE Import As =A= Polygons SHAPEFILE Polygons Use Field for Name =A= Y SHAPEFILE Polygons Name Field =A= TYPE SHAPEFILE Polygons Different Levels =A= N SHAPEFILE Polygons Group by Name =A= N SHAPEFILE Polygons Add Thickness =A= N SHAPEFILE Points One Sample Per Shape =A= Y SHAPEFILE Lines Use as Line Name =A= N SHAPEFILE Lines Line Name Field =A= TYPE SHAPEFILE Lines Header File Name =A= cassure SHAPEFILE Unit =A= Meter SHAPEFILE Isatis Directory Name =A= shape SHAPEFILE Isatis File Name =A= zone SHAPEFILE Use Projection File =A= N SHAPEFILE Projection File Name =A= None SHAPEFILE Keep Geographical Coordinates=A= N # # Import data to be kriged # ******* Bulletin Name ******* =B= Excel Import ***** Bulletin Version ****** =N= 410 Excel File Name =A= K:\CREALP\Projets\51320-LT+GR\BONNARD\01-Deplacements\analyse 2012\data\2008-09.xls Spreadsheet To Load =N= 0 Row To Skip =N= 0 Use First Line As Field Names =A= Y Isatis Directory Name =A= data-new Isatis File Name =A= data_test Field Name =A= Id_points Field Excel Column =N= 1 Field Type =A= alphanumeric Field Unit =A= Field Undef String =A= N/A Field Alpha Length =N= 8 Field Output Format - Type =A= Decimal Field Output Format - Length =A= 19 Field Output Format - Digits =A= 3 Field Name =A= x Field Excel Column =N= 2 Field Type =A= easting(X) Field Unit =A= m Field Undef String =A= N/A Field Alpha Length =N= 0 Field Output Format - Type =A= Decimal Field Output Format - Length =A= 19 Field Output Format - Digits =A= 3 Field Name =A= y Field Excel Column =N= 3 Field Type =A= northing(Y) Field Unit =A= m Field Undef String =A= N/A Field Alpha Length =N= 0 Field Output Format - Type =A= Decimal Field Output Format - Length =A= 19 Field Output Format - Digits =A= 3 Field Name =A= z Field Excel Column =N= 4 Field Type =A= numeric 32 Field Unit =A= m Field Undef String =A= N/A Field Alpha Length =N= 0 Field Output Format - Type =A= Decimal Field Output Format - Length =A= 19 Field Output Format - Digits =A= 3 Field Name =A= D-plan 08-09 Field Excel Column =N= 5 Field Type =A= numeric 32 Field Unit =A= Field Undef String =A= N/A Field Alpha Length =N= 0 Field Output Format - Type =A= Decimal Field Output Format - Length =A= 19 Field Output Format - Digits =A= 3 Field Name =A= D-Z 08-09 Field Excel Column =N= 6 Field Type =A= numeric 32 Field Unit =A= Field Undef String =A= N/A Field Alpha Length =N= 0 Field Output Format - Type =A= Decimal Field Output Format - Length =A= 19 Field Output Format - Digits =A= 3 Field Name =A= CODE Field Excel Column =N= 7 Field Type =A= alphanumeric Field Unit =A= Field Undef String =A= N/A Field Alpha Length =N= 1 Field Output Format - Type =A= Decimal Field Output Format - Length =A= 19 Field Output Format - Digits =A= 3 # # Prepare you background image # # In --> Display --> New page # Load the desired features (points shape as Basemap; isohypse as Isolines; polygone shape as Polygons) # Once ready, the image could be saved in a background directory (Management --> Save as) # During Exploratory analysis, the bacjkgrounded image could be use to better see where you are! # # # Think to save your results as image at each step (Print; then choose Print to file and select the desired format) # The command Save as, allows you only to save images in the Isatis system # # # During the Exploratory analysis, you could "mask" some data (here we have masked all added fix points from the variographic analysis, as well as some points that are far away from the bissectrice line in the H-scatter plot (drawn from Application --> Graphic specific parameter) # # # Explore the data and make adequate selection. On this basis make a variogram map in order to look at the data isotropy and then produce an experimental variogram {directionnal or omnidirectional (by Application --> calculation parameters); that should be saved Application --> Save in parameter file} # ******* Bulletin Name ******* =B= Exploratory Data Analysis ***** Bulletin Version ****** =N= 1100 Directory name =A= data-new File name =A= data_test Selection name =A= without odd Data Variable #1 =A= None Data Variable #2 =A= good for vario Data Variable #3 =A= None Weight Variable =A= None Use the grid organization =A= N Tabulated Representation =A= N Reference Variable =A= good for vario Vertical View Disposal =A= Y Format for Statistics - Type =A= Decimal Format for Statistics - Length =A= 19 Format for Statistics - Digits =A= 3 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= N # Definition of a new Statistical Page Type of representation =A= Base Map List of selected variables =A= good for vario 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 Legend =A= N Section type =N= 2 Section rank =N= 1 Isometric rescale =A= N Viewpoint definition =A= N Downwards orientation =A= N Magnification factor =N= 1 Longitude angle =N= 30 Latitude angle =N= 45 Literal Code Variable =A= None Color Code Variable =A= None Color Scale name =A= Rainbow [READONLY] [AUTO] Grid Highlight Color =A= rgb[0,0,255] Grid Mask Color =A= rgb[255,0,0] Grid High/Mask Opacity =N= 0.75 Page Name =A= Page 14: Base Map (good for vario) 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= Base_Map__good_for_vario_.cgm Page Setup Printer Name =A= Default Page Setup Printer Keep File =A= N Page Setup Printer Keep Name =A= Base_Map__good_for_vario_ Page Setup Save in Batch =A= N Page Setup Save Directory =A= TMP_EDIT Page Setup Save File Name =A= Base Map (good for vario) Page Setup Lock Legend Size =A= N Page Setup Draw Title Block =A= N Page Setup Page Frame =A= N Page Setup Use Gray Scale =A= N Page Setup White BG Print =A= N Page Setup Draw Title =A= Y Page Setup Title On Top =A= Y Page Setup Title Frame =A= N Page Setup Title Text =M= Base Map (good for vario) 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= analyse vitesse GB 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= 0.000106609766 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! good for vario View Label Background Color =A= background View Label Text Color =A= foreground View Label Border Color =A= foreground View Label Text Size =N= 0.2 View Label Position =A= North 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 # # ******* Bulletin Name ******* =B= Variogram Fitting ***** Bulletin Version ****** =N= 1000 Parameter Set Name =A= KO-dir_test-nov-122 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= Y 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= Y Global: Representing variogram names =A= N 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= Minimum number of pairs in dir[ 2] =N= Maximum distance in dir[ 1] =N= Maximum distance in dir[ 2] =N= Graphic Bounds along X-axis =A= N X-bounds for dir[ 1] =N= X-bounds for dir[ 2] =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= Geologist Convention Code =A= A+X-Z 1st Angle =A= 90 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= background Page Setup Save File Name =A= V-Model (D-plan 08-09)(data-new/data_test/good fo Page Setup Lock Legend Size =A= N Page Setup Draw Title Block =A= N Page Setup Page Frame =A= N Page Setup Use Gray Scale =A= N Page Setup White BG Print =A= N 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= analyse vitesse GB 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= 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_-_Global_Window.cgm Page Setup Printer Name =A= Default Page Setup Printer Keep File =A= N Page Setup Printer Keep Name =A= Variogram_Model_-_Global_Window Page Setup Save in Batch =A= N Page Setup Save Directory =A= background Page Setup Save File Name =A= V-Model (data-new/data_test/good for vario) Page Setup Lock Legend Size =A= N Page Setup Draw Title Block =A= N Page Setup Page Frame =A= N Page Setup Use Gray Scale =A= N Page Setup White BG Print =A= N 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 - Global 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= analyse vitesse GB 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 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= KO-dir_test-nov-122 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= No Drift 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= Generalized Cauchy Model Third parameter =N= 1 Switch Anisotropy =A= Y Anisotropic Scale Factor along U =N= 500 m Anisotropic Scale Factor along V =N= 160 m Anisotropic Scale Factor along W =N= 0 m ***** Rotation Version ***** =N= 800 # Local Rotation Convention Name =A= Geologist Convention Code =A= A+X-Z 1st Angle =A= 70.001 degrees 2nd Angle =A= 0 degrees 3rd Angle =A= 0 degrees Coregionalization =N= 0.09 # # Krigeage proprement dit # ******* Bulletin Name ******* =B= Create Grid File ***** Bulletin Version ****** =N= 410 Grid Directory Name =A= grid-2m Grid File Name =A= 2008-09 Grid Type =A= 2D Grid File Samples Directory Name =A= None Samples File Name =A= None Samples Selection Name =A= None Graphical Definition =A= N Graphic Mode =A= Punctual Projection =A= XOY Rescale in the Graphic Page =A= N X unit =A= m Y unit =A= m X format - Type =A= Decimal X format - Length =A= 19 X format - Digits =A= 3 Y format - Type =A= Decimal Y format - Length =A= 19 Y format - Digits =A= 3 X origin =N= 616600 m Y origin =N= 108800 m X Mesh =N= 1 m Y Mesh =N= 1 m X Nodes Number =N= 1000 Y Nodes Number =N= 800 ***** Rotation Version ***** =N= 800 # Grid 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 # ******* Bulletin Name ******* =B= Standard (Co-)Kriging ***** Bulletin Version ****** =N= 1100 Kriging Option =A= Block Number of variables =N= 1 Maximum number of External Drifts =N= 0 Full set of output variables =A= N Rank of the extracted drift =N= 1 Data Directory Name =A= data-new Data File Name =A= data_test Data Selection Name =A= without odd Variable #1 =A= D-plan 08-09 Variance of Measurement error =A= None Kriging Weights =A= None Target Directory Name =A= grid-2m Target File Name =A= 2008-09 Target Selection Name =A= AOI Estimation for Variable #1 =A= KO-dir_test-11.12 St dev for Variable #1 =A= KO-dir_test-11.12 SDEV 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= D-plan 08-09 Collocated CoKriging option =A= N Collocated Variable in Data =A= D-plan 08-09 # Definition of the Model ***** Model Version ***** =N= 903 Modification status =A= N Name of the Model Set =A= KO-dir_test-nov-122 # Definition of the Neighborhood ***** Neighborhood Version ***** =N= 1000 Modification status =A= N Name of the Neighborhood Set =A= neig_1 # 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 # # pour enlever les valeurs négatives # ******* Bulletin Name ******* =B= Calculator ***** Bulletin Version ****** =N= 1000 Directory Name =A= data-displacement File Name =A= pts-2008-09 Selection Name =A= without odds Float Variable Number =N= 15 Selection Variable Number =N= 2 Macro Variable Number =N= 2 Variable Name =A= v1 =D-plan 08-09 Transformation =M= v1=ifelse(v1<0,0,v1) .../... =M=