仅使用OSM地图中心和变焦因子的坐标来计算面积边界

Question

我正在使用 mapscii （由rastapasta）与Python和Nodejs一起进行项目。我需要在MAPSCII输出中通过其位置渲染一些对象。 MAPSCII使用OSM瓷砖层生成ASCII地图。我只知道地图中心的坐标，变焦级别以及ASCII地图的行/列数。

您是否有有关如何计算边界（左上和右下角）的提示，以便我可以将局部坐标系映射到ACSII数据上？

以这些变量为例：

def calculate_boundaries(lat, lng, zoom, width, height) -> tuple:
    ...

lat = 51.5252178
lng = -0.0925642
zoom = 17
width = 80
height = 24
upper_left, lower_right = calculate_boundaries(lat, lng, zoom, width, height)

我偶然发现了 wiki noreferrer“有帮助，因为我不使用瓷砖数字，而是在纬度/经度上。

// 编辑 这甚至可行吗？还是在每个缩放级别的2D MAPSCII数组中移动时更容易注意到多少LAT/LNG更改？

Answer 1

让我们考虑像素的宽度和高度。

计算每个像素的仪表（m/px）
由01（一个）像素表示的距离由：

s = c*cos（y）/2^（z+8）
其中：

c是赤道上地球的圆周；
z是变焦级别；
Y是您想要获得比例的纬度。

来源：

是一个错误，同时转换为Radian。所以我进口了pi和转换。

from math import cos, pi

def calculate_boundaries(lat, lng, zoom, width, height): # -> tuple:
    upper_left, lower_right = {}, {}
    C = 40075 # km - Equator distance around the world
    y = pi * lat / 180 # convert latitude degree to radian
    S = C * cos(y) / 2 ** (zoom + 8) # km distance of 1 px - https://wiki.openstreetmap.org/wiki/Pt:Zoom_levels
    S_deg = S * cos(y) / 100 # convert km (distance of 1 px) to degrees (coordinates)

    upper_left['lat'] = lat + height / 2 * S_deg
    upper_left['lng'] = lng - width / 2 * S_deg

    lower_right['lat'] = lat - height / 2 * S_deg
    lower_right['lng'] = lng + width / 2 * S_deg

    return upper_left, lower_right

# main
lat = 51.5252178
lng = -0.0925642
zoom = 17 # zoom
width = 80 # considered as pixels
height = 24
upper_left, lower_right = calculate_boundaries(lat, lng, zoom, width, height)

print(upper_left, lower_right)