顾及空间关系约束的不同比例尺面目标不一致性同化处理_赵彬彬

2024年2月20日发(作者：狭片鹿角珊瑚)

第４１卷第７期０１６年７月２武汉大学学报·信息科学版ｅｏｍａｔｉｃｓａｎｄＩｎｆｏｒｍａｔｉｏｎＳｃｉｅｎｃｅｏｆＷｕｈａｎＵｎｉｖｅｒｓｉｔＧ　　　　　　ｙＶｏｌ．４１Ｎｏ．７Ｊｕｌ０１６２ｙ　：／ＯＩ１０．１３２０３．ｗｈｕｉｓ２０１４００１１Ｄｊｇ（）文章编号：６７１８６０２０１６０７９１１７１８００－－－顾及空间关系约束的不同比例尺面目标不一致性同化处理赵彬彬１　彭东亮２　张山山２　刘姗姗１　熊旭平１　戴全发１湖南长沙，１　长沙理工大学交通运输工程学院，１０１１４４湖南长沙，２　中南大学地理信息系，１００８３４摘　要：不一致性问题是长期困扰地理信息科学领域的一个难题，亦是衡量空间数据质量的基本指标。不一致性问题直接影响空间数据集成、制图综合和空间数据更新。当前研究较多关注于相同或相近比例尺线目标之间的不一致性问题，而对不同比例尺地图数据或制图综合过程中不一致性处理问题的研究甚少。为此，本文顾及空间目标之间拓扑和距离关系约束，提出了一种基于Ｍｏｒｈｉｎｐｇ变换的不同比例尺地图面目标之间不一致性的同化处理方法，通过对两组不一致性同化处理的对比实验验证了该方法的有效性和实用性。关键词：尺度；空间关系；不一致性；同化ｏｒｈｉｎＭｐｇ变换；中图法分类号：Ｐ２０８；Ｐ２８３　　文献标志码：Ａ　　空间数据不一致性是指在空间数据对象之间存在明显的矛盾或冲突，该研究始于２０世纪９０年代初，９９５年召开的国际空间数据质量会议提１１］，出将一致性作为空间数据质量五大特征之一［对应关系来评价结构不一致性，并通过分析复杂面目标对中两个面目标之间的拓扑关系来评价拓扑不一致性，该方法亦可用于评价不同比例尺地６］。Ｂ图复杂面目标的不一致性［ｒｉｓａｂｏａ等利用数标志着空间数据不一致性问题已成为国际地理信息科学领域的一个基础研究课题。近２在０年来，空间数据之间的不一致性特征描述、分类、表达、度量、探测、评价及处理等方面取得了丰富的成２，３］，但大多数集中在相同（果［或相近）比例尺空据集中目标之间的拓扑不变性定义了一组完整性约束，进而基于约束来度量数据集中每一个点、线和面目标与拓扑约束之间的冲突程度，即不一致７］。上述研究基本是对某一种或两种不一性程度［致性在探测、评价或度量方面进行浅尝辄止的探索性研究，尚未涉及有效的多尺度地图数据不一致性处理等深层次问题。随着制图综合、多尺度地图数据协同应用的深化，对不同比例尺地图空间数据之间不一致性的研究将更具现实意义和应８］。实践中，用价值［许多制图综合算子（如“简间数据之间的不一致性问题上。在多尺度地理空］间数据不一致性方面仅有少量涉及，文献［针对４”操作引起的拓扑不制图综合过程中的“ｏｌｌａｓｅＣｐ一致性，按拓扑属性进行严密分类建立了一套分”操作类规则，提出了基于分类规则的“ｏｌｌａｓｅＣｐ引起的不一致性的评价和探测方法。Ｃｈｅｎ等针对地图更新中河流“爬坡”的冲突问题，提出了基于拓扑、方向和度量关系的集成空间关系模型拓扑链，精细地描述线／线目标空间关系，用于探测５］。Ｄ河流与等高线之间的不一致性问题［ｕ等运、“光滑”和“聚合”等）会改变空间目标几何形化”状、维度和图形结构等细节特征，从而使空间目标６］。的不同比例尺表达发生变化，进而产生不一致［，“（）如图１简化”操作后导致两面目标之间方向ａ（）关系不一致。图１中，对河流的“光滑”操作导ｂ致河流与建筑物之间拓扑不一致，即建筑物“落）入”河流中。本文主要针对图１（的河流与建筑ｂ物之间的不一致性处理问题进行探讨。”两种制图综合算子对和“用“ｅｒｉｎＤｒｏｉｎＭｇｇｐｐｇ”宽边界复杂面目标进行综合，由生成的不同比例尺复杂面目标结构邻域建立复杂面目标之间的收稿日期：１１０１４２２２－－）；；）项目资助：湖南省自然科学基金（国土环境与灾害监测国家测绘地理信息局重点实验室开放国家自然科学基金（４ＪＪ３０８３１４１３０１４０４）。研究基金（ＥＤＭ２０１２Ｂ０２Ｌ第一作者：赵彬彬，博士，副教授，主要从事多源、多尺度地图数据不一致性探测处理理论与方法研究。ｚｂｂｓｉｒ６３．ｃｏｍ＠１

９１２武汉大学学报·信息科学版２０１６年７月法。步骤是：首先，通过叠置、求差和求交等空间１　问题及图形策略（）分析图１可知，由较大比例尺河流图形综ｂ合为较小比例尺河流时，对河流边界的综合涉及细节概括及几何形状，河流边界与建筑物由综合，产生不一致性。前的“相邻”变为综合后的“相交”为了解决此类问题，本文采取的策略是：基于数据同化思想，以Ｍ在河流ｏｒｈｉｎｐｇ变换为技术手段，与建筑物发生不一致性的区域，顾及空间目标之间的拓扑和距离关系约束，对不同比例尺数据中的两条河流进行局部边界的同化，提出一种不同比例尺地图面目标之间的不一致性同化处理方操作获得与建筑物面目标存在不一致性的两个不同比例尺地图数据中河流面目标各自的边界；然后，分别以较大比例尺数据中的河流边界（与建筑物不冲突）和较小比例尺数据中的河流边界（与建筑物冲突）作为变换端点，利用Ｍｏｒｈｉｎｐｇ技术在较大比例尺河流边界和较小比例尺河流边界之间插值，进行河流边界局部同化；进而，根据插值河流边界与建筑物面目标之间的拓扑和距离关系约束决定河流边界局部同化程度，即插值河流边界的位置，以实现对制图综合中建筑物等面目标“落入”河流中这类不一致性进行同化处理的目的，流程如图２。图１　制图综合引起的不一致性ＦＧｅｎｅｒａｌｉｚａｔｉｏｎＯｅｒａｔｉｏｎｉ．１　ＩｎｃｏｎｓｉｓｔｅｎｃｉｅｓＣａｕｓｅｄｂ　　　ｐｇｙ　图２　河流与建筑物间不一致性同化处理框架Ｆｉ．２　ＡＦｒａｍｅｗｏｒｋｆｏｒＡｓｓｉｍｉｌａｔｉｎＩｎｃｏｎｓｉｓｔｅｎｃｉｅｓＢｅｔｗｅｅｎＲｉｖｅｒｓａｎｄＢｕｉｌｄｉｎｓ　　　　　　　ｇｇｇ　

