Geologists seek to understand the relationship between volcanic unrest and eruptions by identifying subtle Volcanic Thermal Features (VTFs) in high-resolution satellite imagery. This analysis requires the careful curation of large databases of relevant volcanic thermal information. However, volcanic unrest is characterized by highly subtle thermal anomalies. Manual identification on a global scale is highly labor- and time-intensive. We propose Hotspotter: an end-to-end system to automatically detect subtle volcanic thermal anomalies in satellite images and derive relevant thermal statistics. Previous solutions for automated VTF detection have limited data size and geographic diversity. To accommodate an unprecedentedly large and diverse volcanic dataset, we propose an automated pipeline combining unsupervised anomaly detection with supervised classification to filter anomalous regions. Hotspotter gives 90% anomaly detection accuracy and robust generalization to new volcanoes. Our automated approach can accelerate scientists’ search for VTFs to help identify relevant thermal precursors and enable more precise forecasts of global volcanic eruptions.
-
Install DORA on your local machine following the steps from the vtf-dora branch in this forked repo -https://github.com/nasaharvest/dora/tree/vtf-dora/dora_exp_pipeline
-
Navigate to directory inference_pipeline with
cd inference_pipeline/ -
Once the results are obtained, run
pip install -r requirements.txtto install all requirements for this app -
Run
vtf_app.pywith your specified input directory and output file for one volcano at a time. For example:python3 vtf_app.py --input="/data/Pacaya/" --output="Pacaya_res.csv"You can runpython3 vtf_app.py --helpif you want to specify any additional arguments
Note: Your directory structure should be of the following format - {Volcano}/"all"/{ASTER file name}, for example: Pacaya/all/AST_08_00302282019042806_20230914075707_28909.SurfaceKineticTemperature.KineticTemperature.tif
- Run
run_all.pywith your specified input directory of a volcanic region. For example:python3 run_all.py --region="Kuril_Islands"You can runpython3 run_all.py --helpif you want help specifying arguments
The VTF app works with ASTER data.
- Volcano_Data_Collected - contains 52 volcanoes and 7941 total TIF files (Box folder)
- Merged_tables_processed.xlsx - latest pruned and processed excel files containing these volcanoes and correct volcano names
