A Python Command Line Interface for classifying fire events from the Collection 6 MODIS Burned Area Product.
This package uses a space-time window to classify individual burn detections from late 2001 to near-present into discrete events and return both a data table and shapefiles of these events. The user is able to specify the spatial and temporal parameters of the window, as well as the spatial and temporal extent, using either a shapefile or a list of MODIS Sinusoidal Projection tile IDs. Shapefiles include full event polygons by default, and the user has the option of having firedpy produce daily-level perimeters, providing a representation of both final and expanding event perimeters.
Any area from the world may be selected. However, in the current version, memory constraints may limit the extent available for a single model run. Equatorial regions have much more fire activity, and may require much more RAM to process than a normal laptop will have.
More methodological information is at:
Balch, J. K., St. Denis, L. A., Mahood, A. L., Mietkiewicz, N. P., Williams, T. P., McGlinchy J, and Cook, M. C. 2020. FIRED (Fire Events Delineation): An open, flexible algorithm & database of U.S. fire events derived from the MODIS burned area product (2001-19). Remote Sensing, 12(21), 3498; https://doi.org/10.3390/rs12213498
Description of the country-level data sets is at:
Mahood, A.L. Lindrooth, E.J., Cook, M.C. and Balch, J.K. Country-level fire perimeter datasets (2001-2021). 2022. Nature Scientific Data, 9(458). https://doi.org/10.1038/s41597-022-01572-3
- 10/14/2024 FIREDpy V2.0
- No longer using setup.py. See new instructions below for running it with Docker or installing it locally.
- Improved fire grouping
- Improved CLI
- Access to MODIS burn area product Version 6.1 with support up to at least 2024
Many of the data products created in Fall 2021 may be shifted by a half pixel, and may lack a coordinate reference system.
The problem is now fixed, so this will not affect new iterations of firedpy. We created a script, R/posthoc_fixes.R that contains a function to fix either or both of these problems.
Sometimes the server (fuoco.geog.umd.edu) that houses the MCD64A1 product used by firedpy is down. If this happens, you just need to wait until it comes back up.
See the issues tab for more bugs, or to report a new bug!
The algorithm and derived data products are under active development. Please take this survey to help us improve firedpy. Or just email [email protected] and Adam will be overjoyed to talk about firedpy.
Already-created products are linked below. They are housed in the CU Scholar data repository in the Earth Lab Data collection, or here.
All of the created products have an event-level shapefile in .gpkg and .shp formats. Many countries also have the daily-level shapefile, but these were not created for most countries in Africa and Asia due to computational restrictions.
- Coterminous USA + Alaska: 2001-2021 2001-2024
- US plus Canada: 2001-2021
- Canada: 2001-2021
- Hawaii: 2001-2021
- Carribean (Barbados, Bahamas, Cayman Islands, Cuba, Dominican Republic, Haiti, Jamaica, Montserrat, Puerto Rico, Saint Kitts And Nevis, Trinidad And Tobago, British Virgin Islands, Guadeloupe, Saint Barthelemy): 2001-2021
- Mexico and Central America (Belize, Guatemala, Honduras, El Salvador, Nicaragua, Costa Rica, Panama): 2001-2021
- Bolivia: 2001-2021
- Argentina: 2001-2021
- Northern South America (Suriname, French Guiana, Guyana): 2001-2021
- Chile: 2001-2021
- Uruguay: 2001-2021
- Brazil: 2001-2021
- Peru: 2001-2021
- Colombia: 2001-2021
- Ecuador: 2001-2021
- Venezuela: 2001-2021
- Paraguay: 2001-2021
- Northern Europe (Iceland, Sweden, Norway, and Denmark): 2001-2021
- Poland: 2001-2021
- Finland: 2001-2021
- Russia: 2001-2021
- Italy: 2001-2021
- Spain & Portugal: 2001-2021
- Western Europe (France, Germany, Poland, Switzerland, Belgium, Netherlands, Luxembourg and Austria): 2001-2021
- Central to Southern Europe (Estonia, Latvia, Lithuania, Belarus, Ukraine, Czech Republic, Slovakia, Hungary, Romania, Bulgaria, Montenegro, Bosnia, Turkey, Republic Of Moldova, Serbia, Albania, Slovenia, and North Macedonia): 2001-2021
- Greece: 2001-2021
- UK and Ireland: 2001-2021
- Angola: 2001-2021
- Benin: 2001-2021
- Botswana: 2001-2021
- Burundi: 2001-2021
- Burkina Faso: 2001-2021
- Cameroon: 2001-2021
- Central North Africa (Libya, Algeria, Tunisia): 2001-2021
- Chad: 2001-2021
- Central African Republic: 2001-2021
- Democratic Republic of the Congo: 2001-2021
- Djibouti: 2001-2021
- Equatorial Guinea: 2001-2021
- Eritrea: 2001-2021
- eSwatini: 2001-2021
- Ethiopia: 2001-2021
- Gabon: 2001-2021
- The Gambia: 2001-2021
- Ghana: 2001-2021
- Guinea: 2001-2021
- Guinea-Bissau: 2001-2021
- Ivory Coast: 2001-2021
- Kenya: 2001-2021
- Lesotho: 2001-2021
- Liberia: 2001-2021
- Madagascar: 2001-2021
- Malawi: 2001-2021
- Mali: 2001-2021
- Mauritania: 2001-2021
- Morocco: 2001-2021
- Mozambique: 2001-2021
- Namibia: 2001-2021
- Niger: 2001-2021
- Nigeria: 2001-2021
- Republic of the Congo: 2001-2021
- Rwanda: 2001-2021
- Senegal: 2001-2021
- Sierra Leone: 2001-2021
- Somalia: 2001-2021
- Somaliland: 2001-2021
- South Africa: 2001-2021
- South Sudan: 2001-2021
- Sudan: 2001-2021
- Tanzania: 2001-2021
- Togo: 2001-2021
- Uganda: 2001-2021
- Zambia: 2001-2021
- Zimbabwe: 2001-2021
- China: 2001-2021
- India: 2001-2021
- Central Asia (Turkmenistan, Kazakhstan, Uzbekistan, Kyrgystan, Tajikistan, Afghanistan, and Pakistan): 2001-2021
- Middle East (Saudi Arabia, Qatar, Oman, Yemen, United Arab Emirates, Iraq, Jordan, Syria, Israel, Palestine, Lebanon, Egypt): 2001-2021
- Mongolia: 2001-2021
- Caucasus (Armenia, Azerbaijan, Georgia): 2001-2021
- Japan: 2001-2021
- South Korea: 2001-2021
- North Korea: 2001-2021
- Taiwan: 2001-2021
- Sri Lanka: 2001-2021
- Nepal: 2001-2021
- Bhutan: 2001-2021
- Bangladesh: 2001-2021
- Vietnam: 2001-2021
- Thailand: 2001-2021
- Laos: 2001-2021
- Myanmar: 2001-2021
- Tasmania: 2001-2021
- Victoria: 2001-2021
- New South Wales + Capital Territory: 2001-2021
- Queensland: 2001-2021
- South Australia: 2001-2021
- Western Australia: 2001-2021
- Northern Territory: 2001-2021
- Philippines: 2001-2021
- Papua New Guinea: 2001-2021
- East Timor: 2001-2021
- New Zealand: 2001-2021
- Malaysia: 2001-2021
- Brunei: 2001-2021
- Indonesia: 2001-2021
There are two ways to install firedpy. Method one is to run it out of a docker container, Method 2 is to install locally.
Note, the docker container has changed from earthlab/firedpy
to earthlabcu/firedpy
-
Run the docker container in a detached state (-d) and bind it to an available port on localhost (-p 127.0.0.1:0:7681)
-
docker run -d -p 127.0.0.1:0:7681 earthlabcu/firedpy:latest
-
Call
docker ps
to get the name of the docker container you just created and the port it is running on. -
Then get into the docker container by either running docker exec:
docker exec -it <silly_name> /bin/bash
-
Or access the CLI from your browser. The output from docker ps will look like this: CONTAINER ID IMAGE COMMAND CREATED STATUS PORTS NAMES 58a8a6ed926a earthlabcu/firedpy:latest "/bin/entry.sh ttyd …" 2 minutes ago Up 2 minutes 127.0.0.1:32768->7681/tcp stupefied_hypatia
In this example the container is running on the host machine at 127.0.0.1:32768. It may be different when you run it. Access this location in your browser by copy and pasting it into your browser's address bar
After creating a new fire product, it might be useful to get it out of the docker container in order to use it.
-
First, exit the docker container by typing
exit
-
Second, copy the file out. Here we will use the example of a container with the name "unruffled_clarke". The
docker cp
command uses the syntaxdocker cp <source> <destination>
. Files inside of a docker container will have a prefix of the docker container name (or container ID) followed by a colon, then with a normal path.Here is an example command using the container name:
docker cp unruffled_clarke:/home/firedpy/proj/outputs/shapefiles/fired_events_s5_t11_2020153.gpkg /home/Documents/fired_events_s5_t11_2020153.gpkg
Another example command using the container ID:
docker cp fa73c6d3e007:/home/firedpy/proj/outputs/shapefiles/fired_events_s5_t11_2020153.gpkg /home/Documents/fired_events_s5_t11_2020153.gpkg
-
Clone this repository to a local folder and change directories into it:
git clone https://github.com/earthlab/firedpy.git
cd firedpy
-
Ensure your anaconda setup has conda-forge, channel_priority set to strict, and update your conda.
conda update conda --yes
conda config --add channels conda-forge
conda config --set channel_priority strict
-
You must have all packages listed in the environment.yaml installed using 'conda install -c conda-forge <package_name>'
-
Create and activate a conda environment:
conda env create -f environment.yml
conda activate firedpy
-
Run firedpy with no options to be prompted with input questions for each option/attribute
python bin/firedpy.py
or if running from Docker container, simplyfiredpy
-
Or use the following commands in your command line to specify the options/attributes you would like:
-
In your terminal use this command to print out the available options and their descriptions:
python bin/firedpy.py --help
-
Run firedpy with the default option to download required data and write a data table of classified fire events to a temporary directory. This uses CONUS as the default area of interest with a spatial parameter of 5 pixels (~2.3 km) and 11 days:
python bin/firedpy.py --default
-
Change the spatial and temporal parameters of the model run:
python bin/firedpy.py -spatial 6 -temporal 10
-
Specify specific tiles and a local project_directory for required data and model outputs:
python bin/firedpy.py -spatial 6 -temporal 10 -aoi h11v09 h12v09 -proj_dir /home/<user>/fired_project
-
Write shapefiles as outputs in addition to the data table:
python bin/firedpy.py -spatial 6 -temporal 10 -aoi h11v09 h12v09 -proj_dir /home/<user>/fired_project --shapefile
-
Add the most common level 3 Ecoregion as an attribute to each event:
python bin/firedpy.py bin/firedpy.py -spatial 6 -temporal 10 -aoi h11v09 h12v09 -proj_dir /home/<user>/fired_project --shapefile -ecoregion_level 3
-
Add landcover information and produce the daily burn file
python bin/firedpy.py -spatial 6 -temporal 10 -aoi h11v09 h12v09 -proj_dir /home/<user>/fired_project --shapefile -ecoregion_level 3 -landcover_type 1 -daily yes
For more information about each parameter, use:
'python bin/firedpy.py --help'
parameter | value(s) | example | description |
---|---|---|---|
-spatial | integer | -spatial 5 | pixel radius for moving window, defaults to 5 |
-temporal | integer | -temporal 11 | day radius for moving window, defaults to 11 |
-aoi | character (MODIS tile) | -aoi h11v09 | which modis tiles should be used |
-aoi | character (shapefile) | -aoi /home/firedpy/individual_countries/canada.gpkg | figures out which modis tiles to download based on the polygon -- polygon must be in the same projection as MODIS MCD64 -- all the polygons in the ref folder are correctly projected and can be used as crs templates to prepare other polygons. |
-proj_dir | character | -proj_dir /home/firedpy/proj | which directory should firedpy operate within? Defaults to a folder called "proj" within the current working directory. |
-ecoregion_type | character | -ecoregion_type na | type of ecoregion, either world or na |
-ecoregion_level | integer | -ecoregion_level 3 | if ecoregion type = na, the level (1-3) of North American ecoregions |
-landcover_type | integer and character | -landcover_type 2:username:password | number (1-3) corresponding with a MODIS/Terra+Aqua Land Cover (MCD12Q1) category. You will need to also make an account at https://urs.earthdata.nasa.gov/home and include your login information within the argument. |
-shp_type | character | -shp_type gpkg | option to build a shapefile for the fired event in gpkg, ESRI shapefile (shp), both, or none |
-file | character | -file fired_colorado | specifies the base of the file name for the tables and shapefile outputs, defaults to "fired", in the format: "(-file aruguement)toYYYYDDD(either events or daily).gpkg", with YYYY being the year, and DDD being the julian day of the last month in the time series. The example would output fired_colorado_to2021031_events.gpkg. |
-daily | character (yes or no) | -daily yes | creates daily polygons, if no just the event-level perimeters will be created. Defaults to no. |
-start_yr | integer | -start_yr 2001 | gets the hdf files from the MODIS tiles starting in this year. The first year avalible is 2001 |
-end_yr | integer | -end_yr 2021 | gets the hdf files from the MODIS tiles ending in this year. The last year avalible is 2021 |
- Country boundaries are in ref/individual_countries
- Continent boundaries are in ref/continents
- United States state boundaries for the United States of America are in ref/us_states
- Australian state boundaries are in ref/australian_states
- For example
python bin/firedpy.py -aoi /home/firedpy/ref/us_states/colorado.gpkg
, and so on. Every space is a '_'. - If using the user input option, when prompted for the name of the continent, country, or state use "_" for spaces.
- Ensure that the input shapefiles are in the modis sinusiodal projection
- step 0.1. install docker (go to the docker website for OS-specific instructions.)
- step 0.2. get a dockerhub account
- step 1. login to docker hub via the command line
docker login
orsudo docker login
- step 2. get the existing docker image set up
- docker run -t -d earthlab/firedpy
- step 3. update from github
- git pull
- step 4. build the docker container
docker build -t earthlab/firedpy:latest .
- step 5. ENSURE THE SOFTWARE STILL WORKS BEFORE PUSHING
firedpy -aoi /home/firedpy/ref/individual_countries/samoa.gpkg
- step 6. push it up to dockerhub
docker push earthlab/firedpy:latest