空间统计分析初步.ppt

甘肃农业大学资源与环境学院变异函数的性质设Z(x)是区域化变量，在满足二阶平稳假设条件下，变异函数式具有如下性质：=0，即在h=0处，变异函数为0；=，即关于直线h=0是对称的，它是一个偶函数；≥0，即只能大于或等于0；甘肃农业大学资源与环境学院|h|→∞时，→c(0)或=c(0),即当空间距离增大时,变异函数接近先验方差[-]必须是一个条件非负定函数，由[-]构成的变异函数矩阵在条件时，为非负定的。甘肃农业大学资源与环境学院变异函数的计算公式设是系统某属性Z在空间位置x处的值，为一区域化随机变量，并满足二阶平稳假设，h为两样本点空间分隔距离，和分别是区域化变量在空间位置和处的实测值[i=1,2,…,N(h)]，那么，变异函数的离散计算公式为(4.2.10)0102甘肃农业大学资源与环境学院这样对不同的空间分隔距离h，计算出相应的和值。如果分别以h为横坐标，或为纵坐标，画出协方差函数和变异函数曲线图，就可以直接展示区域化变量Z(x)的空间变异特点。可见，变异函数能同时描述区域化变量的随机性和结构性，从而在数学上对区域化变量进行严格分析，是空间变异规律分析和空间结构分析的有效工具。甘肃农业大学资源与环境学院甘肃农业大学资源与环境学院例如：假设某地区降水量Z(x)（单位：mm）是二维区域化随机变量，满足二阶平稳假设，其观测值的空间正方形网格数据如图所示（点与点之间的距离为h=1km）。试计算其南北方向及西北和东南方向的变异函数。甘肃农业大学资源与环境学院从图可以看出，空间上有些点，由于某种原因没有采集到。如果没有缺失值，可直接对正方形网格数据结构计算变异函数；在有缺失值的情况下，也可以计算变异函数。只要“跳过”缺失点位置即可（图）。图4.2.1空间正方形网格数据（点间距h=1km）甘肃农业大学资源与环境学院图1空间正方形网格数据（点间距h=100m）甘肃农业大学资源与环境学院甘肃农业大学资源与环境学院图4.2.2缺失值情况下样本数对的组成和计算过程

☉为缺失值首先计算南北方向上的变异函数值，由变异函数的计算公式可得=385/72=5.35甘肃农业大学资源与环境学院方向南北方向西北—东南h12345h1.412.824.245.657.07N(h)362721135N(h)322113825.359.2617.5525.6922.907.0612.9530.8558.1350.00同样计算出最后，得到南北方向和西北—东南方向上的变异函数计算结果见下表。同样可以计算东西方向上的变异函数。甘肃农业大学资源与环境学院变异函数有4个非常重要的参数，即基台值（sill）、变程（range）或称空间依赖范围（rangeofspatialdependence）、块金值（nugget）或称区域不连续性值（localizeddiscontinuity）和分维数（fractaldimension）。前3个参数可以直接从变异函数图中得到。它们决定变异函数的形状与结构。变异函数的形状反映自然现象空间分布结构或空间相关的类型，同时还能给出这种空间相关的范围。变异函数的参数#2022甘肃农业大学资源与环境学院当变异函数随着间隔距离h的增大，从非零值达到一个相对稳定的常数时，该常数称为基台值C0+C。A当间隔距离h=0时，γ(0)=C0，该值称为块金值或块金方差（nuggetvariance）。B基台值是系统或系统属性中最大的变异，变异函数达到基台值时的间隔距离a称为变程。变程表示在h≥a以后，区域化变量Z(x)空间相关性消失。C块金值表示区域化变量在小于抽样尺度时非连续变异，由区域化变量的属性或测量误差决定。D甘肃农业大学资源与环境学院上述3个参数可从变异函数曲线图直接得到，或通过估计曲线回归参数得到。第4个参数，即分维数用于表示变异函数的特性，由变异函数和间隔距离h之间的关系确定分维数D为双对数直线回归方程中的斜率，它是一个无量纲数。分维数D的大小，表示变异函数曲线的曲率，可以作为随机变异的量度。甘肃农业大学资源与环境学院理论变异函数模型实践中，常用的是变异函数图：偏基台值:C(partialsill)