克里格方法内插生成高程曲面课件.ppt

实习8： 克里格方法内 插生成高程曲面 ;1. 背景： 现有某地区一系列高程采样点，需要通过内插生成该地区的高程层面，为后续研究提供合理的数据层面信息。 2. 目的：通过练习熟练掌握并理解每种克里格方法的原理及实现过程，体会他们在具体应用中的适应性。 3. 要求：根据数据特征，至少选用一种克里格方法内插生成高程表面，并分析对于此例，该方法的的适用性。（有能力的同学可选用多种克里格方法进行插值，并比较对于此例不同方法的优劣） 4. 数据：某地区的高程采样点（jyg.shp）；数据存放于Chp10\ex1。 ;操作步骤：（1） 在 ArcMap 中加载 jyg.shp。 （2） 启动地理统计模块 Geostatistical Analyst。 （3） 单击 Geostatistical Analyst 模块的下拉箭头点击 Create Subsets 命令。 （4） 在弹出的对话框中，Input 选项中选择需要生成子集的数??jyg，点击 Next 按钮。（5） 在弹出的对话框中，通过滑块设置训练子集与检验子集的比例，在 Output Personal Geodatabase 中设置子集的输出路径及名称，点击完成按钮（图1）。 ;图1;（6）单击选中 jyg_training 层面，随后点击 Geostatistical Analyst模块的下拉箭头选择 Explore Data 并点击 Histogram 命令，或点击Geostatistical Analyst 模块的下拉箭头选择 Explore Data 并点击Normal QQPlot 命令，生成如下结果（图 2）。由图上可知，数据分布符合正态分布的假设，不需要进行数据变换。 ;图 2;（7） 点击 Geostatistical Analyst 模块的下拉箭头选择 Explore Data 并点击 Trend Analysis 命令，查看数据是否存在趋势，如 图3 所示。;（8）单击 Geostatistical Analyst 模块的下拉箭头点击 Geostatistical Wizard 命令。 （9）在弹出的对话框中，在 Dataset1 选择训练数据 jyg_training 及其属性 STATION，在 Validation 中选择检验数据 jyg_test 及其属性 STATION，在 Methods 中选择Kriging 内插方法，最后点击 Next 按钮。 （见图 4） ;图 4;（10）在弹出的对话框中，展开泛克里格Universal Kriging，在下面的选项中点击预测图（Prediction Map），在 DataSet1 选项卡中的 Transformation 里选择 None变换方式，在 Order of Trend 里选择 First，点击 Next 按钮（见图 5）;图 5;（11） 在弹出的 Detrending 对话框中，点击 Next 按钮。（如图 6） ;（12） 在弹出的Semivariogram/Covariance Modeling 对话框中（图 7），先按照默认参数进行操作，在得到对模型精度评定的结果后，发现结果误差太大，返回更改该对话框中的参数，经比较发现，将分组数 Number of 设为 10 得到的结果较好。需注意的是，在设置分组数时，尽量保证每组中的样点对数大于 10。最后点击 Next 按钮。;图 7;（13） 在弹出的 Searching Neighborhood 对话框中，点击 Next 按钮。 （图 8）;（14） 在弹出的 Cross Validation 对话框中，显示了对模型的精度的评价，如图 9 所示。在对不同参数得到模型的比较中，可参考 Prediction Error 中的几个指标。符合以下标准的模型是最优的：标准平均值（Mean Standardized）最接近于 0，均方根预测误差（Root-Mean-Square）最小，平均标准误（Average Mean Error）最接近于均方根预测误差（Root-Mean-Square），标准均方根预测误差（Root-Mean-Square Standardized）最接近于 1。最后点击 Next 按钮。 ;图 9;（15） 在弹出的 Validation 对话框中，点击 Finish 按钮。泛克立格法内插结果如图 10所示。 ;（16） 双击 jyg 层面，在弹出的属性对话框中，选择 Symbology 选项卡，展开 Quantities列表，选中 Graduated Symbols，在 Value 中选择 STATION，将符号大小的变化范围 Symbol Size f