测绘科学
Sc i ence of Su r veyin g and M app ing
Vo l 133 No 12
第 33卷第 2期 2008年 3月
M a r 1
勇 ① , 霖 ② , 郑新燕 ②
杨 李 ( ①华中农业大学 资源与环境学院 ,武汉
430070; ②武汉大学资源与环境科学学院 , 武汉 430079
) 【摘 要 】当前 , 用成熟的商业数据库管理海量的空间数据已成为当今 G IS 和地图制图软件的主流 , 在作者开发
的地图制图系统中 , 将所使用和处理的数据分为地理数据 、制图辅助数据和地图数据三类 , 本文在分析三类数据 特点的基础上 , 详细介绍了各类数据在数据库中的存储结构 , 特别介绍了如何使用 O rac l e 中的大对象数据类型存 储坐标数据 , 使其该数据库引擎达到功能与效率的统一 。 【关键词 】地图制图 ; 数据库 ; G I S; 存储结构 【中图分类号 】 T P311; P208 【文献标识码 】A 【文章编号 】100922307 ( 2008 ) 02 20137204 DO I: 1013771 / j 1 i ssn 1100922307120081021047
1 引言
经过多年的发展 , G I S 数据或地图数据存储形式依次经 历了以下几个阶段 : 第一阶段 , 用文件存储并管 理所有数 据 ; 第二阶段 , 用文件存储空间坐标数据 , 用数据库存储属 性数据 ; 第三阶段 , 用关系型数据库存储并管理空间坐标数 据和属性数据 , 但空间数据和属性数据在结构上分为多个表 存储 ; 第四阶段 , 基于对象关系模型在数据库中按面向对象 的方法存储空间坐标数据和属性数据 , 显然 , 最
后一种方式 逻辑性和结构性更强 。在地图制图系统中 , 需要处理和存储 三类数据 : 第一类为地理数据 (也称为空间数据 ) , 它是地图 制图的原始数据 , 是对客观显示世界的线划表示 , 它没有宽 度 , 没有样式 , 没有标识 , 它又分为空间坐标数据和属性数 据 ; 第二类为制图辅助数据 , 主要包括符号库 、字体库 、符 号配置规则 、字体配置规则等数据 。第三类为地图数据 , 它 是地理数据经过符号化 , 配置注记 , 解决冲突等制图操作以 后得到的数据 , 即是制图成果数据 , 主要包括空间数据 、制 图图形数据 、符号关联数据及注记数据 , 整饰数据几类 。其 中 , 地图数据是地理数据在制图辅助数据的指导和约束上得 到的 。本文将依次介绍这三种类型数据的存储结构形式 。 2 地理数据的存储
211 地理数据内容分析
地理数据来 源 于 对 实 地 测 绘 数 据 的 数 字 化 跟 踪 结 果 ,
它是按照比例尺和区域分幅存储 , 在单一分幅中 , 又按照 地物类型和几何属性分层存储 , 每一层中存储若干个地物 对象的空间数据和属性数据 。对同一比例尺而言 , 全国地 理数据的分层的形式是统一的 , 表 1 中给出了全国 1 ∶25 万 基础地理信息数据的分层方法 (来源于国家基础地理信息中 心网站 ) 。
表 1 地理数据分层
要素
类型
含几何
类型
层名
主要要素
面状 、线状 河 流 、渠 道 、面 状 湖 泊 、H ydN T 线 、面
水库 、海岸线 、各种滩地等
不依比例水库 、瀑布 、泉 、井 、蓄洪
区范围线 、点状岛屿 、各种礁等
面状居民地 点状居民地 铁路路线及附属物 公路 、其他道路及附属物 境界及行政区划 境界附属物 水系
H ydL K 点 、线 R ESPY R ESPT RA I L K ROAL K
BOUN T BOU PT
面 点 点 、线 点 、线 线 、面 点
居民地 铁路 公路
境界与 政区 沼泽 、盐
田 、水 中 滩 、雪 被 、沙 漠 、TERN T 线 、面 戈壁等
地形
等高线 、等 深
线 、高 程 点 、三 角 点 、TERL K 点 、线 冲沟 、陡石山 、火山 、溶斗等
自然保护区 、长城等 山峰
名 、散列地名中轴线等
其他要素 O T H N T 辅助要素 A T N L K 线 、面 点 、线 从表 1中可以看出 , 地理数据依据要素类型和几何 类 型的不同被划分为 21个图层 。在单一图层中 , 地物的几何 类型和包含的属性项相同 。 212 地理数据存储结构
针对上述的地理数据内容 , 以及 O rac l e 数据库的存 储 方式以及它所提供的 O rac l e S p a t ia l 空间插件对空间数据的 支持 , 我们给出了如下的层次型面向对象的地理数据存储 结构 。
O rac l e 中 , 一个数据库对 象中包含若干个方案 (用户 ) , 每个方案中有表集合 、索引集 合 、视图 集 合 、同 义 词 集 合 、
作者简 介 : 杨 勇 ( 19802) , 男 , 江 西 吉 安人 , 地图学与地理信息系统博士 , 华 中农业大学资源与环境学院讲师 , 现从 事地理信息系统理论与方法研究 。 E 2m a i l: yy04 g is@ 21cn 1com 序列集合 、簇集合 、类型集合 等 , 在本系统中 , 主要用到表 集合来存储地理数据 , 表中又
图 1 地理数据存储结构
包含若干个记录 。这样我们用方案对应于图幅 , 表对应 于 图层 , 记录对应每个地物 , 如图 1所表示 : 213 空间数据的存储方式
按照面向对象的思想 , 每条记录存储一个地物 , 包 括 它的一般属性和空间属性 , 即空间坐标 。可以
用两种方 法 存储空间数据 : 第一种方法是借助 O rac l e 提供的支持空间 数 据 类 型 的 插 件 O rac l e S p a t ia l 中 的 数 据 结 构 类 型
收稿日期 : 2007 201208
基 金 项 目 : 国 家 测 绘 局 基 金 项 目 ( 1469990324231 ) , 数 字 地 图 制 图 知 识 工程
MD SYS 1S DO _ G EOM ETR Y , 其定义如下 :
CR E A T E TYPE S DO _ G EOM ETR Y A S OBJ E CT { S DO _ G TYPE NUMB ER , S DO _
S R I D NUMB ER , S DO _ PO I N T S DO _ PO I N T _ TYPE, S DO _ EL E M _ I N FO MD S YS 1S DO _ E L E M _ I N F O _ ARRA Y ,
S DO _ OR D I N A T ES MD SYS 1S DO _ ORD I N A T E _ AR 2 R A Y } ;
其中 , S DO _ G TYPE 说 明 该 空 间 对 象 的 类 型 ; S DO _ SR I D 记录图形的坐标系统 ; S DO _ PO I N T 存储点对象坐标 ; S DO _ E L E M _ I N FO 为变长数组 , 每三个数为一个单元 , 用 于解释坐标是如何存储在 S DO _ OR D I NA TES 中的 ; S DO _ OR D I N A TES 为变长数组 , 元素类型为 NUMB ER , 用来存储
组成空间对象边界的点坐标 。用这种方式存储空间数据的 优点是结构性和对象性强 , 同时可以利用 O rac le S p a tia l 提 供的空间索引 。缺点是存取速度慢 , 因为访问它的所有内 容需要解析多个数组 , 而 O rac le 提供的访问数据库的接口 对数组的解析需要花费很多的时间 。
第二种方法是利用 O rac le 提供的二进制的大对象类型 BLOB , 它将地物的空间坐标数据和空间对象类型以及其他 一些说明信息按系统设计者规定的规则顺序地存入一个类 型为 BLOB 的 字 段 中 , 当 要 读 取 时 , 再 按 规 则 取 出 即 可 , 这种方法的缺点恰恰是第一种方法的优点 , 而其优点是读 取速度极快 , 这恰恰弥补了第一种方法的缺点 。
综合上述两种方法的优缺点 , 我们将空间数据按两种 方法分 别 存 储 为 两 个 字 段 , 其 字 段 名 分 别 为 S HA P E 和 SHA P E_ B LOB , 当需要读取空间数据时 , 访问大对象类型 字段 SHA P E_ BLOB , 当 要用到空间索引时 , 对 S HA P E 字 段进行操作 。这样做虽然浪费了存储空间 , 但考虑到现在 的硬盘等存储器容量大 , 并且价格便宜 , 牺牲空间换取功
图 2 图元和符号库组织结构
图元代码为内外部关键字进行关联 , 他们的存储关联关系
如图 3所示 。
图 3 符号库存储结构
符号配 置 规 则 用 一 个 表 进 行 存 储 , 其 表 结 构 如 表 2 所示 :
表 2 符号配置规则表结构 符号配置规则表
优先级值 (N um be r ) 符号代码 (N um ber ) 地物国标码 (N um ber ) 当对一个图层配置符号时 , 首先将图层按国标码分 为 若干子图层 , 再按优先级顺序从该表中选取含有该图层中 的国标码的所有记录 , 从记录中获取符号代码 , 依据符 号 代码再从符号库中获取符号参数 , 按优先级从高到低依次
对各类地物配置符号 。从而在符号化的同时解决了各类地 物的绘制顺序问题 。值得一提的是 , 为了实现一些特殊 的 制图要求 (如道路网中不同等级道路相交时 ,交点处融合的 效果 ) , 一个地物可对应多个符号 。 312 字体数据组织
与符号数据类似 , 字体数据部分分为字体参数库和 字 体配置规 则 。相 对 与 符 号 来 说 , 字 体 参 数 描 述 比 较 简 单 , 可以用一个表来存储 , 结构如表 3所示 : 表 3 字体属性表结构 能效率是值 得 的 。这 样 , 每 个 图 层 表 的 字 段 设 置 如 表 所示 。
2
表 2 图层表字段设计
字段名
字段类型
I
D 属性一 NU MB ER (惟一标识 ) N U MB ER 或 V ARCHAR2 属性 N
SH A PE
SH A PE_ BLO B
N U MB ER 或 V ARCHAR2 MD S YS 1S DO _ GEO M ETR Y
BLOB
属性名称 属性类型 备注
字体惟一代码 (表的关键字 ) 单位为 0101mm 单位为 0101mm
如 “宋体 ”, “等线体 ”
依据图式 , 主要有黑 、兰和棕 有 “直立 ”、“左斜 ”和 “耸肩 ”三种
字体代码 字体高度 字体宽度 字体名称 字体颜 字型
数值 数值 数值 字符 数值 数值
3 制图辅助数据存储结构
制图辅助数据包括符号和注记两个部分 。其中符号部
分包括符号库和符号配置规则 , 注记部分包括字体库和字 体配置规则 。
311 符号数据组织
在本系统中 , 符号库有三层 结构 : 符号库 、符号和图 元 。符号分为点线面三种 , 其中面符号层次稍高于点符号 和线符号 , 与点符号和线符号间存在着聚合关系 。点符号 和线符号由点 、线 、面图元组成 , 其具体属性 请参阅参考 文献 [ 2 ]。图元是构成所有符号的图形基础 。符号库和图 元组织结构如图 2所示 :
由此 , 我们设计 8 个表来存储一个符号库 , 表名及含 义分 别 为 PO I N T G U (点 图 元 ) , L I N E G U (线 图 元 ) , AR E 2 A G U (面 图 元 ) , PO I N T S Y MBOL (点 符 号 ) , L I N E S Y MBOL (线符 号 ) , AR E A SY MBOL (面 符 号 ) , M S PEN (画 笔 参 数 表 ) , M S BRU S H (画刷参数表 ) 。各个表中 , 以符号代码或
与符号配置 规 则 中 一 个 条 件 就 可 确 定 一 个 符 号 相 比 ,
字体配置规则比较复杂 , 因为同一类型 (国标 ) 的地物 , 字 体的大小随着该地物的等级而变化 , 即一个字体的约束于 多个条件 , 而且这些条件有时并不等同于该地物某个属性 值的全部 , 而是要从某个属性值的某一位或几位的值来确 定 。例如 , 在 1∶25 万地形图中 , 单线河流分为六级 , 不同 级别所
用字体大小不同 , 其国标码 ( G B ) 都为 21011 , 而级 别信息来自于属性 H y dC 的最后一位的数值 。为此 , 我们用 四个参数描述一个条件 , 即 : 属性名称 , 取值开始位 , 取 值结束位 , 对应值 , 每条字体规则最多设置 5 个条件 。由 此 , 字体规则表结构如表 4所示 :
139 第2期杨勇等地图制图系统数据库存储结构设计
号图形序列(红点在折线和折线多边形中表示节点,在圆
中表示圆心和圆上一点) 。
符号图形的数据分为两类: 第一为图形数据, 即坐标
点数组, 第二类为绘制属性数据, 包括画笔、画刷属性和
圆长短半径。由于图幅中使用相同绘制类型的对象很多,
如果重复存储相同的绘制属性显然会浪费比较多的存储空
间, 为此, 我们以图幅为单位, 维护一个绘制属性结构表,
通过惟一标识符, 确定各个符号图形所用的绘制设备, 该
表(DR AW D E V I CE)结构如表5所示:
表4 字体规则表结构
属性名称属性类型备注
字体代码数值作为字体表的外部关键字
由地物的优先级派生, 解决注记压盖时
作为参考
图层类型标识
由图式规定, 确定该字体所需要的条
件, 《= 5
优先级数值
图层类型
使用条件
个数
条件1 属
性名
条件1 取值
开始位置
条件1 取值
结束位置
数值
数值
表5绘制属性结构表
字符条件一对应与地理数据库表中的字段名
属性名称属性类型备注
D EV I C
E I D
PEN _ D W PEN S TY L E
PEN _ D W PENW I D TH
数值
数值
数值
设备I D(主键)
画笔样式画笔宽
度画笔的其他属
性画刷样式画
刷颜画刷纹理
样式短半径
长半径数值
数值
符合该条件时, 条件一指定字段指定位
置的对应值
以下为另外四个条件的描述, 方式同条
件一
条件1 取值字符
BRU SH _ LB S TY L E
BRU SH _ LBCO LOR
BRU SH _ LBHA TCH
LLO N GRA D I U S
数值
数值
数值
数值
条件2 属性名字符
L S H OR TRA D I U S 数值
当为一个图层配置注记时, 从上表中选择所有与该图
层类型一致的注记配置规则, 再将地物的属性与每条规则
中定义的所有条件进行对比, 如果全部符合, 则为该地物
的注记配置规则所对应的字体。
4 地图数据存储结构
地图数据包括4个部分, 第一部分为地理数据, 第二部
gis地图制作
分地物与符号的对应信息, 第三部分为符号图形数据, 第四
部分为注记数据。其中, 地理数据作为地图数据的基础被带
进地图数据中, 它的存储方式与第1节中介绍的相同。
411 地物符号对应信息存储
一个图层被符号化时,
是分成了若干个子图层,
每个子图层再配上若干个
符号完成的。所以地物与
符号之间通过子图层来联
在一个图层中, 符号图形的坐标数据量非常大, 如果
象存储绘制属性一样在表中每个单独存成一行, 那么当读
取符号图形时, 系统与数据库的I/O 交换量将十分巨大,
速度将十分缓慢。为此, 本系统中将全图的符号图形存储
在一个表中, 以图层名为标识, 用二进制大对象B LOB 存
储具体的坐标数据。
在属性S Y M F IG UR E BLOB 中, 以64位字节为单位, 以
该图层中用到的符号为第一层循环标志, 该符号所对应的
子图层中的地物对象为第二层循环标志, 以地物对象中所
包含的符号图形为第三层循环标志, 以每个符号图形的点
作为第四层循环标志进行存储, 如图6所示:
系的, 如图4所示:图4 地物符号对应图
为此, 为每个图层设计两个表, 假设图层名为L a ye r2
N a m e, 则第一个表名为L a ye r N a m e_ S, 该表记录该图层中
的子图层信息, 第二个表名为L a ye r N a m e_ S P, 记录地物与
子图层的联系信息, 即哪个地物属于哪个子图层。第三个
表名为L a ye r N a m e_ D ,记录该图层中所用到的符号信息。
412 符号图形数据的存储。
符号化的结果数据称为符号图形, 它是地理数据和符
号结合的产物, 是该制图系统的图形基础, 它的存在能加
快系统的绘图速度, 而且施加于其上的编辑操作能改变符
号化后的图形细节, 达到一些特殊的制图效果。针对地图
图形的特点, 将符号图形划分为种基本类型, 即圆、线和
多边形。所有复杂的图形都可以分解为这三种图形。如图5
所示:
图6 符号图形数据存储格式
(每个矩形框代表一个层次的循环)
413 地图注记数据存储
注记是对图形的文字标示, 按所标示的图形几何类型
分为点状注记, 线状注记和面状注记。但无论是何种注记
目的都是将文字字符串摆放在合适的位置, 结果无非是两
种: 一是字符串挨在一起, 此时, 字符串只有一个定位点,
既第一个字符的左上角, 如图7 a; 二是字符串分隔开来,
此时, 字符串有N 个定位点, 分别位于每个字符的左上角,
如图7b和7 c。所以, 注记分为两个层次, 一是注记整体,
表示对一个地物所标示文字的整体信息, 二是注记细胞,
对应注记整体中的有一个定位点定位的字符串对象。注记
整体是由若干个注记细胞组成的。
注记信息分图层存储, 每个图层的注记用两个表存储,
假设图层名为L a ye r N a m e, 则两个表的表名分别为L a ye r
N a m e_ A n no和L a ye r N a m e_ A n noCe l l。其表结果如表6和表
7所示:
图5 符号图形示例
在图5a中用境界符号符号化线状地物以后, 得到一串
圆和线顺序排列的符号图形序列, 在图5b中, 用灯塔符号
符号化点状地物后, 得到由六条线和一个多边形组成的符
5 结束语
本文详细介绍了基于 O rac le 数据库的地图制图系统的 数据库结构设计 , 系统中充分利用了 O rac le 数据库提供的 O rac l e S p a t ia l 空间插件和二进 制 大 对 象 数 据 类 型 , 达 到 了 功能与效率的统一 。同时将系统中数据分为三类存储 , 其 存储结构清晰 , 方便管理 。目前 , 基于该数据库结构的 地 图制图系统已在陕西测绘局和国家基础地理信息中心进行 了实验性生产 , 结果表明系统运行稳定 , 成图质量高 , 速 度快 , 特别是系统与数据库数据交互效率高 。
参考文献
图 7 地物图形和地图注记示意图
表 6 注记整体存储表结构
属性名称
属性类型
备注
AN NO I D A LLANN O STR AN NOCELLCOU N T A NNO TY PE SH A PE I D FON T COD E SER I A LNO
数值
字符
数值
数值 数组 数值
数值 惟一标识 (主键 )
整体字符 含有的注记细胞个数
注记类型 (点 、线或面 ) 注记对应地物的 I D 串 字体代码 优先级
O rac l e 公司 1o r ac l e 817 s p a t ia l g u i de [ S ]. 2000291
李 志 涛 1 地 图 制 图 系 统 中 组 件 技 术 的 应 用 与 实 现
[ D ]. 武汉 :武汉大学 , 20051
毋河海 1地图数据库系统 [M ]. 北京 : 测绘出版社 ,
19912111 祝 国 瑞 1 地 图 学 [ M ]. 武 汉 : 武 汉 大 学 出 版 社 , 2004211 田娇娇 , 等 1动态数据库模型的研究与应用 [ J ]. 测 绘科学 , 2006 , 31 ( 1) 1
田增平 , 周傲英 , 施伯乐 1适合地图数据库应用的可 扩充数据模型研究 [ J ]. 计算机工程 , 1997, 23 ( 4) 1
祝国瑞 , 郭礼珍 , 尹贡白 , 徐永利 1 地图设计与编绘 [M ] 武汉 : 武汉大学出版社 , 2001 2101
于冬梅 , 董罗海 , 张力果 1 数字地图制图理论方法与 应用 [ J ]. 地理信息科学 , 2003, ( 6 ) : 5 271
[ 1 ] [ 2 ] [ 3 ] 表 7 注记细胞存储表结构
[ 4 ] 属性名称
属性类型
备注
对应于注整体的 I D (外键 ) 细胞序号 , 与 AN NO I D 组成主键 注记细胞字符串 定位点 X 坐标值 定位点 Y 坐标值 是否显示背景 背景数值 倾斜角度
是否存在冲突标识
AN NO I D CELL I N D EX A NNOCELL
TEXT STAR T X STAR T Y V I E WBACK COLOR BACK COLOR ESCA
PE M EN T
BH AV ECOLL I S I O N
数值 数值 字符 数值 数值 数值 数值 数值 数值 [ 5 ]
[ 6 ] [ 7 ]
[ 8 ]
T he stora g e 2struc t ure de s i g n of m a p p i n g sy s t em da t a b a se
A b stra c t: Cu rren t ly, u sing m a tu re bu sine ss da taba se to m anag e m a ssive sp a tia l da ta ha s becom e the m a in stream of tod ay ’s G I S
and m app ing system s 1 I n the ca rtograp h ic system deve lop ed by the au tho r , the da ta u sed and p roce ssed is c la ssif ied in to th r ee ca t eg o 2
rie s: geo sp a tia l da ta, m app ing 2a ssisted da ta and m a Pda ta 1 B a sed on ana lyzing the cha rac te r istic s of the se th ree ca t eg o rie s, th i s a r ti c l e in tr oduce s the sto rag e struc tu re of each k ind of da ta in de ta il, e sp ec ia lly how to u se the b ig ob jec t typ e in O rac le to sto re da ta w i th c o o r 2 d i na t e s , wh i ch un i f ie s the func t iona l ity and effic i ency of th i s da t aba s e 1
Key word s : ca r tograp h y; da t aba s e; G IS ; sto r ag e struc t u r e
YAN G Yong ①
, L I L in ②
, ZH E N G X i n 2yan
②
( ①Co l leg e of R e s o u r ce s an d Environm e n t , H u azhong A gricu l tu r a l U n i ve r sity, W u 2
han 430070, Ch i na; ②S choo l O f R e s ou r ce A n d Envir onm e n t a l S c i ence, W uhan U n i ve r sity, W uhan 430079 , Ch i na )
(上接第 184页 )
[ 2 ] 苏洁 , 周东方 , 岳春生 1G P S 车载导航中的实时地图 匹配算法 [ J ]. 测绘学报 , 2001 , 30 ( 3) : 25222561 [ 3 ] 王楠 , 王勇峰 , 刘积仁 1一个基于位置点匹配的地图
匹配算法 [ J ]. 东北大学学报 (自然科学版 ) , 1999 , 20 ( 4) 1 [ 4 ] 谢彩香 , 林宗坚 , 等 1G PS /DR /MM 组合导航中的车 彭飞 , 柳重善 , 张其善 1基于代价函数的组合导航系
统地图匹配算法 [ J ]. 北京航空航天大学学报 , 2002, 28 ( 3 ) : 261 22641 O c h i eng W Y, Q u ddu s M , No l and R B 1MA P 2MA T C H 2 I N G I N C OM P L E X UR BAN R OAD N E T WOR KS [ J ]. B razilian Jo u r na l of Ca r tograp h y ( R e vista B ra s ile i ra de Ca r tograf ia ) , 2003, 55 ( 2 ) : 12181 Q u ddu s M A , O c h i en g W Y, Zhao L e t a l 1 A Gene r a l m a Pm a t ch i ng A lg o r ithm fo r Tran s po r t Te l em a t ic s App lica t ion s [ J ]. G P S S o l u t ion s Jou r na l, 2003, 7 ( 2) : 157 21671
[ 6 ] [ 7 ] 31 ( 1) 1 辆定位精度研究 [ J ]. 测绘科学 , 2006 , [ 5 ] 张其善 1智能车辆定位导航系统及应用 科学出版社 , 20021
[ 8 ] [M ]. 北京 : D e s i g n an d i m p l em e n t a t i on of an i n t egra t ed m a p m a tch i n g m lgor i t hm
A b stra c t: The p ap e r inve stig a ted and rea lized a in tegra ted m ap m a tch ing a lg o rithm acco rd ing to the p rinc i p le of po in t 2b a s ed m a p
m a t ch ing a lg o rithm and we ig h t 2ba sed m ap m a tch ing a lg o rithm 1 I t u sed d iffe ren t m ap m a t ch ing m e thod a t d iffe ren t r oad cond iti o n s , and a lso con side red r oad ne two rk connec t ivity 1 The m a t ch ing accu racy is 86% , and m a tch ing ti m e is ab ou t 0134 s a t ac r o ss 1 S o it can s a t is 2 fy the requ i rem e n t s of veh i c l e navig a t ion syste m 1
Key word s : G PS ; veh i c l e po s iti o n i ng ; m a p m a t ch i ng
L I N N a ① , L I Z h i ① , W A N G B in ②
( ①Chongq i n g Enginee r ing Techno l og y R e s ea r ch Cen t e r fo r Info r m a t i o n M a nag em e n t in D e ve l 2 opm e n t , Chong q ing Techno l og y and B u s ine s s U n i ve r sity, Chong
q i ng 400067, Ch i na; ② C hong q ing Geom a t ic s Cen t e r , Chongq i n g
401121, Ch i na )