Datos de investigación

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)

Cometto, Pablo Marcelo Icon

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

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.

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
 
FS-WEL-2-10_readme_txt.docx
tutorial
8.41 KB
netcdf_0_25_bandpass_2_10y.nc
FIR filtered output file
24.44 MB
clasificacion_kmeans_BP10-2y_robust.nc
Clasification output file
117.58 KB
netcdf_sbai_025_raw.nc
Raw data input file
48.86 MB
welch_fft_2_10y.nc
Welch data output file
1.8 MB
Preparando la descarga 0%