NetCDF Territorial Classification and Python Colab Scripts Based on Welch Power Spectrum SbAI Signals Between 0.1 and 0.5 Cycles per Year in Southern South America (Argentina, Chile, and Uruguay). (FS-WEL-2-10)

Otros responsables
Fecha
Publicador:
Consejo Nacional de Investigaciones Científicas y Técnicas
Fecha de depósito:
27/05/2026
Fecha de creación:
14/5/2026
Clasificación temática
Palabras clave
Resumen
This notebook implements a full pipeline for band-pass spectral analysis and clustering of FIR-filtered SbAI data:
- Loads Welch spectra per pixel after FIR band-pass filtering (10–2 years).
- Applies RobustScaler preprocessing to mitigate outliers and normalize features.
- Performs KMeans++ clustering (n_clusters=4, n_init=10).
- Generates classification maps with geographic boundaries from a geopackage.
- Produces cluster centroids across all frequencies and representative pixels closest to centroids.
- Provides autonomous visualization scripts to replicate maps and spectral plots without recalculating clustering.
Información Técnica
The workflow consists of:
1. Mounting Google Drive and setting the working directory (`base_dir`).
2. Loading FIR band-pass dataset (`welch_fft_2_10y.nc`).
3. Stacking Welch spectra into pixel × frequency matrix.
4. Applying scaling:
- RobustScaler to reduce influence of outliers.
5. Running KMeans++ clustering (4 clusters, n_init=10):
- Computing metrics: explained variance (R²), silhouette index, relative reconstruction error.
- Saving classification results (`clasificacion_kmeans_BP10-2y_robust.nc`).
6. Plotting:
- Geographic cluster map with `canvas.gpkg`.
- Spectral plots of cluster centroids.
- Spectral plots of representative pixels.
7. Autonomous visualization:
- Scripts to reproduce maps and spectral plots directly from saved NetCDFs.
- No need to rerun clustering.
Colecciones
Condiciones de uso
Las buenas prácticas científicas esperan que se otorgue el crédito adecuado mediante una citación. Utilice un formato de citación y aplique estas normas de reutilización:
Licencia Creative Commons
Excepto donde se indique lo contrario, la licencia de este ítem se describe como by-nc-sa
Archivos
ARCHIVO
NOTAS DE USO
TAMAÑO
