3857切片规则

文章目录

目录

前言

一、名词术语定义

1、地球半径(earth radius)：

2、切图范围(extent)：

3、切图原点(origin)：

4、瓦片大小(tile size)：

5、LOD层级(level)：

6、切图比例尺(scale)：

7、地面分辨率(resolution)：

8、屏幕显示分辨率(DPI/PPI)：

二、关系推导

1.分辨率与层级关系推导

2.比例尺、分辨率、DPI关系推导

3.比例尺与层级关系推导

二、瓦片行列号计算

总结

---

前言

提示：。

---

一、名词术语定义

1、地球半径(earth radius)：

WGS84和CGCS2000地球椭球半径， r=6378137m

2、切图范围(extent)：

以赤道和本初子午线交点为中心，赤道长度为边长的正方形范围。

(xmin, ymin, xmax, ymax) = (-πr, -πr, πr, πr)≈

(-20037508.3427892, -20037508.3427892, 20037508.3427892, 20037508.3427892)

（经纬度范围(-180,-85,180,85)）

3、切图原点(origin)：

默认切图原点为切图范围的左上角(X,Y)=(-πr,πr)≈(-20037508.3427892, 20037508.3427892)

4、瓦片大小(tile size)：

单个瓦片大小，用长*宽表示，单位为像素，常见有256*256(默认)、512*512。

5、LOD层级(level)：

瓦片金字塔的层级数，从0级开始，到23级结束，共24级，0级瓦片数为1，其后每一层瓦片数是上一层的4倍（行列数是2倍）,即 瓦片数=4^n，行数=2^n，列数=2^n

6、切图比例尺(scale)：

S=图上距离/实际距离

7、地面分辨率(resolution)：

R=实际距离/像素数

8、屏幕显示分辨率(DPI/PPI)：

P=像素数/图上距离（英寸inch）,1inch = 0.0254m

二、关系推导

1.分辨率与层级关系推导

对于层级n下的整幅地图，其宽度为2πr，切片大小为ts（像素），切片列数为2^n，由此分辨率公式如下：

R = 2πr/(ts*2^n)

在默认情况下（ts=256），R=156543.03392804096153584694438047/2^n，例如：在10层级下（n=10）,R=152.874056570。

2.比例尺、分辨率、DPI关系推导

由术语定义可知：

比例尺：S=图上距离（米）/实际距离

地面分辨率：R=实际距离/像素数

屏幕显示分辨率：P=像素数/图上距离（英寸）

图上距离（米）=0.0254·图上距离（英寸）

则：

S = 0.0254/R*P = 1:(R*P/0.0254)

在切片规则中，习惯上比例尺只取分母表示，即

1/S=R*P/0.0254

例如：R=0.074645535411904163，P=96，则1/S=282.12486。

3.比例尺与层级关系推导

结合1、2公式，有

1/S=2πr·P/(0.0254*ts*2^n)

在默认情况下（ts=256，P=96），1/S=591658710.90913119320635065592619/2^n，例如，在10层级下（n=10），1/S=577791.709872。

从以上关系来看，默认情况下，比例尺和分辨率只与瓦片层级有关，两者之间的比值

S/R=96/0.0254=3779.527559

二、瓦片行列号计算

瓦片的行列号从0开始，依次递增，计算某点(x1,y1)在第n层级的行列号公式如下：

行号：row=|floor((y1-Y)/(ts*R))|

列号：col=|floor((x1-X)/(ts*R))|

结合第二小节推导公式：

行号：row=|floor((y1-Y)*2^n/40075016.685578)|

列号：col=|floor((x1-X)*2^n/40075016.685578)|

---

总结

未完待续。。。