lukasmartinelli’s gists

lukasmartinelli / calc_coords.py

Last active May 4, 2016 15:04

Filter coords whether they are in Switzerland

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

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


	def calculate_center(x, y, zoom):

lukasmartinelli / extract_cities.py

Created April 26, 2016 14:20

Extract city bounding box CSV from Mapzen

	import json

	doc = json.load(open('cities.json'))
	for region, region_obj in doc['regions'].items():
	for city_name, city_obj in region_obj['cities'].items():
	bbox = city_obj['bbox']
	parts = city_name.split('_')

	if len(parts) != 2:
	continue

lukasmartinelli / tile_queries_prior.sql

Created March 11, 2016 08:49

Tiles

	CREATE VIEW poi_label_zoom_level_changes AS (
	SELECT osm_id, timestamp FROM poi_label_z14
	);


	CREATE VIEW road_zoom_level_changes AS (
	SELECT osm_id, timestamp FROM road_z14
	UNION ALL
	SELECT osm_id, timestamp FROM road_z13
	UNION ALL

lukasmartinelli / function.sql

Created March 8, 2016 18:57

Try optimize PostgreSQL function


	CREATE OR REPLACE FUNCTION point_to_tiles(
	_point geometry
	) RETURNS SETOF tile
	AS $$
	DECLARE
	d2r CONSTANT DOUBLE PRECISION := pi() / 180;
	lon CONSTANT DOUBLE PRECISION := st_x(ST_Transform(_point, 4326));
	lat CONSTANT DOUBLE PRECISION := st_y(ST_Transform(_point, 4326));
	_sin CONSTANT DOUBLE PRECISION := sin(lat * d2r);

lukasmartinelli / places.tsv

Created March 4, 2016 14:19

OSMNames Places Export

lukasmartinelli / missing_tiles.csv

Created December 11, 2015 18:26

Missing Tiles

We can make this file beautiful and searchable if this error is corrected: No commas found in this CSV file in line 0.

lukasmartinelli / find_missing_mbtiles.md

Created December 11, 2015 18:22

Find out missing MBTiles at certain zoom level from list of files

Find the real tiles that should exist.

mercantile children "[0,0,0]" --depth 8 > correct_tiles.txt

Find the list of MBTiles from S3

lukasmartinelli / create_extracts.sh

Last active August 9, 2021 16:28

Create coutnry extracts from world MBTiles

	#!/bin/bash

	tilelive-copy --minzoom=0 --maxzoom=14 --bounds="60.403889,29.288333,74.989862,38.5899217" world.mbtiles afghanistan.mbtiles
	tilelive-copy --minzoom=0 --maxzoom=14 --bounds="19.0246095,39.5448625,21.1574335,42.7611669" world.mbtiles albania.mbtiles
	tilelive-copy --minzoom=0 --maxzoom=14 --bounds="-8.7689089,18.868147,12.097337,37.3962055" world.mbtiles algeria.mbtiles
	tilelive-copy --minzoom=0 --maxzoom=14 --bounds="1.3135781,42.3288238,1.8863837,42.7559357" world.mbtiles andorra.mbtiles
	tilelive-copy --minzoom=0 --maxzoom=14 --bounds="11.3609793,-18.1389449,24.18212,-4.2888889" world.mbtiles angola.mbtiles
	tilelive-copy --minzoom=0 --maxzoom=14 --bounds="-63.7391991,17.9609378,-62.6125448,18.8951194" world.mbtiles anguilla.mbtiles
	tilelive-copy --minzoom=0 --maxzoom=14 --bounds="-73.6603073,-55.285076,-53.5374514,-21.6811679" world.mbtiles argentina.mbtiles
	tilelive-copy --minzoom=0 --maxzoom=14 --bounds="43.3471395,38.7404775,46.7333087,41.400712" world.mbtiles armenia.mbtiles

lukasmartinelli / handle_img.py

Last active October 21, 2015 13:27

Handle user uploaded images in Python

	from img_rotate import fix_orientation
	from PIL import Image


	def fix_img(img):
	proper_img = fix_orientation(img)
	proper_img.thumbnail(1500, 1500), Image.ANTIALIAS)
	return proper_img

lukasmartinelli / delete_dups.sql

Last active August 29, 2015 14:27

	DELETE FROM import.filepaths_1
	WHERE data->>'repo' IN
	(SELECT data->>'repo'
	FROM (SELECT data->>'repo', ROW_NUMBER() OVER (partition BY data->>'repo' ORDER BY data->>'repo') AS rnum
	FROM import.filepaths_1) t
	WHERE t.rnum > 1
	);

Lukas Martinelli lukasmartinelli