天地图解析(瓦片层级、比例尺、行列号计算)

2024-01-15 14:38

本文主要是介绍天地图解析(瓦片层级、比例尺、行列号计算),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

目录

瓦片元数据及分级介绍

比例尺如何计算得到

知道一个点坐标,测算某层级下瓦片行列号


常用的天地图瓦片底图有地图和影像两种,每一种又分别有经纬度投影和球面墨卡托投影两种坐标,请求地址中_c表示经纬度投影,_w表示墨卡托投影。

 

                         地图                                 影像

地图是由矢量底图(vec)和矢量注记(cva)叠加显示组合而成。

  =       +    ​​​​​​​

影像是由影像底图(img)和影像注记(cia)叠加显示组合而成。

  =    +  

  影像底图请求的瓦片格式为jpg,影像注记、地图的为png。png格式无数据的地方为透明状态。

下面以经纬度投影的影像底图(img_c)为例,对天地图的瓦片进行分析。

  • 瓦片元数据及分级介绍

元数据获取地址:

https://t0.tianditu.gov.cn/img_c/wmts?request=GetCapabilities&service=wmts

里面记录了天地图服务的相关信息及数据参数。

TileMatrixSet标签下记录了各层级的比例尺,单个瓦片像素大小(256*256),此层级下瓦片的列数(MatrixWidth)、行数(MatrixHeight)。

以下为第1级的相关参数:

通过TileMatrix标签我们可以知道第1级总共有1行、2列,对应的瓦片如下图所示:

左边瓦片对应的地址为:

https://t2.tianditu.gov.cn/img_c/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=img&STYLE=default&TILEMATRIXSET=c&FORMAT=tiles&TILECOL=0&TILEROW=0&TILEMATRIX=1&tk=68d166cfe304fa077ff035bed00edc37

右边瓦片对应的地址为:

https://t2.tianditu.gov.cn/img_c/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=img&STYLE=default&TILEMATRIXSET=c&FORMAT=tiles&TILECOL=1&TILEROW=0&TILEMATRIX=1&tk=68d166cfe304fa077ff035bed00edc37

第2级总共2行、4列,8张瓦片。

 第3级的行数、列数分别是第2级的2倍,瓦片数是4倍。

。。。。

各级行数row=2^{i-1},列数col=2^{i}(i为级数,从1开始起算) 

  • 比例尺如何计算得到

下面我们来通过计算来验证下第一级的比例尺:

两张瓦片本身的宽D=256*2*0.0254/96≈0.13547

注:dpi为96,即1英寸=0.0254米=96像素。

两张瓦片对应的实地宽(即椭球体上赤道面的周长,椭球的长半轴为6378137)L=2*π*6378137≈40,075,016.68559

第1级比例尺分母=\tfrac{2*\pi *6378137}{256*2*0.0254/96}≈295,829,355.454565596≈2.958293554545656E8

  验算出的比例尺与元数据里的比例尺一致。

 实际情况下,同一级的瓦片,位于赤道和靠近两极的比例尺是不一样的。元数据中的比例尺应该是按照赤道线进行计算得到。

推算出各级比例尺分母=\tfrac{2*\pi *6378137}{256*2^{i}*0.0254/96}

各级像素大小=​​​​​​​\tfrac{360}{256*2^{i}}(单位为度)。

  • 知道一个点坐标,测算某层级下瓦片行列号

假定一个物体所在的经度为lg,纬度为la,则在第i级时,该物体所在的瓦片

行号NubRow=[\tfrac{90-la}{180/2^{i-1}}]+1

列号NubCol=[\tfrac{lg-(-180)}{360/2^{i}}]+1

计算得到的行列号在请求瓦片时需减1,因为天地图的瓦片行列号从0开始算,层级从1开始算。

