免费文献传递   相关文献

Some key problems of quantitative estimation in multilevel remote sensing monitoring system of plant medicinal materials

药用植物多级遥感监测体系中定量估计的若干关键问题



全 文 :·综述·
药用植物多级遥感监测体系中定量估计的
若干关键问题
谢彩香1,陈士林1,王 蠫1,周应群2,李 滢1
(1.中国医学科学院 药用植物研究所,北京 100094;2.中国药材集团公司,北京 100094)
[摘要] 多学科技术的发展进步促进了药材资源调查新体系的建立。建立多级遥感平台、粗中细分辨率互补
的药用植物资源监测系统是资源调查技术发展和现实需求的必然结果。作者探讨了多级遥感资源调查的技术思
路,并针对其中存在的关键问题和有关数学理论进行了初步研究,提高药用植物的调查精度,对实际工作具有一定
的指导意义和理论价值。
[关键词] 药用植物;多级遥感;最优估测方程;样方调查
[中图分类号]R282.7 [文献标识码]A [文章编号]10015302(2008)03032604
[收稿日期] 20070828
[基金项目] 科研院所社会公益研究专项(2005DIB3J088);中
央级公益型科研院所基金(YZ102);工业生物技术及产品的研究
(Z0005190043321)
[通讯作者] 陈士林,Tel:(010)62899700,Email:scslchen@
sohu.com
  药用植物资源调查的目的是掌握药用植物的现势情况,
为合理开发利用药用植物、制定宏观政策提供数据支持和依
据[1]。但药用植物资源的分布不同于农业资源和林业资源,
由于药用植物资源品质受生长环境影响很大,因此具有很强
的地域性,分布零散并且地面覆盖范围不规则。传统中药资
源调查是根据收购量推算、主观估计和地面调查相结合,不
仅费时、费力、调查周期长,而且由于主观性很大,因此调查
结果精度较低。遥感技术的不断创新为改进中药资源调查
工作提供了新的可能,积极合理地利用新技术会大大促进中
药资源调查方法的发展,中药资源调查新体系(多级遥
感)[2]的建立正是在不断吸收和利用新技术基础上成长和
发展起来,见图1。利用低空遥感技术的目的是提高边界线
调绘精度,同时根据低空遥感影像可确定地面样方的调查区
域和样方数量。同时,将野外调查样方展绘在校正后的低空
遥感影像上,利用ArcInfo软件自动提取样方对应的遥感信
息及整个研究区域的遥感信息。
多级遥感的资源调查新思路为:高分辨率低空遥感数据
与中分辨率卫星数据相结合进行空中抽样来估计药用植物
面积或蓄积量。在中分辨率影像上布设高分辨率样地,以样
地位置处的高分辨率低空遥感影像为依据布设地面调查样
地,对整个药材分布区进行估计,这样可以大大减少地面工
作量,充分发挥各级遥感平台和多种分辨率卫星数据功能,
  
图1 多级遥感药用植物资源调查示意图
减少经费投入和劳动强度,同时,提高药用植物资源的调查
精度和准确度。
1 药用植物蓄积量遥感定量估计的自变量选择原则
传统药材资源蓄积量的调查重点是分布面积和单位面
积蕴藏量的测算。分布面积采用直接测量法,单位面积药材
蕴藏量的测算采用传统样方测定,对于分布零散不集中的某
些药用植物,测量工作量不仅巨大而且精度难以保证;若用
NDVI单一指标估测生物量,势必丢失大量遥感信息,而且
NDVI对高植被区具有较低的灵敏度。
采用定量遥感技术建立生物量估测模型将改善估测精
度[3,4]。以遥感和地理信息为基础,进行药用植物蓄积量定
量估测,是以地面样地的蓄积观测值为因变量,以样地所对
应的遥感图像各波段灰度值、灰度比值及样地信息为自变
量,通过数量化方法和最小二乘原理(leastsquare,LS)或岭
估计(ridgetrace,RT)建立蓄积和自变量之间的回归模型进
·623·
第33卷第3期
2008年2月
         
    中 国 中 药 杂 志
ChinaJournalofChineseMateriaMedica
       
Vol.33,Issue 3
February,2008
行蓄积量估测。设某观测区域内地面样地总数为 n,样地的
蓄积观测值向量为 V
 
n×10
,在样地对应的遥感和地理信息中,
影响蓄积估测的全部自变量维数为p(含常数项),其观测矩
阵为 X
 
