在服务器上找到正确的传单，并指定经度/纬度

Question

我如何在服务器上找到正确的传单瓦片，以给出经度/纬度。

我遵循此演示文稿的幻灯片19上的公式：https://www.tu-chemnitz.de/urz/stammtisch/rsrc/vortrag-stammtisch.pdf（德语）。

目标： 在给定的经度和纬度下计算（缩放，x和y）

其他信息：

我找到了这个文档https://mapserver.org/mapcache/services.html，其中指出x为column number in zxy naming scheme，y为拖曳编号。

我尝试过的事情：

#London
long = 51.52
lat = -0.18
zoom = 5
lambda = pi/180*long
phi = pi/180*lat

x = 2^zoom*(lambda + pi)/(2*pi)
y = 2^zoom*(pi - log(tan(phi) + 1/(cos(phi))))/(2*pi)
x
y

glue("https://tile.openstreetmap.org/{zoom}/{x}/{y}.png")

Answer 1

您可以将DROP PROCEDURE IF EXISTS insertRandom; CREATE PROCEDURE insertRandom() BEGIN DECLARE mytime timestamp; SET mytime := '2009-01-01 00:00:00' BEGIN test_loop : LOOP while mytime < now() mytime = mytime + interval '8 hours'; insert into tempdata(temp_val, datum) values((select random()*(110)-10), mytime); END LOOP; END; CALL insertRandom; SELECT * FROM `temp_table`; 地图包装在一个非常简单的leaflet应用中，然后使用leaflet shiny map events（在 Inputs / Events 部分）

shiny

在地图上单击时，您将获得图块信息，并且该信息显示在下方

巴黎的结果，请给此链接https://tile.openstreetmap.org/8/129/88.png