Skip to content

Instantly share code, notes, and snippets.

View places.geojson
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View cafe.geojson
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View fraenkische_isochrone.geojson
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View connection_line_and_route.geojson
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View connection_routes_join_with_trips.py
connection_routes_gdf = gpd.sjoin(
routes_segments_gdf, routes_gdf, how='left', op='within').dissolve(by='id', aggfunc='sum')
connection_routes_gdf.to_file(
'connection_routes.geojson', driver='GeoJSON')
View connection_routes_unique_segments.py
routes_segments = routes_gdf.unary_union
routes_segments_gdf = gpd.GeoDataFrame.from_features([{
"type": "Feature",
"properties": {"id": count},
"geometry": mapping(geometry)
} for id, geometry in enumerate(routes_segments)])
View connection_routes_all_routes.py
with open('connection_lines.geojson') as connection_lines_file:
connection_lines = json.loads(connection_lines_file.read())
route_features = []
for line_feature in connection_lines['features']:
origin = {
'type': 'Feature',
'geometry': {
'type': 'Point',
View map.geojson
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View pt_lisbon_lis.geojson
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@ustroetz
ustroetz / mvt.sql
Last active Sep 8, 2017
Get MVT from Postgis
View mvt.sql
SELECT ST_AsMVT('test', 4096, 'geom', q)
FROM (
SELECT
ST_AsMVTGeom(
geom,
TileBBox(z, x, y, 3857), /* https://github.com/mapbox/postgis-vt-util/blob/master/src/TileBBox.sql */
4096,
0,
false
) geom