n×1
,则该观测区域内蓄积线性估测模型[7]可表示为:
V0=Xγ+e,E(e)=0,COV(e)=σ
2I (1)
式中,e为 n×1随机误差向量;σ2为样地蓄积观测方
差;γ为p×1待定参数向量;I为n阶单位阵。
将X表示成分块形式 X=(XqMtXt),相应将 γ表示
成γ'=(γq'┆γt'),则(1)式可变为:
V0=Xqγq+Xtγt+e (2)
式中Xq包含常数项,且Xq和Xt分别为q,t列的列满秩
阵。考虑两种情况:
①蓄积量真实模型为V0=Xγ+e,因丢掉一些必要的自
变量,而被误认为
V0=Xqγq+e (3)
②蓄积量真实模型为 V0=Xqγq+e,因为引进了一些不
必要的自变量,而被误认为:
V0=Xγ+e (4)
设模型(3)和(4)两种模型待定参数的 LS估计及方差
估值分别为:
γq=(Xq′Xq)
-1XqV0
槇σ2q=V′0[I-Xq(Xq′Xq)-1Xq]V′0/(n-q)
γ^=(X′X)-1X′V0=
γ^q
Λ
γ^

1
α^2=V′0[I-X(X′X)
-1X′]V0/(n-p-1)
则有以下结论:
1)若模型(4)正确,则有
E(珘γq)=γq+(X′qXq)
-1X′qXtγt
COV(^γq)≥COV(珘γq)
要使γq为γq的无偏估计,必有(X′qXq)
-1X′qXtγt。表明
γt=0,即后t个自变量与蓄积无关,模型(3)正确;矩阵 X=
(XqMtXt)的两部分正交,即后t个自变量与蓄积估测不起作
用。另外,即使模型(4)正确,丢掉一部分自变量后,总使剩
下的那部分自变量回归系数的LS估计的方差减小。
2)若模型(4)正确,当COV(γ

q)≥γtγ′t时,则有
COV(γ

q)≥MSEM(γ

q)
Var(Δν0)≥MSEP[珓ν0=E(Δν

oq)]
式中,MSEM(珘γq)为待估参数的平均平方误差矩阵;
Δν2oq,Δν0分别为模型(3)和(4)的蓄积预测偏差;MSEP(珓ν0)
为模型(3)的蓄积预测均方误差。
由以上讨论可知:若丢掉的自变量影响确实存在(γt≠
0),但γt难以准确估计(用COV(γ

q)≥γ′tγt刻画),丢掉这样
的变量可使预测精度提高;相反,若模型包含了不必要的自
变量(γt=0),则使预测精度降低。
因此,药用植物资源蓄积预测自变量的选择应遵循少而
精的原则,尽量遗弃影响很小或难以准确估计的变量。
2 药用植物蓄积量定量估计最优方程的选择
在上述确定起主导作用的因子后,随后的问题是要选多
少自变量去建立方程,即怎样选择“最优自变量子集”问题。
据研究表明[5]:选择回归自变量的基本原则是少而精。丢掉
一些变量总会使得剩余变量回归系数LS估计和预测的方差
减少,此时的估计和预测一般是有偏差的。
鉴于其有效性和便于计算机编程,最小二乘(LS)估计
采用平均残差平方和(residualmeansquares,RMS)准则,按
照回归树原理,从所有可能的自变量子集中选择影响蓄积估
测的主要因子子集。当自变量个数为q时,平均残差平方和
的表达式为:
Rq=
VT[I-Xq(X

qXq)
-1XTq]V
n-q
式中,分子 VT[I-Xq(X

qXq)
-1XTq]V为蓄积残差平方
和,且Rq+1≤Rq、即最小二乘(LS)估计蓄积残差平方和随自
变量个数的增加而减少,其在回归方程中包含对蓄积估计有
影响的全部自变量时达到最小值[8]。按照上式选择自变量
时,随着自变量个数q的增加、(n-q)-1将增大,但开始时分
子减少更多,故总起来 Rq将逐渐减少。当自变量增加到一
定程度,重要的自变量基本上都已选上了,此时再增加自变
量,分子减少不多,以至于抵消不了(n-q)-1的增加,最终
还是导致了Rq的增加。到达最小值后,又逐渐增大,见图2。
依据RMS准则,其最小值对应的变量子集即为影响蓄积估
测的主要因子子集。据此可选取最优方程以反映该地区的
药用植物资源概况。
图2 Rq随q的变化曲线
在上述最小二乘方法中,存在的问题是选出的变量间可
能存在一定的复共线性,使观测阵出现病态。当影响资源蓄
积估测的RS和GIS信息间存在复共线性时,用LS估计建立
蓄积估测方程已不再适合,现在许多学者[6]开始研究通过其
他数学方法,如岭估计来提高蓄积量预报的精度和准确度。
3 地面样方的抽样调查设计
抽样调查为科学研究方法中重要技术之一,是指就所要
·723·
第33卷第3期
2008年2月
         
    中 国 中 药 杂 志
ChinaJournalofChineseMateriaMedica
       
