克里金插值法
8页1、克里金插值法克里金插值法又称空间局部插值法,是以变异函数理论和结构分析为基础,在有限区域 内对区域化变量进行无偏最优估计的一种方法,是地统计学的主要内容之一,由南非矿产工 程师D. Mat heron于1951年在寻找金矿时首次提出,法国著名统计学家G. Mat heron随后 将该方法理论化、系统化,并命名为Kriging,即克里金插值法。1克里金插值法原理克里金插值法的适用范围为区域化变量存在空间相关性,即如果变异函数和结构分析的 结果表明区域化变量存在空间相关性,则可以利用克里金插值法进行内插或外推。其实质是 利用区域化变量的原始数据和变异函数的结构特点,对未知样点进行线性无偏、最优估计 无偏是指偏差的数学期望为0,最优是指估计值与实际值之差的平方和最小1。因此,克里 金插值法是根据未知样点有限领域内的若干已知样本点数据,在考虑了样本点的形状、大小 和空间方位,与未知样点的相互空间关系,以及变异函数提供的结构信息之后,对未知样点 进行的一种线性无偏最优估计。假设研究区域a上研究变量Z (x),在点x G A (i=1, 2,n)处属性值为Z(x),II则待插点x G A处的属性值
2、Z(x )的克里金插值结果Z*(x )是已知采样点属性值Z(x) 000i(i=1, 2,,n)的加权和,即:Z*(x ) = 入 Z(x )(1)0 i ii=1式中 i 是待定权重系数。其中Z(x)之间存在一定的相关关系,这种相关性除与距离有关外,还与其相对方向变 i化有关,克里金插值方法将研究的对象称“区域化变量”针对克里金方法无偏、最小方差条件可得到无偏条件可得待定权系数i (匸1,2,,n)满足关系式:才九=1(2)ii=1九以无偏为前提,kriging方差为最小可得到求解待定权系数i的方程组:工九 C(x ,x ) + 卩=C(x ,x )(j = 1,2, n)3)i i j0 j2 X= 1i i=1式中,C(x, x)是Z(x)和Z(x )的协方差函数。i j i j2 方法步骤克里金插值法的应用步骤如下:1、输入原始数据,即采样点,下面以输入三个采样点求待估插值为例来进行说明。如3、数据检验与分析,根据采样值是否合乎实际情况,剔除明显差异点。4、直方图的计算,直方图有助于掌握区域变化的分布规律,以便决定是否对原始数据 进行转换。5、利用变异函数进行变异函数计算,了解
3、变量的空间结构。6、克里金插值估计(1) 待估点权重系数估计 利用多边形估计的方法,首先确定离待估点最近的采样点的权重,根据公式(4)进行 采样点权重估计:1c + 一dw九=i(4)c+dwi=12)根据搜索策略选择合适的参估点,如图2图 2 参估点图示(3)根据已经求出的变异函数以及采样点数量,三个采样点列出三个等式,求出方程组的系数,公式为:_ C (1,1)C (1,2)C (1,3)1-C (0,1)C (2,1) C (2,2) C (2,3)九2=C (0,2)C (3,1) C (3,2) C (3,3)九3C (0,3) _5)4)分析在各向同性条件下改变块金值与在块金值相同条件下改变各向异性对权重值的影响。各向同性条件下改变块金值时对权重值的影响效果如图3 (a),在块金值相同条件下改变各向异性对权重值带来的影响如图3 (b):a)b)图 3 各向同性条件下改变块金值与在块金值相同条件下改变各向异性对权重值的影响(5)根据求出的权重值,代入公式(1),即可求得评估领域内n个采样值的线性组合2。克里金插值法的方法路线图如下:图4 方法路线图3克里金插值法分类及适用类型
4、克里金插值法主要有以下几种类型:普通克里金(Ordinary Kriging)、简单克里金(Simple Kriging)、泛克里金(Universal Kriging)、协同克里金(Co-Kriging)、对数正 态克里金(Logistic Normal Kriging)、指示克里金(Indicator Kriging)、概率克里金(Probability Kriging)和析取克里金(Disjunctive Kriging)等1 克里金插值法可以简单地表达为:Z (s) = P (s) + 8 (s)(6)式中,s为不同位置的点,可以人为是用经纬度表示的空间坐标;Z(s)为s处的变量 值,它可以分解为确定趋势值卩(s)和自相关随机误差(s)。通过对这个公式进行变化,可 以生成克里金插值法的不同类型。首先,对于趋势值卩(s),可以简单地赋予一个常量,即在任何位置s处y(s)二卩,如果卩是未知的,这便是普通克里金基本模型;卩(s)也可表示为空间坐标的线性函数,如:p (s) = B + P x + P y + P x 2 + P y 2 + p xy0 1 2 3 4 5(7)如果趋势
《克里金插值法》由会员公****分享,可在线阅读,更多相关《克里金插值法》请在金锄头文库上搜索。
小熊请客教案.doc
浇剪力墙结构工程模板施工方案-secret【建筑施工资料】.doc
西方国家财政联邦制的比较与借鉴.doc
民用运输机场应急救援规则
苏州大学中国古代史参考文献目录
初中学业水平考试语文模拟试题一.doc
熔铸温度对Al3Nb-铝基复合材料硬度及-摩擦磨损性能影响的研究毕业论文.doc
断指再植的术前准备和术后处理.doc
培训师的工具箱完整版
10分钟懒人鸡翅的做法【健康食谱】
六年级班级德育工作计划范本(2篇).doc
2023年浙江省丽水市遂昌县云峰街道连头村社区工作人员考试模拟试题及答案
室内自动喷水灭火系统安装工艺.doc
产地证企业变更网上操作简介.doc
2023毕业答辩的范文.docx
XXXX年月金融理财师真题解析一
妈妈的爱600字初一作文
试用期工作总结及自我鉴定(3篇).doc
宝宝营养辅食制作大全
世界无烟日科普知识.doc
2023-06-15 19页
2023-05-23 4页
2022-09-03 10页
2023-02-21 3页
2023-02-13 18页
2022-10-29 9页
2023-11-04 1页
2023-01-08 13页
2023-10-10 2页
2024-01-11 22页