GIS设计与实现7.地理模型库设计.ppt

地理建模实例土壤肥料分析模型,1、建模准备,土壤肥力是衡量土壤生产力的重要指标。土壤肥力的好坏主要取决于土壤中N、P、K、有机质含量的高低。不同的土壤单元，各种营养物质的含量不同，用单一的要素来表征土壤的综合肥力状况，往往会产生与实际不符的结果，如下例，土壤单元B中磷（P）、钾（K）、有机质的含量高于土壤单元A，而氮（N）的含量则较低。N(ppm)P(ppm)K(ppm)有机质%土壤单元A 72.91 9.74 174.98 1.4962土壤单元B 70.56 10.24 190.68 1.546 需要建立一个综合的评价模型（指标）来确定到底哪个土壤肥力高些。,返回,2、模型假设,W=（w1 w2 w3 w4）为判断矩阵最大特征值所对应的特征向量，W的各个分量就是各评价因子对土壤肥力的影响权重。,假设一：构造层次分析法的判断矩阵，矩阵各元素取值须有经验专家参与。,假设二 确定各评价因子对土壤肥力的隶属度函数 上面确定的权重只体现了各因子作为一种实体对作物生长及土壤肥力的贡献大小，没有考虑各因子的含量对土壤肥力及作物生长的影响。N、P、K、有机质等因子在土壤中的含量与土壤肥力有很大的相关性，一般来说，含量越高，土壤肥力（作物产量）越大，含量越低，土壤肥力则小。但这种比例关系是有一定的限度的，超出这种限度，则这种比例关系不再成立。根据这种关系及相关专家经验可用指派法来确定各因子的隶属度函数，其计算模型见下页。,E1、E2、E3、E4为N、P、K、有机质四个养分因子对土壤肥力的隶属度函数，e1、e2、e3、e4为各土壤采样点N、P、K、有机质的实际含量。,3、建立模型土壤肥力指数的计算公式如下：S(i)=WE(i)T,S(i)第i个样点的土壤肥力指数；W评价因子对土壤肥力的影响权值矩阵，W=（w1 w2 w3 w4）；E（i）T评价因子隶属度矩阵的转置矩阵，E(i)=（Ei1、Ei2、Ei3、Ei4）；,按该方法求出的土壤单元A的肥力指数是0.936，土壤单元B的肥力指数是0.925，说明A的肥力相对于B来说要稍高一些。,4、模型求解 根据上述模型，求出了所有采样点的土壤肥力指数。,5、模型分析 采用空间数据的内插技术，用GIS 软件绘制出实验区土壤肥力空间分布图，该图按土壤肥力指数S（i）把实验区分成五个不同的等级，反映了土壤肥力的空间差异状况。在进行田间管理时，就可以针对不同的肥力区采取不同的管理措施。,6、模型检验 一般来说，土壤肥力越高，作物的产量越大，下图是实验区的产量分布图，对比产量图和土壤综合肥力分布图，发现：产量高的地区与肥力高的地区相一致，这说明，我们建立的模型是合理的。,返回,