Vol.33,Issue 3
February,2008
研究的某特定现象之母群体中,依据随机原理抽取一部分作
为样本,作为研究母群体的依据。在抽样信赖水平内,根据
样本的抽样研究结果推算母群体的整体特征。
野外样方调查是中药资源调查监测工作中的重要数据
源,也是中药资源调查的一个重要组成部分。为了满足中药
资源调查需要,凡是从遥感影像、地理信息和地图中无法获
取的数据,都必须进行野外样方调查。中药资源调查需要多
方面的样本数据,比如建立数学模型、遥感判读训练、样地动
态变化分析、资源调查监测精度分析评价等。进行野外样方
调查时应注意以下几点:一是样方设置要合理以保证获得正
确信息;二是要省时省力省钱;三是在制定抽样方案时,调查
的抽样方法要与后续的数据分析配套;四是要充分利用已知
信息来提高抽样精度。为充分发掘样地信息的潜在能力和
提高抽样效率,在调查中要特别注意抽样设计,抽样方案的
设计好坏直接影响着调查结果的准确度和精度。
最优样地要能够代表监测区域内药材资源的种类及分
布状况,包括应抽取样地的总数及每类样地应抽取的样地
数。抽样的目的是用所抽样地对应的主要遥感和地理信息
因子建立以样地为单位的蓄积量估测方程,以便利用估测方
程预报与样地类型一致的监测区域的资源蓄积量。
设监测区域已有样地的总数为n,在此n个样地中抽取
q个样地来建立资源蓄积量的估测方程。如何使抽取的样
地能代表监测区域内药材资源的种类及分布状况,根据被抽
样地对应的主要遥感和地理信息建立蓄积估测方程,才能有
效预报监测区域的蓄积量。如何抽取监测区域建立蓄积估
测方程的最优样地,属于多目标优化问题,建立资源蓄积最
优样地的数学模型如下分析。
设监测区域样地对应影响蓄积量估测的主要遥感和地
理信息因子为x=(x1x2…xp)。每个样地对应 p个变量,相
当于p维欧式空间Ep中的一点。设样地的蓄积量 y和 p各
主要影响因子间存在下列关系:
y=β1x1+x2+…+βpxp+ey (1)
式中,β1,β2…βp为待定系数;ey为样地蓄积量观测值的
随机误差。
在监测区域n个样地种抽取用于建立蓄积估测方程的
q个样地,q个样地满足下列关系:
Yq=Xqβ+eYq,E(eYq)=0 (2)
式中,Yq=
y1
y2



q
为q个样地的蓄积观测向量;
Xq=
x11 x12 … x1p
x21 x22 … x2p
… … … …
xq1 xq2 … x

qp
为q个样地p个影响因子的观
测矩阵;eYq=(ey1ey2…eyq)为 q个样地蓄积量观测误差;E
(eYq)=0为 q个样地蓄积量观测误差的均值为 0;cov
(eYq)=σ

2Iq为抽样地蓄积量观测误差等方差且彼此不相
关;Iq为q阶单位矩阵。
作为监测区域建立蓄积量估测方程的最优样地,q个被
抽样地需满足下列条件:
①用q个被抽样地组成方程②,求解药材蓄积量估测方
程;用残差平方和RSSq计算误差方差的估计值 σ

2,其值应
达到最小。②通过式(2)建立蓄积量方程后,对监测区域其
余n-q个未抽样地的蓄积量进行预报。预报偏差的方差应
达到最小。③为确保式(2)的方程组有定解,抽取的样地数
q应满足q>p。
可见,确定监测区域蓄积量估测方程的最优样地的种类
和数量属于多目标优化问题,抽取的q个最优样地应满足的
优化条件有:
min(σ

2),min(σ


n-q)
Yq=Xqβ+eYq,E(eYq)=0,cov(eYq)=σ
2Iq
q>
{

实际计算时可采用线性加权法构造评价函数,优化目标
函数表示为:
f(x)=min(λ1σ^
2+λ2σ^

n-q),x∈E
p (4)
4 小结
信息技术的发展促进了药用植物资源调查技术的进步,
也促使资源调查由定性分析阶段转向定量分析阶段。作者
针对多级遥感体系中的多种实际问题,分析了其中的有关数
学理论,给出了定量分析的数学依据用以指导生产实践。
[参考文献]
[1] 周应群,陈士林,张本刚,等.中药资源调查方法研究[J].世界
科学技术———中医药现代化,2005,7(6):130.
[2] 谢彩香,陈士林,林宗坚,等.无人机遥感技术应用于药用植物
资源调查研究[J].中国现代中药,2007,9(6):4.
[3] 赵宪文.林业遥感定量估测[M].北京:中国林业出版社,
1997.
[4] 李崇贵,赵宪文.以遥感和地理信息为基础的森林蓄积 LS估
计自变量选择研究[J].遥感学报,2001,5(4):277.
[5] 赵宪文,李崇贵.基于“3S”的森林资源定量估测[M].北京:中
国科学技术出版社,2001.
[6] 李崇贵,石 强,赵宪文,等.用岭估计研究以 RS和 GIS为基
础的森林蓄积预报[J].遥感技术与应用,2001,16(3):168.
[7] 唐守正.多元统计分析方法[M].北京:中国林业出版社,
1989.
[8] 陈希孺,王松桂.近代回归分析[M].合肥:安徽教育出版社,
1987.
[9] 孙山泽.抽样调查[M].北京:北京大学出版社,2004.
·823·
第33卷第3期
2008年2月
         
    中 国 中 药 杂 志
ChinaJournalofChineseMateriaMedica
       
Vol.33,Issue 3
February,2008
Somekeyproblemsofquantitativeestimationinmultilevelremote
sensingmonitoringsystemofplantmedicinalmaterials
XIECaixiang1,CHENShilin1,WANGYu1,ZHOUYingqun2,LIYing1
(1.InstituteofMedicinalPlantDevelopment,ChineseAcademyofMedicalSciences,
PekingUnionMedicalColege,Beijing100094,China;
2.ChinaNationalGroupCorp.ofTraditional&HerbalMedicine,Beijing100094,China)
[Abstract] Thedevelopmentandadvancementinmanysubjectsandtechnologypromotethefoundationofnewmonitoringsys
teminplantmedicinalmaterials.Thenewplantmedicinalmaterialsmonitoringsystemwithmultilevelremotesensingpanelanddifer
entresolutionisnecessityoftechnologydevelopmentandneedsofreality.Thepaperdiscussestheideaofthenewmultilevelremote
sensingsystemandstudiestherelatedkeyproblemsandmathematicaltheory,increasingthesurveyprecisionofplantmedicinalmateri
als,providingdefiniteguidanceandtheoryvalue.
[Keywords] plantmedicinalmaterials;multilevelremotesensing;optimalestimationequation;sampleinvestigation
[责任编辑 王亚君]
[收稿日期] 20070521
[通迅作者] 刘晓秋,Tel:(024)23986469,Email:Liuxiao
qiu3388@tom.tom
名贵中成药质量研究现状
林 立1,刘晓秋1,徐 晖2
(1.沈阳药科大学 中药学院,辽宁 沈阳 110016;
2.沈阳药科大学 药学院,辽宁 沈阳 110016)
[摘要] 通过对名贵中成药相关文献及国家标准的统计分析,发现名贵中成药剂型以故老剂型为主,单一落
后;质量研究以化学物质为基础,标准不完善、缺乏整体性,不能完全反映名贵中成药的质量,这些可能是因为中医
理论的“宏观思维”逻辑与现代分析技术“微观研究”的矛盾所导致的。新兴的“代谢组学”和“中药配位化学学说”
2种理论具有“整体性”的特点,既能满足中药化学成分分析需要“微观研究”的要求,又符合中医理论的“宏观思
维”,将是名贵中成药的发展方向。
[关键词] 名贵中成药;剂型;质量研究;代谢组学;中药配位化学学说
[中图分类号]R284.1 [文献标识码]A [文章编号]10015302(2008)03033905
  中成药,即中药成方制剂,系指以中医学基础理论为原
则组方,以传统中药材为原料制备的各类制剂。其处方原则
以“辨证施治”和“君臣佐使”为基础,所以不能等同于通常
所说的中药制剂。名贵中成药系指在民间沿用历史悠久、疗
效独特、口碑良好的中成药。此类中成药处方出自古医名
著、宫廷秘方或者祖传秘方,向上秉承古医精髓,向下融入了
现代制药技术。
本研究选择乌鸡白凤丸、安宫牛黄丸、片仔癀、西黄丸、
龟龄集、八宝丹、云南白药、六神丸、苏合香丸、麝香保心丸、
华佗再造丸、定坤丹、古汉养生精等十余种名贵中成药为对
象,从剂型、质量分析等方面出发,对名贵中成药的质量研究
现状加以分析。
1 传统名贵中成药的剂型现状分析
随着科技的发展,特别是现代制药技术成果广泛运用于
传统制剂中,中成药剂型的发展会越来越快。刘本俊等[1]通
过对286个中成药价格与剂型之间的关系分析得出,技术含
量低的传统剂型,在价位上无法与技术含量高的现代剂型相
比,传统剂型中丸剂、煎膏剂、散剂仍被使用但似有淘汰之趋
·923·
第33卷第3期
2008年2月
         
    中 国 中 药 杂 志
ChinaJournalofChineseMateriaMedica
       
Vol.33,Issue 3
February,2008