CN104750854A - 一种海量三维激光点云压缩存储及快速加载显示方法 - Google Patents

一种海量三维激光点云压缩存储及快速加载显示方法 Download PDF

Info

Publication number
CN104750854A
CN104750854A CN201510178971.XA CN201510178971A CN104750854A CN 104750854 A CN104750854 A CN 104750854A CN 201510178971 A CN201510178971 A CN 201510178971A CN 104750854 A CN104750854 A CN 104750854A
Authority
CN
China
Prior art keywords
data
point cloud
color
block
cloud
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
Application number
CN201510178971.XA
Other languages
English (en)
Other versions
CN104750854B (zh
Inventor
龚书林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
WUHAN HAIDASHU CLOUD TECHNOLOGY Co Ltd
Original Assignee
WUHAN HAIDASHU CLOUD TECHNOLOGY Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by WUHAN HAIDASHU CLOUD TECHNOLOGY Co Ltd filed Critical WUHAN HAIDASHU CLOUD TECHNOLOGY Co Ltd
Priority to CN201510178971.XA priority Critical patent/CN104750854B/zh
Publication of CN104750854A publication Critical patent/CN104750854A/zh
Application granted granted Critical
Publication of CN104750854B publication Critical patent/CN104750854B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

本发明公开了一种海量三维激光点云压缩存储及快速加载显示方法,属于三维激光测量技术领域,通过采集点云数据总体结构、按照不同级别对数据进行分类、压缩,每一级别包括块集、块和包三级索引,实现高度存储和快速加载显示,基于本发明存储格式可解决海量点云占用空间大,不易拷贝问题,编辑分类问题,基于本方法实现的三维激光点云存储、浏览、编辑时具有点云存储空间小,点云浏览处理软件占用内存小,在浏览过程中点云占用内存不超过100Mb;首次加载显示及实时浏览速度快,在浏览点云过程中,不会卡顿现象,实时动态刷新点云,刷新响应时间小于1s;支持编辑分类等优点。

Description

一种海量三维激光点云压缩存储及快速加载显示方法
技术领域
本发明涉及三维激光测量技术领域,具体是一种海量三维激光点云压缩存储及快速加载显示方法。
背景技术
目前三维激光测量技术正在被广泛应用于测绘领域,可以快速采集获取地物表面三维激光点云,三维激光扫描根据载体不同,可以划分为地面站、车载、机载和背包等形式;三维激光扫描仪在移动载体上可以连续扫描测量沿途两侧空间地理要素,获取海量的三维激光点云,目前主流的激光扫描仪如RIEGL、Z+F、Faro及国内的中海达HS300等,可以每秒采集获取30-100万点云;三维激光扫描仪采集获取的点云通常包含三维坐标、强度、时间和颜色等信息,后处理还包括分类信息。点云浏览显示通常包含按高程渲染、按强度渲染、按真彩色渲染、按分类渲染,不同渲染方式对应加载点云属性信息不同;三维激光测量技术的发展与应用的同时,另一方面又伴随着海量点云数据存储与快速显示的问题,目前国内外激光扫描仪厂商配套软件都普遍存在以下问题:占用内存大、首次加载过程缓慢、显示浏览不顺畅、无法连续处理编辑海量点云、数据存储占用空间大;以上问题的存在,在一定程度上制约了三维激光点云预处理及后续加工生产应用。
发明内容
本发明的目的在于提供一种能解决点云编辑问题,实现海量三维激光点云压缩存储及快速加载显示方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:
一种海量三维激光点云压缩存储及快速加载显示方法,包括以下内容:
1)激光点云高度压缩存储方法
1.1)点云数据总体结构:
a.整体由若干层级别构成,从低到高,0级是最详细原始数据,高等级由低等级数据抽稀构成,每一级别对应不同显示比例尺;
b.坐标、强度、时间合并文件存储,颜色、分类单独文件存储;
c.数据文件按级别顺序存储,0,1,2...n;
d.每一级别按坐标、强度、时间顺序存储;
e.颜色数据,通过与点云配准后的影像对点云进行着色后生成;
f.除分类数据无压缩外,其他数据都压缩存储,压缩算法采用FastAC;
g.选择、删除编辑标记信息存放在分类中;
1.2)块集、块、包组织关系:
a.每一级内部包括三级索引:块集、块、包,块集包含块,块包含包;
b.编号从左到右,从下到上;
c.分块包括2×2、4×4、8×8…,默认块集边长64,1块集包含4×4块,1块包含4×4包;
d.每层三级索引分块需均衡,每数据包最大存储1M数据,每个数据包索引包含格网划分规则、点数、范围、坐标、强度、时间、颜色、分类等数据地址信息;
2)数据压缩存储实现方法
2.1)xyz坐标采样网格相对坐标,每个坐标值采用2个字节表达,在保证坐标精度为0.001m情况下,网格边长约设为64m;
2.2)强度属性采用1个字节,表达范围0-255;
2.3)时间属性采用2个字节,精度为1/128秒情况下,表达范围0-512秒;
2.4)颜色属性采用2个字节,红色5bit、绿色6bit、蓝色5bit,即R5G6B5,RGB和R5G6B5转换方法如下:
(a)RGB转R5G6B5 (u16)((color&0x00F80000)>>8|
                     (color&0x0000FC00)>>5|
                     (color&0x000000F8)>>3);
