000000b_谢彩香-中药资源.docVIP

基于地理信息系统的中药材生态适宜性定量化研究 谢彩香1,索风梅1,周应群2,王丽芝1,向丽1,董梁4,贾光林5,孙成忠3,陈士林1 （1.中国医学科学院 药用植物研究所，北京 100193; 2.中国药材集团公司，北京 100094 3.中国测绘科学研究院,北京 100830 4. 河南中医学院，郑州 450003 5.青岛农业大学,青岛 266109） [摘要] 目的：中药材品质与其生长环境密切相关。为应用数学模型对中药材产地适宜性进行定量化研究，需要对影响中药材生长的生态因子和气候数据的插值进行研究。方法：本文首先从生态学和统计分析角度研究对影响中药材生长的生态指标对生产有重要意义的因子，生长发育环境因素是指直接或间接地影响的所有因素生态因素是指环境中直接影响形态、生理和分布的因素，如光、温度、水分、等[3-11][15][18]。光、热、水是构成气候条件的主要因子，它们决定了成土过程的水热条件，是土壤环境系统中的重要因素，而且与中药材生长发育唇齿相关。年平均温度、最冷和最热月平均温度[12]；影响分布的有：影响分布的影响分布的土壤是赖以生存的物质基础，土壤的结构、PH值、肥力、水分等与生长密切相关。一般适宜在有机质含量高，团粒结构，保水、保肥性能好，中性或微酸性的土壤上生长。海拔高度、水域深度等因素属于，因为它们的作用是通过温度、光照、气压等生态因素间接地体现出来的，它们本身对生物没有直接影响。在分析植物的分布时，应区分生态因子间的关系 图1 药用植物生态因子间的层次关系图 因此在进行生态适宜性分析的因子选择时，应该选择具体指标层次的因子。假设从药材分布区布设采样点m个，从数据库中提取出mAi i 1,……,n ，计算每个生态因子的统计指标如下： 平均值Mean： 标准差SD Standard Deviation ： 变异系数CV Coefficient of variation ：C.V ，变异系数可消除单位或平均数不同对多个生态指标变异程度比较的影响[13]。 从上述统计指标数值可知，变异系数较大的生态因子对药材分布影响较弱，而变异系数较小的是影响药材分布的主要生态因子，因此选择变异系数较小的生态因子作为药材生态适宜性分析的区划指标，对区划指标的最适宜范围可用半峰宽公式计算。 区划指标的权重分两种情况，若有每个采样点处的药材质量数据，可采用逐步回归数学方法分析各个因子的权重，否则采用层次分析法（AHP:Analytic Hierarchy Process），AHP是一种定性和定量相结合的系统化和层次化的分析方法[14]。AHP进行权重计算的具体过程如下： 首先根据专家打分建立区划指标间的判断矩阵A，A aij n×n，判断矩阵具有如下特征：aii 1，aji 0，aji ，求解矩阵A的最大特征根λmax。为判别权重分配是否合理，需要对判断矩阵进行一致性检验，一致性检验指标为： CR是判断矩阵的一致性比率，λ为矩阵的最大特征根，n是矩阵A的阶数。RI是随机一致性指标。若判断矩阵A通过一致性检验，即可用特征向量ω作为区划指标的权向量。 3 气象要素插值 气象要素是反映地球表层环境状况的指标参数，也是参与中药区划的关键指标。气象要素如降水、温度、太阳辐射等在空间尺度上是连续分布的数据，但由于自然和社会因素的限制导致气象站点数量有限而且分布不均，需要进行气象要素的空间插值以得到非广大气象站点所在区域内的气象信息数据。空间插值是通过已知的空间数据估算未知空间数据值的数学方法，其前提条件是空间地物具有一定的空间相关性。常用于空间插值的方法有逆距离权重法（Inverse Distance Weighing）、多项式插值法（Interpolating Polynomials）、克吕格插值法（Kriging）、样条插值法（Spline）等。对于应用层面的气象要素空间插值，在选择插值方法时要考虑相对适合同时满足精度要求的方法。 由于气温空间分布与地形小环境密切相关，进行气温要素插值时应考虑气温的水平地带性和垂直地带性分布特征。有研究表明[16-17]，影响气温空间分布的因素主要有经度（主要考虑离海远近）、纬度（主要考虑太阳辐射）、地形条件（山系走向、坡度、坡向）、海拔以及下垫面（土壤和植被状况）等，其中尤以海拔高度和地形影响最显著。而数字高程模型（Digital Elevation Model, DEM）作为地球表面地形的数字描述可容易获取空间位置属性和地形特征，可为研究气温的空间分布提供强大的空间数据库。 基于DEM的多元线性回归插值方法：通过对中国六百多个气象台站三十多年以来的月平均气温与气象站所在的经纬度和海拔高度进行偏相关性分析发现气象站月平均气温与气象站的海拔高度和经纬度有极显著的相关性， （1）