A Python-Based Workflow for Asbestos Roof Mapping and Temporal Monitoring Using Satellite Imagery

The detection and monitoring of asbestos–cement roofing remain a critical public health and environmental challenge, especially in urban and suburban areas where asbestos-containing materials are still widespread due to their extensive use in the 20th century. Although hyperspectral and high-resolution multispectral remote sensing have proven effective for mapping asbestos–cement roofs, many existing approaches rely on proprietary software, limiting transparency, reproducibility, and large-scale adoption. This study presents a fully reproducible, cost-free Python-based workflow for the detection and temporal monitoring of asbestos–cement roofing using high-resolution multispectral WorldView-3 imagery. The workflow integrates atmospheric correction (using the Py6S radiative transfer model), spatial preprocessing, supervised pixel-based classification, postprocessing, and building-level aggregation within an open framework. A Maximum Likelihood Classifier is applied to VNIR and SWIR data using empirically defined roof typologies to enhance class separability. Pixel-level results are aggregated to the building scale through adaptive thresholding enabling the translation of spectral classifications into meaningful building-level information. Tested over the city of Mantua (Italy), the approach achieved reliable classification performance and enabled multi-temporal comparison to identify changes potentially due to roof remediation. Evaluation metrics (precision, recall, and F1-score) highlight the importance of carefully choosing the building-level threshold. By relying exclusively on open-source tools, the workflow enhances transparency, reproducibility, and scalability for long-term monitoring.