１卷第７期　第４赵彬彬等：顾及空间关系约束的不同比例尺面目标不一致性同化处理９１３应用广泛，如动画制作、图像压缩及医学图像重构２　不同比例尺面目标不一致性同化处理相对于较小比例尺地图数据而言，较　　一般地，大比例尺地图数据具有较高的精度和较强的现势性；相对于自然地物而言，人工地物亦具有较高的精度和较强的现势性。为了兼顾较大比例尺地图数据的详细性、精确性和较小比例尺地图数据的概括性与全局性，本文遵循自然地物让位于人工地物的原则，以数据同化思想为指导对较大比例尺河流和较小比例尺河流边界局部进行同化处理。２．１　数据同化数据同化是一种多源数据和模型集成的数据９］，旨在考虑数处理技术，起源于大气与海洋科学［１７］，其最初目的是实现计算机图像的无缝平滑等［１８］。由于采用连续变形思想与多尺度地图渐变［连续综合不谋而合，亦更接近人的视觉感观且符１６］，因而日渐受到制图综合领域学者合人的认知［们的重视。目前，针对不同比例尺地图中道路、河流等线目标的形状内插多采用基于矢量数据的，）在不同比例尺原始ｏｒｈｉｎａＭｐｇ变换。如图４（目标Ｃ和Ｄ之间沿直线移位路径进行Ｍｏｒｈｉｎｐｇ变换插值时，插值结果越移位距离值ｔ越接近０，移位距离值ｔ越接近１，接近原始目标Ｃ；插值结果越接近原始目标Ｄ。当移位距离ｔ＝０时，Ｍｏｒ－ｈｉｎｐｇ插值结果为原始目标Ｃ；当移位距离ｔ＝１时，Ｍｏｒｈｉｎｐｇ插值结果为原始目标Ｄ。图４（）为采用直线移位路径时对应不同移位距离值ｂｔ的插值结果。２．３　空间关系约束由较大比例尺河流综合为较小比例尺河流时，河流与建筑物之间的拓扑关系由综合前的“相，从而导致建筑物“离”变为综合后的“相交”落入”河流中，产生不一致。因此，在对不同比例尺河流边界进行同化处理时，必须顾及河流与建筑物之间的拓扑关系约束，即保持同化后的河流与建筑，，）如图５（河流移位距物之间拓扑关系“合理”ａ离值ｔ其建筑物之间拓＝０．７５对应的同化结果，，显然不满足河流与建筑物间的拓扑关扑“相交”，（）不同的移位距离值对应不同系约束。如图４ｂ），）为了最大程度的河流边界同化结果（如图５（ｂ地契合综合后的较小比例尺河流边界的趋势，则据时空分布及观测场和背景场误差基础上，在数值模型动态运行过程中融合新观测数据的方１０］，其目标是使预报结果与现实趋于一致。目法［１１］、空间数据前，数据同化技术在遥感影像融合［１１２］３］和土地利用变化监测［等地理信息融合更新［科学领域亦取得了初步进展，被认为是解决空间１４］。“同化”数据不一致问题的一条可行思路［的本意是指把不同的事物变得相近或相同的过程。如图３，按两同化对象的同化程度、同化结果、接“整体’近程度，这里指Ａ和Ｂ‘之Ａ和Ｂ同化了（１５］：”———包括五层含义［间的同化）Ｂ①Ａ不变，；变得和Ａ相近或相同，如图３（ｂ）②Ａ基本不；）变，如图３（Ｂ向Ａ同化，ｃ③Ａ和Ｂ以相互靠；拢的方式变得相近或相同，如图３（ｄ）④Ｂ基本；）不变，如图３（Ａ向Ｂ同化，Ａ变得ｅ⑤Ｂ不变，。（）和Ｂ相近或相同，如图３ｆ图４　Ｍｏｒｈｉｎｐｇ变换Ｆｉ．４　ＡｎＥｘａｍｌｅｏｆＭｏｒｈｉｎＴｒａｎｓｆｏｒｍａｔｉｏｎ　　　ｇｐｐｇ　图３　同化的含义Ｆｉ．３　ＩｌｌｕｓｔｒａｔｉｏｎｏｆＡｓｓｉｍｉｌａｔｉｏｎ　　ｇ２．２　Ｍｏｒｈｉｎｐｇ变换为了使同化结果尽可能多地保留不同比例尺河流边界的几何形状等特征，本文以Ｍｏｒｈｉｎｐｇ变换技术为手段进行插值。Ｍｏｒｈｉｎｐｇ变换是一种同时顾及形状和颜色的图像内插技术，亦可以理解为“交叉溶解”过程［６］１图５　顾及空间关系约束的不一致性同化Ｆｉ．５　ＩｎｃｏｎｓｉｓｔｅｎｃＡｓｓｉｍｉｌａｔｉｏｎｗｉｔｈＲｅｓｅｃｔ　　ｇｙｐ　ｔｏＳａｔｉａｌＲｅｌａｔｉｏｎＣｏｎｓｔｒａｉｎｔｓ　　　ｐ，在计算机可视化领域