(b)R5G6B5转RGB 0xFF000000|
             ((color&0xF800)<<8)|
             ((color&0x07E0)<<5)|
             ((color&0x001F)<<3);
2.5)采用1个字节表达分类编辑属性,其中选择标记1bit、删除标记1bit、分类属性5bit,预留1bit,分类表达范围0-31;
其中,坐标、强度、时间属性无需编辑,物理上采用FastAC压缩存储,颜色属性着色时动态生成,并且压缩存储,唯独分类编辑属性不压缩;
3)激光点云快速加载显示方法
a.打开读取点云数据文件,从文件头获取点云范围、数据量等信息,根据屏幕窗口大小及点云范围计算初始全图显示比例尺;
b.获取显示比例尺判断数据加载等级Level;
c.根据窗口范围和三级八叉树分块索引信息,加载分块数据;
d.解压分块数据,并根据坐标编码规则,解算坐标,通知主线程刷新显示分块点云;
e.对当前窗口分块数据进行动态调度管理,内存中保留当前窗口分块点云,实时卸载跨层及窗口范围外分块数据;
4)激光点云编辑分类实现方法
编辑选择点云时将HdPtClass.isSelected属性标记为1,编辑删除点云时将HdPtClass.isDeleted属性标记为1,分类时将HdPtClass.ptclass属性标记为对应类别编码。
与现有技术相比,本发明的有益效果是:基于本方法实现的三维激光点云存储、浏览、编辑时具有以下优点:点云存储空间小,点云浏览处理软件占用内存小,在浏览过程中点云占用内存不超过100Mb;首次加载显示及实时浏览速度快,在浏览点云过程中,不会卡顿现象,实时动态刷新点云,刷新响应时间小于1s;支持编辑分类。
附图说明
图1为点云数据总体结构示意图。
图2为块集、块、包组织关系示意图。
图3为点云快速加载显示流程图。
图4为点云分层分块压缩存储流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1~4,三维激光扫描仪采集获取的点云通常包含三维坐标、强度、时间和颜色等信息,后处理还包括分类信息。点云浏览显示通常包含按高程渲染、按强度渲染、按真彩色渲染、按分类渲染,不同渲染方式对应加载点云属性信息不同;本发明实施例中,一种海量三维激光点云压缩存储及快速加载显示方法,包括以下内容:
1)激光点云高度压缩存储具体实现
(a)根据原始点云范围和点云数据量确定起始划分网格位置,并按每数据包不超过1MB原则,确定划分包数,按每层的块集、块、包三级平衡划分原则,确定块集到块、块到包的格网划分原则。
(b)根据块集划分规则,读取原始数据写入块集数据,如果当前块集数据小于1MB,则直接将当前块集数据压缩写入当前块集,不再划分块;
(c)如果块集数据大于1MB,则根据块集到块的划分规则,将块集数据划分为指定块,如果写入的块数据小于1MB,则直接将当前块数据压缩写入当前块,不再划分包;
(d)如果块数据大于1MB,则根据块到包的划分规则,将块集数据划分为指定包,然后将包数据压缩写入当前包。至此完成0层数据的写入;
(e)根据0层综合抽稀生成1层数据,每次往高一等级抽稀时按1/4比例抽稀;
(f)不断循环抽稀生成高等级数据,直到最高等级数据小于100万;
(g)当为点云着色时,动态生成每数据包的对应颜色数据,并压缩存储,然后更新数据包的颜色数据地址,如果存在颜色数据,则删除重新生成并更新颜色数据地址;
(h)当为点云编辑分类时,动态生成每数据包的对应分类数据,然后更新数据包的分类数据地址;
2)激光点云快速加载显示具体实现
(a)开启OpenGL三维视图,开启GPU硬件加速;
(b)读取获取点云数据范围,设置相机视点、视角等参数;
(c)开启数据读取多线程,线程等待数据读取信号;
(d)读取最高级别层点云数据,按分块建立当前窗口内存MeshBuffer,首次浏览全图俯视显示;
(e)当进行鼠标拖动平移浏览时,将内存中的分块MeshBuffer实时刷新显示。当鼠标弹起松开后,根据窗口地理范围,计算需要新加载的分块,并发送信号给数据读取多线程,同时在内存中卸载窗口范围外分块MeshBuffer;
(f)当多线程每一个分块数据读取解析完毕时,通知主线程进行刷新视图;
(g)当视图进行放大缩小操作时,根据视图比例尺,判读对应显示的数据层级,然后在读取范围分块数据进行显示。
基于本方法,在浏览点云过程中,不会卡顿现象,实时动态刷新点云,在浏览过程中点云占用内存不超过100Mb,刷新响应时间小于2s。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。

