CN116051764A - 一种遥感影像地图绘制方法及系统 - Google Patents
一种遥感影像地图绘制方法及系统 Download PDFInfo
- Publication number
- CN116051764A CN116051764A CN202211721371.XA CN202211721371A CN116051764A CN 116051764 A CN116051764 A CN 116051764A CN 202211721371 A CN202211721371 A CN 202211721371A CN 116051764 A CN116051764 A CN 116051764A
- Authority
- CN
- China
- Prior art keywords
- data
- image
- map
- resolution
- determining
- 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 77
- 230000001965 increasing effect Effects 0.000 claims abstract description 23
- 238000004458 analytical method Methods 0.000 claims abstract description 11
- 238000013507 mapping Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 4
- 102100025444 Gamma-butyrobetaine dioxygenase Human genes 0.000 description 17
- 101000934612 Homo sapiens Gamma-butyrobetaine dioxygenase Proteins 0.000 description 17
- 238000009877 rendering Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 238000013500 data storage Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012952 Resampling Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000012634 fragment Substances 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 241000274965 Cyrestis thyodamas Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/08—Projecting images onto non-planar surfaces, e.g. geodetic screens
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明公开一种遥感影像地图绘制方法及系统,涉及图像处理领域。所述方法包括获取目标层级下地图图像的若干个局部区域;对每个局部区域进行参数解析,得到多个目标层级下的区域参数;确定图像金字塔数据;图像金字塔数据为整体数据或瓦片数据;整体数据包括由下至上层级依次增加、分辨率依次降低的多个地图图像;瓦片数据包括地图图像的多个局部金字塔数据;每个局部金字塔数据包括由下至上层级依次增加、分辨率依次降低的多个局部图像;根据区域参数、图像金字塔数据和调色板进行遥感地图绘制,得到目标层级下的遥感地图;根据区域参数和图像金字塔数据进行图像地图绘制,得到目标层级下的图像地图。本发明能够提高地图的绘制速度。
Description
技术领域
本发明涉及图像处理技术领域,特别是涉及一种遥感影像地图绘制方法及系统。
背景技术
在常见的地图绘制方法中,对于地图图像数据往往是采用分块处理的形式,将整体图像划分为若干个地图块,再根据不同分辨率对每个地图块进行处理及拼接,如早期的谷歌地图采用的方式是提前把每个比例尺下的地图绘制成小块区域,保存在服务器的缓存目录中,在客户端请求地图时,直接获取需要的多个小块区域,拼接成整幅地图,以此来提高访问速度,但这种方式需要将不同分辨率的地图数据文件单独存储,对地图碎片的存储需求很大,而且在根据请求进行地图碎片提取时,需要对不同文件中的数据逐一提取,绘制速度慢,进而导致整体地图的绘制速度不高。
发明内容
本发明的目的是提供一种遥感影像地图绘制方法及系统,能够提高地图的绘制速度。
为实现上述目的,本发明提供了如下方案:
一种遥感影像地图绘制方法,包括:
获取目标层级下地图图像的若干个局部区域;
对每个所述局部区域进行参数解析,得到多个目标层级下的区域参数;所述区域参数包括区域经纬度范围、区域宽度和区域高度;
确定图像金字塔数据;所述图像金字塔数据为整体数据或瓦片数据;所述整体数据包括由下至上层级依次增加、分辨率依次降低的多个地图图像;所述瓦片数据包括地图图像的多个局部金字塔数据;每个所述局部金字塔数据包括由下至上层级依次增加、分辨率依次降低的多个局部图像;
根据所述区域参数、所述图像金字塔数据和调色板进行遥感地图绘制,得到目标层级下的遥感地图;
根据所述区域参数和所述图像金字塔数据进行图像地图绘制,得到目标层级下的图像地图。
可选地,所述整体数据和所述瓦片数据的确定方法为:
获取原始图像;所述原始图像为分辨率最高的地图图像;
根据所述原始图像构建分辨率依次降低的多个地图图像;
按照由下至上分辨率依次降低的排布方式,使得多个所述地图图像由下至上层级依次增加,得到整体数据;
对所述原始图像进行切片处理,得到多个瓦片区域;
根据每个所述瓦片区域构建分辨率依次降低的多个局部图像;
按照由下至上分辨率依次降低的排布方式,使得每个瓦片区域的多个所述局部图像由下至上层级依次增加,得到瓦片数据。
可选地,所述根据所述区域参数、所述图像金字塔数据和调色板进行遥感影像地图绘制,得到目标层级下的遥感地图,具体包括:
构建多维空数组;
根据所述区域参数提取所述图像金字塔数据中的对应数据;
将每个所述对应数据填充至所述多维空数组中,得到填充数组;
利用所述调色板对所述填充数组进行绘制,得到目标层级下的遥感地图。
可选地,所述根据所述区域参数提取所述图像金字塔数据中的对应数据,具体包括:
利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置;
提取所述投影位置处对应的数据,得到对应数据。
可选地,若所述图像金字塔数据为整体数据,则所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置,具体包括:
根据所述区域参数确定所述局部区域的分辨率和尺寸;
利用等经纬投影的方法,根据所述局部区域的分辨率和尺寸确定所述图像金字塔数据的投影位置。
可选地,若所述图像金字塔数据为整体数据,则所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置,具体包括:
根据所述区域参数确定所述局部区域的占图比值和尺寸;
利用标称投影的方法或极射赤面投影的方法,根据所述局部区域的占图比值和尺寸确定所述图像金字塔的投影位置。
可选地,若所述图像金字塔数据为瓦片数据,则所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置,具体包括:
根据所述区域参数确定所述局部区域的分辨率;
利用等经纬投影的方法,根据所述局部区域的分辨率、所述区域经纬度范围和预设经纬度跨度值确定所述图像金字塔数据的投影位置。
可选地,若所述图像金字塔数据为瓦片数据,则所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置,具体包括:
根据所述区域参数确定所述局部区域的分辨率;
利用标称投影的方法或极射赤面投影的方法,根据所述局部区域的分辨率、所述区域经纬度范围和预设区域尺寸确定所述图像金字塔数据的投影位置。
本发明还提出了一种遥感影像地图绘制系统,包括:
数据采集单元,用于获取目标层级下地图图像的若干个局部区域;
参数解析单元,用于对每个所述局部区域进行参数解析,得到多个目标层级下的区域参数;所述区域参数包括区域经纬度范围、区域宽度和区域高度;
数据构建单元,用于确定图像金字塔数据;所述图像金字塔数据为整体数据或瓦片数据;所述整体数据包括由下至上层级依次增加、分辨率依次降低的多个地图图像;所述瓦片数据包括地图图像的多个局部金字塔数据;每个所述局部金字塔数据包括由下至上层级依次增加、分辨率依次降低的多个局部图像;
遥感地图绘制单元,用于根据所述区域参数、所述图像金字塔数据和调色板进行遥感地图绘制,得到目标层级下的遥感地图;
图像地图绘制单元,用于根据所述区域参数和所述图像金字塔数据进行图像地图绘制,得到目标层级下的图像地图。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明公开了一种遥感影像地图绘制方法及系统,所述方法包括获取目标层级下地图图像的若干个局部区域;对每个所述局部区域进行参数解析,得到多个目标层级下的区域参数;确定图像金字塔数据;根据所述区域参数、所述图像金字塔数据和调色板进行遥感地图绘制,得到目标层级下的遥感地图,以及根据所述区域参数和所述图像金字塔数据进行图像地图绘制,得到目标层级下的图像地图。其中,所述图像金字塔数据为整体数据或瓦片数据;所述整体数据包括由下至上层级依次增加、分辨率依次降低的多个地图图像;所述瓦片数据包括地图图像的多个局部金字塔数据;每个所述局部金字塔数据包括由下至上层级依次增加、分辨率依次降低的多个局部图像。本发明通过构建图像金字塔数据,使多个不同分辨率的地图图像数据能够保存在一个数据文件中,能够更加节省数据存储空间,并且能够同时根据多个解析后的区域参数提取目标层级下的对应数据,大大加快了遥感地图和图像地图的绘制速度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明遥感影像地图绘制方法的流程示意图;
图2为本发明遥感影像地图绘制方法的具体应用示意图;
图3为本发明遥感影像地图绘制方法的地图切片层级结构示意图;
图4为本发明遥感影像地图绘制方法的图像金字塔数据结构示意图;
图5为本发明遥感影像地图绘制系统的结构框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种遥感影像地图绘制方法及系统,能够提高地图的绘制速度。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,本发明提供了一种遥感影像地图绘制方法,包括:
步骤100:获取目标层级下地图图像的若干个局部区域。
步骤200:对每个所述局部区域进行参数解析,得到多个目标层级下的区域参数;所述区域参数包括区域经纬度范围、区域宽度和区域高度。
步骤300:确定图像金字塔数据;所述图像金字塔数据为整体数据或瓦片数据;所述整体数据包括由下至上层级依次增加、分辨率依次降低的多个地图图像;所述瓦片数据包括地图图像的多个局部金字塔数据;每个所述局部金字塔数据包括由下至上层级依次增加、分辨率依次降低的多个局部图像。
步骤400:根据所述区域参数、所述图像金字塔数据和调色板进行遥感地图绘制,得到目标层级下的遥感地图。
步骤500:根据所述区域参数和所述图像金字塔数据进行图像地图绘制,得到目标层级下的图像地图。
作为步骤300的一种具体实施方式,所述整体数据和所述瓦片数据的确定方法为:
第一步,获取原始图像;所述原始图像为分辨率最高的地图图像。
第二步,根据所述原始图像构建分辨率依次降低的多个地图图像。
第三步,按照由下至上分辨率依次降低的排布方式,使得多个所述地图图像由下至上层级依次增加,得到整体数据。
第四步,对所述原始图像进行切片处理,得到多个瓦片区域。
第五步,根据每个所述瓦片区域构建分辨率依次降低的多个局部图像。
第六步,按照由下至上分辨率依次降低的排布方式,使得每个瓦片区域的多个所述局部图像由下至上层级依次增加,得到瓦片数据。
作为步骤400的一种具体实施方式,所述根据所述区域参数、所述图像金字塔数据和调色板进行遥感地图绘制,得到目标层级下的遥感地图,具体包括:
第一步,构建多维空数组。
第二步,根据所述区域参数提取所述图像金字塔数据中的对应数据。
本步的具体实施方式包括:利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置;提取所述投影位置处对应的数据,得到对应数据。
第三步,将每个所述对应数据填充至所述多维空数组中,得到填充数组。
第四步,利用所述调色板对所述填充数组进行绘制,得到目标层级下的遥感地图。
对于上述步骤400中第二步的“所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置”,提供第一种具体的实施方式,包括:
若所述图像金字塔数据为整体数据,则:
根据所述区域参数确定所述局部区域的分辨率和尺寸;利用等经纬投影的方法,根据所述局部区域的分辨率和尺寸确定所述图像金字塔数据的投影位置。
对于上述步骤400中第二步的“所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置”,提供第二种具体的实施方式,包括:
若所述图像金字塔数据为整体数据,则:
根据所述区域参数确定所述局部区域的占图比值和尺寸;利用标称投影的方法或极射赤面投影的方法,根据所述局部区域的占图比值和尺寸确定所述图像金字塔数据的投影位置。
对于上述步骤400中第二步的“所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置”,提供第三种具体的实施方式,包括:
若所述图像金字塔数据为瓦片数据,则:
根据所述区域参数确定所述局部区域的分辨率;利用等经纬投影的方法,根据所述局部区域的分辨率、所述区域经纬度范围和预设经纬度跨度值确定所述图像金字塔数据的投影位置。
对于上述步骤400中第二步的“所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置”,提供第四种具体的实施方式,包括:
若所述图像金字塔数据为瓦片数据,则:
根据所述区域参数确定所述局部区域的分辨率;利用标称投影的方法或极射赤面投影的方法,根据所述局部区域的分辨率、所述区域经纬度范围和预设区域尺寸确定所述图像金字塔数据的投影位置。
在上述方法的基础上,提出一种具体的实施例,具体处理过程如下:
本实施例为前后端分离的地图发布系统,在地图进行绘制过程中,需要前后端进行交互,前端全球范按照不同分辨率分为多个层级,每个层级划分为若干个局部区域,并根据目标层级的局部区域提出请求,经过对请求进行解析后,在后端进行绘制并传输展示。其主要分为三个板块,包括后端数据处理与入库,地图发布服务,以及展示服务。
首先是后端数据处理与入库服务,包括两个步骤,生成金字塔图像文件,并入库存档,本文生成金字塔文件的方式选择基本python版本的GDAL库进行实现,生成步骤分为三步:(1)读取图像文件,获取原始地图的数据投影信息;(2)生成金字塔文件;(3)将原始文件的投影赋给对金字塔文件,或者按需定义新的投影。
在图像金字塔文件进行构建过程中,文件仅支持tif、png、jpg等图像格式,如果原始地图的遥感影像为NC或者HDF格式,则需要从该文件中提取对应数据生成图像格式,再生成对应的图像金字塔文件。在对原始地图进行投影采样时,采样的倍率关系如表1所示,采用缩采样,采样方式默认为最近邻插值,连续数据采用双线性插值法或双三次卷积法。第一层的数据为原始文件,随着层级的增加,图像的行列大小均按2倍数减少;图像金字塔文件生成成功后自动存储为图像文件名加.ovr后缀的文件。
表1图像金字塔生成参数说明
采样倍率 | 对应层级 | 图像大小示例 |
1 | 0 | 512*512 |
2 | 1 | 256*256 |
4 | 2 | 128*128 |
8 | 3 | 64*64 |
16 | 4 | 32*32 |
其次是地图发布服务,该部分支持多种遥感地图数据的展示,为了便于管理和查看,需根据不同产品(图像或遥感)、不同需求(服务方式和投影选择)进行后端配置,在本实施例中支持两种地图服务,UTSS(单一瓦片地图)和UTTS(预定义瓦片)地图发布服务,支持图像发布,也支持遥感产品发布。图像地图发布不需要调色板,遥感产品发布需要调色板实时渲染,也支持前端更换调色板。
最后是展示服务,地图展示服务也是前后端分离的服务,前端基于地图框架(openlayers)调用瓦片,调用参数包含层级,瓦片范围(左上角的坐标和右下角的坐标值),以及图像宽高(width,height),瓦片范围计算依据以左上角为起点。瓦片层级分了0到9级,一共十层。层级对应的图像分辨率如表2所示。
表2瓦片分级描述
下面阐述当客户端发出请求时,如何实现实时渲染地图。当前端选择某一个遥感产品展示,前端基于Openlayers框架将多个小块区域向后端发出请求。当后端接收到一个请求,先基于该请求查找远程字典服务数据库(redis)中是否有记录,该数据库用于缓存,如果有记录,从数据库中直接调用结果,并刷新缓存时间;若数据库中没有记录,则调用后台数据处理程序,实时渲染区域。下面根据不同的请求进行算法实现说明。
在本实施例中支持两种图像展示服务方式:UTSS和UTTS,这两种服务方式均支持等经纬投影、极射赤面投影和标称投影。产品绘图是根据前端传入的调色板、以及请求参数,后端服务计算索引位置提取目标区域的原始数据,再根据调色板,实时绘制图像,再经过数据转换和处理传给前端展示;影像绘图的原始数据是一张已经着色好的图像文件,后端服务根据请求参数,计算索引位置,获取目标区域的图像数据,经过数据提取和处理后传给前端展示。
当前端发出请求后对请求的参数进行解析,解析后的请求参数包括图像参数和其他参数。其中,请求参数(URL)如表3所示。
表3图像请求URL参数说明
字段 | 涵义 |
SERVICE | 服务方式 |
VERSION | 系统版本号 |
FORMAT | 请求图片区域返回格式 |
PRODUCT | 产品路径 |
WIDTH | 请求区域的宽度 |
HEIGHT | 请求区域的高度 |
CRS | 投影方式 |
BBOX | 请求区域的范围 |
LEVEL | 请求区域的层级 |
COLORBAR | 调色板 |
如图2所示,当后端收到URL请求参数后,首先判断处理模式,其次判断是否有调色板,最后判断投影方式,根据不同的请求调用不同的处理函数,最后生成相应的产品。
URL请求的中的PRODUCT参数传递的是图像金字塔文件;WIDTH和HEIGHT是请求返回的图像大小;BBOX参数传递的是请求的图像范围的实际经纬度坐标,BBOX参数传递的顺序为(最小经度,最小纬度,最大经度,最大纬度),可结合WIDTH、HEIGHT计算当前图层的分辨率,该分辨率对应的层级是LEVEL(表3所示),也是该URL请求得到的图像的分辨率。
计算公式如下:
横向分辨率=(最大经度-最小经度)/WIDTH;
纵向分辨率=(最大纬度-最小纬度)/HEIGHT;
在本实施例中,主要采用等经纬度投影,即横向分辨率等于纵向分辨率。
作为UTSS服务方式的一种具体实施方式,主要包括如下两种:
第一种为基于上述步骤400中第二步的“所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置”的第一种具体的实施方式。
第一步:根据参数BBOX、参数WIDTH和参数HEIGHT计算出前端服务请求的图像分辨率。
第二步:获取金字塔图像文件的金字塔层数(波段数)、原始分辨率、经纬度范围等信息。
第三步:提取金字塔图像文件和参数BBOX的交集,即有效图像范围,如果请求的参数BBOX的范围超出金字塔文件的范围,即返回一张空白图。否则,进入第四步。
第四步:根据URL参数中的WIDTH、HEIGHT,以及请求图像格式,创建一个多维空数组。
第五步:根据第三步得到的图像交集范围、第二步获得的金字塔文件的经纬度范围以及第一步得到的图像分辨率,计算交集范围在金字塔文件的起始行列位置,以及宽高。
第六步:与上同理,根据交集范围、参数BBOX以及分辨率计算交集在输出图像的起始行列号位置,以及宽高。
第七步:根据图像金字塔文件的层数、分辨率以及URL请求的分辨率,计算该请求的原始数据位于图像金字塔文件的第几个层级。
第八步:根据第七步获得层级数,第五步获得的金字塔文件起始行列位置以及宽高,提取图像数据,将其重采样为第六步的图像宽高大小。如果URL请求的是图像文件,则直接将其转为base64编码返给前端;如果URL请求的是产品文件,则根据调色板,将提取的图像数据绘制成图像,再转为base64编码返给前端。
第二种为基于上述步骤400中第二步的“所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置”的第二种具体的实施方式。
标称投影和极射赤面投影的参数BBOX的数值是以米为单位,固定视角下的视图范围是固定的。处理方式如下:
第一步:根据参数BBOX和固定视图范围,计算当前参数BBOX所占的视图比例;
第二步:获取图像金字塔文件的层级数,数据的宽高。并且根据URL请求层级进行推算对应金字塔文件的层级;
第三步:根据第二步,获取对应层级的图像金字塔文件数据,再根据第一步的比例关系,计算参数BBOX请求的数据范围位于该数据的起始行列,以及宽高;
第四步:根据URL参数中的WIDTH、HEIGHT,以及请求图像格式,创建一个多维空数组;
第五步:根据第二步获得层级数,第三步获得的金字塔文件起始行列位置以及宽高,提取图像数据,将其重采样为URL请求的图像宽高大小。如果URL请求的是图像文件,则直接将其转为base64编码返给前端;如果URL请求的是产品文件,则根据调色板,将提取的图像数据绘制成图像,再转为base64编码返给前端。
作为UTTS服务方式的一种具体实施方式,主要包括如下两种:
第一种为基于上述步骤400中第二步的“所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置”的第三种具体的实施方式。
第一步:根据URL的参数计算请求图像的分辨率;
第二步:根据BBOX、约定的单块图像经纬度跨度值(也即预设经纬度跨度值为11.25°)和分块起算点坐标值,本实施例以左下角为起点。计算BBOX请求范围所对应的瓦片数据编号;
第三步:根据URL参数中的WIDTH、HEIGHT,以及请求图像格式,创建一个多维空数组;
第四步:遍历瓦片数据,获取瓦片数据的经纬度信息,与BBOX计算交集,后面的处理过程同UTSS一样,将瓦片的数据提取出来放进第三步的空数组里。
第五步:如果URL请求的是图像文件,则直接将其转为base64编码返给前端;如果URL请求的是产品文件,则根据调色板,将提取的图像数据绘制成图像,再转为base64编码返给前端。
第二种为基于上述步骤400中第二步的“所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置”的第四种具体的实施方式。
第一步:根据URL的参数计算请求图像的分辨率;
第二步:根据BBOX、预设区域尺寸,以及分块起算点坐标值,本实施例以左下角为起点。计算BBOX请求范围所对应的瓦片数据编号;
第三步:根据URL参数中的WIDTH、HEIGHT,以及请求图像格式,创建一个多维空数组;
第四步:遍历瓦片数据,获取瓦片地图的经纬度信息,与BBOX计算交集,后面的处理过程同UTSS一样,将瓦片的数据提取出来放进第三步的空数组里。
第五步:如果URL请求的是图像文件,则直接将其转为base64编码返给前端;如果URL请求的是产品文件,则根据调色板,将提取的图像数据绘制成图像,再转为base64编码返给前端。
在上述实施例中,UTSS服务和UTTS服务,分别对应Web地图服务(WMS)和Web地图瓦片服务(WMTS)。UTSS/WMS利用具有地理空间位置信息的数据制作地图;WMS针对提供可定制的服务,是一个动态数据或用户定制地图(需结合SLD标准)的理想解决办法。UTTS/WMTS采用预设图块的方法发布数字地图服务,对数据进行切片和缓存,UTTS/WMTS弥补了UTSS/WMS不能提供分块地图的不足,但牺牲了提供定制地图的灵活性,代之以通过提供静态数据(基础地图)来增强伸缩性,这些静态数据的范围框和比例尺被限定在单个图块内。以全球范围的地图展示为例,UTSS服务的每个瓦片请求的文件都是同一个图像文件,不限制瓦片请求的分辨率和范围,而UTTS服务请求的瓦片,需要经过计算来确定请求是某一编号的单块图像,单块图像的经纬度范围约束在11.25度,分辨率为250m。
在上述实施例的图像金字塔数据中,地图切片层级设计为10层,瓦片数据的大小设计为512*512。当第0层时,全球范围被拆分为两个512*512个大小的瓦片,每个瓦片的跨度为180度,瓦片分辨率为0.3515625度;当层级为1级时,全球范围被拆分为4个512*512大小的瓦片,每个瓦片的跨度为90度,瓦片分辨率为0.17578125,以此类推,每一级瓦片数量为2的2n次幂,n表示为当前层级。逐级的分辨率也以2倍的比例缩小,第0级展示的有效区域为全图,随着层级的增加,相同大小的视窗下展示的内容逐渐减少。其变化规律如图3所示。
与瓦片数据的层级相对应的是图像金字塔数据的整体数据,图像金字塔是图像中多尺度表达的一种,最主要用于图像的分割,是一种以多分辨率来解释图像的有效但概念简单的结构。如图4所示,金字塔的底部是高分辨率图像,而顶部是低分辨率图像,连续层级均以2:1的比例进行缩减采样,层级越高,图像越小,分辨率越低。图像金字塔文件包含多个缩采样图层,在数据读取过程中可以通过仅检索使用指定分辨率(取决于显示要求)的数据,可以加快栅格数据的获取速度。利用图像金字塔文件,可以在实时渲染过程中快速绘制较低分辨率的数据,而随着放大操作的进行,更精细的分辨率等级也将快速得到绘制,并且不需要将不同分辨率的数据文件单独存储,可以更加节省数据存储空间。基于nodejs框架,实现并发,加快了图像的实时渲染。
综上,本发明的有益效果是:
构建的图像金字塔数据文件中包含多个缩采样图层,在数据读取过程中可以通过仅检索使用指定分辨率(取决于显示要求)的数据,可以加快栅格数据的获取速度。利用图像金字塔文件,可以在实时渲染过程中快速绘制较低分辨率的数据,而随着放大操作的进行,更精细的分辨率等级也将快速得到绘制,并且不需要将不同分辨率的数据文件单独存储,可以更加节省数据存储空间。基于服务端web应用平台(nodejs)框架,实现并发,加快了图像的实时渲染,缓存采用后端缓存,减少了服务器的缓存压力。
如图5所示,本发明还提供了一种遥感影像地图绘制系统,包括:
数据采集单元,用于获取目标层级下地图图像的若干个局部区域;
参数解析单元,用于对每个所述局部区域进行参数解析,得到多个目标层级下的区域参数;所述区域参数包括区域经纬度范围、区域宽度和区域高度。
数据构建单元,用于确定图像金字塔数据;所述图像金字塔数据为整体数据或瓦片数据;所述整体数据包括由下至上层级依次增加、分辨率依次降低的多个地图图像;所述瓦片数据包括地图图像的多个局部金字塔数据;每个所述局部金字塔数据包括由下至上层级依次增加、分辨率依次降低的多个局部图像。
遥感地图绘制单元,用于根据所述区域参数、所述图像金字塔数据和调色板进行遥感地图绘制,得到目标层级下的遥感地图。
图像地图绘制单元,用于根据所述区域参数和所述图像金字塔数据进行图像地图绘制,得到目标层级下的图像地图。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (9)
1.一种遥感影像地图绘制方法,其特征在于,包括:
获取目标层级下地图图像的若干个局部区域;
对每个所述局部区域进行参数解析,得到多个目标层级下的区域参数;所述区域参数包括区域经纬度范围、区域宽度和区域高度;
确定图像金字塔数据;所述图像金字塔数据为整体数据或瓦片数据;所述整体数据包括由下至上层级依次增加、分辨率依次降低的多个地图图像;所述瓦片数据包括地图图像的多个局部金字塔数据;每个所述局部金字塔数据包括由下至上层级依次增加、分辨率依次降低的多个局部图像;
根据所述区域参数、所述图像金字塔数据和调色板进行遥感地图绘制,得到目标层级下的遥感地图;
根据所述区域参数和所述图像金字塔数据进行图像地图绘制,得到目标层级下的图像地图。
2.根据权利要求1所述的遥感影像地图绘制方法,其特征在于,所述整体数据和所述瓦片数据的确定方法为:
获取原始图像;所述原始图像为分辨率最高的地图图像;
根据所述原始图像构建分辨率依次降低的多个地图图像;
按照由下至上分辨率依次降低的排布方式,使得多个所述地图图像由下至上层级依次增加,得到整体数据;
对所述原始图像进行切片处理,得到多个瓦片区域;
根据每个所述瓦片区域构建分辨率依次降低的多个局部图像;
按照由下至上分辨率依次降低的排布方式,使得每个瓦片区域的多个所述局部图像由下至上层级依次增加,得到瓦片数据。
3.根据权利要求1所述的遥感影像地图绘制方法,其特征在于,所述根据所述区域参数、所述图像金字塔数据和调色板进行遥感影像地图绘制,得到目标层级下的遥感地图,具体包括:
构建多维空数组;
根据所述区域参数提取所述图像金字塔数据中的对应数据;
将每个所述对应数据填充至所述多维空数组中,得到填充数组;
利用所述调色板对所述填充数组进行绘制,得到目标层级下的遥感地图。
4.根据权利要求3所述的遥感影像地图绘制方法,其特征在于,所述根据所述区域参数提取所述图像金字塔数据中的对应数据,具体包括:
利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置;
提取所述投影位置处对应的数据,得到对应数据。
5.根据权利要求4所述的遥感影像地图绘制方法,其特征在于,若所述图像金字塔数据为整体数据,则所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置,具体包括:
根据所述区域参数确定所述局部区域的分辨率和尺寸;
利用等经纬投影的方法,根据所述局部区域的分辨率和尺寸确定所述图像金字塔数据的投影位置。
6.根据权利要求4所述的遥感影像地图绘制方法,其特征在于,若所述图像金字塔数据为整体数据,则所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置,具体包括:
根据所述区域参数确定所述局部区域的占图比值和尺寸;
利用标称投影的方法或极射赤面投影的方法,根据所述局部区域的占图比值和尺寸确定所述图像金字塔的投影位置。
7.根据权利要求4所述的遥感影像地图绘制方法,其特征在于,若所述图像金字塔数据为瓦片数据,则所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置,具体包括:
根据所述区域参数确定所述局部区域的分辨率;
利用等经纬投影的方法,根据所述局部区域的分辨率、所述区域经纬度范围和预设经纬度跨度值确定所述图像金字塔数据的投影位置。
8.根据权利要求4所述的遥感影像地图绘制方法,其特征在于,若所述图像金字塔数据为瓦片数据,则所述利用投影方法确定所述区域参数在所述图像金字塔数据中的投影位置,具体包括:
根据所述区域参数确定所述局部区域的分辨率;
利用标称投影的方法或极射赤面投影的方法,根据所述局部区域的分辨率、所述区域经纬度范围和预设区域尺寸确定所述图像金字塔数据的投影位置。
9.一种遥感影像地图绘制系统,其特征在于,包括:
数据采集单元,用于获取目标层级下地图图像的若干个局部区域;
参数解析单元,用于对每个所述局部区域进行参数解析,得到多个目标层级下的区域参数;所述区域参数包括区域经纬度范围、区域宽度和区域高度;
数据构建单元,用于确定图像金字塔数据;所述图像金字塔数据为整体数据或瓦片数据;所述整体数据包括由下至上层级依次增加、分辨率依次降低的多个地图图像;所述瓦片数据包括地图图像的多个局部金字塔数据;每个所述局部金字塔数据包括由下至上层级依次增加、分辨率依次降低的多个局部图像;
遥感地图绘制单元,用于根据所述区域参数、所述图像金字塔数据和调色板进行遥感地图绘制,得到目标层级下的遥感地图;
图像地图绘制单元,用于根据所述区域参数和所述图像金字塔数据进行图像地图绘制,得到目标层级下的图像地图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211721371.XA CN116051764B (zh) | 2022-12-30 | 2022-12-30 | 一种遥感影像地图绘制方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211721371.XA CN116051764B (zh) | 2022-12-30 | 2022-12-30 | 一种遥感影像地图绘制方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116051764A true CN116051764A (zh) | 2023-05-02 |
CN116051764B CN116051764B (zh) | 2024-01-30 |
Family
ID=86117812
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211721371.XA Active CN116051764B (zh) | 2022-12-30 | 2022-12-30 | 一种遥感影像地图绘制方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116051764B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339570A (zh) * | 2008-08-12 | 2009-01-07 | 北京航空航天大学 | 一种高效的海量遥感数据分布式组织管理方法 |
US20090317020A1 (en) * | 2008-06-24 | 2009-12-24 | Microsoft Corporation | Variable Resolution Images |
JP2011091769A (ja) * | 2009-10-26 | 2011-05-06 | Sony Computer Entertainment Inc | 画像ファイル生成装置、画像処理装置、画像ファイル生成方法、画像処理方法、および画像ファイルのデータ構造 |
WO2016131390A1 (zh) * | 2015-02-16 | 2016-08-25 | 高德软件有限公司 | 电子地图展示方法、装置及电子设备 |
CN108228746A (zh) * | 2017-12-20 | 2018-06-29 | 中国电子科技集团公司电子科学研究院 | 瓦片地图服务的调用方法、装置、存储介质及服务器 |
CN112966134A (zh) * | 2021-03-26 | 2021-06-15 | 北京亿耘科技有限公司 | 一种遥感影像数据处理方法、系统和存储介质 |
CN115170779A (zh) * | 2022-07-18 | 2022-10-11 | 生态环境部环境规划院 | 一种遥感影像渲染和网络发布方法 |
CN115470366A (zh) * | 2022-08-31 | 2022-12-13 | 湖南省第二测绘院 | 一种基于瓦片的遥感图像存储方法及系统 |
-
2022
- 2022-12-30 CN CN202211721371.XA patent/CN116051764B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090317020A1 (en) * | 2008-06-24 | 2009-12-24 | Microsoft Corporation | Variable Resolution Images |
CN101339570A (zh) * | 2008-08-12 | 2009-01-07 | 北京航空航天大学 | 一种高效的海量遥感数据分布式组织管理方法 |
JP2011091769A (ja) * | 2009-10-26 | 2011-05-06 | Sony Computer Entertainment Inc | 画像ファイル生成装置、画像処理装置、画像ファイル生成方法、画像処理方法、および画像ファイルのデータ構造 |
WO2016131390A1 (zh) * | 2015-02-16 | 2016-08-25 | 高德软件有限公司 | 电子地图展示方法、装置及电子设备 |
CN108228746A (zh) * | 2017-12-20 | 2018-06-29 | 中国电子科技集团公司电子科学研究院 | 瓦片地图服务的调用方法、装置、存储介质及服务器 |
CN112966134A (zh) * | 2021-03-26 | 2021-06-15 | 北京亿耘科技有限公司 | 一种遥感影像数据处理方法、系统和存储介质 |
CN115170779A (zh) * | 2022-07-18 | 2022-10-11 | 生态环境部环境规划院 | 一种遥感影像渲染和网络发布方法 |
CN115470366A (zh) * | 2022-08-31 | 2022-12-13 | 湖南省第二测绘院 | 一种基于瓦片的遥感图像存储方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN116051764B (zh) | 2024-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11874855B2 (en) | Parallel data access method and system for massive remote-sensing images | |
CN104134214B (zh) | 一种基于影像金字塔分层的数字切片图像存储及显示方法 | |
CN111723221A (zh) | 基于分布式架构的海量遥感数据处理方法及系统 | |
CN110263116B (zh) | 一种支持动态投影转换的瓦片地图发布方法及装置 | |
CN109885632B (zh) | 一种空间科学与应用数据检索方法、系统、介质及设备 | |
CN113032513B (zh) | 一种针对空间地理数据的切片方法及渲染方法 | |
CN109299298A (zh) | 多尺度影像模型的构建方法、装置、应用方法及系统 | |
CN111354084A (zh) | 一种基于三维模型瓦片的网络地理信息服务系统 | |
CN110597900B (zh) | 一种gdb数据实时按需生成矢量切片的方法 | |
CN112307146B (zh) | Gis路线数据展示系统及方法 | |
CN111949817A (zh) | 基于遥感影像的农作物信息展示系统、方法、设备及介质 | |
CN108228746A (zh) | 瓦片地图服务的调用方法、装置、存储介质及服务器 | |
CN112686997B (zh) | 基于WebGIS的三维模型数据分析展示平台及方法 | |
CN115391699B (zh) | 一种云优化的静态地理信息矢量数据在线服务方法 | |
EP2056217A1 (en) | Geographic XML database management system | |
CN103839479A (zh) | 一种高效电子地图注记交互方法 | |
CN113535867A (zh) | 一种适配多数据源的矢量瓦片生成方法及系统 | |
CN115510337A (zh) | 一种三维立体注记调度显示的方法 | |
CN112380302A (zh) | 基于轨迹数据的热力图生成方法、装置、电子设备以及存储介质 | |
US10013474B2 (en) | System and method for hierarchical synchronization of a dataset of image tiles | |
CN111914042B (zh) | 一种处理电力gis中矢量和栅格切片数据的方法 | |
CN116051764B (zh) | 一种遥感影像地图绘制方法及系统 | |
US20080165207A1 (en) | Method and System for Manipulating Graphical Images | |
CN114661839A (zh) | 一种基于瓦片的轨迹地图的展示方法、装置、介质和设备 | |
CN112115226B (zh) | 地图渲染方法和地图渲染装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |