@pikaa-ai/pikaa 0.3.23 → 0.3.24

package/skills/geomaster/SKILL.md

@@ -1,366 +0,0 @@
----
-name: geomaster
-description: Comprehensive geospatial science skill covering remote sensing, GIS, spatial analysis, machine learning for earth observation, and 30+ scientific domains. Supports satellite imagery processing (Sentinel, Landsat, MODIS, SAR, hyperspectral), vector and raster data operations, spatial statistics, point cloud processing, network analysis, cloud-native workflows (STAC, COG, Planetary Computer), and 8 programming languages (Python, R, Julia, JavaScript, C++, Java, Go, Rust) with 500+ code examples. Use for remote sensing workflows, GIS analysis, spatial ML, Earth observation data processing, terrain analysis, hydrological modeling, marine spatial analysis, atmospheric science, and any geospatial computation task.
-license: MIT License
-metadata:
-  version: "1.1"
-  skill-author: K-Dense Inc.
----
-
-# GeoMaster
-
-Comprehensive geospatial science skill covering GIS, remote sensing, spatial analysis, and ML for Earth observation across 70+ topics with 500+ code examples in 8 programming languages.
-
-## Installation
-
-```bash
-# Core Python stack (conda recommended)
-conda install -c conda-forge gdal rasterio fiona shapely pyproj geopandas
-
-# Remote sensing & ML
-uv pip install rsgislib torchgeo earthengine-api
-uv pip install scikit-learn xgboost torch-geometric
-
-# Network & visualization
-uv pip install osmnx networkx folium keplergl
-uv pip install cartopy contextily mapclassify
-
-# Big data & cloud
-uv pip install xarray rioxarray dask-geopandas
-uv pip install pystac-client planetary-computer
-
-# Point clouds
-uv pip install laspy pylas open3d pdal
-
-# Databases
-conda install -c conda-forge postgis spatialite
-```
-
-## Quick Start
-
-### NDVI from Sentinel-2
-
-```python
-import rasterio
-import numpy as np
-
-with rasterio.open('sentinel2.tif') as src:
-    red = src.read(4).astype(float)   # B04
-    nir = src.read(8).astype(float)   # B08
-    ndvi = (nir - red) / (nir + red + 1e-8)
-    ndvi = np.nan_to_num(ndvi, nan=0)
-
-    profile = src.profile
-    profile.update(count=1, dtype=rasterio.float32)
-
-    with rasterio.open('ndvi.tif', 'w', **profile) as dst:
-        dst.write(ndvi.astype(rasterio.float32), 1)
-```
-
-### Spatial Analysis with GeoPandas
-
-```python
-import geopandas as gpd
-
-# Load and ensure same CRS
-zones = gpd.read_file('zones.geojson')
-points = gpd.read_file('points.geojson')
-
-if zones.crs != points.crs:
-    points = points.to_crs(zones.crs)
-
-# Spatial join and statistics
-joined = gpd.sjoin(points, zones, how='inner', predicate='within')
-stats = joined.groupby('zone_id').agg({
-    'value': ['count', 'mean', 'std', 'min', 'max']
-}).round(2)
-```
-
-### Google Earth Engine Time Series
-
-```python
-import ee
-import pandas as pd
-
-ee.Initialize(project='your-project')
-roi = ee.Geometry.Point([-122.4, 37.7]).buffer(10000)
-
-s2 = (ee.ImageCollection('COPERNICUS/S2_SR_HARMONIZED')
-      .filterBounds(roi)
-      .filterDate('2020-01-01', '2023-12-31')
-      .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20)))
-
-def add_ndvi(img):
-    return img.addBands(img.normalizedDifference(['B8', 'B4']).rename('NDVI'))
-
-s2_ndvi = s2.map(add_ndvi)
-
-def extract_series(image):
-    stats = image.reduceRegion(ee.Reducer.mean(), roi.centroid(), scale=10, maxPixels=1e9)
-    return ee.Feature(None, {'date': image.date().format('YYYY-MM-dd'), 'ndvi': stats.get('NDVI')})
-
-series = s2_ndvi.map(extract_series).getInfo()
-df = pd.DataFrame([f['properties'] for f in series['features']])
-df['date'] = pd.to_datetime(df['date'])
-```
-
-## Core Concepts
-
-### Data Types
-
-| Type | Examples | Libraries |
-|------|----------|-----------|
-| Vector | Shapefile, GeoJSON, GeoPackage | GeoPandas, Fiona, GDAL |
-| Raster | GeoTIFF, NetCDF, COG | Rasterio, Xarray, GDAL |
-| Point Cloud | LAS, LAZ | Laspy, PDAL, Open3D |
-
-### Coordinate Systems
-
-- **EPSG:4326** (WGS 84) - Geographic, lat/lon, use for storage
-- **EPSG:3857** (Web Mercator) - Web maps only (don't use for area/distance!)
-- **EPSG:326xx/327xx** (UTM) - Metric calculations, <1% distortion per zone
-- Use `gdf.estimate_utm_crs()` for automatic UTM detection
-
-```python
-# Always check CRS before operations
-assert gdf1.crs == gdf2.crs, "CRS mismatch!"
-
-# For area/distance calculations, use projected CRS
-gdf_metric = gdf.to_crs(gdf.estimate_utm_crs())
-area_sqm = gdf_metric.geometry.area
-```
-
-### OGC Standards
-
-- **WMS**: Web Map Service - raster maps
-- **WFS**: Web Feature Service - vector data
-- **WCS**: Web Coverage Service - raster coverage
-- **STAC**: Spatiotemporal Asset Catalog - modern metadata
-
-## Common Operations
-
-### Spectral Indices
-
-```python
-def calculate_indices(image_path):
-    """NDVI, EVI, SAVI, NDWI from Sentinel-2."""
-    with rasterio.open(image_path) as src:
-        B02, B03, B04, B08, B11 = [src.read(i).astype(float) for i in [1,2,3,4,5]]
-
-    ndvi = (B08 - B04) / (B08 + B04 + 1e-8)
-    evi = 2.5 * (B08 - B04) / (B08 + 6*B04 - 7.5*B02 + 1)
-    savi = ((B08 - B04) / (B08 + B04 + 0.5)) * 1.5
-    ndwi = (B03 - B08) / (B03 + B08 + 1e-8)
-
-    return {'NDVI': ndvi, 'EVI': evi, 'SAVI': savi, 'NDWI': ndwi}
-```
-
-### Vector Operations
-
-```python
-# Buffer (use projected CRS!)
-gdf_proj = gdf.to_crs(gdf.estimate_utm_crs())
-gdf['buffer_1km'] = gdf_proj.geometry.buffer(1000)
-
-# Spatial relationships
-intersects = gdf[gdf.geometry.intersects(other_geometry)]
-contains = gdf[gdf.geometry.contains(point_geometry)]
-
-# Geometric operations
-gdf['centroid'] = gdf.geometry.centroid
-gdf['simplified'] = gdf.geometry.simplify(tolerance=0.001)
-
-# Overlay operations
-intersection = gpd.overlay(gdf1, gdf2, how='intersection')
-union = gpd.overlay(gdf1, gdf2, how='union')
-```
-
-### Terrain Analysis
-
-```python
-def terrain_metrics(dem_path):
-    """Calculate slope, aspect, hillshade from DEM."""
-    with rasterio.open(dem_path) as src:
-        dem = src.read(1)
-
-    dy, dx = np.gradient(dem)
-    slope = np.arctan(np.sqrt(dx**2 + dy**2)) * 180 / np.pi
-    aspect = (90 - np.arctan2(-dy, dx) * 180 / np.pi) % 360
-
-    # Hillshade
-    az_rad, alt_rad = np.radians(315), np.radians(45)
-    hillshade = (np.sin(alt_rad) * np.sin(np.radians(slope)) +
-                 np.cos(alt_rad) * np.cos(np.radians(slope)) *
-                 np.cos(np.radians(aspect) - az_rad))
-
-    return slope, aspect, hillshade
-```
-
-### Network Analysis
-
-```python
-import osmnx as ox
-import networkx as nx
-
-# Download and analyze street network
-G = ox.graph_from_place('San Francisco, CA', network_type='drive')
-G = ox.add_edge_speeds(G).add_edge_travel_times(G)
-
-# Shortest path
-orig = ox.distance.nearest_nodes(G, -122.4, 37.7)
-dest = ox.distance.nearest_nodes(G, -122.3, 37.8)
-route = nx.shortest_path(G, orig, dest, weight='travel_time')
-```
-
-## Image Classification
-
-```python
-from sklearn.ensemble import RandomForestClassifier
-import rasterio
-from rasterio.features import rasterize
-
-def classify_imagery(raster_path, training_gdf, output_path):
-    """Train RF and classify imagery."""
-    with rasterio.open(raster_path) as src:
-        image = src.read()
-        profile = src.profile
-        transform = src.transform
-
-    # Extract training data
-    X_train, y_train = [], []
-    for _, row in training_gdf.iterrows():
-        mask = rasterize([(row.geometry, 1)],
-                        out_shape=(profile['height'], profile['width']),
-                        transform=transform, fill=0, dtype=np.uint8)
-        pixels = image[:, mask > 0].T
-        X_train.extend(pixels)
-        y_train.extend([row['class_id']] * len(pixels))
-
-    # Train and predict
-    rf = RandomForestClassifier(n_estimators=100, max_depth=20, n_jobs=-1)
-    rf.fit(X_train, y_train)
-
-    prediction = rf.predict(image.reshape(image.shape[0], -1).T)
-    prediction = prediction.reshape(profile['height'], profile['width'])
-
-    profile.update(dtype=rasterio.uint8, count=1)
-    with rasterio.open(output_path, 'w', **profile) as dst:
-        dst.write(prediction.astype(rasterio.uint8), 1)
-
-    return rf
-```
-
-## Modern Cloud-Native Workflows
-
-### STAC + Planetary Computer
-
-```python
-import pystac_client
-import planetary_computer
-import odc.stac
-
-# Search Sentinel-2 via STAC
-catalog = pystac_client.Client.open(
-    "https://planetarycomputer.microsoft.com/api/stac/v1",
-    modifier=planetary_computer.sign_inplace,
-)
-
-search = catalog.search(
-    collections=["sentinel-2-l2a"],
-    bbox=[-122.5, 37.7, -122.3, 37.9],
-    datetime="2023-01-01/2023-12-31",
-    query={"eo:cloud_cover": {"lt": 20}},
-)
-
-# Load as xarray (cloud-native!)
-data = odc.stac.load(
-    list(search.get_items())[:5],
-    bands=["B02", "B03", "B04", "B08"],
-    crs="EPSG:32610",
-    resolution=10,
-)
-
-# Calculate NDVI on xarray
-ndvi = (data.B08 - data.B04) / (data.B08 + data.B04)
-```
-
-### Cloud-Optimized GeoTIFF (COG)
-
-```python
-import rasterio
-from rasterio.session import AWSSession
-
-# Read COG directly from cloud (partial reads)
-session = AWSSession(aws_access_key_id=..., aws_secret_access_key=...)
-with rasterio.open('s3://bucket/path.tif', session=session) as src:
-    # Read only window of interest
-    window = ((1000, 2000), (1000, 2000))
-    subset = src.read(1, window=window)
-
-# Write COG
-with rasterio.open('output.tif', 'w', **profile,
-                   tiled=True, blockxsize=256, blockysize=256,
-                   compress='DEFLATE', predictor=2) as dst:
-    dst.write(data)
-
-# Validate COG
-from rio_cogeo.cogeo import cog_validate
-cog_validate('output.tif')
-```
-
-## Performance Tips
-
-```python
-# 1. Spatial indexing (10-100x faster queries)
-gdf.sindex  # Auto-created by GeoPandas
-
-# 2. Chunk large rasters
-with rasterio.open('large.tif') as src:
-    for i, window in src.block_windows(1):
-        block = src.read(1, window=window)
-
-# 3. Dask for big data
-import dask.array as da
-dask_array = da.from_rasterio('large.tif', chunks=(1, 1024, 1024))
-
-# 4. Use Arrow for I/O
-gdf.to_file('output.gpkg', use_arrow=True)
-
-# 5. GDAL caching
-from osgeo import gdal
-gdal.SetCacheMax(2**30)  # 1GB cache
-
-# 6. Parallel processing
-rf = RandomForestClassifier(n_jobs=-1)  # All cores
-```
-
-## Best Practices
-
-1. **Always check CRS** before spatial operations
-2. **Use projected CRS** for area/distance calculations
-3. **Validate geometries**: `gdf = gdf[gdf.is_valid]`
-4. **Handle missing data**: `gdf['geometry'] = gdf['geometry'].fillna(None)`
-5. **Use efficient formats**: GeoPackage > Shapefile, Parquet for large data
-6. **Apply cloud masking** to optical imagery
-7. **Preserve lineage** for reproducible research
-8. **Use appropriate resolution** for your analysis scale
-
-## Detailed Documentation
-
-- **[Coordinate Systems](references/coordinate-systems.md)** - CRS fundamentals, UTM, transformations
-- **[Core Libraries](references/core-libraries.md)** - GDAL, Rasterio, GeoPandas, Shapely
-- **[Remote Sensing](references/remote-sensing.md)** - Satellite missions, spectral indices, SAR
-- **[Machine Learning](references/machine-learning.md)** - Deep learning, CNNs, GNNs for RS
-- **[GIS Software](references/gis-software.md)** - QGIS, ArcGIS, GRASS integration
-- **[Scientific Domains](references/scientific-domains.md)** - Marine, hydrology, agriculture, forestry
-- **[Advanced GIS](references/advanced-gis.md)** - 3D GIS, spatiotemporal, topology
-- **[Big Data](references/big-data.md)** - Distributed processing, GPU acceleration
-- **[Industry Applications](references/industry-applications.md)** - Urban planning, disaster management
-- **[Programming Languages](references/programming-languages.md)** - Python, R, Julia, JS, C++, Java, Go, Rust
-- **[Data Sources](references/data-sources.md)** - Satellite catalogs, APIs
-- **[Troubleshooting](references/troubleshooting.md)** - Common issues, debugging, error reference
-- **[Code Examples](references/code-examples.md)** - 500+ examples
-
----
-
-**GeoMaster covers everything from basic GIS operations to advanced remote sensing and machine learning.**

package/skills/geopandas/SKILL.md

@@ -1,250 +0,0 @@
----
-name: geopandas
-description: Guidance and local audit tools for Python workflows that directly use GeoPandas GeoSeries, GeoDataFrame, spatial operations, or vector-data I/O.
-license: MIT
-compatibility: Requires Python 3.10+ and uv. Bundled CLIs are local-only; runtime analysis requires the pinned GeoPandas stack below.
-allowed-tools: Read Write Bash Glob Grep
-metadata:
-  version: "1.1"
-  skill-author: K-Dense Inc.
-  last-reviewed: "2026-07-23"
----
-
-# GeoPandas
-
-Use GeoPandas for planar vector data represented as pandas-like `GeoSeries` and
-`GeoDataFrame` objects. This skill targets stable **GeoPandas 1.1.4** (released
-2026-06-26), not the unreleased 1.2 documentation.
-
-## Reproducible environment
-
-GeoPandas 1.1.4 requires Python 3.10+; its tagged source requires NumPy >=1.24,
-pandas >=2.0, Shapely >=2.0, pyproj >=3.5, pyogrio >=0.7.2, and `packaging`.
-This exact Python 3.12 snapshot was smoke-tested on 2026-07-23:
-
-```bash
-uv venv --python 3.12
-uv pip install \
-  "geopandas==1.1.4" \
-  "numpy==2.5.1" \
-  "pandas==3.0.5" \
-  "shapely==2.1.2" \
-  "pyproj==3.7.2" \
-  "pyogrio==0.13.0" \
-  "pyarrow==25.0.0" \
-  "packaging==26.2"
-```
-
-Keep optional plotting and PostGIS packages pinned in the project lock as well.
-Do not mix binary geospatial packages from incompatible package channels.
-
-## Safety and privacy contract
-
-- Treat exact coordinates, addresses, parcel boundaries, trajectories, and
-  small-area joins as sensitive. Default reports to counts, categories, coarse
-  extents, and redacted identifiers. Generalize before publication.
-- Never automatically load a URL, cloud URI, GDAL `/vsi*` path, archive, or
-  geocode an address. Obtain explicit approval, validate provenance and hashes,
-  then stage an unpacked local file in an isolated workspace.
-- GDAL/OGR drivers, GEOS, PROJ, pyogrio, Shapely, pyproj, and their wheels are a
-  native-code trust boundary. Prefer official wheels/conda-forge, record native
-  versions, restrict drivers, and process untrusted data in a sandbox.
-- Do not open macro-enabled office files or nested archives through permissive
-  GDAL drivers. The bundled CLIs use an extension allowlist and reject archives.
-- Read only named database secrets such as `GEOPANDAS_POSTGIS_PASSWORD`; use a
-  secret manager or scoped environment variable. Never embed a password in a
-  URL or source, print an engine/URL, or dump the environment.
-- Every derived artifact needs source hashes/versions, CRS, operation parameters,
-  predicate, join cardinality, precision/repair choices, and row-count checks.
-
-## Correctness gates
-
-Apply these gates before trusting a result:
-
-1. **Identity and provenance** — identify the source layer, stable feature key,
-   duplicate IDs, row count, geometry column, parser/driver, and content hash.
-2. **Geometry state** — count null, empty, invalid, mixed, Z/M, and collapsed
-   geometries separately. `None` is missing; an empty Shapely geometry is real.
-3. **CRS semantics** — require CRS metadata. `set_crs()` assigns metadata;
-   `to_crs()` transforms coordinates. Never guess a CRS from coordinate ranges.
-4. **Units and operation** — GeoPandas is planar. Geographic coordinates are
-   angular; do not use them directly for buffer, distance, area, nearest joins,
-   precision grids, or tolerances. Choose a fit-for-purpose local/equal-area CRS
-   or a geodesic method.
-5. **Transform quality** — inspect axis order, area of use, datum pipeline,
-   expected accuracy, ballpark status, and missing grids. Keep PROJ network
-   disabled unless the user explicitly approves grid retrieval.
-6. **Topology and precision** — validate before and after repair/overlay. Pick a
-   precision grid from source accuracy and CRS units; arbitrary snapping can
-   collapse features or create bias.
-7. **Cardinality** — state expected one-to-one, one-to-many, or many-to-many
-   behavior before `merge`, `sjoin`, or `sjoin_nearest`; audit unmatched and
-   multiplied rows afterward.
-8. **Output contract** — use a new output path, preserve a stable feature ID,
-   document schema/CRS/encoding, reopen the artifact, and compare counts/types.
-
-## CRS and antimeridian rules
-
-GeoPandas stores CRS as `pyproj.CRS`. Coordinate arrays use traditional GIS
-`(x, y)` order, while authority definitions can advertise latitude-first axes.
-Use `Transformer(..., always_xy=True)` for explicit coordinate-array pipelines,
-and record that choice.
-
-`to_crs()` transforms vertices and assumes each segment is straight in the
-source CRS; it does not transform geodesic arcs. Geometries crossing ±180° or a
-projection boundary can be badly wrapped. Detect crossings, split/unwrap and
-densify in a documented geographic representation, transform parts, then
-validate. Do not use Web Mercator as a general measurement CRS.
-
-```python
-crs = gdf.crs  # a pyproj.CRS when present
-if crs is None or crs.is_geographic:
-    raise ValueError("Choose a justified projected CRS before planar measurement")
-
-unit_names = [axis.unit_name for axis in crs.axis_info]
-areas = gdf.geometry.area  # square CRS units, not automatically square metres
-```
-
-See [CRS management](references/crs-management.md).
-
-## Core API decisions
-
-### Data structures
-
-- A `GeoDataFrame` can hold multiple geometry columns, each with CRS metadata,
-  but only `active_geometry_name` drives frame-level spatial operations.
-- Binary `GeoSeries` methods are row-wise and align by index by default. Use
-  `align=False` only when positional pairing is explicitly intended and lengths
-  and order were verified.
-- Duplicate column names and duplicate feature IDs are ambiguous; reject or
-  resolve them before joins and exports.
-
-See [data structures](references/data-structures.md).
-
-### Geometry validity, precision, and union
-
-Use `is_valid` and redacted `is_valid_reason()` categories before
-`make_valid(method="linework"|"structure", keep_collapsed=...)`. Repair can
-change geometry type or dimension; retain the original and compare counts,
-area, types, empties, and collapsed parts.
-
-`set_precision(grid_size, mode=...)` uses **CRS units** and may remove duplicate
-vertices or collapse features. `union_all(method="unary", grid_size=...)` is the
-robust default. Use `coverage` only after `is_valid_coverage()` proves
-non-overlap and edge matching; use `disjoint_subset` with Shapely >=2.1 when its
-partitioning assumption is useful.
-
-See [geometric operations](references/geometric-operations.md).
-
-### Joins, overlay, clip, and dissolve
-
-- `sjoin` predicates are directional: `left.within(right)` is not
-  `left.contains(right)`. `intersects` includes boundary contact; `contains`
-  excludes boundary-only points, while `covers` includes boundary points.
-- `predicate="dwithin"` requires `distance`; scalar or per-left-row distances
-  are in CRS units. `sjoin_nearest` returns all equidistant nearest matches and
-  does **not** implement a `k=` parameter.
-- `overlay(..., make_valid=True)` repairs invalid input but can change types;
-  `keep_geom_type=None` drops other types with a warning. Precision mismatch can
-  create slivers; quantify them rather than silently deleting them.
-- `clip` dissolves the mask. Rectangle clipping is fast but possibly dirty and
-  may omit a line collapsed to a point; validate its output.
-- `dissolve` combines `groupby.agg` with `union_all`; choose explicit attribute
-  aggregations and audit null group keys.
-
-See [spatial analysis](references/spatial-analysis.md).
-
-### I/O, Arrow, and PostGIS
-
-GeoPandas 1.x defaults to pyogrio. Driver availability and semantics come from
-the installed GDAL, not GeoPandas alone. Prefer local GeoPackage for general
-interchange and WKB GeoParquet for columnar interoperability.
-
-GeoParquet defaults to stable schema 1.0.0. Native GeoArrow encodings and bbox
-covering require schema 1.1.0 and remain less interoperable. A missing GeoParquet
-`crs` key means `OGC:CRS84`; explicit `crs: null` means unknown—do not conflate
-them. Reopen and validate every export.
-
-Use parameterized SQL and a SQLAlchemy `Engine`/`Connection` for PostGIS.
-`if_exists="replace"` is destructive; default to `"fail"` and use a transaction.
-
-See [data I/O](references/data-io.md).
-
-## Migration checklist
-
-For code moving from GeoPandas 0.14 or earlier:
-
-- GeoPandas 1.0 supports Shapely >=2 only; PyGEOS, Shapely <2, and the rtree
-  spatial-index backend were removed.
-- pyogrio replaced Fiona as the installed/default I/O engine. Set `engine=`
-  explicitly and test schema, empty, datetime, encoding, and append behavior.
-- Replace `sjoin(op=...)` with `predicate=`, `sindex.query_bulk()` with
-  `sindex.query()`, `unary_union` with `union_all()`, and
-  `GeometryArray.data` with `to_numpy()`/`np.asarray`.
-- Replace `read_file(include_fields=...|ignore_fields=...)` with `columns=`.
-  Use `schema_version=`, not the removed GeoParquet `version=` compatibility.
-- Do not use removed `geopandas.datasets`, internal `geopandas.io.*` entry
-  points, plot `axes`/`colormap`, or set-operation operators.
-- `explode()` now defaults `index_parts=False`; a named Series passed to
-  `set_geometry()` supplies the new active-column name; a named right index can
-  replace `index_right` in `sjoin` output.
-- Do not assign `.crs` to override metadata or rely on deprecated
-  `set_geometry(drop=...)`; use explicit `set_crs()` and rename/drop steps.
-- GeoPandas 1.1 requires Python >=3.10, pandas >=2.0, NumPy >=1.24, and pyproj
-  >=3.5. Version 1.1.2 fixed SQL injection through a PostGIS geometry-column
-  name; the pinned 1.1.4 includes that fix.
-
-### Plotting and exploration
-
-Maps are analytical outputs: label units, classification method, missing data,
-normalization denominator, and date. `explore()` can expose every attribute in
-tooltips/popups and contact tile/CDN servers; generalize first and use
-`tiles=None`, `tooltip=False`, and `popup=False` for a local draft.
-
-See [visualization](references/visualization.md).
-
-## Bundled local CLIs
-
-All helpers are deterministic, reject network/archive paths, bound input bytes
-and feature counts, keep imports lazy so `--help` is dependency-free, and emit
-JSON without coordinates or record identifiers.
-
-| CLI | Purpose |
-|---|---|
-| `scripts/vector_inventory.py` | Redacted local vector/GeoParquet technical inventory |
-| `scripts/crs_reprojection_plan.py` | CRS units, axes, candidate transform and antimeridian plan |
-| `scripts/geometry_validity_report.py` | Dry-run validity audit; optional repair to a new GeoPackage |
-| `scripts/spatial_join_audit.py` | Predicate semantics, duplicate IDs and join cardinality |
-| `scripts/export_plan.py` | Non-executing vector/GeoParquet export contract |
-| `scripts/sensitive_coordinates_checklist.py` | Privacy/generalization release gate |
-
-```bash
-python skills/geopandas/scripts/vector_inventory.py --help
-python skills/geopandas/scripts/crs_reprojection_plan.py \
-  --source-crs EPSG:4326 --target-crs EPSG:32631
-python skills/geopandas/scripts/geometry_validity_report.py data.gpkg
-python skills/geopandas/scripts/spatial_join_audit.py points.gpkg zones.gpkg \
-  --predicate within --left-id point_id --right-id zone_id
-python skills/geopandas/scripts/export_plan.py data.gpkg result.parquet \
-  --format geoparquet --schema-version 1.0.0 \
-  --stable-id-column feature_id --id-unique-verified
-python skills/geopandas/scripts/sensitive_coordinates_checklist.py \
-  --public-output --precise-points --contains-addresses
-```
-
-## Reference index
-
-- [Data structures](references/data-structures.md)
-- [CRS management](references/crs-management.md)
-- [Geometric operations](references/geometric-operations.md)
-- [Spatial analysis](references/spatial-analysis.md)
-- [Data I/O](references/data-io.md)
-- [Visualization](references/visualization.md)
-
-## Sources (verified 2026-07-23)
-
-- [GeoPandas 1.1.4 on PyPI](https://pypi.org/project/geopandas/1.1.4/) — released 2026-06-26.
-- [GeoPandas 1.1.4 release](https://github.com/geopandas/geopandas/releases/tag/v1.1.4) — bug-fix release.
-- [GeoPandas 1.1.4 tagged dependencies](https://github.com/geopandas/geopandas/blob/v1.1.4/pyproject.toml).
-- [Stable GeoPandas documentation](https://geopandas.org/en/stable/).
-- [GeoPandas 1.0 migration release](https://github.com/geopandas/geopandas/releases/tag/v1.0.0).