DATA=$(cat <<'END_HEREDOC'
{"name": "Calwer/Herrenberger",
 "entry": {"type":"Feature","properties":{"type": "entry"},
   "geometry":{"type":"Polygon",
   "coordinates":[[[9.003386,48.683486],[9.003185,48.68335],[9.003226,48.683275],
     [9.003769,48.683492],[9.003679,48.68359],[9.003386,48.683486]]]}},
 "exit": {"type":"Feature","properties":{"type": "exit"},
   "geometry":{"type":"Polygon",
   "coordinates":[[[9.002798,48.683225],[9.002868,48.683157],[9.002577,48.683058],
     [9.002519,48.683127],[9.002798,48.683225]]]}}}
END_HEREDOC
)
echo $DATA | sqlite-utils insert crossings.db crossings -

var example = [
    // entry
    {"type":"Feature","properties":{"type":"entry"},"geometry":{"type":"Polygon","coordinates":[[[9.0033,48.683546],[9.003185,48.68335],[9.003226,48.683275],[9.003769,48.683492],[9.003608,48.683654],[9.0033,48.683546]]]}},
    // exit
    {"type":"Feature","properties":{"type":"exit"},"geometry":{"type":"Polygon","coordinates":[[[9.002745,48.683301],[9.002868,48.683157],[9.002577,48.683058],[9.002447,48.6832],[9.002592,48.683245],[9.002745,48.683301]]]}},
    // track
    {"type": "Feature", "geometry": {"type": "LineString",
            "coordinates": [...<see above>...]}}
        ]

  L.geoJSON(example, {
    style: function(feature) {
        if (feature.properties){
            switch (feature.properties.type) {
                   case 'entry': return {color: "green"}
                   case 'exit':   return {color: "violet"}
            }}}}).addTo(map);

import json
from shapely.geometry import shape
from sqlite_utils.db import Database

# create a new database
db = Database("crossings_demo.db")

# add the spatialize tables
db.init_spatialite("/usr/lib/mod_spatialite.so")  # path to library file

# table for polygons
table = db["areas"]
table.insert_all([
        {"id": 1, "type": "entry",
            "geojson": {
                "type": "Polygon",
                "coordinates": [
                    [
                        [9.003386, 48.683486], [9.003185, 48.68335], [9.003226, 48.683275],
                        [9.003769, 48.683492], [9.003679, 48.68359], [9.003386, 48.683486],
                    ]]}},
        {"id": 2, "type": "exit",
            "geojson": {
                "type": "Polygon",
                "coordinates": [
                    [
                        [9.002798, 48.683225], [9.002868, 48.683157], [9.002577, 48.683058],
                        [9.002519, 48.683127], [9.002798, 48.683225],
                    ]]}},
    ], pk="id"
)

# add a geometry column
table.add_geometry_column(f"geometry", "POLYGON")

# convert geojson to geometry field
for row in table.rows:
    wkt = shape(json.loads(row["geojson"])).wkt
    db.execute(
        f"update areas set geometry = GeomFromText(?, 4326) where id=?",
        [wkt, row["id"]],
    )

# test if points are in areas
for x, y in (
    (9.003473, 48.683439),  # in entry
    (9.002692, 48.683142),  # in exit
    (9, 48),  # outside of both
):
    r = db.query(
        "select :x, :y, type, WITHIN(MakePoint(:x, :y), areas.geometry) as inside from areas",
        {"x": x, "y": y}
    )
    print(list(r))

# returns
# [{':x': 9.003473, ':y': 48.683439, 'type': 'entry', 'inside': 1},
#  {':x': 9.003473, ':y': 48.683439, 'type': 'exit', 'inside': 0}]
# [{':x': 9.002692, ':y': 48.683142, 'type': 'entry', 'inside': 0},
#  {':x': 9.002692, ':y': 48.683142, 'type': 'exit', 'inside': 1}]
# [{':x': 9, ':y': 48, 'type': 'entry', 'inside': 0}, {':x': 9, ':y': 48, 'type': 'exit', 'inside': 0}]

# query for a track intersection with the entry/exit areas
track = {"type": "LineString", "coordinates": [
  [9.003845, 48.683697], [9.003795, 48.683676], [9.003744, 48.683656],
  [9.003694, 48.683636], [9.003652, 48.683621], [9.003619, 48.683609],
  [9.003592, 48.683598], [9.003592, 48.683598], [9.003567, 48.683579],
  [9.003337, 48.683441], [9.003328, 48.683438], [9.003304, 48.683428],
  [9.00271, 48.683207], [9.002634, 48.683182], [9.002555, 48.68316],
]}

wkt = shape(track).wkt
r = db.query(
  "select type, Intersects(GeomFromText(?), Envelope(areas.geometry)) as crosses from areas",
  [wkt],
)
print(list(r))

# returns
# [{'type': 'entry', 'crosses': 1}, {'type': 'exit', 'crosses': 1}]

import httpx

r = httpx.get(
    "https://api-eu.stadiamaps.com/geocoding/v1/reverse",
    params={
        "point.lat": "48.68330",
        "point.lon": "9.00311",
        "boundary.circle.radius": 1,  # 1 km radius
        "layers": "street",
        "sources": "openstreetmap",
        "api_key": "your-api-key-goes-here",
    },
)

for index, street in enumerate(r.json()["features"]):
    print(
        {
            k: v
            for k, v in street.get("properties").items()
            if k in ["distance", "label", "gid"]
        }
    )
    if index == 4:
        # stop after max 5 streets
        break

r = httpx.get(
    "https://api.opencagedata.com/geocode/v1/json",
    params={
        "q": "48.68330,9.00311",
        "language": "en",
        "key": "your-api-keys-goes-here",
        "pretty": 1,
    },
)

result = r.json()["results"][0]
print(result.get("formatted"))