Documenter le panneau unique de la carte

Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>

AGENTS.md

@@ -19,7 +19,8 @@
 
 ## Conventions
 
-- **Un seul serveur web : `mapserve.py`** (image `lidar-maps`, port 8975 léger / 8973 worker). L'ancienne webapp (`webapp.py`, `export.py`, index.html/`_APP_JS`) a été supprimée : la génération de tuiles (portée de la webapp historique — `/api/preview`, `/api/generate`, `/api/status`, `/api/stop`, `/api/queue/clear`, `/api/cell`) vit dans `mapserve.py`, l'interface dans `mapui.py` (constantes `_MAP_HTML`/`_MAP_CSS`/`_MAP_JS`, écrites par `write_map_assets()` et bâchées dans les images). Sur l'image légère sans `LIDAR_GENERATION_URL`, `/api/status` répond `available: false` et l'interface masque les boutons.
+- **Un seul serveur web : `mapserve.py`** (image `lidar-maps`, port 8975 léger / 8973 worker). L'ancienne webapp (`webapp.py`, `export.py`, index.html/`_APP_JS`) a été supprimée : la génération de tuiles (portée de la webapp historique — `/api/preview`, `/api/generate`, `/api/status`, `/api/stop`, `/api/queue/clear`, `/api/cell`) vit dans `mapserve.py`, l'interface dans `web/map.{html,css,js}` (relus par `mapui.py`, constantes `_MAP_*` conservées, écrites par `write_map_assets()` et bâchées dans les images). Sur l'image légère sans `LIDAR_GENERATION_URL`, `/api/status` répond `available: false` et l'interface masque l'onglet Génération.
+- **Interface en panneau unique** (`web/map.{html,css,js}`) : un seul panneau à onglets Affichage / Export PDF / Génération (masqué si le générateur est indisponible ou non autorisé) / Partager, remplaçant l'ancienne pile de blocs empilés. Sous 720 px de large, le panneau devient un volet en bas d'écran à trois hauteurs (`closed`/`half`/`full`, poignée `#sheetGrip` glissée ou simplement touchée, ou onglet actif retouché). Un clic sur une dalle ouvre une bulle (date, densité) → « Détails » bascule une fiche complète dans le panneau (« ← Retour » rouvre l'onglet précédent) ; Échap ferme la bulle, puis la fiche, puis replie le panneau (bande d'icônes sur ordinateur, volet fermé sur téléphone — la bande reste utilisable, un clic sur un onglet redéplie le panneau). Raccourcis clavier : 1–4 (onglets, sans effet si l'onglet est masqué), P (mode d'affichage suivant), Échap. Deux thèmes clair/sombre (`lidar-theme`, `auto` par défaut, suit le système) et tous les réglages communs posés sur `:root` en variables CSS (aucune couleur en dur dans les composants). Clés `localStorage` (via `lsGet`/`lsSet`, silencieux en navigation privée) : `lidarMapView_v2` (vue/couches), `lidar-print` (réglages d'export), `lidar-panel` (onglet/hauteur du volet), `lidar-theme`.
 - **`index.py` = catalogue + registres partagés** (plus d'interface) : `VIZ_LABELS`/`VIZ_LEGENDS`, défauts d'affichage (`DEFAULT_VIZ`/`PRECISION_VIZ`/`VIEW_MODES`), `PANEL_VIZ`/`KEYWORD_TO_STEP`, `scan_tiles`/`cells_with_all_viz`, vignettes + sous-tuiles + inventaire `index_tiles.json` (`build_index`). L'inventaire est servi par `/api/tiles` de mapserve aux machines légères (`LIDAR_SOURCE_URL`).
 - **Generation is 0.2 m only** (policy): `/api/generate` (`GENERATE_RESOLUTIONS` in `mapserve.py`), the compose `process` command and the CLI `-r` default all produce 0.2 m exclusively; 0.5 m stays available via explicit `-r 0.5`. Completeness detection (`complete_cells`) requires the viz at 0.2 m only.
 - **Génération du nord au sud** : les tuiles sont traitées par ligne décroissante (row = nord en km), colonnes croissantes — `find_laz_files` (pipeline.py) pour les passes batch et `_resolve_request` (mapserve.py) pour les runs lancés depuis la carte. Les workers prennent les fichiers dans l'ordre de soumission : la carte se remplit de haut en bas pendant un run (`--file` explicite au CLI = ordre utilisateur préservé). Parallélisme de génération : `LIDAR_WORKERS` (10 dans les compose ; `auto` sinon).
@@ -50,7 +51,7 @@
 - **Calage vertical des faisceaux de vol** : chaque tuile mélange plusieurs passes (1-2 `PointSourceId` par passe) parfois biaisées verticalement de quelques cm (±2,5 cm mesurés sur 1000_6882). `create_dtm_fast` mesure l'offset robuste de chaque faisceau (points sol, maille 1 m, surface médiane itérée 3×) et retranche les offsets ≥ 0,5 cm (`STRIP_ALIGN_THRESHOLD` dans `dtm.py`) avant rastérisation. Offsets calculés **par tuile** (ils dérivent le long d'une ligne de vol : jamais de table globale), mémoïsés par LAS sol, consignés dans `DTM/*_dtm*_stripalign.json` (sidecar de cache : absent, ou version/seuil/paramètres différents ⇒ régénération du DTM). Désactivable : `--no-strip-align`.
 - **Gigue intra-faisceau (2ᵉ passe du calage)** : les lignes de balayage successives d'une MÊME passe peuvent être décalées verticalement de façon aléatoire (vibration capteur / bruit haute fréquence de trajectoire) — un offset constant par faisceau n'y suffit pas. `_strip_jitter_offsets` découpe chaque faisceau en fenêtres de temps GPS (`STRIP_JITTER_BIN` = 0,1 s, origine de temps propre à chaque faisceau), mesure l'offset robuste de chaque fenêtre contre la surface médiane des AUTRES faisceaux (maille 1 m partagée, ≥ `STRIP_JITTER_MIN_CELLS` = 40 cellules), lisse la série (médiane glissante `STRIP_JITTER_SMOOTH` = 5 fenêtres), la borne à ± `STRIP_JITTER_MAX` (10 cm) puis l'interpole au temps GPS de chaque point (`_apply_strip_jitter`) ; en recouvrement à deux faisceaux, chacun reçoit une série (chacun absorbe sa part). Requiert la dimension `gps_time` (silencieusement ignorée sinon). Sidecar version 2 (séries dans `jitter`), couverte par `--no-strip-align`.
 - **Couches produites et servies = `PANEL_VIZ` (`index.py`, aujourd'hui `('relief_oriente', 'densite_sol')`)** : pipeline sans `--only`/`--skip` (`panel_steps()`), génération lancée depuis la carte (`_panel_viz_steps` dans `mapserve.py`) et couches servies (`tiles.available_layers` filtre : panneau, `/tiles/…`, TileJSON, WMTS, JOSM). Les autres visualisations restent calculables avec `--only`. Les tests de la carte lèvent la restriction via `_setup(..., panel=None)`.
-- **Précision (`densite_sol`) et affichage de la carte** : `create_dtm_fast` écrit la densité des points sol retenus (bande de raccord comprise) dans `DTM/*_dtm*_density.tif` (mailles 1 m, moyenne 3 × 3, pts/m²) ; `generate_densite_sol` la quantifie en 16 niveaux (`density_levels` : échelle log fixe, niveau k dès 0,25 × 2^(k/2) pts/m²) et la garde à 1 m (dalle 1000², ~200 Ko en WebP sans perte, `rendering.LOSSLESS_GRAY_KEYWORDS`) ; tuiles au plus proche voisin (`tiles.NEAREST_LAYERS`) pour garder 16 gris nets aux zooms 18–19. Interface (`mapui.py`) : plus de pile de couches — une couche principale (`DEFAULT_VIZ`) et trois modes relief / précision / les deux (`VIEW_MODES`, précision en `multiply` à `DEFAULT_PRECISION_OPACITY` 0,6), touche P, lien `&M=…&P=mode:opacité`, défauts figés `{main, mode, precision_opacity, base}`.
+- **Précision (`densite_sol`) et affichage de la carte** : `create_dtm_fast` écrit la densité des points sol retenus (bande de raccord comprise) dans `DTM/*_dtm*_density.tif` (mailles 1 m, moyenne 3 × 3, pts/m²) ; `generate_densite_sol` la quantifie en 16 niveaux (`density_levels` : échelle log fixe, niveau k dès 0,25 × 2^(k/2) pts/m²) et la garde à 1 m (dalle 1000², ~200 Ko en WebP sans perte, `rendering.LOSSLESS_GRAY_KEYWORDS`) ; tuiles au plus proche voisin (`tiles.NEAREST_LAYERS`) pour garder 16 gris nets aux zooms 18–19. Interface (`web/map.js`, onglet Affichage) : plus de pile de couches — une couche principale (`DEFAULT_VIZ`) et trois modes relief / précision / les deux (`VIEW_MODES`, précision en `multiply` à `DEFAULT_PRECISION_OPACITY` 0,6), touche P, lien `&M=…&P=mode:opacité`, défauts figés `{main, mode, precision_opacity, base}`.
 - **Relief orienté (`relief_oriente`, seule couche de la carte)** : image RGB unique (GeoTIFF uint8 3 bandes, rendue telle quelle comme ortho/topo via `RGB_KEYWORDS` dans `rendering.py`) — clarté CIELAB = openness positive **locale** (MNT − gaussienne `RELIEF_DETREND_M` = 10 m, rayons `RELIEF_RADII_M` = 5/10/20 m, 16 directions) 65 % + ombrage 35 % ; teinte = aspect, chroma fixe (`RELIEF_CHROMA`). Échelle log fixe `RELIEF_OPEN_RANGE` (pas de statistique par dalle) et support total 40 m < bande de raccord 100 m : dalles jointives. Rapide : détendance + rayons sur grille décimée ~0,8 m (`RELIEF_GRID_M`), noyau dédié qui n'accumule que la moyenne des angles (`_mean_horizon_*` : CuPy `RawKernel` sur GPU, numba parallèle sur CPU, numpy en repli), colorisation fusionnée (numba) ou vectorisée sans trigonométrie (CuPy) via une table L* × teinte (`_relief_lut`). ~5 s de calcul hors préparation sur CPU 12 cœurs. Tout changement de constante change le rendu : régénérer les dalles (`--only relief_oriente --force`).
 - **Openness sous-échantillonnée** : `generate_openness` calcule le lancé de rayons (l'étape la plus coûteuse : 532 s/tuile à 0,2 m sur CPU) sur une grille décimée par blocs (`OPENNESS_DOWNSAMPLE = 2` : max par bloc en positive, min en négative — préserve les reliefs qui bornent l'horizon) puis rééchantillonne en bilinéaire. Coût ÷ facteur³ : 532 s → 40 s (×13). Signal archéologique préservé (corr. 0,93 après lissage) ; la texture de bruit sub-métrique disparaît. `--openness-downsample 1` = pleine résolution. SVF et openness anisotrope ne sont PAS concernées.
 - **Raccord des bords entre tuiles** : les rendus à grand noyau (openness/SVF : rayons 100 m ; LRM : 15 m) tronquent leur fenêtre au bord de dalle — bandes d'artefacts à chaque changement de tuile. `--edge-buffer N` (défaut 100 ; toujours appliqué par la génération depuis la carte, `EDGE_BUFFER_METERS` = 100 m dans `mapserve.py`, plus de case à cocher ; 0 = off en ligne de commande) fait rastériser le DTM sur la **dalle nominale 1 km alignée sur la grille** plus une bande de N m remplie avec les points sol des 8 LAZ voisines (`_neighbor_ground_points` dans `dtm.py` : PDAL en flux, découpe + filtre de classes IGN ; voisine absente = téléchargement automatique depuis le catalogue IGN avant le run, **isolée dans `input/edge_neighbors/`** pour ne pas gonfler le corpus des passes globales (dédupliqué sur tout le lot, `_fetch_edge_neighbors` dans `pipeline.py`) ; introuvable ou échec = bande vide). Les visualisations calculent sur l'emprise étendue puis `rendering.py` (`_core_tile_window`, via `tif_to_crop`/`tif_to_png`) recadre les sorties sur la dalle 1 km exacte lue dans le nom LHD — les AVIF restent des carrés 1 km alignés dans la mosaïque. Tampon consigné dans le tag GeoTIFF `LIDAR_EDGE_BUFFER` du DTM : changer `--edge-buffer` invalide le cache DTM automatiquement (tag absent = 0). Bandes voisines non calées par faisceaux (contexte seul, recadrée hors image finale). Coût : ~7 s de lecture par voisine + ~44 % de pixels en plus à 100 m/0,2 m. Nom hors pattern LHD : option ignorée (bornes d'en-tête, pas de recadrage).

docs/MAPS.md

@@ -16,16 +16,40 @@ docker compose -f docker-compose.maps.yml logs -f maps
 ./run.sh --serve-maps            # équivalent en conteneur au premier plan
 ```
 
+## Interface
+
+Un **panneau unique à onglets** (`web/map.{html,css,js}`) regroupe tous les
+réglages, à la place de l'ancienne pile de blocs empilés : **Affichage**
+(couche principale, mode relief/précision, fond de carte), **Export PDF**,
+**Génération** (masqué si le générateur est indisponible ou non autorisé
+pour ce navigateur) et **Partager**. Sur ordinateur, le panneau occupe une
+colonne fixe, repliable en bande d'icônes (**‹**, toujours utilisable : un
+clic sur un onglet redéplie le panneau). Sous 720 px de large, il devient un
+**volet en bas d'écran** à trois hauteurs — fermé, mi-hauteur, plein — que
+l'on change en glissant la poignée, en la touchant simplement (un cran), ou
+en retouchant l'onglet déjà actif.
+
+Un clic sur la carte ouvre une **bulle** sur la dalle LiDAR HD sous le
+curseur (date d'acquisition, densité de points sol) ; son bouton **Détails**
+bascule une **fiche complète** dans le panneau (emprise, IGN, recalage des
+passes) — **← Retour** rouvre l'onglet précédent. Raccourcis clavier :
+**1–4** (onglets du panneau, sans effet si l'onglet est masqué), **P**
+(mode d'affichage suivant), **Échap** (ferme la bulle, puis la fiche, puis
+replie le panneau). Deux thèmes clair/sombre (bouton ◐/☀/☾, `auto` par
+défaut = suit le système) ; réglages et état du panneau retenus dans
+`localStorage` (`lidarMapView_v2`, `lidar-print`, `lidar-panel`,
+`lidar-theme` — silencieusement ignorés en navigation privée).
+
 ## Générer des tuiles depuis la carte
 
-- **＋ Zone** (barre d'outils) — dessiner un rectangle : les dalles LHD de 1 km
+- **＋ Zone** (onglet Génération) — dessiner un rectangle : les dalles LHD de
-  qui l'intersectent sont téléchargées depuis la géoplateforme IGN puis
+  1 km qui l'intersectent sont téléchargées depuis la géoplateforme IGN puis
   traitées (0,2 m, options ci-dessous) ; les zones et clics successifs
   s'additionnent ;
 - **⤒ Compléter** — toutes les dalles déjà présentes dans `input/` qui
   manquent au moins une des couches demandées ;
-- **↻ Générer cette dalle** (fiche d'infos au clic) — une dalle précise,
+- **↻ Générer/Régénérer cette dalle** (bouton de la fiche de dalle, bulle →
-  même sans données existantes.
+  Détails) — une dalle précise, même sans données existantes.
 
 Options du run : couches visées (défaut : les couches du panneau) et
 régénération forcée. La classification du sol (IGN, sol seul) et le raccord
@@ -106,15 +130,18 @@ source manquée pendant qu'il était éteint est reprise au scan suivant.
 ## Fiche de dalle et rose des vents
 
 Un clic sur la carte sélectionne la dalle LiDAR HD sous le curseur (cadre
-jaune en pointillés), qu'elle soit générée ou non. La fiche s'affiche tout de
+jaune en pointillés) et ouvre une **bulle** — nom, date d'acquisition,
-suite avec le nom, l'emprise Lambert 93, puis, si la dalle est rendue, la
+densité de points sol — qu'elle soit générée ou non. Son bouton **Détails**
-résolution, la date de génération et le recalage vertical des passes
+bascule une fiche complète dans le panneau : nom, emprise Lambert 93, puis,
-(faisceaux, décalages, correction des lignes). Les informations IGN arrivent
+si la dalle est rendue, la résolution, la date de génération et le recalage
-à part (`GET /api/map/ign?col&row`, catalogue STAC mis en cache dans
+vertical des passes (faisceaux, décalages, correction des lignes). Les
-`output/ign_meta/`) : date et heure du scan LiDAR, capteurs, mission,
+informations IGN arrivent à part (`GET /api/map/ign?col&row`, catalogue STAC
-opérateur, date d'édition, procédé de classement, nombre de points et lien de
+mis en cache dans `output/ign_meta/`) : date et heure du scan LiDAR,
-téléchargement du nuage `.copc.laz` sur la géoplateforme. Un catalogue lent
+capteurs, mission, opérateur, date d'édition, procédé de classement, nombre
-ou injoignable n'empêche jamais la sélection.
+de points et lien de téléchargement du nuage `.copc.laz` sur la
+géoplateforme. Un catalogue lent ou injoignable n'empêche jamais la
+sélection. **← Retour** (ou Échap) ferme la fiche et rouvre l'onglet
+précédent.
 
 Quand le relief orienté est affiché, une rose des vents donne la couleur de
 chaque orientation de pente (même formule CIELAB que le rendu) ; la clarté
@@ -169,17 +196,19 @@ principale, mode validé, opacités bornées à 0–1. Un fichier de l'ancienne
 
 ## Export PDF (planche d'impression terrain)
 
-Le bouton **⎙** ouvre une carte de réglages (format A4/A3, paysage/portrait,
+L'onglet **Export PDF** affiche les réglages (format A4/A3, paysage/portrait,
-échelle 1:1 000 à 1:10 000, titre optionnel) et affiche un **cadre jaune en
+échelle 1:1 000 à 1:10 000, titre optionnel) et, tant qu'il est ouvert, un
-pointillés** qui montre la zone qui sera imprimée. Le cadre est **posé sur le
+**cadre jaune en pointillés** qui montre la zone qui sera imprimée (il
+disparaît en changeant d'onglet). Le cadre est **posé sur le
 terrain** : il se place au centre de la vue à l'ouverture (ou garde sa position
 précédente si elle est encore visible), la carte zoome pour le montrer en
-entier au-dessus des panneaux, puis on navigue librement sans qu'il bouge. On le
+entier au-dessus du panneau, puis on navigue librement sans qu'il bouge. On le
 déplace en faisant glisser sa **poignée ✥** (souris ou doigt) ; au relâchement,
-la géométrie Lambert 93 exacte est recalculée (`GET /api/export/frame`).
+la géométrie Lambert 93 exacte est recalculée (`GET /api/export/frame`) — un
-« ⌖ Centrer ici » le ramène au centre de la vue, « ⤢ Voir le cadre » zoome
+relâchement suivi d'un clic immédiat n'ouvre pas la bulle de dalle (garde de
-dessus ; changer de format, d'orientation ou d'échelle recadre la vue. Les
+300 ms). « ⌖ Centrer ici » le ramène au centre de la vue, « ⤢ Voir le cadre »
-réglages et la position du cadre sont mémorisés dans `localStorage` du
+zoome dessus ; changer de format, d'orientation ou d'échelle recadre la vue.
+Les réglages et la position du cadre sont mémorisés dans `localStorage` du
 navigateur. **Exporter le PDF** télécharge la planche (`GET
 /api/export/pdf`), nommée `relief_{x_km}_{y_km}_1-{échelle}.pdf` (centre
 Lambert 93 en km, à trois décimales).
@@ -420,5 +449,5 @@ consultable même worker éteint.
 
 - `lidar_pipeline/tiles.py` — grille, reprojection, cache, pré-chauffage.
 - `lidar_pipeline/mapserve.py` — routes tuiles/TileJSON/WMTS/JOSM et API carte.
-- `lidar_pipeline/mapui.py` — interface (Leaflet `L.tileLayer`, LOD natif).
+- `lidar_pipeline/web/map.{html,css,js}` — interface (panneau à onglets, Leaflet `L.tileLayer`, LOD natif), relue par `lidar_pipeline/mapui.py`.
 - `Dockerfile.maps`, `docker-compose.maps.yml`, `./run.sh --serve-maps`.