９１４武汉大学学报·信息科学版２０１６年７月必须顾及同化后河流边界与建筑物之间的距离约束，即满足不同比例尺地图相关规范以及人眼对１９，２０］（本文取０．图上最小距离分辨率的要求［１用两组不同比例尺的河流和居民地数据进行实验，）（并与Ｄ图６ｅｌａｕｎａｙ三角网骨架线中轴化法进行对比。３．１　不一致性同化处理用两个实验说明不一致性同化处理方法。实验一　采用比例尺１∶２０００数据的河流和　建筑物及综合后的比例尺１∶１万数据的河流面目。由较大比例尺河））和６（标进行实验，如图６（ｂａ流综合为较小比例尺河流时，因河流边界“简化”导。（）致河流与建筑物之间发生拓扑不一致，如图６ｃ（）。图５中，移位距离值ｔ＝０．２５所对应的ｍｍ）ｂ河流边界同化结果与建筑物之间的图上距离ｄＭａｐ约为６．亦不满足距离关系约束。２ｍｍ，３　实验分析１９，２０］，采按照§２提出的策略，参照相关规范［图６　河流与建筑物之间的不一致性ｉ．６　ＩｎｃｏｎｓｉｓｔｅｎｃＢｅｔｗｅｅｎａＲｉｖｅｒａｎｄａＢｕｉｌｄｉｎＦ　　　　　ｇｙｇ　（）根据图６得到两不同比例尺河流之ａ　　首先，），）进一步获得两不差，边界不一致区域（图７（ａ）。进而，同比例尺河流边界（图７（通过Ｍｏｒｂ）－ｈｉｎｐｇ变换选取不同的移位距离值ｔ获得不同程），）显然，度的同化结果（图７（当ｔ∈ｃ）ｆ～７（时，同化结果都不满足两面状地物“相０．６，０］１．［时，同化结果离”的拓扑关系约束；当ｔ０，６）０．∈［满足拓扑“相离”的约束。继而，通过距离关系约同化结果与束确定移位距离值，当ｔ＝０．３６３时，（）如图７建筑物之间的图上距离ｄＭａ．１ｍｍ，ｇｐ为０。从插值河流边界与较大比例尺河流（和图７ｈ）边界的形态接近程度来看，同化处理结果较好地）。（）体现了同化本质的第三种含义（如图３ｄ实验二　采用比例尺１∶１万数据的河流和建筑物及综合后的比例尺１∶５万数据的河流面。目标进行实验，如图８按两不同比例尺河流面目），）进标之差，获得两河流边界不一致区域（图８（ａ）。进）一步获得两不同比例尺河流边界（图８（ｂ而，通过Ｍｏｒｈｉｎｐｇ变换选取不同的移位距离值ｔ）。显然，（）（）获得不同程度的同化结果（图８ｃｆ～８，时，同化结果都不满足两面目标当ｔ０．４．０］１∈［，“时，同化结相离”的拓扑关系约束；当ｔ０．４）０∈［果满足拓扑“相离”的约束。继而，通过距离关系约同化结束进一步确定移位距离值，当ｔ＝０．０７９时，果与建筑物之间的图上距离ｄＭａ．１ｍｍｐ为０（（））。从插值河流边界与较大比例（）图８和图８ｈｇ尺河流边界的形态接近程度来看，同化处理结果较）。（）好地体现了同化本质的第二种含义（图３ｃ图７　河流与建筑物间不一致性同化处理实验一Ｆｉ．７　ＴｈｅＦｉｒｓｔＥｘｅｒｉｍｅｎｔｏｆＡｓｓｉｍｉｌａｔｉｎＩｎｃｏｎｓｉｓｔｅｎｃＢｅｔｗｅｅｎａＲｉｖｅｒａｎｄａＢｕｉｌｄｉｎ　　　　　　　　　ｇｐｇｙｇ　　

１卷第７期　第４赵彬彬等：顾及空间关系约束的不同比例尺面目标不一致性同化处理９１５图８　河流与建筑物不一致性同化处理实验二Ｆｉ．８　ＴｈｅＳｅｃｏｎｄＥｘｅｒｉｍｅｎｔｏｆＡｓｓｉｍｉｌａｔｉｎＩｎｃｏｎｓｉｓｔｅｎｃＢｅｔｗｅｅｎａＲｉｖｅｒａｎｄａＢｕｉｌｄｉｎ　　　　　　　　　ｇｐｇｙｇ　　３．２　Ｄｅｌａｕｎａｙ三角网骨架线中轴化法（）（图９分别为Ｄｅｌａｕｎａａｈ）～９ｙ三角网提取骨架线法对上述实验一和实验二中河流与建筑物之间的不一致区域进行处理，通过提取Ｄｅｌａｕｎａｙ２１］。图９（）中，对于实验三角网骨架中轴线实现［ｄ扑关系约束，亦不满足数据规范的最小距离关系约束。由此可知，ｅｌａｕｎａＤｙ三角网骨架线中轴化法存在明显不足，骨架中轴线与两个不同比例尺河流边界在几何形态上相似度低，亦无法有效地处理此类不一致性。通过对不同比例尺地图河流面目标边界与建筑物之间的不一致性同化处理结果的对比可以看出：同化处理后的河流边界与两个不同比例尺河流边界相似度高，在几何形态上较好地兼顾了两种不同比例尺河流边界的特征，即较大比例尺河流边界的弯曲细节和较小比例尺河流边界的综合，处理结果更合理。同时，通过同化处理后的河流边界特征亦较明显地体现了同化本质的具体含义。一，中轴线与建筑物之间满足拓扑“相离”条件，但红圈所两者之间的图上距离ｄＭａ．０３９ｍｍ（ｐ为０，小于人眼对图上最小距离分辨率０．示）１ｍｍ，，从人眼角度来看两者之间为拓扑“即不满相接”）。对足数据规范的最小距离关系约束（图９（ｈ）于实验二，ｅｌａｕｎａＤｙ三角网骨架线中轴与建筑物（，红圈所示）不满足原之间的拓扑关系为“相交”较大比例尺地图中河流与建筑物之间“相离”的拓图９　Ｄｅｌａｕｎａｙ三角网骨架线中轴化处理结果Ｆｉ．９　Ｒｅｓｕｌｔｓｆｒｏｍ　ＡｘｉａｌＳｋｅｌｅｔｏｎｏｆＤｅｌａｕｎａＴｒｉａｎｕｌａｔｉｏｎ　　　　ｇｙｇ　调一致的地理空间数据的巨大需求，随着制图综４　结　语高速发展的国民经济催生了对实时准确、协合技术在空间数据服务中的应用广泛深入，由此产生的不一致性问题降低了空间数据质量，极大地影响了地理信息的协同应用。为此，本文分析

