Recipes

import json, os, statistics, urllib.parse, urllib.request
 
BASE = "https://api.gimlabs.io/v1"
KEY = os.environ["GIMLABS_API_KEY"]
 
def get(path, **params):
    url = BASE + path + ("?" + urllib.parse.urlencode(params) if params else "")
    req = urllib.request.Request(url, headers={"Authorization": f"Bearer {KEY}"})
    with urllib.request.urlopen(req) as resp:
        return json.load(resp)
 
def paged(path, **params):
    """Follow meta.nextCursor until it comes back null."""
    cursor = None
    while True:
        body = get(path, limit=200, cursor=cursor, **params)
        yield from body["data"]
        cursor = body["meta"].get("nextCursor")
        if not cursor:
            return
 
def ground_model(project_id):
    strata, unnamed = {}, 0
 
    for row in paged(f"/projects/{project_id}/tables/GEOL"):
        d = row["data"]
        try:
            top, base = float(d["GEOL_TOP"]), float(d["GEOL_BASE"])
        except (KeyError, ValueError):
            continue
        if base <= top:
            continue
 
        # Correlate on the named formation only.
        name = (d.get("GEOL_GEOL") or "").strip()
        if not name:
            unnamed += 1
            continue
 
        s = strata.setdefault(name, {"holes": set(), "tops": [], "thick": []})
        s["holes"].add(row["locationId"])
        s["tops"].append(top)
        s["thick"].append(base - top)
 
    # Down the page in the order you would meet them.
    for name, s in sorted(strata.items(), key=lambda kv: statistics.median(kv[1]["tops"])):
        print(f"{name:<30}{len(s['holes']):>4} holes"
              f"  top {min(s['tops']):>6.2f}-{max(s['tops']):<6.2f}"
              f"  thickness {min(s['thick']):>5.2f}-{max(s['thick']):<5.2f}")
 
    if unnamed:
        print(f"{unnamed} beds had no stratum named and are not included.")