I've been playing with #OpenStreetMaps via #osmnx, which is awesome , but I struggle with simple stuff like adding a bunch of places as markers. Everything looks a bit like the owl drawing meme, either showing something too easy and useless, or something too advanced and also useless or beyond my comprehension. Maybe some other Python tools?
(I know about Marcelo's fabulous PrettyMaps but it is not exactly a viz tool)
@villares GeoDesk for Python (https://docs.geodesk.com/python) covers a lot of what you're looking for. You can easily query and analyze OpenStreetMap features and display them on a Leaflet map.
@villares - Example: To visualize all primary roads in São Paulo:
import geodesk
brazil = geodesk.Features("brazil.gol") # The file created by the GOL Tool
sao_paulo = brazil("a[boundary=administrative][admin_level=8][name='São Paulo']").one
roads = brazil("w[highway=primary]")
roads(sao_paulo).map.show()
I hope it won't be an abuse of your friendly helpfulness to ask some follow up questions:
1a. How would you geocode a list of 10 addresses to show on that map...
1b. If the former is cumbersome, it could be a list of lat-long coordinates of points of interest.
How would you show a translucent circle showing a 100 meter radius around the points of interest.
If you just need the city, you can do a contains query:
cities = brazil("a[boundary=administrative][admin_level=8]")
cities.contains(point)
print(f"Point is in {cities.first.name}")
To get the lon/lat of features:
for poi in features:
print(f"{poi.name}: {poi.lon}, {poi.lat}")
@geodesk@opencage great, I've been using Nominatim via osmnx, I guess. I was wondering if you wrapped some other gecoding feature... I'll look at @opencage!!!
@villares Yes, to_mercator() and from_mercator() are used to convert between latitude/longitude and Mercator units. Shapely simply assumes a Euclidean plane without any particular unit of measure. to_mercator(meters=100, lat=-23) gives you the number of Shapely units that equate to 100 meters at -23 degrees latitude (In the Mercator projection the scale depends on latitude -- there are more units per meter the closer you move towards the polar regions). https://docs.geodesk.com/python/mercator
@villares In theory, Leaflet can work with other projections, but virtually all the base map tiles are in Web Mercator (EPSG:3857). GeoDesk uses a slightly modified projection based on EPSG:3857, here is more info: https://docs.geodesk.com/core-concepts#coordinate-system
@villares I have used #Leaflet before and it is pretty easy to have a map running with custom markers. I didn't do anything big, but for smaller maps (maybe bigger ones too), it is a great suggestion from @slackline
Add comment