９１６武汉大学学报·信息科学版（）：ｎｃｅ，０１４，２８１５６２ｅ２８－２０１６年７月了较大比例尺地图综合为较小比例尺地图过程中制图综合算子对空间目标形态、几何形状等的影响。针对建筑物“落入”河流中这类不一致性，依据数据同化原理，以Ｍｏｒｈｉｎｐｇ变换为技术手段，提出了一种顾及河流和建筑物之间拓扑和距离关系约束的不一致性同化处理方法，该方法能够有效地继承较大比例尺地图空间数据的丰富性，同时较好地兼顾较小比例尺地图空间目标概括性。两组不同比例尺的河流、建筑物目标之间的不一致性进行同化处理实验，结果验证了该方法的有效性。与Ｄｅｌａｕｎａｙ三角网骨架线中轴化法对比，顾及空间关系约束的不一致性同化处理方法优势明显，亦为类似不一致性问题的处理提供了方法参考。后续工作将进一步深入研究不同类型空间目标之间的不一致性及其处理方法，以服务于地图数据库更新和维护，以及多源多尺度数据集成和协同应用。参　考　文　献［］Ｊ］．ＥａｉｎｚＷ．ＬｏｉｃａｌＣｏｎｓｉｓｔｅｎｃｌｅｍｅｎｔｓｏ１　　　　Ｋｇｙ［ｆ（）：ＤＱ１１ｔｉａｌａｔａ　ｕａｌｉｔ９９１５１０９３７Ｓａ　－ｙ，ｐ［］２ＩｓｓｕｅｓｉｎＳａｔｉａｌＤａｏｄｒíｕｅｚＭ　Ａ．Ｉｎｃｏｎｓｉｓｔｅｎｃ－　　　　Ｒ　ｐｇｙ　［］ａｂａｓｅｓＪ．ｎｃｏｎｓｉｓｔｅｎｅｏｌｅｒａｎｃｅｅｃｔｕｒｅ　ｏｔｅｓｔＩＴＬＮ　ｙ　２ＣＳ２ｉｎ　ｏｍｕｔｅｒｃｉｅｎｃｅ，００５，３３００：２３７６９　－　ｐ，［］，，ＤＭｉｎｅｔａｌ．Ｄｅ３ＺｈａｏＢｉｎｂｉｎｅｎｉａｎＣａｎｌｉａｎ－　　Ｊ　　ｇｇ　ｔｓａＤｅａｉｌｅｄＲｅｖｉｅｗｏｎＧｅｏｔｉａｌＤａｔａＩｎｃｏｎｓｉｓｔｅｎｃ－－　　　　　ｐｙ　ｔＥＪ］．ＣａｎｄＨａｎｄｌｉｎｅｃｔｉｎｏｍｕｔｅｒｎｉｎｅｅｒｉｎ　　ｇｇ［ｐｇｇ　（）（：简灿良，ｌｉｃａｔｉｏｎ，０１３，４９１０１５０５５Ａ１２ａｎｄ　－ｐｐ赵彬彬，邓敏，等．地理空间数据不一致性探测处］理方法研究［０１３，４９Ｊ．计算机工程与应用，２：）（）５０５５１１１０－［］４Ｈ　Ｋ，ＬｉＫＪ．ＡｓｓｅｓｓｉｎＴｏｏｌｏｉｃａｌＣｏｎｓｉｓａｎ　　　　Ｋ－ｐｇｇｇ　　ｔｆｏｒＣｏｌｌａｓｅＯｅｒａｔｉｏｎｉｎＧｅｎｅｒａｌｉｚａｔｉｏｎｏｆｅｎｃ　　　　　ｐｐｙ　［］，Ｚ８ｉＤｅｒｅｎｈａｎＪＸ，ＷｕＨｕａｉ．ＳａｔｉａｌＤａｔａ　Ｌ　　　　ｇｙｐ　］ＱｕａｌｉｔａｎｄＢｅｏｎｄ［Ｊ．ｎｔｅｒｎａｔｉｏｎａｌｏｕｒｎａｌＩＪｏ　　　ｙｙｆ　：ｅｏｒａｈｉｃａｌｎｏｒｍａｔｉｏｎｃｉｅｎｃｅ，０１２，２６（１２）ＩＳ２Ｇ　　ｇｐｆ２２２２７７９０－　　［］９ａｂｉｅｒＦ．ＯｖｅｒｖｉｅｗｏｆＧｌｏｂａｌＤａｔａＡｓｓｉｍｉｌａｔｉｏｎＤｅ－　Ｒ　　　　　　ｒｅｌｏｍｅｎｔｓｉｎＮｕｍｅｒｉｃａｌＷｅａｔｈｅｒｅｄｉｃｔｉｏｎＣｅｎｔｒｅｓｖ　　　　－ｐｐ［］ＱＪｏｔＲＭｕａｒｔｅｒｌｏｕｒｎａｌｈｅｏａｌ　ｅｔｅｏｒｏｌｏｉＪ．　　－ｆｙ　ｙｇ　ａｌｏｃｉｅｔ００５，１３１：３２１５３３２ｃＳ３２　　－　ｙ，［１０］ＭｃｌａｕｈｌｉｎＤ．ＲｅｃｅｎｔＤｅｖｅｌｏｍｅｎｔｉｎＨｄｒｏｌｏｉｃ　　　　ｇｐｙｇＤａｔａＡｓｓｉｍｉｌａｔｉｏｎ［Ｊ］．Ｒｅｖｉｅｗｓｅｏｈｓｉｃｓ，ｏＧ　　ｐｙｆ　（）：９９５ｓｕ９７７８４９１－ｐｐ［］Ｃ，，Ｗ，１１ｈｅｎＲｏｎｕａｎＬｉｕＧｕｏｉｎａｎＬｅｉｕａｎｅｔ　　ｇｙｙｇｇｇｇ　ａｌ．ＦｕｓｉｏｎＡｌｏｒｉｔｈｍｏｆＭｕｌｔｉｓｅｃｔｒａｌａｎｄＰａｎｃｈｒｏ－　　　　　ｇｐ］ａｔｉｃＩｍａｅｓＢａｓｅｄｏｎｔｈｅＤａｔａＡｓｓｉｍｉｌａｔｉｏｎ［Ｊ．ｍ　　　　　　ｇＧａＩＳｏＷＵｅｏｍａｔｉｃｓｎｄｎｏｒｍａｔｉｏｎｃｉｅｎｃｅｕｈａｎ　ｎｉ　　　　－ｆ　ｆ：陈荣元，刘国英，ｅｒｓｉｔ００９，３４（８）９１９２２（ｖ９２－ｙ，王雷光，等．基于数据同化的全色和多光谱遥感影］像融合［００９，３４Ｊ．武汉大学学报·信息科学版，２（）：）８９１９２２９－，［］Ａ，，Ｘｉａｏｍｉｎｅｔａｌ．Ａｃ１２ｎＸｉａｏａＳｕｎＣｈｕｎＺｈａｎ　　　－ｇｙｇ　ａｎｄＡｌｉｃａｔｉｏｎＡｎａｌｓｉｓｏｆＭｕｌｔｉＵｄａｔｉｎｉｖｅｌｔ　　　　－ｐｇｐｐｙｙ　　］ＧＪｏＧｅｏｓａｔｉａｌＤａｔａＡｓｓｉｍｉｌａｔｉｏｎ［Ｊ．ｏｕｒｎａｌｅｏ　－　　ｐｆ　：安晓ｎｏｒｍａｔｉｏｎｃｉｅｎｃｅ，０１０，１２（４）５４１４８（２ＩＳ５　－ｆ亚，孙群，张小明，等．多源地理空间数据同化的主动更新与应用分析［Ｊ］．地球信息科学学报，（）：）０１０，１２４５４１４８５２－，［］Ｌ，，ｅｔａｌ．ＡｓｓｉｍｉｌａＹｉｈａｎＬｉｕＸｉａｏｉｎ１３ｉＸｉａＺｈａｎ　　　－ｇｐｇ　ｔＰｒｏｃｅｓｓＣｏｎｔｅｘｔＩｎｆｏｒｍａｔｉｏｎｏｆＣｅｌｌｕｌａｒＡｕｔｏｍａｔａｉｎ　　　　　ｇ　ＬａｎｄＵｓｅＣｈａｎｉｎｔｏＣｈａｎｅＤｅｔｅｃｔｉｏｎｆｏｒＭｏｎｉｔｏｒｉｎ－　　　　　　ｇｇ　［］ｅＩＪｏＧＩｓＪ．ｎｔｅｒｎａｔｉｏｎａｌｏｕｒｎａｌｅｏｒａｈｉｃａｌｎｏｒ－　　　ｇｆ　ｇｐｆ，（）：ａｔｉｏｎｃｉｅｎｃｅ，０１２２６９１６６７８７２ｍＳ１６　　　－：［］Ｌ，Ｈ１４ｉＸｉｎｕａｎＣｈｕｎｌｉｎ．ＤａｔａＡｓｓｉｍｉｌａｔｉｏｎＡ　Ｎｅｗ　　ｇ　ｓＭｅａｎｓｆｏｒＭｕｌｔｉｏｕｒｃｅＧｅｏｓａｔｉａｌＤａｔａＩｎｔｅｒａｔｉｏｎ－　　　　　ｐｇ［ｃｉｅｎｃｅｎｄ　ｅｃｈｎｏｌｏｅｖｉｅｗ，２００４，２５ａＴＲＪ］．Ｓ　ｇｙ　（）：（李新，黄春林．数据同化：一种集成多１２１３６１－］源地理空间数据的新思路［００４，２５Ｊ．科技导报，２：）（）１３６１１２－［—，１５］ＷａｎＹｕｅｓｈａｎ．ＤａｔａＡｓｓｉｍｉｌａｔｉｏｎＩｔｓＣａｕｓｅＩｔｓ　　ｇ　ｅｃｔｕｒｅ　ｏｔｅｓｎ　ｏｍｕｔｅｒＮｉＣＳａｔｉａｌＤａｔａｂａｓｅｓ［Ｊ］．Ｌ　　ｐｐＳｃｉｅｎｃｅ，００５，３７７０，２４９５８２２－［］，Ｌ，Ｌ，５ｈｅｎＪｕｎｉｕ　ＷｅｉｉＺｈｉｌｉｎｅｔａｌ．Ｄｅｔｅｃｔｉｏｎｏｆ　Ｃ　　　　ＳａｔｉａｌＣｏｎｆｌｉｃｔｓＢｅｔｗｅｅｎＲｉｖｅｒｓａｎｄＣｏｎｔｏｕｒｓｉｎ　　　　　　ｐ［］ＩＪｏｎｔｅｒｎａｔｉｏｎａｌｏｕｒｎａｌＤｉｉｔａｌＭａＵｄａｔｉｎＪ．　　　ｇｐｐｇｆ　：ＩＳ２Ｇｅｏｒａｈｉｃａｌｎｏｒｍａｔｉｏｎｃｉｅｎｃｅ，００７，２１（１０）　　ｇｐｆ９３１４１１１０－　　［］，，Ｗ，６ｕＳｈｉｈｏｎＱｉｎＱｉｍｉｎａｎＱｉａｏｅｔａｌ．Ｅｖａｌｕａ　Ｄ　　　－ｇｇ　ｔｉｎＳｔｒｕｃｔｕｒａｌａｎｄＴｏｏｌｏｉｃａｌＣｏｎｓｉｓｔｅｎｃｏｆＣｏｍ－　　　　ｇｐｇｙ　　ｅｘＲｅｉｏｎｓｗｉｔｈＢｒｏａｄＢｏｕｎｄａｒｉｅｓｉｎＭｕｌｔｉｅｓｏｌｕｌｒ　　　　　　－－ｇｐｉｏｎＳａｔｉａｌＤａｔａｂａｓｅｓ［Ｊ］．ｎｏｒｍａｔｉｏｎｃｉｅｎｃｅｓ，ｔＩＳ　　　ｐｆ（）：６００８，１７８１５２８２－［］７ｒｉｓａｂｏａＮ　Ｒ，ＬｕａｃｅｓＭ　Ｒ，ＲｏｄｒíｕｅｚＭ　Ａ，ｅｔａｌ．　Ｂ　　　　ｇＡｎＩｎｃｏｎｓｉｓｔｅｎｃＭｅａｓｕｒｅｏｆＳａｔｉａｌＤａｔａＳｅｔｓｗｉｔｈ　　　　　　ｙｐ　ＦＭｅａｎｉｎａｎｄＭａｉｎＰｒｏｃｅｄｕｒｅｓ［Ｊ］．Ｍａｒｉｎｅｏｒｅ　－　　ｇ　（：——王跃山．数据同化—ａｓｔｓ，９９９，１８１）１１０（ｃ２１－］它的缘起、含义和主要方法［１Ｊ．海洋预报，９９９，：）（）１１０２１８１－［］Ｌ，Ｗ１６ｉＺｈｉｌｉｎｏｎＭｅｉｌｉｎ．ＡｎｉｍａｔｉｎＢａｓｉｃＯｅｒａ－　　ｇｇｇｐ　　ｔｉｏｎｓｆｏｒＤｉｉｔａｌＭａＧｅｎｅｒａｌｉｚａｔｉｏｎｗｉｔｈＭｏｒｈｉｎ　　　　　ｇｐｐｇ　［］ｈｅｎｔｅｒｎａｔｉｏｎａｌｒｃｈｉｖｅｓｈｅＩＡｏｔＴｅｃｈｎｉｕｅｓＪ．Ｔ　　　ｑｆ　ＲＳａＳａＩＰｈｏｔｏｒａｍｍｅｔｒｅｍｏｔｅｅｎｓｉｎｎｄ　ｔｉａｌｎ－　　ｐｇｙ，ｇ　）：ｏＳ２６ｒｍａｔｉｏｎｃｉｅｎｃｅ，００８，（Ｂ２６３７４２　－ｆ［］Ｐ，Ｄ，Ｚ１７ｅｎＤｏｎｌｉａｎｅｎＭｉｎｈａｏＢｉｎｂｉｎ．Ｍｕｌｔｉ－　ｇｇｇｇ　　ｓｃａｌｅＴｒａｎｓｆｏｒｍａｔｉｏｎｏｆＲｉｖｅｒＮｅｔｗｏｒｋｓＢａｓｅｄｏｎ　　　　　　ＭｏｒｈｉｎＴｅｃｈｎｏｌｏＪ］．ＪｏｕｒｎａｌｅｍｏｔｅｏＲ　ｐｇｇｙ［ｆ　　ｎｔｅｒｎａ－ＲｅｓｅｃｔｔｏＴｏｏｌｏｉｃａｌＣｏｎｓｔｒａｉｎｔｓ［Ｊ］．Ｉ　　　ｐｐｇｉｏｎａｌｏｕｒｎａｌｅｏｒａｈｉｃａｌｎｏｒｍａｔｉｏｎｃｉｔＪｏＧＩＳ　　　　－ｇｐｆｆ　

