Hermes Agent
c0c9441930
- Geocoding via PDOK Locatieserver - Isochroon (auto 20 min) via OpenRouteService - Rapport generatie in markdown formaat - Voorbeeld adres: Griendvelden 14, Best
131 lines
4.0 KiB
Python
131 lines
4.0 KiB
Python
#!/usr/bin/env python3
|
|
"""
|
|
Adres Analyse - Rapport generator
|
|
Gebruikt PDOK voor geocoding en ORS voor isochronen
|
|
"""
|
|
|
|
import json
|
|
import sys
|
|
import urllib.request
|
|
import urllib.parse
|
|
from datetime import datetime
|
|
|
|
# Config
|
|
ORS_BASE = "http://192.168.1.71:9080/ors/v2"
|
|
PDOK_BASE = "https://api.pdok.nl/bzk/locatieserver/search/v3_1"
|
|
|
|
def geocode_adres(adres):
|
|
"""Geocode adres naar coördinaten via PDOK Locatieserver"""
|
|
# Gebruik suggest endpoint voor betere adres matching
|
|
params = urllib.parse.urlencode({
|
|
'q': adres,
|
|
'rows': 5,
|
|
'fl': 'centroide_ll,weergavenaam,postcode,woonplaatsnaam,type,huisnummer'
|
|
})
|
|
|
|
url = f"{PDOK_BASE}/free?{params}"
|
|
with urllib.request.urlopen(url, timeout=10) as resp:
|
|
data = json.loads(resp.read())
|
|
|
|
if not data['response']['docs']:
|
|
raise ValueError(f"Adres niet gevonden: {adres}")
|
|
|
|
# Zoek naar exact adres match (type=adres)
|
|
doc = None
|
|
for d in data['response']['docs']:
|
|
if d.get('type') == 'adres':
|
|
doc = d
|
|
break
|
|
|
|
if not doc:
|
|
doc = data['response']['docs'][0] # Fallback naar eerste resultaat
|
|
|
|
# Parse POINT(lon lat) format
|
|
point = doc['centroide_ll']
|
|
coords = point.replace('POINT(', '').replace(')', '').split()
|
|
lon, lat = float(coords[0]), float(coords[1])
|
|
|
|
return {
|
|
'adres': doc.get('weergavenaam', adres),
|
|
'postcode': doc.get('postcode', ''),
|
|
'woonplaats': doc.get('woonplaatsnaam', ''),
|
|
'lat': lat,
|
|
'lon': lon
|
|
}
|
|
|
|
def haal_isochroon(lon, lat, profiel='driving-car', minuten=20):
|
|
"""Haal isochroon op van ORS"""
|
|
url = f"{ORS_BASE}/isochrones/{profiel}"
|
|
payload = {
|
|
"locations": [[lon, lat]],
|
|
"range": [minuten * 60], # seconden
|
|
"attributes": ["area", "reachfactor", "total_pop"]
|
|
}
|
|
|
|
req = urllib.request.Request(
|
|
url,
|
|
data=json.dumps(payload).encode(),
|
|
headers={'Content-Type': 'application/json'},
|
|
method='POST'
|
|
)
|
|
|
|
with urllib.request.urlopen(req, timeout=30) as resp:
|
|
return json.loads(resp.read())
|
|
|
|
def format_oppervlakte(m2):
|
|
"""Converteer m² naar leesbaar formaat"""
|
|
if m2 >= 1_000_000:
|
|
return f"{m2/1_000_000:.1f} km²"
|
|
elif m2 >= 10_000:
|
|
return f"{m2/10_000:.1f} ha"
|
|
else:
|
|
return f"{m2:.0f} m²"
|
|
|
|
def genereer_rapport(adres_input):
|
|
"""Genereer volledig adres analyse rapport"""
|
|
print(f"📍 **Adres Analyse Rapport**")
|
|
print(f"Gegenereerd: {datetime.now().strftime('%d-%m-%Y %H:%M')}")
|
|
print()
|
|
|
|
# Geocoding
|
|
print("🔍 Geocoding...")
|
|
locatie = geocode_adres(adres_input)
|
|
print(f"**Adres:** {locatie['adres']}")
|
|
print(f"**Coördinaten:** {locatie['lat']:.6f}, {locatie['lon']:.6f}")
|
|
print()
|
|
|
|
# Isochroon
|
|
print("🚗 **Bereikbaarheid (Auto 20 min)**")
|
|
iso_data = haal_isochroon(locatie['lon'], locatie['lat'], 'driving-car', 20)
|
|
|
|
if 'features' in iso_data and iso_data['features']:
|
|
props = iso_data['features'][0]['properties']
|
|
area_m2 = props.get('area', 0)
|
|
reach = props.get('reachfactor', 0)
|
|
|
|
print(f"**Oppervlakte:** {format_oppervlakte(area_m2)}")
|
|
print(f"**Reach factor:** {reach:.3f}")
|
|
print()
|
|
|
|
# Bereken bbox
|
|
bbox = iso_data.get('bbox', [])
|
|
if len(bbox) == 4:
|
|
print(f"**Bereik:** {bbox[1]:.3f}°-{bbox[3]:.3f}°N, {bbox[0]:.3f}°-{bbox[2]:.3f}°E")
|
|
print()
|
|
|
|
# Save GeoJSON voor kaart generatie
|
|
output_file = f"/tmp/isochroon_{datetime.now().strftime('%Y%m%d_%H%M%S')}.geojson"
|
|
with open(output_file, 'w') as f:
|
|
json.dump(iso_data, f)
|
|
print(f"💾 Data opgeslagen: `{output_file}`")
|
|
else:
|
|
print("❌ Geen isochroon data ontvangen")
|
|
|
|
print()
|
|
print("---")
|
|
print("_Onderdeel 1/?: Isochroon bereikbaarheid_")
|
|
|
|
if __name__ == '__main__':
|
|
adres = sys.argv[1] if len(sys.argv) > 1 else "Griendvelden 14 Best"
|
|
genereer_rapport(adres)
|