3pp-mirror/map-machine
1
0
Fork 0
mirror of https://github.com/enzet/map-machine.git synced 2025-08-03 16:49:10 +02:00
Code Issues Projects Releases Packages Wiki Activity

Issue #76: make equator length a parameter.

Browse source
This commit is contained in:
Sergey Vartanov 2021-09-01 03:04:40 +03:00
parent 902956b1b3
commit 7fcfc85d44
1 changed files with 10 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

15
roentgen/flinger.py
View file

@ -10,8 +10,6 @@ from roentgen.boundary_box import BoundaryBox
__author__ = "Sergey Vartanov"
__email__ = "me@enzet.ru"
EQUATOR_LENGTH: float = 40_075_017 # meters
def pseudo_mercator(coordinates: np.ndarray) -> np.ndarray:
"""
@ -27,15 +25,19 @@ def pseudo_mercator(coordinates: np.ndarray) -> np.ndarray:
return np.array((coordinates[1], y))
def osm_zoom_level_to_pixels_per_meter(zoom_level: float) -> float:
def osm_zoom_level_to_pixels_per_meter(
zoom_level: float, equator_length: float = 40_075_017.0
) -> float:
"""
Convert OSM zoom level to pixels per meter on Equator. See
https://wiki.openstreetmap.org/wiki/Zoom_levels
:param zoom_level: integer number usually not bigger than 20, but this
function allows any non-negative float value
:param equator_length: celestial body equator length in meters, default
value is set for Earth
"""
return 2 ** zoom_level / EQUATOR_LENGTH * 256
return 2 ** zoom_level / equator_length * 256
class Flinger:
@ -48,17 +50,20 @@ class Flinger:
geo_boundaries: BoundaryBox,
zoom_level: float = 18,
border: np.ndarray = np.array((0, 0)),
equator_length: float = 40_075_017.0,
) -> None:
"""
:param geo_boundaries: minimum and maximum latitude and longitude
:param zoom_level: zoom level in OpenStreetMap terminology
:param border: size of padding in pixels
:param equator_length: celestial body equator length in meters, default
value is set for Earth
"""
self.geo_boundaries: BoundaryBox = geo_boundaries
self.border: np.ndarray = border
self.ratio: float = (
osm_zoom_level_to_pixels_per_meter(zoom_level)
* EQUATOR_LENGTH
* equator_length
/ 360
)
self.size: np.ndarray = border * 2 + self.ratio * (