１卷第７期　第４赵彬彬等：顾及空间关系约束的不同比例尺面目标不一致性同化处理９１７（）：彭东亮，邓敏，赵ｅｎｓｉｎ０１２，１６５９５３６８（Ｓ２９－ｇ，］彬彬．河网多尺度ＭｏｒｈｉｎＪ．ｐｇ的变换方法研究［（）：）遥感学报，０１２，１６５９５３６８２９－［］Ｗ：Ａ］１８ｏｌｂｅｒＧ．ＩｍａｅＭｏｒｈｉｎＳｕｒｖｅＪ．Ｔｈｅ　　ｇｇｐｇｙ［　［］Ｓ／２０ｔａｔｅＢｕｒｅａｕｏｆＳｕｒｖｅｉｎａｎｄＭａｉｎ．ＧＢＴ　　　　ｙｇｐｐｇ　２０２５７．２—２００６．ＣａｒｔｏｒａｈｉｃＳｍｂｏｌｓｆｏｒＮａｔｉｏｎａｌ　　　ｇｐｙ：ＦｕｎｄａｍｅｎｔａｌＳｃａｌｅＭａｓａｒｔ２ＳｅｃｉｆｉｃａｔｉｏｎｓｆｏｒＰ　　　　－－ｐｐＣａｒｔｏｒａｈｉｃＳｍｂｏｌｓ１∶５０００１∶１００００Ｔｏｏ　　　－ｇｐｙｐ：Ｃ，ｒａｈｉｃＭａｓ［Ｓ］．ＢｅｉｉｎｈｉｎａＳｔａｎｄａｒｄＰｒｅｓｓ　　　ｇｐｐｊｇ（中华人民共和国国家质量监督检验检疫总２００６／局．ＧＢＴ２０２５７．２—２００６．国家基本比例尺地形图　图式第２部分：１∶５０００１∶１００００地形图图式［Ｓ］．　）北京：中国标准出版社，２００６［］Ａ，Ｗ２１ｉＴｉｎｈｕａｕＨｅｈａｉ．ＣｏｎｓｉｓｔｅｎｃＣｏｒｒｅｃｔｉｏｎｏｆ　　　ｇｙ　ＳｈａｒｅｄＢｏｕｎｄａｒＢｅｔｗｅｅｎＡｄａｃｅｎｔＰｏｌｏｎｓ［Ｊ］．　　　ｙｊｙｇ　Ｖｉｓｕａｌｏｍｕｔｅｒ，９９８，１４：３６０７２Ｃ３１　－ｐ［］Ｓ／１９ｔａｔｅＢｕｒｅａｕｏｆＳｕｒｖｅｉｎａｎｄＭａｉｎ．ＧＢＴ　　　　ｙｇｐｐｇ　１２３４３．１—２００８．ＣｏｍｉｌａｔｉｏｎＳｅｃｉｆｉｃａｔｉｏｎｓｆｏｒＮａ　　　－ｐｐ—ＰｉｏｎａｌＦｕｎｄａｍｅｎｔａｌＳｃａｌｅＭａｓａｒｔ１：Ｃｏｍｉｌａｔ　　　　－ｐｐｔｉｏｎＳｅｃｉｆｉｃａｔｉｏｎｆｏｒ１∶２５０００００　　　　０　１∶５　０ｐ：Ｃ１∶１０００００ＴｏｏｒａｈｉｃＭａｓ［Ｓ］．Ｂｅｉｉｎｈｉｎａ　　ｐｇｐｐｊｇ，中华人民共和国国家质量监ＳｔａｎｄａｒｄＰｒｅｓｓ２００８（　／督检验检疫总局．ＧＢＴ１２３４３．１—２００８国家基本　比例尺地图编绘规范第１部分：１∶２５０００　１∶５００００１∶１０００００地形图编绘规范［Ｓ］．北京：　　　）中国标准出版社，００８２ｏｕｒｎａｌｕｈａｎ　ｅｃｈｎｉｃａｌｎｉｖｅｒｓｉｔｕｒｖｅｉｎＪｏＷＴＵｏＳ　　ｙｙｇｆ　ｆ　　，（）：（艾廷华，毋河ａｉｎ０００２５５２６３１４２Ｍ４ａｎｄ　－ｐｐｇ，］海．相邻多边形共享边界的一致化改正［Ｊ．武汉测，（）：）绘科技大学学报，２０００２５５４２６３１４－ＡｎＡｓｓｉｍｉｌａｔｉｏｎＭｅｔｈｏｄｏｆＩｎｃｏｎｓｉｓｔｅｎｃＢｅｔｗｅｅｎＡｒｅａＯｂｅｃｔｓａｔ　　　　　　　ｙｊ　ＤｉｆｆｅｒｅｎｔＳｃａｌｅｓｗｉｔｈＲｅｓｅｃｔｔｏＳａｔｉａｌＲｅｌａｔｉｏｎＣｏｎｓｔｒａｉｎｔｓ　　　　　　　ｐｐ１２ＺＨＡＢＤＳＯ　ｉｎｂｉｎＥＮＧ　ｏｎｌｉａｎｈａｎｓｈａｎ２　Ｐｇｇ　ＺＨＡＮＧ　１１ＬＳＸＱＩＵ　ｈａｎｓｈａｎＩＯＮＧ　ｕｉｎＩ　ｕａｎａ１ｇｇ　Ｘｐｇ　ＤＡｆ，１　ＳｃｈｏｏｌｏｆＴｒａｆｆｉｃａｎｄＴｒａｎｓｏｒｔａｔｉｏｎＥｎｉｎｅｅｒｉｎＣｈａｎｓｈａＵｎｉｖｅｒｓｉｔｏｆＳｃｉｅｎｃｅａｎｄ　　　　　　　ｐｇｇｇｙ　，Ｃｈａｎｓｈａ４１０１１４，ＣｈｉｎａＴｅｃｈｎｏｌｏ　ｇｇｙ，，Ｃｈａｎｓｈａ４１００８３，Ｃｈｉｎａ２　ＤｅａｒｔｍｅｎｔｏｆＧｅｏｎｆｏｒｍａｔｉｃｓＣｅｎｔｒａｌＳｏｕｔｈＵｎｉｖｅｒｓｉｔｉ　　　　　－ｙｇｐ：Ａｕａｌｉｔａｓｓｅｓｓｍｅｎｔ．ＡｓａｎｏｏｓｉｔｅｏｆｃｏｎｓｉｓｂｓｔｒａｃｔＣｏｎｓｉｓｔｅｎｃｉｓａｂａｓｉｃｉｎｄｉｃａｔｏｒｏｆｓａｔｉａｌｄａｔａ－　　　　　　　　　　　ｑｙｐｐｙｐ　　，ｔｅｎｃｉｎｃｏｎｓｉｓｔｅｎｃｏｆｓａｔｉａｌｄａｔａｒｅｆｅｒｓｔｏｔｈｅｃｏｎｆｌｉｃｔｓｏｒｃｏｎｔｒａｄｉｃｔｉｏｎｓｂｅｔｗｅｅｎｓａｔｉａｌｏｂｅｃｔｓａｔ　　　　　　　　　　　　ｙｙｐｐｊ　ｔｈｅｓａｍｅｓｃａｌｅｏｒｄｉｆｆｅｒｅｎｔｓｃａｌｅｓ．Ｓａｔｉａｌｄａｔａｉｎｃｏｎｓｉｓｔｅｎｃｉｓａｃｏｍｍｏｎｃｏｎｃｅｒｎｉｎｔｈｅｉｎｔｅｒｎａｔｉｏｎａｌ　　　　　　　　　　　　　ｐｙ　’，ｅｏｒａｈｉｃａｌｉｎｆｏｒｍａｔｉｏｎｃｏｍｍｕｎｉｔｓｉｎｃｅ１９９０ｓａｎｄｗｉｌｌｈａｖｅａｄｉｒｅｃｔｅｆｆｅｃｔｏｎｓａｔｉａｌｄａｔａｉｎｔｅｒａ－　　　　　　　　　　　　ｇｇｐｙｐｇ　，ｔｅｎｅｒａｌｉｚａｔｉｏｎａｎｄｓａｔｉａｌｄａｔａｕｄａｔｉｎ．Ｐｒｅｖｉｏｕｓｓｔｕｄｉｅｓｍａｉｎｌｆｏｃｕｓｏｎｈａｎｄｌｉｎｉｏｎｃａｒｔｏｒａｈｉｃ　　　　　　　　　ｇｐｐｇｙｇｇｐ　ｒｏｂｌｅｍｏｆｈａｎａｉｄｔｏｔｈｅｉｎｃｏｎｓｉｓｔｅｎｃｉｅｓｂｅｔｗｅｅｎｌｉｎｅｏｂｅｃｔｓｗｉｔｈｔｈｅｓａｍｅｓｃａｌｅ．Ｌｅｓｓａｔｔｅｎｔｉｏｎｉｓ－　　　　　　　　　　　　　　　ｐｐｊｄｌｉｎｉｎｃｏｎｓｉｓｔｅｎｃｉｅｓｂｅｔｗｅｅｎｏｂｅｃｔｓｆｒｏｍ　ｍａｓｗｉｔｈｄｉｆｆｅｒｅｎｔｓｃａｌｅｓｏｒｒｅｌａｔｅｄｉｓｓｕｅｓａｂｏｕｔ“ａｂｕｉｌｄ－　　　　　　　　　　　ｇｊｐ　”ｉｒｏｃｅｓｓｏｆｃａｒｔｏｒａｈｉｃｅｎｅｒａｌｉｚａｔｉｏｎｗｈｉｌｅｄｅｒｉｖｉｎａｓｍａｌｌｅｒｓｃａｌｅｍａｎｄｒｏｅｄｉｎａｒｉｖｅｒｉｎｔｈｅ　　　　　　　　　　　　　ｐｇｐｇｇｐｇｐｐ　　，ｆｒｏｍａｌａｒｅｒｓｃａｌｅｍａｅｎｃｅｔｈｉｓａｅｒｒｏｏｓｅｓａｎａｒｏａｃｈｂａｓｅｄｏｎｄａｔａａｓｓｉｍｉｌａｔｉｏｎｒｉｎｃｉｌｅｓ　　　　　　　　　　　　　ｇｐ．Ｈｐｐｐｐｐｐｐｐｆｏｒｈａｎｄｌｉｎｉｎｃｏｎｓｉｓｔｅｎｃｂｅｔｗｅｅｎｒｉｖｅｒｓａｎｄｂｕｉｌｄｉｎｓ．Ｔｈｅｒｏｏｓｅｄｍｅｔｈｏｄｆｉｒｓｔｅｘｔｒａｃｔｓｂｏｕｎｄａ－　　　　　　　　　ｇｙｇｐｐ　　，），ｒｉｅｓｏｆｂｏｔｈｒｉｖｅｒｓａｔｄｉｆｆｅｒｅｎｔｓｃａｌｅｓｂｍｅａｎｓｏｆｓａｔｉａｌｏｅｒａｔｉｏｎｓ（ｅ．．ｓｌｉｔｉｎｔｅｒｓｅｃｔｉｏｎ．Ｔｈｅｎ　　　　　　　　　　ｙｐｐｇｐ　，ｓｃａｌｅｒｉｖｅｒｂｏｕｎｄａｒｔｏｔｈｅｏｔｈｅｒｏｎｅ．Ｆｉｎａｌｌｍｏｒｈｉｎｔｒａｎｓｆｏｒｍａｔｉｏｎｓａｒｅｅｒｆｏｒｍｅｄｆｒｏｍａｌａｒｅｒ－　　　　　　　　　　ｐｇｙｙｐｇ　　，ａｎｏｔｉｍａｌｗｉｔｈｒｅｓｅｃｔｔｏｔｏｏｌｏｉｃａｌａｎｄｄｉｓｔａｎｃｅｃｏｎｓｔｒａｉｎｔｓｂｅｔｗｅｅｎｒｉｖｅｒｂｏｕｎｄａｒａｎｄａｂｕｉｌｄｉｎ　　　　　　　　　　　　ｐｐｇｙｇｐ　ａｓｓｉｍｉｌａｔｅｄｒｉｖｅｒｂｏｕｎｄａｒｉｓｄｅｔｅｒｍｉｎｅｄ．Ｅｘｅｒｉｍｅｎｔｓｄｅｍｏｎｓｔｒａｔｅｔｈａｔｔｈｅｒｏｏｓｅｄｍｅｔｈｏｄｉｓａｂｌｅ　　　　　　　　　　ｙｐｐｐ　，ｓｔｏｈａｎｄｌｅｉｎｃｏｎｓｉｓｔｅｎｃｂｅｔｗｅｅｎｅｎｅｒａｌｉｚｅｄｒｉｖｅｒａｎｄｂｕｉｌｄｉｎｓａｔａｌａｒｅｒｃａｌｅｅｆｆｅｃｔｉｖｅｌｒｏｖｉｄｉｎ－　　　　　　　　　　ｙｇｇｇｙｐｇ　ｓａｐｒｏｍｉｓｉｎｓｏｌｕｔｉｏｎｆｏｒｅｌｉｍｉｎａｔｉｎｔｏｏｌｏｉｃａｌｃｏｎｆｌｉｃｔｓｗｈｉｌｅｄｅｒｉｖｉｎｃｏｎｓｉｓｔｅｎｔｌａｓｍａｌｌｅｒｃａｌｅ－　　　　　　ｇｇｐｇｇｙ　　　　ｓｍａｆｒｏｍａｌａｒｅｒｃａｌｅｍａ－　　　ｐｇｐ．　：；；；ｍ；Ｋｅｗｏｒｄｓｓｃａｌｅｓａｔｉａｌｒｅｌａｔｉｏｎｉｎｃｏｎｓｉｓｔｅｎｃｏｒｈｉｎｔｒａｎｓｆｏｒｍａｔｉｏｎａｓｓｉｍｉｌａｔｉｏｎ　ｐｙｐｇｙ　　：，，ＦｉｒｓｔａｕｔｈｏｒＺＨＡＯＢｉｎｂｉｎＰｈＤ，ａｓｓｏｃｉａｔｅｒｏｆｅｓｓｏｒｓｅｃｉａｌｉｚｅｓｉｎｔｈｅｔｈｅｏｒｉｅｓａｎｄｍｅｔｈｏｄｓｏｆｉｎｃｏｎｓｉｓｔｅｎｃｄｅｔｅｃｔｉｏｎａｎｄｈａｎｄｌｉｎｏｆ　　　　　　　　　　　　ｐｐｙｇ　　：ｍｕｌｔｉｏｕｒｃｅａｎｄｍｕｌｔｉｃａｌｅｍａｄａｔａ．Ｅ－ｍａｉｌｚｂｂｓｉｒ６３．ｃｏｍｓｓ　　　＠１－－ｐ　：，ＦｏｕｎｄａｔｉｏｎｓｕｏｒｔＴｈｅＮａｔｉｏｎａｌＮａｔｕｒａｌＳｃｉｅｎｃｅＦｏｕｎｄａｔｉｏｎｏｆＣｈｉｎａＮｏ．４１３０１４０４；ｔｈｅＨｕｎａｎＰｒｏｖｉｎｃｉａｌＮａｔｕｒａｌＳｃｉｅｎｃｅＦｏｕｎｄａｔｉｏｎ　　　　　　　　　　　　ｐｐ，Ｎ，ｏｆＣｈｉｎａｏ．１４ＪＪ３０８３；ｔｈｅＳｃｉｅｎｔｉｆｉｃＲｅｓｅａｒｃｈＦｏｕｎｄａｔｉｏｎｆｏｒＮＡＳＧ　ＫｅＬａｂｏｒａｔｏｒｏｆＬａｎｄＥｎｖｉｒｏｎｍｅｎｔａｎｄＤｉｓａｓｔｅｒＭｏｎｉｔｏｒｉｎ　　　　　　　　　　　ｙｙｇ　　Ｎｏ．ＬＥＤＭ２０１２Ｂ０２．