Skip to content

Instantly share code, notes, and snippets.

Uli Stroetz ustroetz

Block or report user

Report or block ustroetz

Hide content and notifications from this user.

Learn more about blocking users

Contact Support about this user’s behavior.

Learn more about reporting abuse

Report abuse
View GitHub Profile
View places.geojson
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View cafe.geojson
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
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
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
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
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
You can’t perform that action at this time.