-
摘要: 月球地形制图是了解月球形貌构造，开展月球工程探测及科学研究的基础性工作。自从发射嫦娥一号月球探测卫星以来，中国正式步入月球与行星制图研究的行列。首先，回顾了国内外月球探测与地形制图的基本情况；然后，结合中国嫦娥工程中在月球地形制图领域的研究工作，从月球地形制图标准体系建设、月球地形制图产品研制、月球地形制图理论及技术研究3个方面介绍了中国在月球地形制图领域开展的研究进展和取得的初步成果；最后，总结了当前月球地形制图研究领域中存在的关键理论和技术难题，并探讨了下一步亟需开展的研究工作，可为中国未来的月球地形制图研究提供借鉴。Abstract:Objectives Lunar topographic mapping is one of the fundamental research work for lunar topographic study and lunar exploration program. Since China launched the Chang'E-1 probe in 2007, we have officially joined the group of lunar and planetary mapping research. With the successive success of the following lunar exploration missions, China has successfully obtained a large number of lunar topographic data. Using these data, scholars have carried out preliminary research on lunar topographic mapping, the main purpose of this paper is to review and sort out these research results and understand the research progress of Chang'E lunar topographic mapping in China.Methods Lunar topographic mapping is a subject integrating science, engineering and application. Therefore, in order to understand the progress of China's lunar topographic mapping in the Chang'E project, this paper conducts research and combing from the following three aspects. Firstly, from the aspect of engineering, mainly starting from the standard development of lunar topographic mapping, this paper reviews how domestic scholars carry out the standard development work according to the needs of Chang'E project and related scientific research, starting from the blank, and gradually establish the preliminary research situation of the relevant standard system of lunar topographic mapping, and summarizes the existing standard research results. Secondly, from the aspect of product application, starting from three more commonly used map products, such as lunar atlas product, multimedia map product and lunar globe product, this paper surveys and reviews in detail how domestic scholars design and use reasonable data processing and mapping methods to develop lunar topographic map product according to the quality characteristics of Chang'E lunar data, and summarizes the lunar map products available. Finally, based on the experience gained from Chang'E lunar topographic mapping and its application, this paper reviews and combs the research on the theory and technology of lunar topographic mapping in China, starting from the elements and symbolic expression of lunar topographic mapping, multi-scale topographic mapping, and visualization of lunar topography.Results To meet the needs of Chang'E project and lunar scientific research in China, Chinese scholars have successfully developed 10 national standards related to lunar topographic mapping, and initially established a basic lunar topographic mapping standard system in China. Based on Chang'E lunar data, domestic scholars have compiled five atlas products, developed multimedia lunar map products covering map websites, APP and other forms, developed and produced various lunar globes, and these mapping products have won many awards in the field of mapping at home and abroad. At the same time, domestic scholars have preliminarily studied and established the classification system of lunar topography and geomorphology, and have made some progress in the research of cartography theory and visualization technology of multi-scale cartography for the special lunar topography.Conclusions Lunar and planetary cartography could not be carried out without the need and promotion of the corresponding lunar and deep space exploration programs. Based on the Chang'E project, China has obtained a large number of lunar topographic data. Using these data, China has carried out the research on the theory, technology and application of lunar topographic mapping, and made some progress. The corresponding mapping products have been widely recognized by the international cartographic community, which indicates that in the world, China has entered the forefront of research in the field of lunar topographic mapping. In the future, under the guidance of follow-up lunar exploration program and scientific research, Chinese scholars will continue to carry out further research from the difficult problems in the field of lunar topographic mapping, such as multi-source lunar exploration data fusion, quantitative indicators system construction and automatic extraction of lunar topography features, and automatic multi-scale mapping of massive lunar topographic data.
-
Keywords:
- lunar and planetary cartography /
- topographic map /
- lunar image /
- Chang'E project
-
尺度是空间数据的主要特性之一, 也是多尺度空间数据建模与分析中的重要内容。在地理学、地理信息科学等研究领域, 不同的认知过程, 其尺度含义不同[1-4], 如可划分为空间尺度、时间尺度和语义尺度等[4]。尺度变化可能导致空间数据的维数、几何形态、属性信息以及空间关系的变化[5-6]。多尺度数据处理与分析的一个主要内容是根据尺度变化对空间信息进行概化和综合, 通过对详细尺度的空间数据进行空间综合(选取、合并、化简)和属性综合, 达到自动或者半自动处理粗略尺度下的空间数据的目的。由尺度变化引起的拓扑关系、方向关系等空间关系的变化称为尺度依赖性建模[7], 主要应用于拓扑关系、方向关系的一致性分析[8-11]、多尺度空间数据查询[12-13]等领域。语义尺度影响下的拓扑关系多尺度计算是空间关系尺度依赖性建模问题之一[7, 14]。Tryfona等[8]、杜世宏等[7, 10]指出该问题的语义尺度主要表现为多尺度属性划分, 通过属性归纳从详细尺度向粗略尺度转换; 尺度变化后的空间对象通过区域合并进行概化综合, 此过程不改变区域对象的维数和几何形态, 空间关系由推理得到。
目前, 拓扑关系的多尺度计算主要依赖组合推理的方法, 通过简单对象间的基本拓扑关系(简单区域间的8种基本关系, 线面间的19种关系等), 构建出推理组合表进行查表运算, 比如区域合并组合推理[8]、简单区域间的组合推理[15]等。这种方法简单直观, 易于理解, 但其不足主要表现为:①由于基于空间对象间的基本关系, 其推理结果的值域仍是基本关系; ②只能适用于简单空间对象的处理, 复杂空间对象间的推理组合不易实现; ③推理结果存在多解性。
在拓扑关系的描述模型中, 9交模型[16-17]仍是广泛使用的一种方法, 被OpenGIS采用[18]。9交模型通过两个对象间内部(°)、边界(∂)和外部(+)的交集形成一个9交矩阵描述拓扑关系, 已被应用于描述复杂对象的拓扑关系, 如复杂区域对象间存在33种拓扑关系[19-20]。基于9交模型可以提供多种层次的空间查询, 包括谓词查询和9交矩阵查询[18, 21]。如果直接基于9交矩阵进行拓扑关系的多尺度计算, 可弥补现有方法的不足, 一方面可以不限于简单对象, 另一方面计算结果的值域可扩展为9交模型表达的所有拓扑关系。针对语义尺度影响下的拓扑关系多尺度计算, 本文提出了基于9交矩阵的计算方法, 通过定义9交关系矩阵操作算子, 直接计算出9交矩阵。
1 语义尺度与空间区域划分
区域对象在OpenGIS[18]和文献[19, 22]中进行了详细描述。设区域对象A⊂R2, 满足$A = \overline {A^\circ }, \overline {A^\circ } $是包含A°的最小正则闭包; -A是A的补集, 即-A= R2-A°。R2=A°∪ ∂A ∪ A+, -A=∂A ∪ A+, -(-A) = A。简单区域是一种特殊的区域对象, 同胚于一个圆盘, 具有完全连通的内部、边界和外部。
对于R2上的两个简单区域A1和A2, 如图 1所示, 将两个相离区域的并A1∪ A2记为A1+A2, 其合并结果是一个复杂区域, 称为组合区域; 两个相邻区域(存在公共边, 拓扑关系为边相接)的并A1∪ A2记为A1⊕A2, 其合并结果是一个去除公共边内部的简单区域。
语义尺度主要表现为多尺度属性划分[7], 比如:①概念划分, 即部分和整体的关系, 如行政区划, 村庄是县的一部分, 县是市的一部分; ②多尺度分类系统。属性是表征空间对象可被人类认知理解的概念化语义、符号或者名称。语义尺度的变化引起空间区域的多尺度划分, 表现为不同尺度下平面空间被划分成有限个区域对象的集合[7]。
定义1:语义尺度属性集S是一组属性的集合, 记为S ={ s1…sm} (m>1)。语义尺度属性集S2比语义尺度属性集S1详细, 则存在映射f:S2→ S1, f表示部分和整体的关系或者分类关系映射。
定义2:A是语义尺度S1下的一个区域对象, P ={ p1…pn} (n>1)为一个简单区域集合, 称P为A在语义尺度S2下的一个区域划分, 当且仅当下列条件满足:
1) 语义尺度S2比语义尺度S1详细。
2) p1…pn均为简单区域。
3) ∀ i, j = 1…n且i ≠ j, pi和pj的拓扑关系只能是边相接和相离情况之一。
4) $A = \bigcup\nolimits_{i = 1}^n {{p_i}} $。
5) pi属于语义尺度S2下的某一个属性。
6) 同一属性的区域间只能是相离关系。
以上区域划分是一个无缝划分, 不允许由粗到细的划分过程中出现空白区域。边相接关系说明相邻区域须存在公共边。区域合并是区域划分的逆过程, 即从一个详细的语义尺度S2向粗略语义尺度S1转换的过程。
2 基于9交模型的拓扑关系计算
2.1 9交模型
9交模型[16-17]是通过两个区域对象A、B的内部(°)、边界(∂)和外部(+)的交集来描述拓扑关系。由一个3×3的0/1型9交矩阵R(A, B)表示:
$$ \begin{array}{*{20}{c}} {\left[ {\begin{array}{*{20}{c}} {r_{11}^{AB}}&{r_{12}^{AB}}&{r_{13}^{AB}}\\ {r_{21}^{AB}}&{r_{22}^{AB}}&{r_{23}^{AB}}\\ {r_{31}^{AB}}&{r_{32}^{AB}}&{r_{33}^{AB}} \end{array}} \right] = }\\ {\left[ {\begin{array}{*{20}{c}} {{A^ \circ } \cap {B^ \circ }}&{{A^ \circ } \cap \partial B}&{{A^ \circ } \cap {B^ + }}\\ {\partial A \cap {B^ \circ }}&{\partial A \cap \partial B}&{\partial A \cap {B^ + }}\\ {{A^ + } \cap {B^ \circ }}&{{A^ + } \cap \partial B}&{{A^ + } \cap {B^ + }} \end{array}} \right]} \end{array} $$ (1) R(A, B)所能表达的拓扑关系总数为33种[19]。若A和B均为简单区域, 仅存在8种关系:相离、相接、相交、包含、覆盖、包含于、覆盖于和相等。
2.2 相离区域合并与拓扑关系尺度计算
如果一个粗尺度下的区域对象B由两个相离的简单区域B1、B2构成, 即B=B1+B2, 则简单区域A和B的9交拓扑关系可通过A和B1、B2的拓扑关系分解计算得到[23]:
$$ \begin{array}{*{20}{c}} {\mathit{\boldsymbol{R}}\left( {A,B} \right) = \mathit{\boldsymbol{R}}\left( {A,{B_1} + {B_2}} \right) = }\\ {\mathit{\boldsymbol{R}}\left( {A,{B_1}} \right)\underline \vee \mathit{\boldsymbol{R}}\left( {A,{B_2}} \right) = }\\ {\left[ {\begin{array}{*{20}{c}} {r_{11}^{A{B_1}} \vee r_{11}^{A{B_2}}}&{r_{12}^{A{B_1}} \vee r_{12}^{A{B_2}}}&{r_{13}^{A{B_1}} \wedge r_{13}^{A{B_2}}}\\ {r_{21}^{A{B_1}} \vee r_{21}^{A{B_2}}}&{r_{22}^{A{B_1}} \vee r_{22}^{A{B_2}}}&{r_{23}^{A{B_1}} \wedge r_{23}^{A{B_2}}}\\ {r_{31}^{A{B_1}} \vee r_{31}^{A{B_2}}}&{r_{32}^{A{B_1}} \vee r_{32}^{A{B_2}}}&{r_{33}^{A{B_1}} \wedge r_{33}^{A{B_2}}} \end{array}} \right]} \end{array} $$ (2) 式中, 运算符⊻称为相离区域间拓扑关系加算子。⊻满足交换律R(A, B1)⊻ R(A, B2)= R(A, B2)⊻ R(A, B1), 并且, R (B1+B2, A)=[R(A, B1+B2)]T=[R(A, B1)⊻ R(A, B2)]T, T表示矩阵转置。
2.3 相邻区域合并与拓扑关系尺度计算
如果一个粗尺度下的区域对象B由两个边相接的简单区域B1、B2构成, 记为B=B1⊕B2, 此时利用运算符⊻不能计算出正确的拓扑关系。文献[8]通过推理组合表定性处理, 将计算结果限定在某几个关系中, 再根据实际交集情况逐一排查。如图 2所示, A与B1和B2的拓扑关系均为相交, 利用式(1)计算出的结果是相交, 显然与图示A包含于B不相符。利用文献[8]中的推理表, 图 2中R(A, B)可以是相交、覆盖于和包含于3种取值情况之一, 这导致了推理的多解性。
![]()
图 2 具有公共边界的区域合并与拓扑关系计算Figure 2. Regions Merging with Common Boundaries and Topological Relations Computing由于9交矩阵以及8个基本拓扑关系是一种定性方法, 这类拓扑关系的描述模型通过0和1概化处理交集关系时, 已经丢失了大量的信息, 推理表只能给出取值的多种可能性。要消除这一多解问题, 显然需要补充信息, 弥补概化造成的信息损失。
通过分析图 2拓扑关系计算失败的原因, 可以发现:①区域对象B1(或B2)的外部与B2(或B1)的内部相交, 通过9交矩阵R(A, B1)无法区分哪些部分是B2的内部; ②区域对象B1和B2的公共边界属于B的内部, 从R(A, B1)和R(A, B2)无法获知A是否与公共边界相交。因此, 相邻区域合并后, 从9交矩阵无法有效判断出A与B的外部和边界的真实交集信息。
因此增加辅助信息, 使得拓扑关系可直接计算。当区域对象B由两个边相接的简单区域B1、B2构成时, 增加辅助区域对象b1…bn (n>1), 用以区分B的外部。添加规则如下:
1) b1…bn均为简单区域。
2) 在{ b1…bn, B1, B2}中, 其任意两个区域的关系只能是边相接和相离情况之一。
3) 任意bi与B1或者B2存在边相接关系, 其中i=1…n。
4) B的边界是辅助区域对象边界的一部分, 即$\left( {\partial {B_{1}} \cup \partial {B_2}} \right)-\left( {\partial {B_{1}} \cap \partial {B_2}} \right) \subset \bigcup\nolimits_{i = 1}^n {\partial {b_i}} $。
5) B1和B2公共边界的两个端点位于辅助区域对象边界, 即$\partial {B_{1}} \cap \partial {B_{2}} \subset \bigcup\nolimits_{i = 1}^n {\partial {b_i}} $。
如图 3所示, 增加辅助区域对象, 图 3(b)和图 3(c)为合理的添加规则; 而图 3(d)为不合理的添加规则, 因为b1与b2不为边相接关系。图 4为区域合并后形成洞的情况, 洞内也需要增加辅助区域。
![]()
图 3 添加辅助区域消除计算错误(当B1和B2为边相接关系)Figure 3. Extra Regions Added to Solve Ambiguities When B1 Meets B2 with Common Boundaries![]()
图 4 添加带洞辅助区域处理具有公共边界的B1和B2区域合并Figure 4. Adding Extra Regions with Holes to Merge B1 and B2 with Common Boundaries当添加辅助区域后, 区域对象B的内部、边界和外部将分别由式(3)至式(5)计算得到:
$$ {B^ \circ } = \left( {B_1^ \circ \cup B_2^ \circ } \right) \cap \left( {\bigcap\nolimits_{i = 1}^n {b_i^ + } } \right) $$ (3) $$ \partial B = \left( {\partial {B_1} \cup \partial {B_2}} \right) \cap \left( {\bigcup\nolimits_{i = 1}^n {\partial {b_i}} } \right) $$ (4) $$ {B^ + } = \left( {B_1^ + \cup B_2^ + } \right) \cap \left( {\bigcup\nolimits_{i = 1}^n {b_i^ \circ } } \right) $$ (5) 当区域对象B由两个相邻的简单区域B1、B2构成时, 则简单区域A和B的9交拓扑关系可通过A和B1、B2以及辅助区域对象b1…bn (n>1)间的拓扑关系计算得到:
$$ \begin{array}{*{20}{c}} {\mathit{\boldsymbol{R}}\left( {A,B} \right) = \mathit{\boldsymbol{R}}\left( {A,{B_1} + {B_2}} \right) \wedge }\\ {\left[ { \bot \mathit{\boldsymbol{R}}\left( {A,{b_1} + {b_2} + \cdots + {b_n}} \right)} \right]} \end{array} $$ (6) 式中, ⊥表示互换9交矩阵的第1列和第3列; ∧表示矩阵对应元素进行逻辑与运算。
证明:设Ma={A°, ∂A, A+}, ∀a ∈ Ma, 可得:
$$ \begin{array}{l} a \cap {B^ \circ } = a \cap \left( {\left( {B_1^ \circ \cup B_2^ \circ } \right) \cap \left( {\bigcap\nolimits_{i = 1}^n {b_i^ + } } \right)} \right) = \\ \left[ {\left( {a \cap B_1^ \circ } \right) \cup \left( {a \cap B_2^ \circ } \right)} \right] \cap \left[ {a \cap \left( {\bigcap\nolimits_{i = 1}^n {b_i^ + } } \right)} \right]\\ \;\;\;\;\;\;\;\;\;\; \to \left( {r_{i1}^{A{B_1}} \vee r_{i1}^{A{B_2}}} \right) \wedge \left( {\mathop \wedge \limits_{i = 1}^n r_{i3}^{A{b_i}}} \right) \end{array} $$ $$ \begin{array}{l} a \cap \partial B = a \cap \left( {\left( {\partial {B_1} \cup \partial {B_2}} \right) \cap \left( {\bigcup\nolimits_{i = 1}^n {\partial {b_i}} } \right)} \right) = \\ \left[ {\left( {a \cap \partial {B_1}} \right) \cup \left( {a \cap \partial {B_2}} \right)} \right] \cap \left[ {a \cap \left( {\bigcup\nolimits_{i = 1}^n {\partial {b_i}} } \right)} \right]\\ \;\;\;\;\;\;\;\;\;\; \to \left( {r_{i2}^{A{B_1}} \vee r_{i2}^{A{B_2}}} \right) \wedge \left( {\mathop \vee \limits_{i = 1}^n r_{i2}^{A{b_i}}} \right) \end{array} $$ $$ \begin{array}{l} a \cap {B^ + } = a \cap \left( {\left( {B_1^ + \cap B_2^ + } \right) \cap \left( {\bigcup\nolimits_{i = 1}^n {b_i^ \circ } } \right)} \right) = \\ \left[ {\left( {a \cap B_1^ + } \right) \cap \left( {a \cap B_2^ + } \right)} \right] \cap \left[ {a \cap \left( {\bigcup\nolimits_{i = 1}^n {b_i^ \circ } } \right)} \right]\\ \;\;\;\;\;\;\;\;\;\; \to \left( {r_{i3}^{A{B_1}} \vee r_{i3}^{A{B_2}}} \right) \wedge \left( {\mathop \vee \limits_{i = 1}^n r_{i1}^{A{b_i}}} \right) \end{array} $$ 式中, $\mathop \wedge \limits_{i = 1}^n $表示n-1个逻辑与运算; $\mathop \vee \limits_{i = 1}^n $表示n-1个逻辑或运算。由于A、B1和B2以及辅助区域对象b1…bn (n>1)均为简单区域, A与B1、B2间的计算可以记为R(A, B1+ B2), 而A与辅助区域对象b1…bn的计算可以记为R(A, b1 + b2 +…+ bn)。R(A, b1 + … + bn)的第1列与第3列调换, 即得R(A, B)= R(A, B1+ B2) ∧ [⊥R(A, b1 + b2 +…+ bn)]。
3 拓扑关系计算应用实例
图 5所示的气象信息系统空气污染条件中, 图 5(a)为华北4省市(北京市、天津市、河北省和山西省)气象信息系统某一时段的空气污染状况。其中一项查询服务为查询行政区域组合与污染条件分布区的拓扑关系。气象信息系统中, 最详细的尺度为县区级行政区(见图 5(b)), 系统只存储县区级区域和空气污染分布区域间的拓扑关系。县区级行政区采用简单区域对象进行表达, 省、市或者特定地区的区域对象由县区级行政区合并而成, 系统内仅存储了其对县区级行政区域的构成关系。由于存在相邻区域合并情况, 故在整个地图(华北4省市)外部增加了两个附加的县区级区域, 如图 6所示。
![]()
图 5 华北4省气象信息系统空气污染条件(2016年10月5日6时)Figure 5. Meteorological Conditions of Air Pollution in Four Provinces of North China (At 6:00 on October 5, 2016)![]()
图 6 河北省区域构成与附加区域Figure 6. Administrative Regions of Hebei Province and the Added Extra Regions记A为某一时段某一等级空气污染条件分布区域(简单区域), 多个县区级行政区域合并与A之间的拓扑关系计算步骤为:
1) 获取查询区域内所有县级区域对象的集合MB={B1 …Bm} (m>1)。
2) 将MB中的区域对象根据相邻和相离关系拆分成多个相离部分Ci, 即MB= ∪Ci。其中, Ci满足:Ci中任意区域对象与非Ci中其他区域对象是相离关系, 且Ci中任意区域对象必与Ci中另一区域对象是相接关系, 否则Ci中只有一个区域对象。
3) 对任意Ci, 存在与Ci中的对象为相接关系但不属于Ci的区域对象, 构成辅助区域对象集合Mb={b1…bn} (n>1)。
4) 针对任意Ci, 按式(6)计算A与Ci的拓扑关系。
5) 按式(2)计算A与B的拓扑关系。
查询实例:计算2016年10月5日6时河北省与五级空气污染分布区的拓扑关系。步骤如下:
1) 记简单区域A为五级空气污染分布区。
2) 记区域B为河北省区域, 其由兴隆县等172个县级区域构成, 即:
MB={B1…B172}={兴隆县…宽城}
从图 6(a)可知, 河北由两个相离部分构成, 划分为两个区域, 即区域1(C1)和区域2(C2)。
3) C1由169个县级区域构成, 即C1={B1…B169}={兴隆县…}。记Mb= {b1…b32, b33}={延庆…附加区域1, 附加区域2}, 按式(6)计算A与C1的拓扑关系为:
$$ \begin{array}{*{20}{c}} {\mathit{\boldsymbol{R}}\left( {A,{C_1}} \right) = \mathit{\boldsymbol{R}}\left( {A,{B_1} + \cdots + {B_{169}}} \right) \wedge }\\ {\left[ { \bot \mathit{\boldsymbol{R}}\left( {A,{b_1} + \cdots + {b_{33}}} \right)} \right] = \left[ {\begin{array}{*{20}{c}} 1&0&0\\ 1&0&0\\ 1&1&1 \end{array}} \right]} \end{array} $$ 4) C2由3个县级区域构成, 即C2={B1, B2, B3}={三河, 大厂, 香河}。记Mb={b1…b6}={通州, 顺义, 平谷, 武清, 蓟县, 宝坻}, 计算A与C2的拓扑关系为:
$$ \begin{array}{*{20}{c}} {\mathit{\boldsymbol{R}}\left( {A,{C_2}} \right) = \mathit{\boldsymbol{R}}\left( {A,{B_1} + {B_2} + {B_3}} \right) \wedge }\\ {\left[ { \bot \mathit{\boldsymbol{R}}\left( {A,{b_1} + \cdots + {b_6}} \right)} \right] = \left[ {\begin{array}{*{20}{c}} 0&0&1\\ 0&0&1\\ 1&1&1 \end{array}} \right]} \end{array} $$ 5) 按式(2)计算A与B的拓扑关系为:
$$ \begin{array}{*{20}{c}} {\mathit{\boldsymbol{R}}\left( {A,B} \right) = \mathit{\boldsymbol{R}}\left( {A,{C_1} + {C_2}} \right) = }\\ {\mathit{\boldsymbol{R}}\left( {A,{C_1}} \right)\underline \vee \mathit{\boldsymbol{R}}\left( {A,{C_2}} \right) = \left[ {\begin{array}{*{20}{c}} 1&0&0\\ 1&0&0\\ 1&1&1 \end{array}} \right]} \end{array} $$ 通过计算可得9交矩阵R(A, B)的拓扑语义为“包含于”。
4 结语
空间数据间的拓扑关系往往随着尺度的变化而变化, 快速有效地计算拓扑关系有助于多尺度空间数据查询以及拓扑一致性分析。针对区域合并导致的拓扑关系变化与多尺度计算, 区别于目前组合推理的方法[8-10], 本文提出了直接基于9交矩阵的拓扑关系多尺度计算方法, 并进行了应用分析。不同语义尺度下的复杂区域由有限个简单区域合并而成, 包括相离区域合并和相邻区域合并, 通过定义两个9交矩阵操作算子, 可利用详细语义尺度下的拓扑关系计算出粗略语义尺度下的拓扑关系。该方法的特点为:计算结果值域为复杂区域间所有可能的拓扑关系, 可以适用于不同层次的查询分析; 主要针对数据变化不大的情况, 通过增加辅助信息消除了相邻区域合并引起的拓扑关系计算多解性; 目前只考虑了多个相离或相邻的简单区域合并情况。
由于布尔型9交矩阵对真实拓扑关系进行了简化, 需进一步研究采取增加信息的办法, 扩展本方法适用于尺度变化引起的拓扑关系计算问题, 并结合栅格领域、方向关系消除计算歧义性; 同时后续研究还可考虑将该方法扩展到带洞区域或者其他类型的空间对象(比如点、线)的多尺度分析中等。
-
[1] Schenk P M, Bussey D B J. Galileo Stereo Topography of the Lunar North Polar Region[J]. Geophysical Research Letters, 2004, 31(23): L23701
[2] Heiken G H, Vaniman D T, French B M. Lunar Sourcebook: A User's Guide to the Moon[M]. Cambridge : Cambridge University Press, 1991
[3] Bowker D E, Hughes J K. Lunar Orbiter Photographic Atlas of the Moon[M]. Washington D C: Government Printing Office, 1970
[4] 欧阳自远, 邹永廖, 李春来, 等. 月球——人类走进深空的前哨站[M]. 北京: 科学出版社, 2002 Ouyang Ziyuan, Zou Yongliao, Li Chunlai, et al. The Moon—An Outpost for Human Beings to Enter Deep Space[M]. Beijing: Science Press, 2002
[5] 欧阳自远. 月球科学概论[M]. 北京: 中国宇航出版社, 2005 Ouyang Ziyuan. Introduction to Lunar Science[M]. Beijing: China Astronautic Publishing House, 2005
[6] Binder A B. Lunar Prospector: Overview[J]. Science, 1998, 281(5382): 1475-1476 doi: 10.1126/science.281.5382.1475
[7] Foing B H, Racca G D, Marini A, et al. SMART-1 Mission to the Moon: Technology and Science Goals[J]. Advances in Space Research, 2003, 31(11): 2323-2333 doi: 10.1016/S0273-1177(03)00541-6
[8] Araki H, Tazawa S, Noda H, et al. Lunar Global Shape and Polar Topography Derived from KaguyaLALT Laser Altimetry[J]. Science, 2009, 323(5916): 897-900 doi: 10.1126/science.1164146
[9] 李春来, 刘建军, 任鑫, 等. 嫦娥一号图像数据处理与全月球影像制图[J]. 中国科学: 地球科学, 2010, 40(3): 294-306 doi: 10.3969/j.issn.0253-2778.2010.03.015 Li Chunlai, Liu Jianjun, Ren Xin, et al. The Global Image of the Moon by the Chang'E-1: Data Processing and Lunar Cartography[J]. Scientia Sinica Terrae, 2010, 40(3): 294-306 doi: 10.3969/j.issn.0253-2778.2010.03.015
[10] 李春来, 任鑫, 刘建军, 等. 嫦娥一号激光测距数据及全月球DEM模型[J]. 中国科学: 地球科学, 2010, 40(3): 281-293 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201003002.htm Li Chunlai, Ren Xin, Liu Jianjun, et al. Laser Altimetry Data of Chang'E-1 and the Global Lunar DEM Model[J]. Scientia Sinica Terrae, 2010, 40(3): 281-293 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201003002.htm
[11] 李春来, 刘建军, 任鑫, 等. 基于嫦娥二号立体影像的全月高精度地形重建[J]. 武汉大学学报·信息科学版, 2018, 43(4): 485-495 doi: 10.13203/j.whugis20170400 Li Chunlai, Liu Jianjun, Ren Xin, et al. Lunar Global High-Precision Terrain Reconstruction Based on Chang'E-2 Stereo Images[J]. Geomatics and Information Science of Wuhan University, 2018, 43(4): 485-495 doi: 10.13203/j.whugis20170400
[12] 欧阳自远. 嫦娥二号的初步成果[J]. 自然杂志, 2013, 35(6): 391-395 https://www.cnki.com.cn/Article/CJFDTOTAL-ZRZZ201306002.htm Ouyang Ziyuan. Chang'E-2 Preliminary Results[J]. Chinese Journal of Nature, 2013, 35(6): 391-395 https://www.cnki.com.cn/Article/CJFDTOTAL-ZRZZ201306002.htm
[13] 欧阳自远, 李春来, 邹永廖, 等. 绕月探测工程的初步科学成果[J]. 中国科学: 地球科学, 2010, 40(3): 261-280 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201003001.htm Ouyang Ziyuan, Li Chunlai, Zou Yongliao, et al. The Primary Science Result from the Chang'E-1 Probe[J]. Scientia Sinica Terrae, 2010, 40(3): 261-280 https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201003001.htm
[14] Li C L, Mu L L, Zou X D, et al. Analysis of the Geomorphology Surrounding the Chang'E-3 Landing Site[J]. Research in Astronomy and Astrophysics, 2014, 14(12): 1514-1529 doi: 10.1088/1674-4527/14/12/002
[15] Liu J J, Ren X, Yan W, et al. Descent Trajectory Reconstruction and Landing Site Positioning of Chang'E-4 on the Lunar Farside[J]. Nature Communications, 2019, 10: 4229 doi: 10.1038/s41467-019-12278-3
[16] Goswami J N, Annadurai M. Chandrayaan-1 Mission to the Moon[J]. Acta Astronautica, 2008, 63(11): 1215-1220
[17] Sinha R K, Sivaprahasam V, Bhatt M, et al. Geological Characterization of Chandrayaan-2 Landing Site in the Southern High Latitudes of the Moon[J]. Icarus, 2020, 337: 113449 doi: 10.1016/j.icarus.2019.113449
[18] Vondrak R, Keller J, Chin G, et al. Lunar Reconnaissance Orbiter(LRO): Observations for Lunar Exploration and Science[J]. Space Science Reviews, 2010, 150(1): 7-22
[19] Konopliv A S, Park R S, Yuan D N, et al. HighResolution Lunar Gravity Fields from the GRAIL Primary and Extended Missions[J]. Geophysical Research Letters, 2014, 41(5): 1452-1458 doi: 10.1002/2013GL059066
[20] Isbell C, Gaddis L, Garcia P, et al. Kaguya Terrain Camera Mosaics[C]//Lunar and Planetary Science Conference, Woodlands, Texas, USA, 2014
[21] 刘峰, 燕琴, 高武俊. 月球探测的发展[J]. 测绘科学, 2004, 29(4): 63-66 doi: 10.3771/j.issn.1009-2307.2004.04.021 Liu Feng, Yan Qin, Gao Wujun. The Development of Lunar Exploration[J]. Science of Surveying and Mapping, 2004, 29(4): 63-66 doi: 10.3771/j.issn.1009-2307.2004.04.021
[22] 陈俊勇. 月球地形测绘和月球大地测量(1)[J]. 测绘科学, 2004, 29(2): 1-5 doi: 10.3771/j.issn.1009-2307.2004.02.001 Chen Junyong. Topography Surveying and Geodesy on the Moon[J]. Science of Surveying and Mapping, 2004, 29(2): 1-5 doi: 10.3771/j.issn.1009-2307.2004.02.001
[23] 陈俊勇, 章传银, 党亚民. 月球航天探测和月球测绘[J]. 测绘学报, 2005, 34(3): 189-195 doi: 10.3321/j.issn:1001-1595.2005.03.001 Chen Junyong, Zhang Chuanyin, Dang Yamin. Lunar Spatial Exploration, Surveying and Mapping[J]. Acta Geodaetica et Cartographic Sinica, 2005, 34(3): 189-195 doi: 10.3321/j.issn:1001-1595.2005.03.001
[24] 葛之江, 殷礼明, 林宗坚. 月球测图技术发展及其关键技术分析[J]. 航天器环境工程, 2006, 23(2): 67-74 doi: 10.3969/j.issn.1673-1379.2006.02.002 Ge Zhijiang, Yin Liming, Lin Zongjian. Development of Lunar Topography Surveying and Its Key Techniques[J]. Spacecraft Environment Engineering, 2006, 23(2): 67-74 doi: 10.3969/j.issn.1673-1379.2006.02.002
[25] Li C L, Liu J J, Ren X, et al. The Chang'E 3 Mission Overview[J]. Space Science Reviews, 2015, 190(1): 85-101
[26] Bowker D E, Hughes J K. Lunar Orbiter Photographic Atlas of the Moon[M]. Houston: Lunar and Planetary Institute, 1971
[27] Tiernan M, Roth L, Thompson T W, et al. Lunar Cartography with the Apollo 17 Alse Radar Imagery[J]. The Moon, 1976, 15(1): 155-163
[28] Kirk R L, Archinal B A, Gaddis L R, et al. Lunar Cartography: Progress in the 2000s and Prospects for the 2010s[J]. The International Archives of the Photogrammetry, Remote Sensing and Spatial Infor- mation Sciences, 2012, DOI: 10.5194/isprsarchives-XXXIX-B4-489-2012
[29] Bussey B, Spudis P. The Clementine Atlas of the Moon[M]. Cambridge: Cambridge University Press, 2004
[30] Araki H, Tazawa S, Noda H, et al. Preliminary Results of Lunar Topography by KAGUYA-LALT Mission[J]. Transactions of the Japan Society for Aeronautical and Space Sciences, Space Technology Japan, 2009, 7(26): 23-25
[31] Jeffrey J G, Rueb D, Cohen J, et al The Lunar Orbiter Photographic Atlas Digital Archive[C]// Lunar and Planetary Science XXXI, Lunar and Planetary Institute, Houston, Texas, USA, 2000
[32] Byrne C J. Lunar Orbiter Photographic Atlas of the near Side of the Moon[M]. London: Springer London, 2005
[33] Archinal B A, Rosiek M R, Kirk R L, et al. A Report on the Unified Lunar Control Network 2005 and Lunar Topographic Model[C]// European Planetary Science Congress 2006, Berlin, Germany, 2006
[34] 孔祥元, 张松林, 魏二虎, 等. "嫦娥"工程与月球测绘学[J]. 大地测量与地球动力学, 2004, 24(1): 117-120 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB200401027.htm Kong Xiangyuan, Zhang Songlin, Wei Erhu, et al. The Goddess of the Moon' Project with Lunar Geodesy(Selenodesy)[J]. Journal of Geodesy and Geodynamics, 2004, 24(1): 117-120 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB200401027.htm
[35] 全国空间科学及其应用标准化技术委员会. GB/T 30114. 4-2014空间科学及其应用术语第4部分: 月球与行星科学[S]. 北京: 中国标准出版社, 2014 SAC/TC312. GB/T 30114. 4-2014 Terminology for Space Science and Application-Part 4: Lunar and Planetary Science[S]. Beijing: Standards Press of China, 2014
[36] 全国空间科学及其应用标准化技术委员会. GB/T 37541-2019月球地形图要素分类、代码与图式[S]. 北京: 中国标准出版社, 2019 SAC/TC312. GB/T 37541-2019 Classification, Codes and Symbols for the Features of Lunar Topographic Map[S]. Beijing: Standards Press of China, 2019
[37] 全国空间科学及其应用标准化技术委员会. GB/T 30112-2013月球空间坐标系[S]. 北京: 中国标准出版社, 2013 SAC/TC312. GB/T 30112-2013 The Lunar Coordinate System[S]. Beijing: Standards Press of China, 2013
[38] 全国空间科学及其应用标准化技术委员会. GB/T 32521-2016月球基本比例尺分幅和编号[S]. 北京: 中国标准出版社, 2016 SAC/TC312. GB/T 32521-2016 Subdivision and Numbering for the Lunar Primary Scale Topographic Maps[S]. Beijing: Standards Press of China, 2016
[39] 全国空间科学及其应用标准化技术委员会. GB/T 33990-2017月球数字地图质量检查与验收[S]. 北京: 中国标准出版社, 2017 SAC/TC312. GB/T 33990-2017 Inspection and Acceptance for Lunar Digital Map Quality[S]. Beijing: Standards Press of China, 2017
[40] 全国空间科学及其应用标准化技术委员会. GB/T 33991-2017月球信息要素数据字典[S]. 北京: 中国标准出版社, 2017 SAC/TC312. GB/T 33991-2017 Data Dictionary for Lunar Information Essentials[S]. Beijing: Standards Press of China, 2017
[41] 全国空间科学及其应用标准化技术委员会. GB/T 33996-2017月球探测数据产品分级与命名[S]. 北京: 中国标准出版社, 2017 SAC/TC312. GB/T 33996-2017 Hierarchy and Nomenclature of Lunar Exploration Data Products[S]. Beijing: Standards Press of China, 2017
[42] 全国空间科学及其应用标准化技术委员会. GB/T 33997-2017月球与行星数据产品格式规范[S]. 北京: 中国标准出版社, 2017 SAC/TC312. GB/T 33997-2017 Specification of Lunar and Planetary Data Products Format[S]. Beijing: Standards Press of China, 2017
[43] 全国空间科学及其应用标准化技术委员会. GB/T 33998-2017月球信息元数据[S]. 北京: 中国标准出版社, 2017 SAC/TC312. GB/T 33998-2017 Lunar Information Metadata[S]. Beijing: Standards Press of China, 2017
[44] 全国空间科学及其应用标准化技术委员会. GB/T 34054-2017月球影像平面图制作规范[S]. 北京: 中国标准出版社, 2017 SAC/TC312. GB/T 34054-2017 Specifications for Lunar Image Plan Making[S]. Beijing: Standards Press of China, 2017
[45] 全国空间科学及其应用标准化技术委员会. GB/T 34055-2017月球数字高程模型数据制作规范[S]. 北京: 中国标准出版社, 2017 SAC/TC312. GB/T 34055-2017 Production Specifications of Lunar Digital Elevation Model[S]. Beijing: Standards Press of China, 2017
[46] 李春来, 刘建军, 牟伶俐, 等, 嫦娥一号全月球影像图集[M]. 北京: 中国地图出版社, 2010 Li Chunlai, Liu Jianjun, Mu Linli, et al. The Chang'e-1 Image Atlas of the Moon[M]. Beijing: Sinomap Press, 2010
[47] 李春来, 刘建军, 牟伶俐, 等, 嫦娥一号全月球地形图集[M]. 北京: 中国地图出版社, 2012 Li Chunlai, Liu Jianjun, Mu Linli, et al. The Chang'E-1 Topographic Atlas of the Moon[M]. Beijing: Sinomap Press, 2012
[48] 李春来, 刘建军, 牟伶俐, 等. 嫦娥二号高分辨率月球影像图集[M]. 北京: 测绘出版社, 2012 Li Chunlai, Liu Jianjun, Mu Linli, et al. The Chang'E-2 High Resolution Image Atlas of the Moon[M]. Beijing: Surveying and Mapping Press, 2012
[49] 李春来, 刘建军, 牟伶俐, 等. 嫦娥二号高分辨率月球虹湾影像图集[M]. 北京: 中国地图出版社, 2014 Li Chunlai, Liu Jianjun, Mu Linli, et al. The Chang'E-2 High Resolution Image Atlas of Lunar Sinus Iridum[M]. Beijing: Sinomap Press, 2014
[50] 李春来, 刘建军, 任鑫, 等. 嫦娥三号着陆区地形地貌[M]. 北京: 测绘出版社, 2018 Li Chunlai, Liu Jianjun, Ren Xin, et al. Topography of Chang'E-3 Landing Site[M]. Beijing: Surveying and Mapping Press, 2018
[51] 曾兴国, 左维, 牟伶俐, 等. 面向移动端月球形貌展示的快速三维地图平台搭建[J]. 测绘通报, 2017(3): 42-45 https://www.cnki.com.cn/Article/CJFDTOTAL-CHTB201703011.htm Zeng Xingguo, Zuo Wei, Mu Lingli, et al. Research on Fast Development of 3D Lunar Topography Representation Platform Based on Mobile Devices[J]. Bulletin of Surveying and Mapping, 2017(3): 42-45 https://www.cnki.com.cn/Article/CJFDTOTAL-CHTB201703011.htm
[52] 王任享. 月球卫星三线阵CCD影像EFP光束法空中三角测量[J]. 测绘科学, 2008, 33(4): 5-7 doi: 10.3771/j.issn.1009-2307.2008.04.001 Wang Renxiang. EFP Bundle Triangulation Using Lunar Imagery Obtained from Satellite Three-LineArray Camera[J]. Science of Surveying and Mapping, 2008, 33(4): 5-7 doi: 10.3771/j.issn.1009-2307.2008.04.001
[53] Di K C, Hu W M, Liu Y L, et al. Co-registration of Chang'E-1 Stereo Images and Laser Altimeter Data with Crossover Adjustment and Image Sensor Model Refinement[J]. Advances in Space Research, 2012, 50(12): 1615-1628 doi: 10.1016/j.asr.2012.06.037
[54] 王涛, 项琳, 曹锋. 月球卫星CE-1三线阵影像数据的解算试验[J]. 遥感信息, 2010, 25(3): 18-20 doi: 10.3969/j.issn.1000-3177.2010.03.004 Wang Tao, Xiang Lin, Cao Feng. Study on Calculating Three-Line-Array Imagery Obtained from Satellite Chang'E-1[J]. Remote Sensing Information, 2010, 25(3): 18-20 doi: 10.3969/j.issn.1000-3177.2010.03.004
[55] 平劲松, 黄倩, 鄢建国, 等. 基于嫦娥一号卫星激光测高观测的月球地形模型CLTM-s01[J]. 中国科学: 物理学力学天文学, 2008, 38(11): 1601-1612 https://www.cnki.com.cn/Article/CJFDTOTAL-JGXK200811020.htm Ping Jinsong, Huang Qian, Yan Jianguo, et al. Lunar Topographic Model Based on the Laser Altimetry Observation of Chang 'E-1 Satellite[J]. Science in China : Physics, Mechanics & Astronomy, 2008, 38(11): 1601-1612 https://www.cnki.com.cn/Article/CJFDTOTAL-JGXK200811020.htm
[56] 胡文敏, 邸凯昌, 岳宗玉, 等. 嫦娥一号激光高度计数据交叉点分析与平差处理[J]. 测绘学报, 2013, 42(2): 218-224 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201302011.htm Hu Wenmin, Di Kaichang, Yue Zongyu, et al. Crossover Analysis and Adjustment for Chang'E-1 Laser Altimeter Data[J]. Acta Geodaetica et Carto-graphica Sinica, 2013, 42(2): 218-224 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201302011.htm
[57] Liu Z Q, Di K C, Peng M, et al. High Precision Landing Site Mapping and Rover Localization for Chang'E-3 Mission[J]. Science China Physics, Mechanics & Astronomy, 2015, 58(1): 1-11
[58] 王奋飞, 任鑫, 刘建军, 等. 嫦娥二号影像密集匹配及DEM制作方法[J]. 武汉大学学报· 信息科学版, 2017, 42(7): 918-923 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH201707008.htm Wang Fenfei, Ren Xin, Liu Jianjun, et al. The Methods of Chang'E-2 Image Dense Matching and Lunar DEM Extraction[J]. Geomatics and Information Science of Wuhan University, 2017, 42(7): 918-923 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH201707008.htm
[59] 彭嫚, 万文辉, 吴凯, 等. 嫦娥三号导航相机测图能力分析及地形重建[J]. 遥感学报, 2014, 18(5): 995-1002 https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB201405004.htm Peng Man, Wan Wenhui, Wu Kai, et al. Topographic Mapping Capbility Analysis of Chang'E-3 Navcam Stereo Images and Three-Dimensional Terrain Reconstruction for Mission Operations[J]. Journal of Remote Sensing, 2014, 18(5): 995-1002 https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB201405004.htm
[60] Li F, Zhu C, Hao W F, et al. An Improved Digital Elevation Model of the Lunar Mons Rümker Region Based on Multisource Altimeter Data[J]. Remote Sensing, 2018, 10(9): 1442 doi: 10.3390/rs10091442
[61] Di K C, Zhu M H, Yue Z Y, et al. Topographic Evolution of von Karman Crater Revealed by the Lunar Rover Yutu-2[J]. Geophysical Research Letters, 2019, 46(22): 12764-12770 doi: 10.1029/2019GL085252
[62] 邸凯昌, 刘斌, 刘召芹, 等. 月球遥感制图回顾与展望[J]. 遥感学报, 2016, 20(5): 1230-1242 https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB201605044.htm Di Kaichang, Liu Bin, Liu Zhaoqin, et al. Review and Prospect of Lunar Mapping Using Remote Sensing Data[J]. Journal of Remote Sensing, 2016, 20(5): 1230-1242 https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB201605044.htm
[63] Zou X D, Liu J J, Li C L. Applications of the Visualization of Lunar Surface Based on Chang'E-1 Observations[J]. Chinese Journal of Space Science, 2011, 31(1): 100-105
[64] Zeng X G, Mu L L, Liu J J, et al. Setting Diverging Colors for a Large-Scale Hypsometric Lunar Map Based on Entropy[J]. Entropy, 2015, 17(12): 5133-5144 doi: 10.3390/e17075133
[65] 周增坡, 程维明, 周成虎, 等. 基于"嫦娥一号"的月表形貌特征分析与自动提取[J]. 科学通报, 2011, 56(1): 18-26 https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201101005.htm Zhou Zengpo, Cheng Weiming, Zhou Chenghu, et al. Characteristic Analysis of the Lunar Surface and Automatically Extracting of the Lunar Morphology Based on Chang'E-1[J]. Chinese Science Bulletin, 2011, 56(1): 18-26 https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201101005.htm
-
期刊类型引用(2)
1. 卢万杰，刘伟，牛朝阳，谢丽敏，吕亮. 服务于空天协同对地观测的语义信息模型构建. 地球信息科学学报. 2022(08): 1421-1431 . 
百度学术
2. 李雯静，张馨心，焦宇豪. 基于精细拓扑的矿井排水系统多尺度建模. 金属矿山. 2022(11): 216-221 . 百度学术
其他类型引用(1)
下载:
计量
- 文章访问数: 1118
- HTML全文浏览量: 403
- PDF下载量: 176
- 被引次数: 3
