Nathaniel V. KELSO nvkelso

Geocoding with Mapzen Search and using geographic prioritization

The Mapzen Search API is built from open-source data and technology and allows up to 30,000 geocoding requests per day. This gist shows how to use the Prioritize within a circular region feature to get more accurate results.

Here's the results of geocoding the Starbucks locations in New York City, using the priorization feature:

The data files in this gist:

The contents of this gist now reside in their own repo at https://github.com/urschrei/router_comparison!

MapRoulette Challenge Tutorial (ßeta)

This tutorial is for MapRoulette 1 which is no longer active. Please see the MapRoulette 2 wiki for up to date documentation.

You have played MapRoulette. You have seen some of the fun challenges. If you are reading this, you are probably thinking: 'I have a great idea for the next MapRoulette challenge!'

Great! That is exactly what I am here to explain step by step. So let's get started!

Finding the Simulator

You can test with the iOS Simulator that comes with Xcode. Navigate to the Xcode app in the Finder, right click and select "Show Package Contents":

Then navigate to Contents > Applications, and open the shortcut to "iPhone Simulator" (it may be called "iOS Simulator" depending on which version of Xcode you're running):

This script processes VIIRS Nighttime Lights 2012 GeoTIFFs from the Earth Observation Group at NOAA National Geophysical Data Center to prepare them for rendering in TileMill and uploading to MapBox Hosting.

Read Lights of the Night on MapBox to learn more about NPP the functions of this script.

	#-------------------------------------------------------------------------------
	# Name: combine_mbtiles.py
	# Purpose: Processes multiple sqlite .mbtiles files into single file (or by level? or by degree grids?)
	#
	# Author: joshn. Modified by ache015
	#
	# Created: 18/11/2014
	#
	# Last Modified: 10/11/2018
	#-------------------------------------------------------------------------------

	WITH RECURSIVE
	x(i) AS ( VALUES (0)
	UNION ALL SELECT i + 1
	FROM x
	WHERE i < 101),
	Z(Ix, Iy, Cx, Cy, X, Y, I) AS (
	SELECT
	Ix,
	Iy,
	X :: FLOAT,

	# coding=utf-8
	''' Guess at ISO-3166-2 admin 1 codes for Natural Earth data.

	Uses subdivision names from Ola Holmström’s iso-3166-2:
	https://github.com/olahol/iso-3166-2.js/blob/master/data.csv

	ne_10m_admin_1_states_provinces.csv is a direct CSV
	export of Natural Earth admin-1 data as of Nov 8, 2012:
	https://github.com/nvkelso/natural-earth-vector/tree/c97860331/10m_cultural

	————————————————————————————————————————————————————————————————————————————————————————————————————
	BBEdit / BBEdit-Lite / TextWrangler Regular Expression Guide Modified: 2018/08/10 01:19
	————————————————————————————————————————————————————————————————————————————————————————————————————

	NOTES:

	The PCRE engine (Perl Compatible Regular Expressions) is what BBEdit and TextWrangler use.

	Items I'm unsure of are marked '# PCRE?'. The list while fairly comprehensive is not complete.

	from sys import stderr
	from osgeo import ogr
	from time import time
	from math import sqrt

	def not_big(feature):
	highway = feature.GetFieldAsString(feature.GetFieldIndex('type'))

	return highway in ('residential', 'unclassified', 'road', 'minor')