Claims (1)

1.一种海量三维激光点云压缩存储及快速加载显示方法,其特征在于,包括以下内容:
1)激光点云高度压缩存储方法
1.1)点云数据总体结构:
a.整体由若干层级别构成,从低到高,0级是最详细原始数据,高等级由低等级数据抽稀构成,每一级别对应不同显示比例尺;
b.坐标、强度、时间合并文件存储,颜色、分类单独文件存储;
c.数据文件按级别顺序存储,0,1,2...n;
d.每一级别按坐标、强度、时间顺序存储;
e.颜色数据,通过与点云配准后的影像对点云进行着色后生成;
f.除分类数据无压缩外,其他数据都压缩存储,压缩算法采用FastAC;
g.选择、删除编辑标记信息存放在分类中;
1.2)块集、块、包组织关系:
a.每一级内部包括三级索引:块集、块、包,块集包含块,块包含包;
b.编号从左到右,从下到上;
c.分块包括2×2、4×4、8×8…,默认块集边长64,1块集包含4×4块,1块包含4×4包;
d.每层三级索引分块需均衡,每数据包最大存储1M数据,每个数据包索引包含格网划分规则、点数、范围、坐标、强度、时间、颜色、分类等数据地址信息;
2)数据压缩存储实现方法
2.1)xyz坐标采样网格相对坐标,每个坐标值采用2个字节表达,在保证坐标精度为0.001m情况下,网格边长约设为64m;
2.2)强度属性采用1个字节,表达范围0-255;
2.3)时间属性采用2个字节,精度为1/128秒情况下,表达范围0-512秒;
2.4)颜色属性采用2个字节,红色5bit、绿色6bit、蓝色5bit,即R5G6B5,RGB和R5G6B5转换方法如下:
(a)RGB转R5G6B5(u16)((color&0x00F80000)>>8|
                      (color&0x0000FC00)>>5|
                      (color&0x000000F8)>>3);
(b)R5G6B5转RGB 0xFF000000|
                 ((color&0xF800)<<8)|
                    ((color&0x07E0)<<5)|
                    ((color&0x001F)<<3);
2.5)采用1个字节表达分类编辑属性,其中选择标记1bit、删除标记1bit、分类属性5bit,预留1bit,分类表达范围0-31;
其中,坐标、强度、时间属性无需编辑,物理上采用FastAC压缩存储,颜色属性着色时动态生成,并且压缩存储,唯独分类编辑属性不压缩;
3)激光点云快速加载显示具体实现
a.开启OpenGL三维视图,开启GPU硬件加速;
b.读取获取点云数据范围,设置相机视点、视角等参数;
c.开启数据读取多线程,线程等待数据读取信号;
d.读取最高级别层点云数据,按分块建立当前窗口内存MeshBuffer,首次浏览全图俯视显示;
e.当进行鼠标拖动平移浏览时,将内存中的分块MeshBuffer实时刷新显示,当鼠标弹起松开后,根据窗口地理范围,计算需要新加载的分块,并发送信号给数据读取多线程,同时在内存中卸载窗口范围外分块MeshBuffer;
f.当多线程每一个分块数据读取解析完毕时,通知主线程进行刷新视图;
g.当视图进行放大缩小操作时,根据视图比例尺,判读对应显示的数据层级,然后在读取范围分块数据进行显示。
CN201510178971.XA 2015-04-16 2015-04-16 一种海量三维激光点云压缩存储及快速加载显示方法 Active CN104750854B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510178971.XA CN104750854B (zh) 2015-04-16 2015-04-16 一种海量三维激光点云压缩存储及快速加载显示方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510178971.XA CN104750854B (zh) 2015-04-16 2015-04-16 一种海量三维激光点云压缩存储及快速加载显示方法

Publications (2)

Publication Number Publication Date
CN104750854A true CN104750854A (zh) 2015-07-01
CN104750854B CN104750854B (zh) 2018-03-02

Family

ID=53590538

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510178971.XA Active CN104750854B (zh) 2015-04-16 2015-04-16 一种海量三维激光点云压缩存储及快速加载显示方法

Country Status (1)

Country Link
CN (1) CN104750854B (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106844610A (zh) * 2017-01-18 2017-06-13 上海交通大学 一种分布式结构化三维点云图像处理方法及系统
CN107346215A (zh) * 2017-05-10 2017-11-14 国网浙江省电力公司宁波供电公司 点云显示方法和设备
CN107886561A (zh) * 2017-11-30 2018-04-06 公安部物证鉴定中心 基于犯罪现场点云数据的存储调度与显示方法
CN109410118A (zh) * 2018-09-12 2019-03-01 广州南方卫星导航仪器有限公司 一种用于海量的点云数据快速加载显示方法
CN109444878A (zh) * 2018-09-28 2019-03-08 中国科学院电子学研究所苏州研究院 一种基于sar影像的微波三维数据集构建方法
CN110363822A (zh) * 2018-04-11 2019-10-22 上海交通大学 一种3d点云压缩方法
CN110413808A (zh) * 2019-07-17 2019-11-05 深兰科技(上海)有限公司 一种三维点云数据存储、查询方法及装置
WO2020107137A1 (en) * 2018-11-26 2020-06-04 Beijing Didi Infinity Technology And Development Co., Ltd. Systems and methods for point cloud rendering using video memory pool
CN111352861A (zh) * 2020-02-19 2020-06-30 Oppo广东移动通信有限公司 内存压缩方法、装置及电子设备
WO2020154975A1 (zh) * 2019-01-30 2020-08-06 深圳市大疆创新科技有限公司 用于点云数据的存储方法、回放方法及计算机可读介质
WO2021046817A1 (zh) * 2019-09-12 2021-03-18 深圳市大疆创新科技有限公司 点云数据的编解码方法、系统和存储介质
CN113297200A (zh) * 2020-06-24 2021-08-24 阿里巴巴集团控股有限公司 数据存储和获取方法、装置、电子设备及可读存储介质
CN115695393A (zh) * 2022-12-28 2023-02-03 山东矩阵软件工程股份有限公司 一种雷达点云数据的格式转换方法、系统及存储介质
CN116095181A (zh) * 2022-12-30 2023-05-09 天翼云科技有限公司 一种基于对象存储的点云压缩存储方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200506A (zh) * 2014-08-04 2014-12-10 广东威创视讯科技股份有限公司 三维gis海量矢量数据渲染方法及装置
CN104391906A (zh) * 2014-11-18 2015-03-04 武汉海达数云技术有限公司 车载海量点云数据动态浏览方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200506A (zh) * 2014-08-04 2014-12-10 广东威创视讯科技股份有限公司 三维gis海量矢量数据渲染方法及装置
CN104391906A (zh) * 2014-11-18 2015-03-04 武汉海达数云技术有限公司 车载海量点云数据动态浏览方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
郝刚等: ""一种三维激光扫描点云数据的预处理方法"", 《测绘科学》 *
黄承亮等: ""三维激光扫描点云数据压缩方法"", 《测绘科学》 *

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106844610A (zh) * 2017-01-18 2017-06-13 上海交通大学 一种分布式结构化三维点云图像处理方法及系统
CN106844610B (zh) * 2017-01-18 2020-03-24 上海交通大学 一种分布式结构化三维点云图像处理方法及系统
CN107346215A (zh) * 2017-05-10 2017-11-14 国网浙江省电力公司宁波供电公司 点云显示方法和设备
CN107886561A (zh) * 2017-11-30 2018-04-06 公安部物证鉴定中心 基于犯罪现场点云数据的存储调度与显示方法
CN107886561B (zh) * 2017-11-30 2021-11-09 公安部物证鉴定中心 基于犯罪现场点云数据的存储调度与显示方法
CN110363822A (zh) * 2018-04-11 2019-10-22 上海交通大学 一种3d点云压缩方法
CN109410118A (zh) * 2018-09-12 2019-03-01 广州南方卫星导航仪器有限公司 一种用于海量的点云数据快速加载显示方法
CN109444878A (zh) * 2018-09-28 2019-03-08 中国科学院电子学研究所苏州研究院 一种基于sar影像的微波三维数据集构建方法
CN112005275B (zh) * 2018-11-26 2023-04-21 北京嘀嘀无限科技发展有限公司 使用视频内存池进行点云渲染的系统和方法
CN112005275A (zh) * 2018-11-26 2020-11-27 北京嘀嘀无限科技发展有限公司 使用视频内存池进行点云渲染的系统和方法
WO2020107137A1 (en) * 2018-11-26 2020-06-04 Beijing Didi Infinity Technology And Development Co., Ltd. Systems and methods for point cloud rendering using video memory pool
WO2020154975A1 (zh) * 2019-01-30 2020-08-06 深圳市大疆创新科技有限公司 用于点云数据的存储方法、回放方法及计算机可读介质
CN110413808B (zh) * 2019-07-17 2021-09-24 深兰科技(上海)有限公司 一种三维点云数据存储、查询方法及装置
CN110413808A (zh) * 2019-07-17 2019-11-05 深兰科技(上海)有限公司 一种三维点云数据存储、查询方法及装置
WO2021046817A1 (zh) * 2019-09-12 2021-03-18 深圳市大疆创新科技有限公司 点云数据的编解码方法、系统和存储介质
CN113906681A (zh) * 2019-09-12 2022-01-07 深圳市大疆创新科技有限公司 点云数据的编解码方法、系统和存储介质
CN113906681B (zh) * 2019-09-12 2022-10-18 深圳市大疆创新科技有限公司 点云数据的编解码方法、系统和存储介质
CN111352861A (zh) * 2020-02-19 2020-06-30 Oppo广东移动通信有限公司 内存压缩方法、装置及电子设备
CN111352861B (zh) * 2020-02-19 2023-09-29 Oppo广东移动通信有限公司 内存压缩方法、装置及电子设备
CN113297200A (zh) * 2020-06-24 2021-08-24 阿里巴巴集团控股有限公司 数据存储和获取方法、装置、电子设备及可读存储介质
CN115695393A (zh) * 2022-12-28 2023-02-03 山东矩阵软件工程股份有限公司 一种雷达点云数据的格式转换方法、系统及存储介质
CN115695393B (zh) * 2022-12-28 2023-03-21 山东矩阵软件工程股份有限公司 一种雷达点云数据的格式转换方法、系统及存储介质
CN116095181A (zh) * 2022-12-30 2023-05-09 天翼云科技有限公司 一种基于对象存储的点云压缩存储方法及装置
CN116095181B (zh) * 2022-12-30 2024-06-07 天翼云科技有限公司 一种基于对象存储的点云压缩存储方法及装置

Also Published As

Publication number Publication date
CN104750854B (zh) 2018-03-02

Similar Documents

Publication Publication Date Title
CN104750854A (zh) 一种海量三维激光点云压缩存储及快速加载显示方法
CN103093499B (zh) 一种适用于网络传输的城市三维模型数据组织方法
US8194071B2 (en) Tile based graphics rendering
US8570322B2 (en) Method, system, and computer program product for efficient ray tracing of micropolygon geometry
CN102890829B (zh) 一种基于gpu的渲染地形的方法
CN109979008A (zh) 一种基于属性的点云条带划分方法
CN110706341A (zh) 一种城市信息模型的高性能渲染方法、装置及存储介质
CN105261066A (zh) 一种三维地理信息系统实时绘制多线程分配与控制方法
CN104391906A (zh) 车载海量点云数据动态浏览方法
CN106023297A (zh) 一种精细三维模型的纹理动态组织方法
CN116127115B (zh) 面向时序更新的倾斜摄影三维模型瓦片集可视化方法
CN108765538A (zh) 基于cad平台的osgb数据分级渲染的方法
CN112785696A (zh) 基于游戏引擎和倾斜摄影数据的三维实景模型生成方法
CN115630588A (zh) 一种新风流场仿真和动态可视化系统
CN103136782B (zh) 一种三维模型地图动态渲染方法及装置
Bao et al. Realistic real-time rendering for large-scale forest scenes
CN103927395A (zh) 结构简化图像三维几何及材质纹理的数据结构及渲染方法
Andújar et al. Inexpensive reconstruction and rendering of realistic roadside landscapes
CN112843700A (zh) 地形图像的生成方法、装置、计算机设备和存储介质
CN104601950B (zh) 一种视频监控方法
CN116758206A (zh) 矢量数据的融合渲染方法、装置、计算机设备、存储介质
CN108269301B (zh) 一种三维气象雷达显示方法
CN114037792B (zh) 一种Unreal Engine平台下的百万量级地下管线模型动态生成方法
CN115964455A (zh) 一种基于MapVGL开源GIS引擎的数据可视化处理方法
CN109712237B (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
GR01 Patent grant
GR01 Patent grant