重庆朝天门的经纬度坐标大概为(106.58828259,29.56782092),通过公式计算在15级时所在的瓦片行号为5500,列号为26085,对应的url为https://t2.tianditu.gov.cn/img_c/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=img&STYLE=default&TILEMATRIXSET=c&FORMAT=tiles&TILECOL=26085&TILEROW=5500&TILEMATRIX=15&tk=68d166cfe304fa077ff035bed00edc37

这篇关于天地图解析(瓦片层级、比例尺、行列号计算)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/609236

相关文章

在C#中合并和解析相对路径方式

《在C#中合并和解析相对路径方式》Path类提供了几个用于操作文件路径的静态方法,其中包括Combine方法和GetFullPath方法,Combine方法将两个路径合并在一起,但不会解析包含相对元素... 目录C#合并和解析相对路径System.IO.Path类幸运的是总结C#合并和解析相对路径对于 C

如何用Java结合经纬度位置计算目标点的日出日落时间详解

《如何用Java结合经纬度位置计算目标点的日出日落时间详解》这篇文章主详细讲解了如何基于目标点的经纬度计算日出日落时间,提供了在线API和Java库两种计算方法,并通过实际案例展示了其应用,需要的朋友... 目录前言一、应用示例1、天安门升旗时间2、湖南省日出日落信息二、Java日出日落计算1、在线API2

Java解析JSON的六种方案

《Java解析JSON的六种方案》这篇文章介绍了6种JSON解析方案,包括Jackson、Gson、FastJSON、JsonPath、、手动解析,分别阐述了它们的功能特点、代码示例、高级功能、优缺点... 目录前言1. 使用 Jackson:业界标配功能特点代码示例高级功能优缺点2. 使用 Gson:轻量

Java如何接收并解析HL7协议数据

《Java如何接收并解析HL7协议数据》文章主要介绍了HL7协议及其在医疗行业中的应用,详细描述了如何配置环境、接收和解析数据,以及与前端进行交互的实现方法,文章还分享了使用7Edit工具进行调试的经... 目录一、前言二、正文1、环境配置2、数据接收:HL7Monitor3、数据解析:HL7Busines

python解析HTML并提取span标签中的文本

《python解析HTML并提取span标签中的文本》在网页开发和数据抓取过程中,我们经常需要从HTML页面中提取信息,尤其是span元素中的文本,span标签是一个行内元素,通常用于包装一小段文本或... 目录一、安装相关依赖二、html 页面结构三、使用 BeautifulSoup javascript

网页解析 lxml 库--实战

lxml库使用流程 lxml 是 Python 的第三方解析库,完全使用 Python 语言编写,它对 XPath表达式提供了良好的支 持,因此能够了高效地解析 HTML/XML 文档。本节讲解如何通过 lxml 库解析 HTML 文档。 pip install lxml lxm| 库提供了一个 etree 模块,该模块专门用来解析 HTML/XML 文档,下面来介绍一下 lxml 库

无人叉车3d激光slam多房间建图定位异常处理方案-墙体画线地图切分方案

墙体画线地图切分方案 针对问题:墙体两侧特征混淆误匹配,导致建图和定位偏差,表现为过门跳变、外月台走歪等 ·解决思路:预期的根治方案IGICP需要较长时间完成上线,先使用切分地图的工程化方案,即墙体两侧切分为不同地图,在某一侧只使用该侧地图进行定位 方案思路 切分原理:切分地图基于关键帧位置,而非点云。 理论基础:光照是直线的,一帧点云必定只能照射到墙的一侧,无法同时照到两侧实践考虑:关

【C++】_list常用方法解析及模拟实现

相信自己的力量,只要对自己始终保持信心,尽自己最大努力去完成任何事,就算事情最终结果是失败了,努力了也不留遗憾。💓💓💓 目录   ✨说在前面 🍋知识点一:什么是list? •🌰1.list的定义 •🌰2.list的基本特性 •🌰3.常用接口介绍 🍋知识点二:list常用接口 •🌰1.默认成员函数 🔥构造函数(⭐) 🔥析构函数 •🌰2.list对象

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc