多源信息融合数字模型
科技前沿与学术评论
多源信息融合数字模型
研究员、博导 岳天祥研究员、博导 刘纪远
(中国科学院地理学与资源研究所, 北京100101)
摘 要:研究结果表明,在目前基础条件下,多源信息融合数字模型的实现需要解决现行数字地面模型和空间插值模型的误差问题、点—面信息有效融合问题、多尺度转换问题和多维GIS面临的理论问题。建立多源信息融合数字模型的基本步骤可归纳为:(a)建立基于曲面论数字模型的基本方
程,(b)运用遥感数据反演数字模型的首次近似表达形式,(c)如果有更新信息,重复以上过程,直至理论模型与实际需求完全相符。
关键词:曲面论 遥感反演 多源信息融合 数字模型
ADigitalModelforMulti-SourcesInformationFusion
Professor YUETianxiang
Professor LIUJiyuan
(InstituteofGeographicalSciencesandNaturalResourcesResearch,CAS,Beijing100101)
Abstract:Ourresearchresultshowsthatrealizationofthedigitalmodelformulti-sourcesinformationfusionneedstosolveproblemsoferrorsofexistingdigitalterrainmodelandspatialinterpolationmodel,virtualfusion
ofpointandsurfaceinformation,informationtransformationatvariousscales,andmulti-dimensionGIS.Thebasicstepsofconstructingthedigitalmodelinclude,(a)establishingbasicequationsofthedigitalmodelbymeansofsurfacetheory,(b)retrievingfirstapproximateformulationusingremotesensingdata,(c)ifthere
aremoreavailableinformation,thestepaboveisrepeateduntilrequirementisreached.Keywords:surfacetheory,remotesensingretrieval,informationfusion,digitalmodel
海洋监视系统和小型战术系统的第一代信息融合系统,它们包括军用分析系统(TCAC)、多平台多传感器跟踪信息相关处理系统(INCA)、全员分析系统(PAAS)、海军战争状态分析显示系统(TOP)、辅助空中作战命令分析专家系统(DAGR)、空中目标确定和截击武器选择专家系统(TATR)、自动多传感器部队识别系统(AMSUI)和目标获取与武器输送系统(TR-WDS)。九十年代研制的主要数据融合系统包括全源信息分
1 引言
七十年代初,美国研究机构发现,利用计算机技术对多个独立的连续声纳信号进行融合后,可以自动检测出敌方潜艇的位置[1]。这一发现使信息融合作为一门独立的技术首先在军事应用中得到青睐,美国相继研究开发了几十个军事融合系统。进入八十年代,研制出了应用于大型战略系统、
中国科学院知识创新工程项目(No.kzcx2-308-02)
23卷5期1
科技前沿与学术评论
可表达为与均值有关的结构成分、与空间有关的随机成分和随机噪音三部分之和;(2)所研究区域是均质的;(3)对景观的不同部分应使用不同的半变图(semivariogram)。然而,GIS软件的使用者往往忽视了这三个重要的前提条件,因此得出的结论可信度很差。其根源在于Kriging模型在理论上不够完善,有待进一步研究;与此同时,到目前为止,在非线性插值和非Gaussian过程方面还没有任何可用的研究成果,这将是我们要深入探讨的主要问题。
2.2 数字地面模型
数字地面模型(DTM)被定义为描述地面属性空间分布的有序数据阵列[4][5]。数字高程模型(DEM)是DTM最重要的基础模型,它被定义为描述高程空间分布的有序数字阵列。DEM由通过离散点上抽样获取的高程数据组成。构建DEM的方法主要有三种:(a)基于不规则三角网(TIN)的DEM常常通过地面特定点的抽样,形成以为坐标的不规则点网,其网格是由三个相邻点构成的不规则多边形;(b)基于规则栅格网(RGN)的DEM与所研究区域的尺度密切相关,其网格是由三个或四个相邻点构成的正三角形或正方形;(c)基于等高线网(CBN)的DEM由数字化等高线组成,其格网是由相邻等高线围成的不规则多边形。
目前,在GIS中应用较广泛的DEM是基于RGN的DEM和基于TIN的DEM。然而,RGN的运算结果缺乏精确性:(1)不易处理高程的陡峭变化,(2)网格的大小影响运算结果和计算效率,(3)在水文分析中,其运算得出的斜坡流路为不符合现实的锯齿形。TIN与RGN相比,使用起来比较灵活,但由于其网格的不规则性,使得实际操作和属性运算都比较困难。也就是说,在上述三种方法中,基于CBN的DEM是最有优势的模型,但它对数字高程数据的数量和质量都有比较高的要求,没有考虑插值点组合的任何曲面性质。
2.3 误差问题
空间插值模型和数字地面模型的误差问题在二十世纪六十年代末就引起了有关专家的重视[6],在近30年的研究中,产生了许多误差检测方法和评价指标
[7][8]
析系统(ASAS)、战术陆军与空军指挥员自动情报保障系统(LENSCE)和敌态势分析系统(ENSCS)。随着计算机技术和地理信息系统的发展,许多算法的复杂性不再成为开发新型信息融合模型的障碍。新的融合方法将不仅仅限于军事问题,它有更普遍的作用和意义。然而,要实现多源信息融合模型,需要解决三方面问题。
2 从根本上解决数字地面模型和空间插值模
型的误差问题
2.1 空间插值模型
在观测点覆盖地区内,估计非观测点或非抽样点属性值的过程,称为空间插值法。目前,主要的空间插值法包括边界线勾画法、趋势面分析法、样条函数法、滑动平均法和Krig-ing插值法
[2]
。
根据景观特征勾画景观单元轮廓的边界线勾画法是最简单的插值法,例如,土壤、地质、植被和土地利用等专题地图以及影像分析中的边界搜寻算法等。
通过多元回归模拟空间渐变的趋势面分析法也是一种便于理解的简单方法,它可以用来处理线性趋势面上的概括性特征,但很难用于二次以上的高次多项式所表达的趋势面。因此,趋势面分析的主要用途通常不是插值,而是在使用其它插值方法之前,提取数据的整体性特征。
样条函数(splinefunctions)是活动曲线规(flexibleruler)的数学表达,它们是少量数据点精确拟合的函数片段,并假定各曲线片段之间是连续的。这就意味着不需要重新运算整条曲线就可以用样条函数修正曲线的一个片段,这对趋势面分析法是不可能的。但样条插值存在着以下缺点:(1)选择已知观测点作为样条函数的重合点得出的结论与选择已知观测点作为样条函数的间隔点得出的结论大相径庭;(2)样条插值的误差不能直接估计;(3)在三维插值中,如何把斑块缝合在一起形成一个完整的表面仍是亟待解决的实际问题。
滑动平均法是一种非常简单的插值法,它通过非观测点邻域中已知点的平均值来估计非观测点的值。其中加权滑动平均法是较常用的一种滑动平均法,已知数据点距插值点的距离越近,其权重越大。滑动平均法的缺点包括:(1)估计值由其估计函数及其参数主导;(2)插值点的正确性取决于已知数据点的数据质量;(3)插值点随所取邻域面积的大小而发生变化,尤其在已知点分布不规则的情况下,每个插值点都可能有不同的大小和形状,因此,插值点的准确性是很值得怀疑的;(4)在许多情况下,已知数据点很不规则,滑动平均所得出的曲面很可能不通过一些已知点,这些点常常被直接输出,因此导致了插值面中的奇点。为了克服滑动平均法的缺点,法国地质数学家G.Matheron和南非采矿工程师D.G.Krige为采矿业发展了一种最优插值法,它是随机表面的最佳线性无偏估计(BestLinearUnbiasedEstimate)。后来成为GIS中大众化的空间数据内插法,人们称其为Kriging插值法[3]。Kriging插值法有三个重要的前提条件:(1)区域变量
。例如,使用
滤波方法提高DEM数据的质量、基于趋势面的误差检测、基于坡度信息的误差检测、基于等高线拓扑关系的误差检测和DEM原始数据误差检测等,量度误差的均方根指数(RMSE)、标准差(SD)、误差直方图和空间自相关系数等。有关误差研究方法可概括为三类:(1)评价和比较各种内插算法的误差,(2)分析和评价原始数据误差,(3)研究数字地面模型主要变量和参数的不完备性所引起的误差。
我们认为,现有的空间插值模型和数字地面模型获取内插点的数值时,将采样点或观测点视为相互独立的点,使得已知点之间的相互关系造到破坏,最终导致比较大的误差甚至错误的结论[9]。虽然在以往的研究,Evans已认识到坡面、坡向和曲率是反映局部规律的重要变量[10],但他并没有将它们用于表达数字地面模型。事实上,根据微分几何学理论,坡面、坡向和曲率只是地面剖面线的决定要素,数字地面
世界科技研究与发展
科技前沿与学术评论
的决定要素是曲面的第一类基本量和第二类基本量。 3 实现点-面信息的有效融合
目前,探测地表过程和生态过程的主要方式有两种:一是单点定位观测或空间采样;二是利用遥感数据进行反演。单点定位观测或空间采样能较为准确地获取高时间分辨率的点上观测数据。但定位观测或空间采样只局限于离散、稀疏的有限空间点,不能直接得到区域尺度的有关参数,只能根据点数据进行估计,或以点代面,或利用简单的插值方法估算。而由于地理系统和生态系统在空间属性上的非均匀性,利用点数据直接描述区域属性不能真实反映空间分布规律。
遥感方法能频繁和持久地提供地理过程和生态过程的面状信息,相对于单纯依靠地面站点定位观测是一个革命性的变化。但遥感描述并不能直接获得过程参数,需要通过遥感信息(光谱信号和雷达信号)结合地表特性进行反演。传统的解决方法是利用简单内插把离散数据空间化,但是这样处理的结果精度较差,大大影响了遥感反演的精度。提高离散参数空间化精度的一个有效途径是引入新的、有力的数学工具,设计全新的数字模型,建立高精度的非线性内插模型。也就是说,运用曲面论基本定理使模型的内插点能够准确反映局部规律;引入遥感参数,使模型能够包容整体规律和采样点或观测点之间的有机联系。基于遥感反演和曲面论基本定理的离散点数据空间化数字模型是实现点面信息有效融合的有效途径。 4 解决多维GIS
5 结论
从理论上讲,曲面论和遥感反演方法的结合可以彻底解决多源空间信息融合原型系统所面临的上述问题。
微分几何的出发点是微积分,微积分在几何上的应用演变为曲线论和曲面论。Gauss是微分几何学的创始人,他将欧氏几何推广到了曲面几何上,建立了曲面论。曲面论[11]是微分几何学的核心内容之一,也是自然要素数字模型的理论基础。曲面论基本定理可概述为:则除运动不计外,空间曲面由第一基本形式和第二基本形式完全确定。
在卫星遥感中,反演是指利用卫星探测数据对含有特定信息微分方程求解、确定特定参数的过程。在反演过程中,必须合理地添加有用的信息,同时应尽量避免浪费有用的信息[12]。目前,在遥感反演方法和模型方面已有许多研究成果。例如,线性混合像元模型(inversionoflinearmixturemod-[13][14]
el),迭代反演法(iterativeinversionapproach),经验散射[15]模型(empiricalscatteringmodel),基于人工神经网络的反演[16]方法(inversionapproachbasedonartificialneuralnetwork),基
于分形的相干散射模型(fractal-basedcoherentscatteringmod-[17]el),基于数理统计的反演方法(inversionapproachbasedon
statistics)
[18]
和基于半经验后向散射模型的反演方法(inversion
[19]methodbasedonsemiempiricalbackscattermodel)。
建立多源信息融合数字模型的基本步骤可归纳为:运用微分几何学曲面论基本定理,建立数字模型的基本方程;运用遥感数据,反演数字模型的一次近似表达形式;根据点上的采样数据或观测数据,修正遥感反演模型,实现数字模型
面临的理论问题和多尺度转换
的第二次近似表达形式;如果有更新的遥感数据和点上的采样数据或观测数据,重复以上过程,直至理论模型与实际情况完全相符。二十一世纪信息处理技术的首要任务是应用现代信息技术提高信息的生产力,解决决策中的复杂性问题。信息融合是为了实现上述任务目标的一种重要方法。虽然信息融合方法的初衷是军事问题,但是它在非军事领域的应用是它的必然趋势。作者欲通过研究开发多源空间信息融合数字模型原型系统,弥合军事应用和非军事应用之间的带沟,促进军事信息技术和非军事信息技术之间的相互转换。
问题
地理信息系统软件平台发展到今天仍是面向二维的空间数据,对三维空间数据和时序数据缺乏有效的数字模型,更难以组织和处理包括空间数据在内的多维数据。目前,虽然有三维地理信息系统问世,但理论完善的数字模型一直没有解决,其功能远远不能满足实际问题的需要。
对地球表面上被定义为生态系统的任何大小的斑块,都有一系列外部因素影响其功能行为,进而产生了生态系统与其外部环境的能流和物流。在每个尺度,为了刻划生态系统类型的范围、压力、状况、价值和趋势以及土地利用格局和基本的生态系统结构特征,需要确定一套空间尺度明确的指标。在较大尺度进行的研究包含了这些能流和物流的较详细信息,然而却损失了政策制定者常常所需要的地域特性。也就是说,不管在那个单一尺度,都不可能系统地认识土地利用和土地覆被动态的所有细节。因此,多尺度转换模型研究势在必行。在某种意义上,尺度转换问题可以看作是点-面信息融合的一个特例。例如,一个AVHRR像元可以描述为由269个MSS像元组成的曲面;一个MSS像元可描述为5个TM像元组成的曲面;一个TM像元可描述为由36个IRS像元组成的曲面。23卷5期参考文献
[1]康耀红.数据融合理论与应用.西安:西安电子科学大学出版
社,1997
[2]岳天祥,刘纪远.第四代地理信息系统研究中的尺度转换数字模
型.中国图象图形学报,2001,6A(9):907-911
[3]SteinML.InterpolationofSpatialData.NewYork:Springer-Verlag
NewYork,Inc.,1999
[4]MooreID,GraysonRB,LadsonAR.Digitalterrainmodelling:are-viewofhydrological,geomorphological,andbiologicalapplication.In:BevenKJandMooreID(eds).TerrainAnalysisandDistributedmod-ellinginHydrology.Chichester:JohnWiley&Sons,1992,7-34
3
科技前沿与学术评论
[14]PolatinPF,SarabandiK.Aniterativeinversionalgorithmwithapplica-tiontothepolarimetricradarresponseofvegetationcanopies.IEEETransactionsonGeoscienceandRemoteSensing,1994,32:62-71[15]OhY,SarabandiK,UlabyFT.Anempiricalmodelandaninversion
techniqueforradarscatteringfrombaresoilsurfaces.IEEETransactionsonGeoscienceandRemoteSensing,1992,30:370-381
[16]PierceL,SarabandiK,UlabyFT.Applicationofartificialneuralnet-workincanopyscatteringinversion.InternationalJournalofRemoteSensing,1994,15:3263-3270
[17]LinYC,SarabandiK.AMonteCarlocoherentscatteringmodelforfor-estcanopiesusingfractal-generatedtrees.IEEETransactionsonGeo-scienceandRemoteSensing,1999,37:440-451
[18]LinYC,SarabandiK.Retrievalofforestparametersusingafractal-basedcoherentscatteringmodelandageneticalgorithm.IEEETransac-tionsonGeoscienceandRemoteSensing,1999,37:1415-1424[19]KurvonenL,PulliainenJ,HallikainenM.Retrievalofbiomassinbore-alforestsfrommultitempotalERS-1andJERS-1SARimages.IEEETransactionsonGeoscienceandRemoteSensing,1999,37:198-205
社,2000
[5]陈述彭,鲁学军,周成虎.地理信息系统导论.北京:科学出版[6]CrainIK.Digitalrepresentationoftopographicsurfaces.Photogrammetric
EngineeringandRemoteSensing,1970,54(11):1577-1580[7]TangG.AResearchontheAccuracyofDigitalElevationModels.Bei-jing:SciencePress,2000
[8]李志林,朱庆.数字高程模型.武汉:武汉测绘科技大学,2000[9]WiseS.GISdatamodelling-lessonsfromtheanalysisofDTMs.Interna-tionalJournalofGeographicalInformationScience,2000,14(4):313-318
[10]EvansIS.Anintegratedsystemofterrainanalysisandslopemapping.
ZeitschriftfuerGeomorphologie,Suppl-Bd,1980,36:274-295[11]苏步青,华宣积.实用微分几何学.北京:科学出版社,1998[12]陈述彭.遥感大词典.北京:科学出版社,1990,41
[13]CherchaliS,AmramO,FlouzatG.Retrievaloftemporalprofilesofre-flectancesfromsimulatedandrealNOAA-AVHRRdataoverheteroge-neouslandscapes.InternationalJournalofRemoteSensing,2000,21(4):753-775
作者简介
岳天祥(1963-),男,研究员,博士生导师,长期从事资源环境、数字模型与地理信息系统研究。
刘纪远(1947-),男,研究员,博士生导师,中国科学院地理科学与资源研究所所长,长期从事资源环境信息分析研究。
(责任编辑:曙光)
国外新闻
日本大幅度修改航天器发射计划
日本宇宙开发事业团决定大幅度修改卫星等航天器的发射计划。据悉,主要原因是经费短缺和火箭技术存在问题。
据该机构公布的材料说,它在前不久举行的会议上决定把国际空间站上日本实验舱的发射时间推迟1年,即从原定的2003年度(到2004年3月31日结束)延期到2004年度,为这座实验舱运送物资的货运飞船“HTV”也从2004年度向后推迟1年发射。
日本参加了正在建设的国际空间站的部分工作,预备在它的实验舱里进行细胞培养和放射线检测等生命科学和生
物技术实验。
该会议决定推迟发射计划的还有原定在2003年度发射的“H2A”加强式国产大型火箭,它安装有助推器,将推迟2年发射;原定在2003年度和2004年度发射的技术试验卫星“ETS—8”、陆地观测卫星“ALOS”和绕月卫星“SELENE”等也均延期1年;新型数据中转技术卫星“DRTS—E”和“J—1”小型火箭停止开发。
日本文部科学省的官员说,航天计划的大幅度修改主要原因是经费削减。但据认为,国产大型火箭“H2A”的研究开发进展得不顺利也是一个重要原因。(新华社供本刊稿)
世界科技研究与发展