### Year-round habitat suitability of Antarctic blue whales (Balaenoptera musculus intermedia) in the Southern Ocean
### Year-round habitat suitability of Antarctic minke whales (Balaenoptera bonaerensis) in the Southern Ocean
### Year-round habitat suitability of fin whales (Balaenoptera physalus) in the Southern Ocean
### Year-round habitat suitability of humpback whales (Megaptera novaeangliae) in the Southern Ocean

### Data downloaded on 2025-05-21 from: https://doi.org/10.1594/PANGAEA.946772; https://doi.org/10.1594/PANGAEA.946769; https://doi.org/10.1594/PANGAEA.946773; https://doi.org/10.1594/PANGAEA.946774
    El-Gabbas A, Van Opzeeland I, Burkhardt E, Boebel O (2022). Year-round habitat suitability of Antarctic blue whales (Balaenoptera musculus intermedia) in the Southern Ocean [dataset]. PANGAEA, https://doi.org/10.1594/PANGAEA.946772
    El-Gabbas A, Van Opzeeland I, Burkhardt E, Boebel O (2022). Year-round habitat suitability of Antarctic minke whales (Balaenoptera bonaerensis) in the Southern Ocean [dataset]. PANGAEA, https://doi.org/10.1594/PANGAEA.946769
    El-Gabbas A, Van Opzeeland I, Burkhardt E, Boebel O (2022). Year-round habitat suitability of fin whales (Balaenoptera physalus) in the Southern Ocean [dataset]. PANGAEA, https://doi.org/10.1594/PANGAEA.946773
    El-Gabbas A, Van Opzeeland I, Burkhardt E, Boebel O (2022). Year-round habitat suitability of humpback whales (Megaptera novaeangliae) in the Southern Ocean [dataset]. PANGAEA, https://doi.org/10.1594/PANGAEA.946774

### Description
    Year-round habitat suitability of four whale species in the Southern Ocean. The data extracted for this analysis is the MEAN across the full study period by month, and quantiles 10% and 90% for each month.

### Processing steps

    1.	Download Monthly_Habitat_suitability_Summary files from https://doi.org/10.1594/PANGAEA.946772
        Download Monthly_Habitat_suitability_Summary files from https://doi.org/10.1594/PANGAEA.946769
        Download Monthly_Habitat_suitability_Summary files from https://doi.org/10.1594/PANGAEA.946773
        Download Monthly_Habitat_suitability_Summary files from https://doi.org/10.1594/PANGAEA.946774
        ...file all in species_folder/input/

    2.  Run 02_process_tifs.py2.  Run 02_process_tifs.py
        * Read monthly Mean, SD, 10th percentile and 90th percentile rasters from each species_folder/input/
        * Reproject all rasters to EPSG:3031 (bilinear resampling)
        * Create a 48-band monthly habitat-suitability raster for each species, ordered as:
            months 1-12 mean
            months 1-12 SD
            months 1-12 10th percentile
            months 1-12 90th percentile
        * Calculate annual mean habitat suitability from the 12 monthly means
        * Calculate annual SD by combining within-month interannual variation with variation between monthly means
        * Save DEFLATE-compressed COGs in each species_folder/output/:
            <prefix>_Habitat_Monthly.tif
            <prefix>_Habitat_Annual.tif