Filter coords whether they are in Switzerland

calc_coords.py

#!/usr/bin/env python
import csv
import sys
import os
import mercantile

DELIMITER = os.getenv('DELIMITER', ' ')


def calculate_center(x, y, zoom):
    bounds = mercantile.bounds(x, y, zoom)
    height = bounds.north - bounds.south
    width = bounds.east - bounds.west
    center = (bounds.north + height / 2, bounds.west + width / 2)

    return center


if __name__ == '__main__':
    reader = csv.reader(sys.stdin, delimiter=DELIMITER)
    writer = csv.writer(sys.stdout, delimiter=DELIMITER)

    for row in reader:
        zoom = int(row[0])
        x = int(row[1])
        y = int(row[2])
        requests = int(row[3])

        center = calculate_center(x, y, zoom)

        writer.writerow([zoom, x, y, requests] + list(center))

filter_switzerland.py

#!/usr/bin/env python
"""
Filter out all points in Switzerland
"""
import csv
import sys
import os
from decimal import *

DELIMITER = os.getenv('DELIMITER', ' ')
NORTH = Decimal('47.9922193487799')
WEST = Decimal('5.99235534667969')
EAST = Decimal('11.1243438720703')
SOUTH = Decimal('45.6769214851596')


def in_switzerland(coords):
    lat, lng = coords
    return lat < NORTH and lat > SOUTH and lng > WEST and lng < EAST


if __name__ == '__main__':
    reader = csv.reader(sys.stdin, delimiter=DELIMITER)
    writer = csv.writer(sys.stdout, delimiter=DELIMITER)

    for row in reader:
        requests = int(row[0])
        lat = Decimal(row[1])
        lng = Decimal(row[2])

        if in_switzerland((lat, lng)):
            writer.writerow([requests, lat, lng])