CN102938230A - 电子地图上二维空间标量场的动态分层渲染表达方法 - Google Patents
电子地图上二维空间标量场的动态分层渲染表达方法 Download PDFInfo
- Publication number
- CN102938230A CN102938230A CN201210460475XA CN201210460475A CN102938230A CN 102938230 A CN102938230 A CN 102938230A CN 201210460475X A CN201210460475X A CN 201210460475XA CN 201210460475 A CN201210460475 A CN 201210460475A CN 102938230 A CN102938230 A CN 102938230A
- Authority
- CN
- China
- Prior art keywords
- data
- boundary
- property value
- electronic chart
- interval
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000009877 rendering Methods 0.000 title abstract description 6
- 238000006243 chemical reaction Methods 0.000 claims abstract description 8
- 239000000284 extract Substances 0.000 claims abstract description 6
- 239000003086 colorant Substances 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 5
- 239000007787 solid Substances 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 claims description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 6
- 238000000605 extraction Methods 0.000 abstract description 2
- 238000004040 coloring Methods 0.000 abstract 1
- 239000010410 layer Substances 0.000 abstract 1
- 239000002344 surface layer Substances 0.000 abstract 1
- 239000004744 fabric Substances 0.000 description 2
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Processing Or Creating Images (AREA)
- Image Processing (AREA)
Abstract
本发明涉及电子地图上二维空间标量场的动态分层渲染表达方法，其针对存在复杂的空间嵌套关系的边界多边形集合，利用空心环是由双边界组成的思想，提出根据包含多边形的奇偶性判断嵌套关系，提取矢量面的方法。该方法能够一次性的，以较快的速度完成场数据（栅格数据）的分层矢量面的提取，并返回给电子地图前台程序，进行不同区域面状图层的渲染。与传统方法相比，该方法计算速度较快，可以一次性提取出场数据中不同区间的所有矢量面，解决了传统电子地图或WebGIS中对网格数据的渲染因数据量大、分层转换慢而造成的实时性、动态性差的问题。在与其他地图图层进行叠置显示时，能够较好地展现“透明”或“半透明”设色效果。
Description
技术领域
本发明涉及电子地图上二维空间标量场的动态分层渲染表达方法，属于地图制图、地理信息工程技术领域。
背景技术
在电子地图、大气/水文等过程模拟时，常需对实时计算的场数据（栅格数据）进行表达；目前，主要的表达方式是采用基于二维空间标量场的分层显示手段。即将场数据依据不同的值区间划分为多个区域，以不同的色彩渲染出来。
然而，在电子地图或Web GIS中，对网格数据渲染存在一些技术难点：a. 数据量大，转换、传输、渲染慢，难以满足实时性要求；b. 网格边界效果不好；c. 如果要与其他地图图层进行叠置显示，则很难使用“透明”或“半透明”效果。
发明内容
本发明要解决技术问题是：克服现有技术的上述不足，提供一种电子地图上二维空间标量场的动态分层渲染表达方法，其能够以较快的速度动态、实时对场数据（栅格数据）进行渲染，且能得到较好的表达效果。
为了解决以上技术问题，本发明提供的一种电子地图上二维空间标量场的动态分层渲染表达方法，包括以下步骤：
步骤1、电子地图后台程序周期性的接收有关于大气或水文的场数据，所述场数据为格式统一的栅格数据；
步骤2、针对场数据中栅格的属性值预先设置有若干属性值区间，电子地图的后台程序依次根据各属性值区间从场数据中提取相应的区域边界，并将每个属性值区间所对应的区域边界转换为相应属性值区间的矢量面集合，与同一属性值区间对应的矢量面赋予同一属性值；
步骤3、电子地图的后台程序将所有属性值区间对应的矢量面集合传回给电子地图的前台程序；
步骤4、电子地图的前台程序根据矢量面数据的属性值，对矢量面一一进行渲染，相同属性值的矢量面渲染相同颜色，不同属性值的矢量面渲染不同的颜色。
进一步的所述步骤2中，获得与属性值区间对应的矢量面集合的方法如下：
步骤a、对整个场数据进行二值化，若栅格的属性值在该属性值区间内，则将该栅格的属性值赋予1，否则赋予0；
步骤b、对整个场数据的1、0交接处进行边界追踪，得到该属性值区间在场数据中的所有区域边界；
步骤c、采用内插法对得到的区域边界形做平滑处理，消除区域边界的锯齿；
步骤d、将区域边界组成边界多边形集合，在边界多边形集合中，区域边界根据面积的从小到大排序；
步骤e、对边界多边形集合中面积最小的区域边界进行判断，如果包含它的区域边界是偶数个，则该区域边界是一个单独的实心面，把该区域边界转化为面数据添加至矢量面集合中，并将该区域边界从边界多边形集合中删去；如果包含它的区域边界数目是奇数个，则与包含它的面积最小的区域边界组成一个环形矢量面，将该环形矢量面添加至矢量面集合中，并将构成该环形矢量面的两个区域边界一并从边界多边形集合中删去;
步骤f、重复步骤e，直到边界多边形集合为空，此时已获得与该属性值区间对应的矢量面集合。
本发明将栅格数据根据其属性值进行划分，并对属于同一区域的栅格进行边界追踪，获得有关于该区域的若干个边界，由于边界只是线数据，因此要将其转化为面数据。本步骤利用多边形嵌套原理进行面数据的转化，该方法构思巧妙，简单易行，处理效率高。
有时传来的场数据中，目标区域的栅格顶在场数据的边缘，从而导致经步骤C平滑处理后，场数据的边缘信息会有一定的损失，因此为了更好地还原数据真实性，取得更好的直观性，本发明执行步骤a之前，首先对场数据外扩5*5个栅格作为该场数据的外边界，并对外扩的栅格赋予不属于任一预设属性值区间的值。这样一来，数据场边缘进行了一定程度的外扩，使目标能更准确清晰的体现。
有时传来的场数据数据格式不同一，那么就需要有个中转站对数据进行规范化；于是本发明设计了一个与电子地图后台程序进行数据传输的数据转化单元，所述数据转化单元读取有关于大气或水文的原始数据，并将该原始数据转化为能够被电子地图后台程序读取的场数据。该数据化单元良好的解决了该问题。
本发明针对存在复杂的空间嵌套关系的边界多边形集合，利用空心环是由双边界组成的思想，提出根据包含多边形的奇偶性判断嵌套关系，提取矢量面的方法。该方法能够一次性的，以较快的速度完成场数据（栅格数据）的分层矢量面的提取。与同传统方法相比，该方法计算速度较快，可以一次性提取出场数据中不同区间的所有矢量面，解决了传统电子地图或Web GIS中对网格数据的渲染因数据量大、分层转换慢而造成的实时性、动态性差的问题。在与其他地图图层进行叠置显示时，能够较好地展现“透明”或“半透明”设色效果。
附图说明
下面结合附图对本发明作进一步的说明。
图1是一幅简单的栅格数据的局部区域。
图2是以区间(0.5,1.5)对该区域栅格数据进行二值化。
图3是根据二值化结果进行边界追踪。
图4是对图3得到的边界多边形做平滑处理。
图5是图1在区间 (0.5,1.5)下得到的渲染结果。
具体实施方式
下面根据附图详细说明本发明，本发明的目的和效果将变得更加明显。
电子地图上二维空间标量场的动态分层渲染表达方法，包括以下步骤：
步骤1、电子地图后台程序周期性的接收有关于大气或水文的场数据，所述场数据为格式统一的栅格数据；
本步中，电子地图后台程序对传来的场数据判断是否满足格式要求，是否为有效数据。如果格式不符合要求，则由数据转化单元将其转化为能够被电子地图后台程序读取的场数据。然后，给出执行命令和区间集合，令后台对场数据（栅格数据） 开始计算处理。
步骤2、针对场数据中栅格的属性值预先设置有若干属性值区间，电子地图的后台程序依次根据各属性值区间从场数据中提取相应的区域边界，并将每个属性值区间所对应的区域边界转换为相应属性值区间的矢量面集合，与同一属性值区间对应的矢量面赋予同一属性值。具体步骤如下：
首先对场数据（栅格数据）向外扩5*5个栅格，并赋予不属于区间集合内任何区间的值（比如很大的值9999），作为该场数据（栅格数据）的外边界。如图1所示，为场数据（栅格数据）中某一区域的栅格值，为了说明方便，这里使用较简单的数字表示该区域内的栅格值。
对整个场数据进行二值化，若栅格的属性值在该属性值区间内，则将该栅格的属性值赋予1，否则赋予0；本例中，以 (0.5，1.5) 区间进行二值化为例，图2是对图1二值化后的结果。
对整个场数据的1、0交接处进行边界追踪，得到该属性值区间在场数据中的所有区域边界；图3是边界追踪的结果，从图中可以看出，共追踪出3条闭合的区域边界线。由于场数据（栅格数据）存在外边界，所以最终得到的所有区域边界都是闭合的多边形。
采用内插法对得到的区域边界形做平滑处理，消除区域边界的锯齿；处理结果见图4。
将区域边界组成边界多边形集合，在边界多边形集合中，区域边界根据面积的从小到大排序；
由上述几幅图可以看出，该区域在区间（0.5,1.5）上得到的实际矢量面包括一个空心环和一个实心面，目前得到的是一个存在层层嵌套关系的三个边界多边形的集合（见图4，从外向内依次是0-1交接线、1-0交接线、0-1交接线）。
因此接下来需要对这三个边界多边形做嵌套处理，方法如下：
如图4，先从内部最小的多边形开始，因为包含它的多边形有两个，说明它是一个单独的实心面，将其移出边界多边形集合，移入矢量面集合，然后判断面积次小的多边形，即中间那个多边形，因为包含它的多边形只有一个，说明它与这个多边形会组成一个空心环，将这两个多边形做个拓扑嵌套，移入矢量面集合，此时边界多边形集合为空，证明所有边界多边形都已处理，判断结束。
当所有边界多边形判断结束后，最终得到的矢量面集合即为该区间在场数据（栅格数据）下的所有矢量面，将每个矢量面赋予该区间的相应属性值，如图4中各矢量面可以赋予属性值1（0.5与1.5的中间值）。
步骤3、电子地图的后台程序将所有属性值区间对应的矢量面集合传回给电子地图的前台程序；并以WKT文的格式将矢量面集合返回给电子地图。
步骤4、电子地图的前台程序根据矢量面数据的属性值，对矢量面一一进行渲染，相同属性值的矢量面渲染相同颜色，不同属性值的矢量面渲染不同的颜色。图5是对区间（0.5,1.5）得到的矢量面渲染后的结果。
除上述实施例外，本发明还可以有其他实施方式。凡采用等同替换或等效变换形成的技术方案，均落在本发明要求的保护范围。
Claims (4)
1. 电子地图上二维空间标量场的动态分层渲染表达方法，包括以下步骤：
步骤1、电子地图后台程序周期性的接收有关于大气或水文的场数据，所述场数据为格式统一的栅格数据；
步骤2、针对场数据中栅格的属性值预先设置有若干属性值区间，电子地图的后台程序依次根据各属性值区间从场数据中提取相应的区域边界，并将每个属性值区间所对应的区域边界转换为相应属性值区间的矢量面集合，与同一属性值区间对应的矢量面赋予同一属性值；
步骤3、电子地图的后台程序将所有属性值区间对应的矢量面集合传回给电子地图的前台程序；
步骤4、电子地图的前台程序根据矢量面数据的属性值，对矢量面一一进行渲染，相同属性值的矢量面渲染相同颜色，不同属性值的矢量面渲染不同的颜色。
2. 根据权利要求1所述的电子地图上二维空间标量场的动态分层渲染表达方法，其特征在于所述步骤2中，获得与属性值区间对应的矢量面集合的方法如下：
步骤a、对整个场数据进行二值化，若栅格的属性值在该属性值区间内，则将该栅格的属性值赋予1，否则赋予0；
步骤b、对整个场数据的1、0交接处进行边界追踪，得到该属性值区间在场数据中的所有区域边界；
步骤c、采用内插法对得到的区域边界形做平滑处理，消除区域边界的锯齿；
步骤d、将区域边界组成边界多边形集合，在边界多边形集合中，区域边界根据面积的从小到大排序；
步骤e、对边界多边形集合中面积最小的区域边界进行判断，如果包含它的区域边界是偶数个，则该区域边界是一个单独的实心面，把该区域边界转化为面数据添加至矢量面集合中，并将该区域边界从边界多边形集合中删去；如果包含它的区域边界数目是奇数个，则与包含它的面积最小的区域边界组成一个环形矢量面，将该环形矢量面添加至矢量面集合中，并将构成该环形矢量面的两个区域边界一并从边界多边形集合中删去;
步骤f、重复步骤e，直到边界多边形集合为空，此时已获得与该属性值区间对应的矢量面集合。
3. 根据权利要求3所述的电子地图上二维空间标量场的动态分层渲染表达方法，其特征在于：执行步骤a之前，首先对场数据外扩5*5个栅格作为该场数据的外边界，并对外扩的栅格赋予不属于任一预设属性值区间的值。
4. 根据权利要求3所述的电子地图上二维空间标量场的动态分层渲染表达方法，其特征在于：具有一个与电子地图后台程序进行数据传输的数据转化单元，所述数据转化单元读取有关于大气或水文的原始数据，并将该原始数据转化为能够被电子地图后台程序读取的场数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210460475.XA CN102938230B (zh) | 2012-11-15 | 2012-11-15 | 电子地图上二维空间标量场的动态分层渲染表达方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210460475.XA CN102938230B (zh) | 2012-11-15 | 2012-11-15 | 电子地图上二维空间标量场的动态分层渲染表达方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102938230A true CN102938230A (zh) | 2013-02-20 |
| CN102938230B CN102938230B (zh) | 2014-12-24 |
Family
ID=47697122
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210460475.XA Expired - Fee Related CN102938230B (zh) | 2012-11-15 | 2012-11-15 | 电子地图上二维空间标量场的动态分层渲染表达方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102938230B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105760581A (zh) * | 2016-01-29 | 2016-07-13 | 中国科学院地理科学与资源研究所 | 一种基于osg的沟道流域整治规划仿真方法及系统 |
| CN104602271B (zh) * | 2013-11-01 | 2018-11-02 | 中国移动通信集团设计院有限公司 | 一种覆盖优化中问题区域的生成方法及装置 |
| CN113485411A (zh) * | 2021-06-21 | 2021-10-08 | 安徽农业大学 | 一种航空精准施药三维航线规划方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1584932A (zh) * | 2004-06-10 | 2005-02-23 | 西安煤航遥感信息有限责任公司 | 用于图像转图形中边界线追踪的优化方法 |
| CN101488221A (zh) * | 2009-01-22 | 2009-07-22 | 浙江工商大学 | 遥感影像图中道路矢量化方法 |
-
2012
- 2012-11-15 CN CN201210460475.XA patent/CN102938230B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1584932A (zh) * | 2004-06-10 | 2005-02-23 | 西安煤航遥感信息有限责任公司 | 用于图像转图形中边界线追踪的优化方法 |
| CN101488221A (zh) * | 2009-01-22 | 2009-07-22 | 浙江工商大学 | 遥感影像图中道路矢量化方法 |
Non-Patent Citations (1)
| Title |
|---|
| 尹潘等: "《综合数学形态矢量化在城市规划中的应用》", 《计算机工程》, vol. 34, no. 1, 31 January 2008 (2008-01-31) * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104602271B (zh) * | 2013-11-01 | 2018-11-02 | 中国移动通信集团设计院有限公司 | 一种覆盖优化中问题区域的生成方法及装置 |
| CN105760581A (zh) * | 2016-01-29 | 2016-07-13 | 中国科学院地理科学与资源研究所 | 一种基于osg的沟道流域整治规划仿真方法及系统 |
| CN113485411A (zh) * | 2021-06-21 | 2021-10-08 | 安徽农业大学 | 一种航空精准施药三维航线规划方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102938230B (zh) | 2014-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105069020B (zh) | 自然资源数据三维可视化方法和系统 | |
| CN103093499B (zh) | 一种适用于网络传输的城市三维模型数据组织方法 | |
| CN103559374B (zh) | 一种多子网格模型上进行面分裂型曲面细分的方法 | |
| CN101393649A (zh) | 多视点图像渲染装置及方法 | |
| CN101630419A (zh) | 一种用于城市综合管网三维可视化系统的架构方法 | |
| CN101976468B (zh) | 一种多分辨率动态地形可视化方法及系统 | |
| CN103927424B (zh) | 利用体网格简化建筑三维模型结构的方法 | |
| Zhao | Application of 3D CAD in landscape architecture design and optimization of hierarchical details | |
| CN102663800A (zh) | 一种顾及城市意象的城市建筑综合与渲染的方法 | |
| CN101833780A (zh) | 一种基于游程表达和运算的地图成图方法 | |
| CN103927779B (zh) | 一种基于配置的二维动画生成方法 | |
| CN108022273A (zh) | 一种图数可拆合制图方法及系统 | |
| CN103514625A (zh) | 基于多视角自适应特征配准的三维重建方法 | |
| CN103927781A (zh) | 一种三维场景下视觉相关的模型渐进渲染方法 | |
| CN102938230B (zh) | 电子地图上二维空间标量场的动态分层渲染表达方法 | |
| CN105488840B (zh) | 一种信息处理方法及电子设备 | |
| CN102881028A (zh) | 一种栅格数字图像快速矢量化方法 | |
| CN103345534A (zh) | 一种动态图处理方法及装置 | |
| CN103473984A (zh) | 网络环境下基于模板的动态地图获取方法 | |
| CN101789030A (zh) | Vge符号模型及基于该模型的地图符号共享系统和方法 | |
| CN100407231C (zh) | 用于图像转图形中边界线追踪的优化方法 | |
| CN103927395A (zh) | 结构简化图像三维几何及材质纹理的数据结构及渲染方法 | |
| CN106649776A (zh) | 一种半自动化综合矢量多边形的方法 | |
| Franke et al. | Procedural generation of 3D karst caves with speleothems | |
| CN105512194A (zh) | 一种游戏场景的管理方法和装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141224 Termination date: 20151115 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |