pygeospace 0.6.1__tar.gz

pygeospace-0.6.1/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 PyGeoSpace contributors
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

pygeospace-0.6.1/MANIFEST.in

@@ -0,0 +1,4 @@
+include README.md ROADMAP.md LICENSE pyproject.toml
+recursive-include pygeospace *.py
+recursive-include examples *.py
+recursive-include tests *.py

pygeospace-0.6.1/PKG-INFO

@@ -0,0 +1,249 @@
+Metadata-Version: 2.4
+Name: pygeospace
+Version: 0.6.1
+Summary: Powerful, honest, open-source geospatial visualization with true 3D.
+Author: PyGeoSpace contributors
+License: MIT
+Project-URL: Homepage, https://github.com/pygeospace/pygeospace
+Project-URL: Documentation, https://github.com/pygeospace/pygeospace#readme
+Project-URL: Roadmap, https://github.com/pygeospace/pygeospace/blob/main/ROADMAP.md
+Keywords: gis,geospatial,visualization,deck.gl,mapping
+Classifier: Development Status :: 4 - Beta
+Classifier: Intended Audience :: Science/Research
+Classifier: Intended Audience :: Developers
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Topic :: Scientific/Engineering :: GIS
+Classifier: Topic :: Scientific/Engineering :: Visualization
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: geopandas>=1.0
+Requires-Dist: shapely>=2.0
+Requires-Dist: pyproj>=3.5
+Requires-Dist: pydeck>=0.9
+Requires-Dist: matplotlib>=3.7
+Requires-Dist: pillow>=10.0
+Requires-Dist: scikit-learn>=1.3
+Requires-Dist: fastapi>=0.110
+Requires-Dist: uvicorn>=0.27
+Requires-Dist: click>=8.1
+Requires-Dist: pydantic>=2.0
+Requires-Dist: numpy>=1.24
+Requires-Dist: pandas>=2.0
+Requires-Dist: scipy>=1.10
+Provides-Extra: raster
+Requires-Dist: rasterio>=1.3; extra == "raster"
+Provides-Extra: 3d
+Requires-Dist: pyvista>=0.45; extra == "3d"
+Requires-Dist: meshio>=5.3; extra == "3d"
+Requires-Dist: trame>=3.0; extra == "3d"
+Requires-Dist: trame-vtk>=2.8; extra == "3d"
+Requires-Dist: trame-vuetify>=2.7; extra == "3d"
+Provides-Extra: pointcloud
+Requires-Dist: laspy[lazrs]>=2.5; extra == "pointcloud"
+Provides-Extra: streaming
+Requires-Dist: websockets>=12.0; extra == "streaming"
+Requires-Dist: paho-mqtt>=2.0; extra == "streaming"
+Provides-Extra: clustering-h3
+Requires-Dist: h3>=4.0; extra == "clustering-h3"
+Provides-Extra: network
+Requires-Dist: osmnx>=2.0; extra == "network"
+Provides-Extra: threed
+Requires-Dist: pyvista>=0.45; extra == "threed"
+Requires-Dist: meshio>=5.3; extra == "threed"
+Requires-Dist: trame>=3.0; extra == "threed"
+Requires-Dist: trame-vtk>=2.8; extra == "threed"
+Requires-Dist: trame-vuetify>=2.7; extra == "threed"
+Provides-Extra: distributed
+Requires-Dist: dask[complete]>=2024.1; extra == "distributed"
+Provides-Extra: all
+Requires-Dist: rasterio>=1.3; extra == "all"
+Requires-Dist: laspy[lazrs]>=2.5; extra == "all"
+Requires-Dist: websockets>=12.0; extra == "all"
+Requires-Dist: paho-mqtt>=2.0; extra == "all"
+Requires-Dist: h3>=4.0; extra == "all"
+Requires-Dist: pyvista>=0.45; extra == "all"
+Requires-Dist: meshio>=5.3; extra == "all"
+Requires-Dist: trame>=3.0; extra == "all"
+Requires-Dist: trame-vtk>=2.8; extra == "all"
+Requires-Dist: trame-vuetify>=2.7; extra == "all"
+Provides-Extra: dev
+Requires-Dist: pytest>=8.0; extra == "dev"
+Requires-Dist: pytest-cov>=4.1; extra == "dev"
+Requires-Dist: ruff>=0.4; extra == "dev"
+Dynamic: license-file
+
+# PyGeoSpace
+
+**Powerful, honest, open-source geospatial visualization for Python.**
+
+PyGeoSpace is a 2D / true-3D mapping library built on a real geospatial stack
+(GeoPandas, Shapely, pyproj, deck.gl, PyVista). It reads the common vector and
+raster formats, renders interactive 2D maps and **true 3D terrain, point clouds,
+and extruded geometry**, runs spatial + spectral analytics, and exports
+interactive HTML, static images, or 3D model files (glTF/STL/PLY/OBJ).
+
+This is the **0.6.0 (Beta)** release — "Advanced Visualization & True 3D". It is
+deliberately honest about its boundaries: every capability under "What works
+today" is implemented and covered by **81 passing tests**. Capabilities under
+"Roadmap" are **not silently stubbed** — calling them raises a clear error that
+names exactly what to install or that the feature is planned.
+
+---
+
+## Install
+
+```bash
+pip install pygeospace            # core: vector IO, analytics, deck.gl, web, CLI
+pip install "pygeospace[raster]"  # + GeoTIFF/COG/JPEG2000 reading (rasterio)
+pip install "pygeospace[3d]"      # + true 3D (PyVista, trame, meshio)
+pip install "pygeospace[streaming]"  # + WebSocket & MQTT clients
+pip install "pygeospace[all]"     # everything optional
+```
+
+Python 3.10+.
+
+## True 3D quickstart (0.6.0)
+
+```python
+import pygeospace as pgs
+
+# Terrain straight from an elevation GeoTIFF — interactive in the browser.
+m = pgs.Map(mode="3d")
+m.add_terrain("srtm_everest.tif", exaggeration=2.0, cmap="gist_earth")
+m.export_3d("terrain.html")     # interactive vtk.js page
+m.render_3d("terrain.png")      # static render
+m.export_3d("terrain.gltf")     # 3D model for Blender/Unity/printing
+
+# Classified LIDAR point cloud in 3D.
+pgs.Map(mode="3d").add_pointcloud_3d("forest.las").render_3d("canopy.png")
+
+# Extrude building footprints by a height attribute.
+pgs.Map(mode="3d").add_extruded("buildings.shp", "height", factor=1.0).export_3d("city.gltf")
+
+# Scripted camera flythrough.
+cam = pgs.Camera3D(position=(0, 0, 5000), focal_point=(0, 0, 0))
+cam.fly_to(4000, 4000, 2000, duration=3.0).orbit(45).tilt(20)
+m.render_3d("flyover.png", camera=cam)
+```
+
+```bash
+pygeospace 3d terrain elevation.tif --exaggeration 2.0 -o terrain.html
+pygeospace 3d pointcloud data.las --classify -o cloud.html
+pygeospace spectral sentinel.tif --indices ndvi,ndwi -o indices/
+pygeospace flow od_data.csv -o flow.html
+pygeospace export-gltf elevation.tif -o scene.gltf
+```
+
+## Quickstart (under 5 minutes)
+
+```python
+import pygeospace as pgs
+
+m = pgs.Map()
+
+# Read any vector format; style fluently; results chain.
+cities = m.add_layer("cities.geojson")
+cities.style(get_fill_color=[255, 90, 0, 200], get_radius=3000)
+
+# Classify an attribute into a choropleth (quantiles / equal_interval / jenks).
+from pygeospace.analytics import choropleth
+districts = m.add_layer("districts.shp")
+choropleth(districts, "population", method="jenks", k=5)
+
+# A geometry pipeline, returning new layers at each step.
+buffered = m.add_layer("rivers.gpkg").buffer(250)  # 250 m, in true meters
+
+m.fit().save("map.html")     # interactive, offline-capable deck.gl page
+m.save("map.png", dpi=300)   # static export
+```
+
+Command line:
+
+```bash
+pygeospace visualize cities.geojson -o map.html --style choropleth --attribute pop --method jenks
+pygeospace export cities.geojson -o map.png --dpi 300
+pygeospace serve cities.geojson --port 8000      # live preview in the browser
+pygeospace serve --api --port 8000               # REST CRUD API
+pygeospace info data.gpkg
+```
+
+## What works today (0.5.0)
+
+| Area | Capability | Status |
+|------|-----------|--------|
+| **IO** | Shapefile, GeoJSON, GeoPackage, KML, GPX, GML (via GeoPandas/OGR) | ✅ |
+| | CSV with coordinate columns | ✅ |
+| | PostGIS query → layer | ✅ |
+| | Format autodetect (extension + magic bytes) | ✅ |
+| | GeoTIFF / COG / JPEG2000 read | ✅ with `[raster]` |
+| | LAS/LAZ point cloud read | ✅ with `[pointcloud]` |
+| **Vector viz** | deck.gl scatter / GeoJSON / heatmap / hexagon-bin | ✅ |
+| | Choropleth: quantiles, equal-interval, **Jenks natural breaks** | ✅ |
+| | pandas-query attribute filtering (instant, no reload) | ✅ |
+| | CRS auto-reprojection to EPSG:4326 | ✅ |
+| **Raster** | Multi-band arrays, on-the-fly NDVI, hillshade | ✅ |
+| **Analytics** | Buffer (true meters via UTM), intersection, difference, dissolve | ✅ |
+| | KMeans / DBSCAN clustering, hex/grid binning (H3 if installed) | ✅ |
+| **3D (0.6.0)** | True 3D terrain from rasters (exaggeration, elevation colormap) | ✅ `[3d]` |
+| | 3D point clouds with LAS-classification coloring; polygon extrusion; cut planes | ✅ `[3d]` |
+| | Render to PNG + interactive HTML; export glTF/GLB/STL/PLY/OBJ/VTK | ✅ `[3d]` |
+| | `Camera3D`: fly_to / orbit / tilt / path recording | ✅ `[3d]` |
+| | Tilted (pitch/bearing) pseudo-3D deck.gl view | ✅ |
+| **Raster (0.6.0)** | Spectral indices NDVI/NDWI/NDBI/SAVI/EVI; band compositing; slope | ✅ |
+| | On-the-fly reprojection + multi-raster mosaicing | ✅ `[raster]` |
+| **Vector (0.6.0)** | Contour lines from scattered points; O-D flow maps (arcs); H3 hexbins | ✅ |
+| **Streaming** | WebSocket & MQTT clients; trail buffer; threshold alerts | ✅ with `[streaming]`¹ |
+| **Publishing** | Interactive HTML, static PNG/PDF, Jupyter inline (`_repr_html_`) | ✅ |
+| | FastAPI dev server + REST CRUD API | ✅ |
+| **Extensibility** | Decorator plugin system (data sources, renderers, tools) | ✅ |
+
+¹ The clients are real and connect to any reachable endpoint. Connecting to a
+specific public broker requires network access to that host.
+
+## Roadmap — what is *not* in 0.5.0
+
+These are tracked in [`ROADMAP.md`](ROADMAP.md). Where the public API touches
+them, it raises a clear, actionable error rather than pretending.
+
+| Planned for | Feature |
+|-------------|---------|
+| **0.6.0** | True PyVista 3D globe + terrain extrusion; raster↔vector opacity blending in the static renderer; vector-tile generation |
+| **0.7.0** | Getis-Ord Gi\* hotspot maps with significance shading; zonal statistics charts; OSMnx shortest-path routing |
+| **1.0.0** | WebGPU renderer for multi-million-point scenes; COG tile server with LRU disk cache; Dask-parallel distributed loading; flow-map & time-slider widgets as first-class UI |
+
+## Design notes
+
+- **Honesty over surface area.** A smaller set of things that genuinely work
+  beats a large set that breaks on first use.
+- **Real units.** Buffers reproject to the local UTM zone so "250 m" means 250
+  meters, not 250 degrees.
+- **Lazy optional deps.** Heavy/optional libraries (rasterio, laspy, paho-mqtt,
+  h3, pyvista) are imported only when used, with install hints on failure.
+
+## Development
+
+```bash
+pip install -e ".[dev]"
+pytest            # 43 tests
+ruff check .
+```
+
+## License
+
+MIT.
+
+## 3D & 2D Gallery
+
+Rendered headlessly by the test/example suite (see `examples/`):
+
+| | |
+|---|---|
+| ![terrain](docs/gallery/terrain.png) 3D terrain (elevation colormap) | ![mountain](docs/gallery/terrain_mountain.png) Terrain from GeoTIFF |
+| ![pointcloud](docs/gallery/pointcloud_classified.png) Classified point cloud | ![ndvi](docs/gallery/ndvi.png) NDVI spectral index |
+| ![true_color](docs/imagery/true_color.png) True color satellite imagery | ![false_color_ir](docs/imagery/false_color_ir.png) False color infrared (FCIR) |
+| ![agriculture](docs/imagery/agriculture.png) Agriculture/Vegetation analysis | ![nyc_map](docs/imagery/nyc_map.png) 2D NYC interactive map |
+

pygeospace-0.6.1/README.md

@@ -0,0 +1,171 @@
+# PyGeoSpace
+
+**Powerful, honest, open-source geospatial visualization for Python.**
+
+PyGeoSpace is a 2D / true-3D mapping library built on a real geospatial stack
+(GeoPandas, Shapely, pyproj, deck.gl, PyVista). It reads the common vector and
+raster formats, renders interactive 2D maps and **true 3D terrain, point clouds,
+and extruded geometry**, runs spatial + spectral analytics, and exports
+interactive HTML, static images, or 3D model files (glTF/STL/PLY/OBJ).
+
+This is the **0.6.0 (Beta)** release — "Advanced Visualization & True 3D". It is
+deliberately honest about its boundaries: every capability under "What works
+today" is implemented and covered by **81 passing tests**. Capabilities under
+"Roadmap" are **not silently stubbed** — calling them raises a clear error that
+names exactly what to install or that the feature is planned.
+
+---
+
+## Install
+
+```bash
+pip install pygeospace            # core: vector IO, analytics, deck.gl, web, CLI
+pip install "pygeospace[raster]"  # + GeoTIFF/COG/JPEG2000 reading (rasterio)
+pip install "pygeospace[3d]"      # + true 3D (PyVista, trame, meshio)
+pip install "pygeospace[streaming]"  # + WebSocket & MQTT clients
+pip install "pygeospace[all]"     # everything optional
+```
+
+Python 3.10+.
+
+## True 3D quickstart (0.6.0)
+
+```python
+import pygeospace as pgs
+
+# Terrain straight from an elevation GeoTIFF — interactive in the browser.
+m = pgs.Map(mode="3d")
+m.add_terrain("srtm_everest.tif", exaggeration=2.0, cmap="gist_earth")
+m.export_3d("terrain.html")     # interactive vtk.js page
+m.render_3d("terrain.png")      # static render
+m.export_3d("terrain.gltf")     # 3D model for Blender/Unity/printing
+
+# Classified LIDAR point cloud in 3D.
+pgs.Map(mode="3d").add_pointcloud_3d("forest.las").render_3d("canopy.png")
+
+# Extrude building footprints by a height attribute.
+pgs.Map(mode="3d").add_extruded("buildings.shp", "height", factor=1.0).export_3d("city.gltf")
+
+# Scripted camera flythrough.
+cam = pgs.Camera3D(position=(0, 0, 5000), focal_point=(0, 0, 0))
+cam.fly_to(4000, 4000, 2000, duration=3.0).orbit(45).tilt(20)
+m.render_3d("flyover.png", camera=cam)
+```
+
+```bash
+pygeospace 3d terrain elevation.tif --exaggeration 2.0 -o terrain.html
+pygeospace 3d pointcloud data.las --classify -o cloud.html
+pygeospace spectral sentinel.tif --indices ndvi,ndwi -o indices/
+pygeospace flow od_data.csv -o flow.html
+pygeospace export-gltf elevation.tif -o scene.gltf
+```
+
+## Quickstart (under 5 minutes)
+
+```python
+import pygeospace as pgs
+
+m = pgs.Map()
+
+# Read any vector format; style fluently; results chain.
+cities = m.add_layer("cities.geojson")
+cities.style(get_fill_color=[255, 90, 0, 200], get_radius=3000)
+
+# Classify an attribute into a choropleth (quantiles / equal_interval / jenks).
+from pygeospace.analytics import choropleth
+districts = m.add_layer("districts.shp")
+choropleth(districts, "population", method="jenks", k=5)
+
+# A geometry pipeline, returning new layers at each step.
+buffered = m.add_layer("rivers.gpkg").buffer(250)  # 250 m, in true meters
+
+m.fit().save("map.html")     # interactive, offline-capable deck.gl page
+m.save("map.png", dpi=300)   # static export
+```
+
+Command line:
+
+```bash
+pygeospace visualize cities.geojson -o map.html --style choropleth --attribute pop --method jenks
+pygeospace export cities.geojson -o map.png --dpi 300
+pygeospace serve cities.geojson --port 8000      # live preview in the browser
+pygeospace serve --api --port 8000               # REST CRUD API
+pygeospace info data.gpkg
+```
+
+## What works today (0.5.0)
+
+| Area | Capability | Status |
+|------|-----------|--------|
+| **IO** | Shapefile, GeoJSON, GeoPackage, KML, GPX, GML (via GeoPandas/OGR) | ✅ |
+| | CSV with coordinate columns | ✅ |
+| | PostGIS query → layer | ✅ |
+| | Format autodetect (extension + magic bytes) | ✅ |
+| | GeoTIFF / COG / JPEG2000 read | ✅ with `[raster]` |
+| | LAS/LAZ point cloud read | ✅ with `[pointcloud]` |
+| **Vector viz** | deck.gl scatter / GeoJSON / heatmap / hexagon-bin | ✅ |
+| | Choropleth: quantiles, equal-interval, **Jenks natural breaks** | ✅ |
+| | pandas-query attribute filtering (instant, no reload) | ✅ |
+| | CRS auto-reprojection to EPSG:4326 | ✅ |
+| **Raster** | Multi-band arrays, on-the-fly NDVI, hillshade | ✅ |
+| **Analytics** | Buffer (true meters via UTM), intersection, difference, dissolve | ✅ |
+| | KMeans / DBSCAN clustering, hex/grid binning (H3 if installed) | ✅ |
+| **3D (0.6.0)** | True 3D terrain from rasters (exaggeration, elevation colormap) | ✅ `[3d]` |
+| | 3D point clouds with LAS-classification coloring; polygon extrusion; cut planes | ✅ `[3d]` |
+| | Render to PNG + interactive HTML; export glTF/GLB/STL/PLY/OBJ/VTK | ✅ `[3d]` |
+| | `Camera3D`: fly_to / orbit / tilt / path recording | ✅ `[3d]` |
+| | Tilted (pitch/bearing) pseudo-3D deck.gl view | ✅ |
+| **Raster (0.6.0)** | Spectral indices NDVI/NDWI/NDBI/SAVI/EVI; band compositing; slope | ✅ |
+| | On-the-fly reprojection + multi-raster mosaicing | ✅ `[raster]` |
+| **Vector (0.6.0)** | Contour lines from scattered points; O-D flow maps (arcs); H3 hexbins | ✅ |
+| **Streaming** | WebSocket & MQTT clients; trail buffer; threshold alerts | ✅ with `[streaming]`¹ |
+| **Publishing** | Interactive HTML, static PNG/PDF, Jupyter inline (`_repr_html_`) | ✅ |
+| | FastAPI dev server + REST CRUD API | ✅ |
+| **Extensibility** | Decorator plugin system (data sources, renderers, tools) | ✅ |
+
+¹ The clients are real and connect to any reachable endpoint. Connecting to a
+specific public broker requires network access to that host.
+
+## Roadmap — what is *not* in 0.5.0
+
+These are tracked in [`ROADMAP.md`](ROADMAP.md). Where the public API touches
+them, it raises a clear, actionable error rather than pretending.
+
+| Planned for | Feature |
+|-------------|---------|
+| **0.6.0** | True PyVista 3D globe + terrain extrusion; raster↔vector opacity blending in the static renderer; vector-tile generation |
+| **0.7.0** | Getis-Ord Gi\* hotspot maps with significance shading; zonal statistics charts; OSMnx shortest-path routing |
+| **1.0.0** | WebGPU renderer for multi-million-point scenes; COG tile server with LRU disk cache; Dask-parallel distributed loading; flow-map & time-slider widgets as first-class UI |
+
+## Design notes
+
+- **Honesty over surface area.** A smaller set of things that genuinely work
+  beats a large set that breaks on first use.
+- **Real units.** Buffers reproject to the local UTM zone so "250 m" means 250
+  meters, not 250 degrees.
+- **Lazy optional deps.** Heavy/optional libraries (rasterio, laspy, paho-mqtt,
+  h3, pyvista) are imported only when used, with install hints on failure.
+
+## Development
+
+```bash
+pip install -e ".[dev]"
+pytest            # 43 tests
+ruff check .
+```
+
+## License
+
+MIT.
+
+## 3D & 2D Gallery
+
+Rendered headlessly by the test/example suite (see `examples/`):
+
+| | |
+|---|---|
+| ![terrain](docs/gallery/terrain.png) 3D terrain (elevation colormap) | ![mountain](docs/gallery/terrain_mountain.png) Terrain from GeoTIFF |
+| ![pointcloud](docs/gallery/pointcloud_classified.png) Classified point cloud | ![ndvi](docs/gallery/ndvi.png) NDVI spectral index |
+| ![true_color](docs/imagery/true_color.png) True color satellite imagery | ![false_color_ir](docs/imagery/false_color_ir.png) False color infrared (FCIR) |
+| ![agriculture](docs/imagery/agriculture.png) Agriculture/Vegetation analysis | ![nyc_map](docs/imagery/nyc_map.png) 2D NYC interactive map |
+

pygeospace-0.6.1/ROADMAP.md

@@ -0,0 +1,83 @@
+# PyGeoSpace Roadmap
+
+PyGeoSpace follows one rule: **ship what works, name what doesn't.** This file
+is the single source of truth for which features are real in each release. If a
+feature is listed as planned, the corresponding code path raises a clear error
+(often an `ImportError` naming an extra, or a message pointing here) instead of
+returning a fake result.
+
+## 0.5.0 — current (Beta)
+
+Implemented and tested:
+
+- Vector IO for Shapefile, GeoJSON, GeoPackage, KML, GPX, GML, CSV-with-coords,
+  PostGIS; format autodetection by extension and magic bytes.
+- Raster IO (GeoTIFF/COG/JPEG2000) via the `[raster]` extra.
+- LAS/LAZ point-cloud reading via the `[pointcloud]` extra.
+- deck.gl rendering: scatter, GeoJSON, heatmap, hexagon binning.
+- Choropleth classification: equal-interval, quantiles, Fisher-Jenks natural
+  breaks (implemented from scratch, no optional dependency).
+- pandas-query attribute filtering.
+- Spatial analytics: metric buffer (UTM reprojection), intersection, difference,
+  dissolve — all returning chainable layers.
+- Clustering: KMeans, DBSCAN, hex/grid binning (H3 via `[clustering-h3]`).
+- Raster analytics: NDVI, hillshade.
+- Streaming: WebSocket and MQTT clients with rolling trail buffers and threshold
+  alerts; message handling fully unit-tested without a network (`[streaming]`).
+- Publishing: interactive offline HTML, static PNG/PDF, Jupyter inline repr,
+  FastAPI dev server, REST CRUD API.
+- Plugin system: decorator registration of custom data sources, renderers, tools.
+- CLI: `visualize`, `export`, `serve`, `stream`, `info`.
+
+## 0.6.0 — current (Beta): Advanced Visualization & True 3D
+
+**Delivered and tested** (headless via PyVista's offscreen GL context):
+
+- True 3D engine (PyVista): terrain extrusion from elevation rasters with
+  vertical exaggeration and elevation colormapping; 3D point clouds with LAS
+  classification coloring; polygon extrusion to 3D blocks; 3D polylines; cut
+  planes; a graticuled globe (`mode="3d_globe"`).
+- 3D rendering to PNG, interactive HTML (vtk.js/trame), and export to
+  glTF / GLB / STL / PLY / OBJ / VTK.
+- `Camera3D` API: `fly_to` (with flythrough keyframes), `orbit`, `tilt`,
+  `get_position`, `save_path`.
+- Advanced raster: spectral indices (NDVI, NDWI, NDBI, SAVI, EVI), multi-band
+  true/false-color compositing, slope analysis, on-the-fly reprojection, and
+  multi-raster mosaicing (the last two via the `[raster]` extra).
+- Advanced vector: contour lines from scattered points (SciPy interpolation +
+  isolines), origin-destination flow maps (deck.gl arcs), H3 hex binning.
+- CLI: `3d terrain`, `3d pointcloud`, `spectral`, `flow`, `export-gltf`.
+- Test suite grown from 43 to 81.
+
+**Honestly deferred from 0.6.0** (these were in the original 0.6.0 wishlist but
+are not yet real, so they are not shipped as stubs):
+
+- Satellite-textured globe (needs bundled basemap imagery).
+- Animated flow particles / time-series **GIF/MP4** export (needs an offscreen
+  animation+encoding pipeline that isn't validated yet) → 0.7.0.
+- Interactive browser widgets: drag-drop layer panel, colormap editor, measure
+  tool, draw/annotate, lasso brushing (need a full JS frontend) → 0.7.0.
+- Synchronized split 2D/3D view → 0.7.0.
+- LOD / tile-based 3D streaming / octree billion-point clouds / explicit
+  GPU-FPS guarantees → 1.0.0 (research-grade systems work).
+- Proportional pie/bar chart markers, curved-label placement engine → 0.7.0.
+
+## 0.7.0
+
+- Getis-Ord Gi\* hotspot detection with significance visualization.
+- Zonal statistics (raster summarized under vector polygons) with chart output.
+- Network analysis: shortest path on OSM data via OSMnx.
+- Animated flow maps and origin-destination rendering.
+
+## 1.0.0
+
+- WebGPU renderer targeting multi-million-point interactivity.
+- COG-native tiling with an LRU disk tile cache.
+- Dask-parallel distributed data loading.
+- First-class interactive widgets: time slider, draggable legends, color-scale
+  editors, layer toggles.
+
+## Explicitly out of scope (for now)
+
+AR/VR, a full Google Earth Engine equivalent, and a hosted SaaS. These are not
+on a dated milestone and will not be stubbed.

pygeospace-0.6.1/examples/crop_health.py

@@ -0,0 +1,58 @@
+"""Crop Health (spectral) — runnable with synthetic Sentinel-like bands.
+
+    python examples/crop_health.py
+
+Computes NDVI, NDWI, and a false-color composite from a synthetic 5-band raster
+and saves colorized PNGs. Swap the synthetic array for
+`read_file("sentinel.tif")` to analyze real imagery.
+"""
+
+from __future__ import annotations
+
+import numpy as np
+import matplotlib
+
+matplotlib.use("Agg")
+import matplotlib.pyplot as plt
+
+from pygeospace.core.layer import RasterLayer
+from pygeospace.analytics.spectral import spectral_index, composite
+
+
+def synthetic_scene() -> RasterLayer:
+    """5 bands (blue, green, red, nir, swir). A vegetated patch + a water body."""
+    rng = np.random.default_rng(11)
+    h = w = 128
+    veg = np.zeros((h, w), bool)
+    veg[20:90, 30:110] = True
+    water = np.zeros((h, w), bool)
+    water[95:120, 10:50] = True
+
+    blue = rng.uniform(0.05, 0.15, (h, w))
+    green = rng.uniform(0.08, 0.2, (h, w))
+    red = np.where(veg, 0.08, 0.25) + rng.uniform(0, 0.03, (h, w))
+    nir = np.where(veg, 0.55, 0.2)
+    nir = np.where(water, 0.05, nir) + rng.uniform(0, 0.03, (h, w))
+    swir = np.where(water, 0.04, 0.3) + rng.uniform(0, 0.03, (h, w))
+    arr = np.stack([blue, green, red, nir, swir]).astype("float32")
+    return RasterLayer(arr, transform=None, crs="EPSG:4326", name="field")
+
+
+def main() -> None:
+    scene = synthetic_scene()
+
+    ndvi = spectral_index(scene, "ndvi")
+    ndwi = spectral_index(scene, "ndwi")
+    plt.imsave("crop_ndvi.png", ndvi.array[0], cmap="RdYlGn", vmin=-1, vmax=1)
+    plt.imsave("crop_ndwi.png", ndwi.array[0], cmap="Blues", vmin=-1, vmax=1)
+
+    false_color = composite(scene, [4, 3, 2])  # NIR-R-G highlights vegetation red
+    plt.imsave("crop_falsecolor.png", np.transpose(false_color.array, (1, 2, 0)))
+
+    veg_frac = float((ndvi.array[0] > 0.3).mean())
+    print(f"NDVI mean {ndvi.array.mean():.3f}; vegetated fraction {veg_frac:.0%}")
+    print("Wrote crop_ndvi.png, crop_ndwi.png, crop_falsecolor.png")
+
+
+if __name__ == "__main__":
+    main()

pygeospace-0.6.1/examples/quickstart.py

@@ -0,0 +1,64 @@
+"""Runnable end-to-end example. No external files needed.
+
+    python examples/quickstart.py
+
+Produces quickstart_map.html (interactive) and quickstart_map.png (static) in
+the current directory, and demonstrates the chaining API, choropleth
+classification, clustering, and offline streaming logic.
+"""
+
+from __future__ import annotations
+
+import numpy as np
+import geopandas as gpd
+
+import pygeospace as pgs
+from pygeospace.analytics import choropleth, kmeans
+from pygeospace.core.layer import StreamLayer
+from pygeospace.streaming.source import BaseSource
+
+
+def make_data() -> gpd.GeoDataFrame:
+    rng = np.random.default_rng(7)
+    n = 500
+    lon = rng.normal(-0.18, 0.05, n)  # around Accra
+    lat = rng.normal(5.60, 0.04, n)
+    return gpd.GeoDataFrame(
+        {"incidents": rng.integers(1, 100, n)},
+        geometry=gpd.points_from_xy(lon, lat),
+        crs="EPSG:4326",
+    )
+
+
+def main() -> None:
+    gdf = make_data()
+
+    # 1) Build a map, classify a choropleth, and cluster.
+    m = pgs.Map(title="Accra incidents")
+    layer = m.add_vector(gdf, name="incidents")
+    choropleth(layer, "incidents", method="jenks", k=5)
+
+    clustered = kmeans(layer, n_clusters=6)
+    print("cluster counts:\n", clustered.data["cluster"].value_counts())
+
+    # 2) Chaining: buffer hotspots and add as a new red layer.
+    hot = layer.filter("incidents > 80").buffer(300).style(get_fill_color=[220, 30, 30, 90])
+    m.add_layer(hot)
+
+    m.fit().save("quickstart_map.html")
+    m.save("quickstart_map.png", dpi=150)
+    print("Wrote quickstart_map.html and quickstart_map.png")
+
+    # 3) Streaming logic (no network): feed synthetic messages.
+    stream = StreamLayer(name="buses", trail_length=10)
+    stream.add_alert("speed", 80, ">")
+    src = BaseSource(stream)
+    for speed in (30, 50, 95, 40):
+        alerts = src.handle_raw({"id": "bus_1", "lon": -0.18, "lat": 5.6, "speed": speed})
+        for a in alerts:
+            print(a)
+    print(f"stream received {src.received}, buffered {len(stream)}")
+
+
+if __name__ == "__main__":
+    main()