基于ArcGIS_10.0的地表高程序贯高斯模拟和地质统计建模
1. 数据需求：
总共需要四个文件夹：即两个输入型文件夹和两个输出型文件夹，其中输入型文件夹包括一个存放原始数据的文件夹和一个存放shp格式数据的文件夹；而输出型文件夹包含一个存放简单克里金估计结果（dbf格式）的文件夹和一个存放序贯高斯模拟结果（dbf格式）的文件夹。
以及三个gdb数据库：即一个存放简单克里金输出结果（Table格式）的数据库，一个存放100次序贯高斯模拟的直接云图结果的数据库以及一个存放序贯高斯模拟输出结果（Table格式）的数据库。（之所以需要给出简单克里金插值结果，是因为序贯高斯模拟只能在简单克里金插值结果的基础上进行计算和分析）
2. 数据关系
1）Import_Surface_Data&Background_Grid_xlsx&txt文件夹：
① Import_Surface_Data.xlsx（原始数据，包含位置坐标和属性值，供建立shp用）
② Import_Surface_Data.txt（原始数据，包含位置坐标和属性值，供Surfer插值生成背景 格用以及后期和ArcGIS插值结果作对比）
③ Background_Grid.xlsx（只包含待插值点的位置坐标，一般来说插值高程的话为x，y；插值属性的话可以为x，y，z；）
④ Import_Surface_Data.grd
⑤ Output_grid_of_kriging_standard_deviation_from_surfer.grd
⑥ Out_import_Surface_Data.dat
⑦ Out_Transform.dat
⑧ Out_Convert.dat
2）Import_Map_shp文件夹：
① Import_Surface_Data.dbf
② Background_Grid.dbf
③ Kriging_Estimation_Result_Second.dbf
④ Kriging_Estimation_Result_Third.dbf
3）Export_Kriging_Result_shp文件夹：
① Export_Output_Second.dbf
② Export_Output_Third.dbf
4）Extract_Simulation_Results_shp文件夹：
1～100次的模拟结果（dbf格式），在此不详细列出
5）Sequential_Gaussian_Simulation.gdb数据库：
1～100次的模拟结果，在此不详细列出
6）Export_Kriging_Result.gdb数据库：
① Export_Output_Second
② Export_Output_Third
7）Extract_Simulation_Results.gdb数据库：
1～100次的模拟结果（Table格式），在此不详细列出
3. 实施方案
A. 首先需要建立基于龙头山表面高程原始数据的shp文件和利用Surfer软件生成的待插值背景 格节点的shp文件，其中需要通过Surfer软件先生成10m*10m的插值结果grd文件，在把其转换成dat格式文件,方便得到背景 格的shp文件；
B. 然后通过简单克里金对原始数据进行插值，在得到克里金插值结果的同时结合背景 格文件生成一个基于验证和预测的克里金插值结果shp文件，从中可以看到每个 格节点上的预测结果值与预测标准差；
C. 最后利用上述克里金插值结果来获得100次序贯高斯随机模拟的结果，其中包含每个 格节点上的最大值、最小值、均值和标准差等等，以及每一次模拟所生成的云图，并通过提取功能得到模拟的输出结果以便于后期外部处理和利用。
4. 具体技术路线
启动ArcCatalog生成shp文件，如下图所示：
启动ArcMap导入之前生成的shp文件，如下图所示：
修改标记，方便查看和论文成图，如下图所示：
如果第一次安装ArcGIS，地统计学工具条并未可用，则需要点击菜单栏中的Customize 下的Extensions按钮，勾选其中的Geostatistical Analyst前置方框，然后Close，地统计学工具条就可以用了，如下图所示：
点击已经变绿的地质统计学分析工具条，对Import_Surface_Data进行简单克里金统计分析（这里尽量使用简单克里金估计类型，为了方便后期模拟使用），如下图所示：
注意选择正态变换和去趋势，因为一般数据通常不符合插值所需的标准条件：
注意漂移阶数，一定要选择合适的漂移阶数使得去趋势和转换后的数据尽可能的呈现出正态分布的特征，如下图所示（例如本次模拟选择2阶多项式漂移）：
协方差函数和变异函数模型的拟合曲线如下图所示：
选择合适的模型，并计算其基台值，偏基台值，块金值以及变程等参数：
General
Examine bivariate distribution False
Variable Semivariogram
Optimize model
Model Nugget
Enable True
Calculate Nugget True
Nugget 0.0071928488230167
Measurement Error 100
Model #1
Type Gaussian
Major Range 255.65093131342041
Anisotropy False
Calculate Partial Sill True
Partial Sill 1.2536437163695078
Model #2
Type
Model #3
Type
Lag
Lag Size 29.37115822112131
Number of Lags 12
得到最终结果，如下图所示：
预测值和测量值有很高的相关性，如下图所示：
Kriging_Estimation_Result_Second.dbf
Kriging_Estimation_Result_Third.dbf
然后通过对简单克里金估计结果进行验证和预测，导出克里金估计结果shp文件，其中需要用到背景 格文件，如下图所示：
注意下面这个对话框中的第二项Input pointobservation locations应该选择背景 格（背景 格是之前建立用来展示克里金估计结果）
最后开始准备执行序贯高斯模拟，如下图所示：
在ArcCatalog里面新建立一个Sequential_Gaussian_Simulation.gdb数据库文件，如下图所示：
谨慎进行序贯高斯模拟参数的选择，序贯高斯模拟需要之前做的简单克里金估计结果作为模板来进行条件随机模拟（注意这里需要选择条件数据，就是原始测量数据），如下图所示：
再在ArcCatalog里面重新建立两个gdb数据库来存放克里金估计和序贯高斯模拟所导出的Table格式文件，即：
Export_Kriging_Result.gdb；
Extract_Simulation_Results.gdb
再在ArcCatalog里面重新建立两个文件夹来存放克里金估计和序贯高斯模拟所导出的shp格式的dbf文件，方便后期的excel处理，即：
Export_Kriging_Result_shp
Extract_Simulation_Results_shp
导出克里金估计结果到Export_Kriging_Result.gdb数据库，并得到Table格式数据，如下图所示：
再从之前建立的Table中导出克里金估计结果到Export_Kriging_Result_shp文件夹，并得到dbf格式数据，方便excel直接打开处理，如下图所示：
注意这里文件格式应选为dBASE Table：
再从之前的Extract_Simulation_Results.gdb数据库中的各个Table表中提取数据到Extract_Simulation_Results_shp文件夹，并得到dbf格式数据，方便后期的excel直接打开处理，如下图所示：
注意这里文件格式应选为dBASE Table：
声明：本站部分文章及图片源自用户投稿，如本站任何资料有侵权请您尽早请联系jinwei@zod.com.cn进行处理,非常感谢！