BS_Mapper_simulations_202309.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import ee\n",
    "ee.Initialize(project='pyregence-ee')\n",
    "import datetime\n",
    "import scripts.analysis_functions as af\n",
    "import scripts.utils as utils\n",
    "#individual fire export uses geetools\n",
    "import geetools\n",
    "%load_ext autoreload\n",
    "%autoreload 2"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Run BS Mapper for a particular year"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Simulation runs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Simulated run date: 2022-09-30\n",
      "Fire filter date: 2022-09-20\n",
      "Total Fires in FeatureCollection: 24\n",
      "Fire count post simulated date filtering: 24\n",
      "export task started: fire window diagnostics bs_simrun2022_20220930_nifctest2_recentfixed_windows_20230926 to Google drive of user.\n",
      "Skipping individual fire export\n",
      "export task started: projects/pyregence-ee/assets/fires_bs_tool/dmn_testing/bs_simrun2022_20220930_nifctest2_recentfixed_20230926\n"
     ]
    }
   ],
   "source": [
    "# For running a year of fires with a SIMULATED run date (as if run on this date)\n",
    "# REQUIRES USER INPUT/SELECTION, see starred section below. \n",
    "#  Uses specialized subfunctions that are variants of original functions\n",
    "#  Developed using MTBS fires, would need modification for NIFC (most likely, unless field names matched)\n",
    "\n",
    "# Last edited 2023-09-27 by dnekorchuk\n",
    "\n",
    "#***\n",
    "# USER: \n",
    "# 1. Set data import, align year and import dataset\n",
    "# 2. Scroll down and set export flags\n",
    "# 3. Select simulation date (by index)\n",
    "# 4. Set export flags for various results\n",
    "# 5. Change out folder if needed\n",
    "# 6. Note that fire filter date is set to the same for all sims (in a year), change if needed.\n",
    "#  Developed with MTBS, will likely need modifications for NIFC\n",
    "#***\n",
    "\n",
    "#### IMPORT DATA ******************************************\n",
    "#One set of yearly fires, no need for remove_recent() or year/ac filtering, that was already done upon data fetch\n",
    "\n",
    "# ## OPTION 1a: MTBS option 2020\n",
    "# #align year with fires import below\n",
    "# yr = '2020' \n",
    "# fires = ee.FeatureCollection(\"projects/pyregence-ee/assets/conus/mtbs/mtbs_fires_2020_20230406\")\n",
    "# data_origin = 'mtbs'\n",
    "# name_field = 'Incid_Name'\n",
    "# id_field = 'irwinID'\n",
    "# size_field = 'BurnBndAc'\n",
    "# asset_folder = 'projects/pyregence-ee/assets/fires_bs_tool/dmn_testing'\n",
    "\n",
    "# ## OPTION 1b: MTBS option 2021\n",
    "# #align year with fires import below\n",
    "# yr = '2021' \n",
    "# fires = ee.FeatureCollection(\"projects/pyregence-ee/assets/conus/mtbs/mtbs_fires_2021_20230810\")\n",
    "# data_origin = 'mtbs'\n",
    "# name_field = 'Incid_Name'\n",
    "# id_field = 'irwinID'\n",
    "# size_field = 'BurnBndAc'\n",
    "# asset_folder = 'projects/pyregence-ee/assets/fires_bs_tool/dmn_testing'\n",
    "\n",
    "# ## OPTION 1c: MTBS option 2019\n",
    "# #align year with fires import below\n",
    "# yr = '2019' \n",
    "# fires = ee.FeatureCollection(\"projects/pyregence-ee/assets/conus/mtbs/mtbs_fires_2019_20230920\") \n",
    "# data_origin = 'mtbs'\n",
    "# name_field = 'Incid_Name'\n",
    "# id_field = 'irwinID'\n",
    "# size_field = 'BurnBndAc'\n",
    "# asset_folder = 'projects/pyregence-ee/assets/fires_bs_tool/dmn_testing'\n",
    "\n",
    "# ## OPTION 1d: MTBS option 2017\n",
    "# #align year with fires import below\n",
    "# yr = '2017' \n",
    "# fires = ee.FeatureCollection(\"projects/pyregence-ee/assets/conus/mtbs/mtbs_fires_2017_20230920\") \n",
    "# data_origin = 'mtbs'\n",
    "# name_field = 'Incid_Name'\n",
    "# id_field = 'irwinID'\n",
    "# size_field = 'BurnBndAc'\n",
    "# asset_folder = 'projects/pyregence-ee/assets/fires_bs_tool/dmn_testing'\n",
    "\n",
    "\n",
    "# ## OPTION 2a: Test / development using filtered MTBS, limited numbers of fires\n",
    "# yr = '2020'\n",
    "# fires_raw = ee.FeatureCollection('projects/pyregence-ee/assets/conus/mtbs/mtbs_fires_2020_20230406')\n",
    "# irwin_list = [\n",
    "#     '56C3E0AE-98EF-4433-863D-23A5DFF3D0C2', '60146AC0-3AC8-462B-AC80-60D31F1B4D96', \n",
    "#     '4CE00061-B897-49DC-8839-D08E040FF65D', 'BDE833A4-B278-4418-8526-67F2BFC7D7CD'] \n",
    "#     #BOBCAT = 56C3E0AE-98EF-4433-863D-23A5DFF3D0C2, DOME = 60146AC0-3AC8-462B-AC80-60D31F1B4D96, \n",
    "#     #CASTLE = 4CE00061-B897-49DC-8839-D08E040FF65D (note multipolygon), SILVERADO = BDE833A4-B278-4418-8526-67F2BFC7D7CD (late fire 2020-10-26, should be filtered out)\n",
    "# #irwin_list_failure = ['6F0ECD19-02D6-4987-A2A4-FEEDA94B35BD'] #09-28 fire\n",
    "# fires = fires_raw.filter(ee.Filter.inList('irwinID',irwin_list))\n",
    "# data_origin = 'mtbs4'\n",
    "# name_field = 'Incid_Name'\n",
    "# id_field = 'irwinID'\n",
    "# size_field = 'BurnBndAc'\n",
    "# #output folder\n",
    "# asset_folder = 'projects/pyregence-ee/assets/fires_bs_tool/dmn_testing'\n",
    "\n",
    "# # ## OPTION 2b: Test / development using filtered MTBS, limited numbers of fires\n",
    "# # Hunting 2019 Expanding window issue - pick 'expanding' and Sept run below\n",
    "# yr = '2019'\n",
    "# fires_raw = ee.FeatureCollection(\"projects/pyregence-ee/assets/conus/mtbs/mtbs_fires_2019_20230920\") \n",
    "# #bad fire is AR3481609418020190917?\n",
    "# fireBad = fires_raw.filter(ee.Filter.eq('Event_ID', 'AR3481609418020190917'))\n",
    "# fires = fireBad\n",
    "# data_origin = 'mtbsfireBad'\n",
    "# name_field = 'Incid_Name'\n",
    "# id_field = 'irwinID'\n",
    "# size_field = 'BurnBndAc'\n",
    "# #output folder\n",
    "# asset_folder = 'projects/pyregence-ee/assets/fires_bs_tool/dmn_testing'\n",
    "\n",
    "\n",
    "## OPTION 3: NIFC option\n",
    "yr = '2022'\n",
    "fires_raw = ee.FeatureCollection('projects/pyregence-ee/assets/conus/nifc/nifc_fires_2022_gte100ac_20230406')\n",
    "#relevant field names: Acres, Discovery, GlobalID, Name\n",
    "name_field = 'Name'\n",
    "id_field = 'GlobalID'\n",
    "size_field = 'Acres'\n",
    "# testing with limited set\n",
    "id_list = ['dd114282-a54d-414b-bdb1-4e180fc1ea7d', '221b33e5-594c-4f68-9131-9d412e057dfd']\n",
    "fires_known = fires_raw.filter(ee.Filter.inList(id_field,id_list))\n",
    "#plus some other fires\n",
    "fires_add = fires_raw.limit(22, 'Discovery', True) \n",
    "fires = fires_known.merge(fires_add)\n",
    "data_origin = 'nifctest2'\n",
    "asset_folder = 'projects/pyregence-ee/assets/fires_bs_tool/dmn_testing' #'projects/pyregence-ee/assets/conus/nifc'\n",
    "\n",
    "\n",
    "## USER SETTINGS **********************************\n",
    "\n",
    "export_flag = True #severity classes, composite/mosiac conus\n",
    "per_fire_export = False #indiv fire results for dNBR and severity. #trying fitoprincipe batch export, failing\n",
    "csv_export_flag = True #csv log/diagnostic dates and info\n",
    "#Variants of recent fire windows\n",
    "# 'fixed' - fixed 90 days post fire, but won't go past sim date\n",
    "# 'expanding' - variable days expanding from fire to current/sim date\n",
    "# 'sliding' - fixed 90 day window, sliding along to most recent 90 days to current/sim date.  \n",
    "recent_type = 'fixed' \n",
    "\n",
    "#Set simulation date (by index)\n",
    "sim_index = 0 #pick 0 - 3 for sept, dec, march, june simulation\n",
    "sim_date_list = [\n",
    "    datetime.datetime(int(yr), 9, 30), \n",
    "    datetime.datetime(int(yr), 12, 31), \n",
    "    datetime.datetime((int(yr)+1), 3, 31), \n",
    "    datetime.datetime((int(yr)+1), 6, 30)]\n",
    "\n",
    "# NOTE: Fires are filtered to the SAME date across all simulations\n",
    "fire_filter_date = datetime.datetime(int(yr), 9, 20) \n",
    "\n",
    "\n",
    "##### SET (SIMULATED) RUN DATE #####\n",
    "\n",
    "sim_date = sim_date_list[sim_index] \n",
    "sim_date_ee = ee.Date(sim_date)\n",
    "#For consistency, foramt the sim date the same as Discovery date as from pull_mtbs_fires.py\n",
    "sim_date_formatted = str(sim_date)[0:10]\n",
    "print(f'Simulated run date: {sim_date_formatted}')\n",
    "\n",
    "# NOTE: Fires are filtered to the SAME date across all simulations\n",
    "filter_date = fire_filter_date\n",
    "filter_date_ee = ee.Date(filter_date)\n",
    "filter_date_formatted = str(filter_date)[0:10]\n",
    "print(f'Fire filter date: {filter_date_formatted}')\n",
    "\n",
    "# capture the code run date (today's / current date) as well\n",
    "code_date = datetime.datetime.now()\n",
    "code_date_formatted = str(code_date)[0:10]\n",
    "\n",
    "# Map over the fires, adding the simulated run date \n",
    "#  note: must use subfunction set_windows_sim() in bs calc function\n",
    "#  Also add id field for labeling bands in exported intermediate products\n",
    "def set_sim_date(feat: ee.Feature):\n",
    "    fire_date = ee.Date(feat.getString('Discovery'))\n",
    "    #find out if fire is after filter date (negative is before filter date)\n",
    "    days_from_filter = fire_date.difference(filter_date_ee,'day')\n",
    "    #Also want to know discovery in reference to simulation date\n",
    "    #negative days mean Discovery was BEFORE simulated date\n",
    "    days_from_sim = fire_date.difference(sim_date_ee,'day')\n",
    "    #Create a name for the individual fire export. Not currently used until indiv export fixed.\n",
    "    #concatenate name pieces. Note: cannot start with number, if doing band-method. \n",
    "    fire_id_prefix = ee.String('Fire_')\n",
    "    fire_name_raw = ee.String(feat.get(name_field))\n",
    "    # cannot have () and probably other special characters in band names, also removes all whitespace\n",
    "    fire_name = fire_name_raw.replace('[^a-zA-Z0-9]', \"\", 'g')\n",
    "    sep = ee.String('_')\n",
    "    fire_event_id = ee.String(feat.get(id_field))\n",
    "    fire_id = fire_id_prefix.cat(fire_name).cat(sep).cat(fire_event_id)\n",
    "    return feat.set('sim_date',sim_date_formatted, \n",
    "                    'days_from_sim',days_from_sim,\n",
    "                    'filter_date',filter_date_formatted,\n",
    "                    'days_from_filter',days_from_filter,\n",
    "                    'fire_id',fire_id,\n",
    "                    'code_run_date',code_date_formatted,\n",
    "                    'recent_type',recent_type)\n",
    "fires_dt = fires.map(set_sim_date)\n",
    "\n",
    "# id_list = fires_dt.aggregate_array('fire_id')\n",
    "# print(id_list.getInfo())\n",
    "\n",
    "# Filter fires to only those that were discovered PRIOR to the filter date (plus a day padding for time ranges)\n",
    "fires_sim = fires_dt.filter(ee.Filter.lt('days_from_filter', -1))\n",
    "\n",
    "##### METADATA COLLECTION 1 #####\n",
    "# Metadata collection and print\n",
    "fires_count = fires.size().getInfo()\n",
    "print(f'Total Fires in FeatureCollection: {fires_count}')\n",
    "fires_sim_count = fires_sim.size().getInfo()\n",
    "print(f'Fire count post simulated date filtering: {fires_sim_count}')\n",
    "\n",
    "metadata = {\n",
    "    'simulated_run_date':sim_date.strftime('%Y%m%d'),\n",
    "    'code_run_date':code_date.strftime('%Y%m%d'),\n",
    "    'fire_count_original':fires_count,\n",
    "    'fire_filter_date':fire_filter_date.strftime('%Y%m%d'),\n",
    "    'fire_count_filtered':fires_sim_count,\n",
    "    'code_origin':'BS_Mapper_testing202308',\n",
    "    'recent_type':recent_type\n",
    "}\n",
    "#print(metadata)\n",
    "\n",
    "\n",
    "###### GET FIRE WINDOW DIAGNOSTICS ######\n",
    "\n",
    "bs_windows = ee.FeatureCollection(fires_sim.map(af.bs_get_windows))\n",
    "\n",
    "desc_windows = f'bs_simrun{yr}_{sim_date.strftime(\"%Y%m%d\")}_{data_origin}_recent{recent_type}_windows_{code_date.strftime(\"%Y%m%d\")}'\n",
    "#wanted properties. Event_ID only in MTBS, will be empty in NIFC runs\n",
    "wanted_cols = [\n",
    "    id_field, name_field, 'Event_ID', size_field,\n",
    "    'filter_date', 'days_from_filter', 'Discovery', \n",
    "    'sim_date', 'days_from_sim', \n",
    "    'mode', 'recent_type', 'pre_start', 'pre_end', 'post_start', 'post_end',\n",
    "    'code_run_date']\n",
    "task_csv = ee.batch.Export.table.toDrive(\n",
    "    collection=bs_windows.set('code_run_date', code_date),\n",
    "    description=desc_windows,\n",
    "    selectors=wanted_cols,\n",
    "    folder='Fire',\n",
    "    fileFormat='CSV')\n",
    "\n",
    "if csv_export_flag == True:\n",
    "    print(f'export task started: fire window diagnostics {desc_windows} to Google drive of user.')\n",
    "    task_csv.start()\n",
    "\n",
    "\n",
    "##### RUN BURN SEVERITY #####\n",
    "\n",
    "# Returns RdNBR & BS for each fire ee.Feature\n",
    "bs_coll_combined = ee.FeatureCollection(fires_sim).map(af.bs_calc_v2309)\n",
    "\n",
    "# Separate out RdNBR and Miller Threshold severity classes\n",
    "#bs_rdnbr = ee.ImageCollection(bs_coll_combined).select('RdNBR')\n",
    "bs_coll = ee.ImageCollection(bs_coll_combined).select('MillersThresholds')\n",
    "\n",
    "\n",
    "##### Export intermediate: individual features #####\n",
    "\n",
    "if per_fire_export == True:\n",
    "\n",
    "    #print('INDIVIDUAL FIRES NOT RUN: Need to write export.')\n",
    "\n",
    "    desc_indiv = f'bs_simrun{yr}_{sim_date.strftime(\"%Y%m%d\")}_{data_origin}_recent{recent_type}_perfire_{code_date.strftime(\"%Y%m%d\")}'\n",
    "    \n",
    "    # TODO: batch export isn't picking up the path correctly??\n",
    "    # Parameter \"name\" value \"projects/pyregence-ee/assets\" does not match the pattern \"^projects/[^/]+/assets/.*$\"\n",
    "    \n",
    "    path_coll = asset_folder+\"/\"+desc_indiv #f'{asset_folder}/{desc_indiv}'\n",
    "    print(f'individual fire image collection: {path_coll}') \n",
    "    \n",
    "    #hard code test\n",
    "    asset_path = 'projects/pyregence-ee/assets/bs_simrun2020_20200930_mtbs4_recentfixed_perfire_20230925'\n",
    "\n",
    "    # TODO: Doesn't seem necessary in python geetools? \n",
    "    # use ! to submit a Command-Line call to the `earthengine` Command-Line Tool (https://developers.google.com/earth-engine/guides/command_line)\n",
    "    #!earthengine create collection $path_coll\n",
    "\n",
    "    #use geetools to batch export to an image collection\n",
    "    geetools.batch.Export.imagecollection.toAsset(\n",
    "        collection=bs_coll_combined, \n",
    "        assetPath=asset_path, #path_coll, \n",
    "        namePattern='{fire_id}', \n",
    "        scale=30,\n",
    "        dataType=\"float\", \n",
    "        region=None, \n",
    "        datePattern=None,\n",
    "        extra=None, \n",
    "        verbose=False\n",
    "    )\n",
    "\n",
    "else:\n",
    "    print('Skipping individual fire export')\n",
    "\n",
    "\n",
    "###### METADATA COLLECTION 3 ######\n",
    "\n",
    "#fire_calc_count = ee.ImageCollection(bs_coll_combined).size() # always an image, but some are no data\n",
    "# Instead add up flag (0/1) counts in 'post_fire_calc'\n",
    "fire_calc_count = ee.ImageCollection(bs_coll_combined).aggregate_sum('post_fire_calc')\n",
    "metadata['fire_count_results'] = fire_calc_count\n",
    "\n",
    "##### COMPOSITE BURN SEVERITY & FINAL EXPORT #####\n",
    "\n",
    "# composite that ee.ImageCollection with a max() reducer, add metadata\n",
    "bs_composite = bs_coll.max().rename('SEVERITY').set('Year',int(yr)).set(metadata)\n",
    "\n",
    "# To Asset\n",
    "desc = f'bs_simrun{yr}_{sim_date.strftime(\"%Y%m%d\")}_{data_origin}_recent{recent_type}_{code_date.strftime(\"%Y%m%d\")}'\n",
    "\n",
    "utils.exportImgtoAsset(bs_composite, \n",
    "                    desc=desc,\n",
    "                    region=None,\n",
    "                    asset_folder=asset_folder, \n",
    "                    export_type='conus',\n",
    "                    export=export_flag)\n",
    "\n"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# For one fire \n",
    "### Select a fire from a pre-existing fire featurecollection or provide your own fire feature asset (must have 'Discovery' property with value of format ee.String('yyyy-mm-dd') )\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Running for a particular fire\n",
    "# OUTDATED, using functions before simulation code & various fixes in there. \n",
    "\n",
    "# # Select the fire based on irwinID (or name and edit filter below)\n",
    "# #fire_name = 'KNP COMPLEX' \n",
    "# #fire_irwin = '499D6545-2A19-458A-8C7C-554932C6B86E'\n",
    "# #fire_name = 'WINDY' #there are 2 WINDY fires, using irwinID # fire_irwin = 'ACD5C207-BB83-4419-977D-3F64778081A3'\n",
    "# fire_irwin = 'ACD5C207-BB83-4419-977D-3F64778081A3'\n",
    "\n",
    "# # Read in fire data\n",
    "# fires = ee.FeatureCollection('projects/pyregence-ee/assets/conus/mtbs/mtbs_fires_2021_20230810')\n",
    "# print(f'Total Fires in FeatureCollection: {fires.size().getInfo()}')\n",
    "# #print(fires.first().getInfo()['properties'])\n",
    "# #print(fires.aggregate_array('Incid_Name').distinct().sort().getInfo())\n",
    "\n",
    "# # Filter to desired fire\n",
    "# #fires_f = fires.filter(ee.Filter.eq('Incid_Name',fire_name))\n",
    "# fires_f = fires.filter(ee.Filter.eq('irwinID',fire_irwin))\n",
    "# print(f'Filtered fire size: {fires_f.size().getInfo()}')\n",
    "\n",
    "# fire_name = fires_f.first().get('Incid_Name').getInfo().lower() #brittle, make sure only 1 feature\n",
    "# print(fire_name)\n",
    "\n",
    "# # Run burn severity tool \n",
    "# bs_coll = ee.FeatureCollection(fires_f).map(af.bs_calc_new)\n",
    "# #bs_fire = af.bs_calc_new(fires_f.first()) #testing direct function call\n",
    "# #print(bs_fire.propertyNames().getInfo())\n",
    "# # while max() isn't really needed for 1 fire, keeping the pattern of the multi-fire runs\n",
    "# bs_composite = ee.ImageCollection(bs_coll).max().add(1).rename('SEVERITY') #shift severity values +1 so 0 can be nodata\n",
    "\n",
    "# # To Asset\n",
    "# # get current date to put into asset name\n",
    "# cur_date = datetime.datetime.now().strftime(\"%Y%m%d\")\n",
    "# desc = f'mtbs_bs_{str(fire_name).replace(\" \", \"\")}_{str(cur_date)}'\n",
    "# print(desc)\n",
    "# #fails with Image.clipToBoundsAndScale: The geometry for image clipping must be bounded.\n",
    "# utils.exportImgtoAsset(bs_composite,  \n",
    "#                     desc=desc,\n",
    "#                     region=bs_coll.geometry(),\n",
    "#                     asset_folder='projects/pyregence-ee/assets/fires_bs_tool/dmn_testing', \n",
    "#                     export_type='single_fire', \n",
    "#                     export=False) \n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## RDNBR calculation only, as a test"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# # Outdated, check before use\n",
    "\n",
    "# # Running for a particular fire - ONLY CALC RDNBR\n",
    "\n",
    "# # Select the fire based on irwinID \n",
    "#     #fire_name = 'KNP COMPLEX' \n",
    "# #fire_irwin = '499D6545-2A19-458A-8C7C-554932C6B86E'\n",
    "#     #fire_name = 'WINDY' #there are 2 WINDY fires, using irwinID # fire_irwin = 'ACD5C207-BB83-4419-977D-3F64778081A3'\n",
    "# fire_irwin = 'ACD5C207-BB83-4419-977D-3F64778081A3'\n",
    "\n",
    "# # Read in fire data\n",
    "# fires = ee.FeatureCollection('projects/pyregence-ee/assets/conus/mtbs/mtbs_fires_2021_20230810')\n",
    "# print(f'Total Fires in FeatureCollection: {fires.size().getInfo()}')\n",
    "# #print(fires.first().getInfo()['properties'])\n",
    "# #print(fires.aggregate_array('Incid_Name').distinct().sort().getInfo())\n",
    "\n",
    "# # Filter to desired fire\n",
    "# #fires_f = fires.filter(ee.Filter.eq('Incid_Name',fire_name))\n",
    "# fires_f = fires.filter(ee.Filter.eq('irwinID',fire_irwin))\n",
    "# print(f'Filtered fire size: {fires_f.size().getInfo()}')\n",
    "\n",
    "# fire_name = fires_f.first().get('Incid_Name').getInfo().lower() #brittle, make sure only 1 feature\n",
    "# print(fire_name)\n",
    "\n",
    "# # CalC RDNBR\n",
    "# rdnbr_coll = ee.FeatureCollection(fires_f).map(af.rdnbr_only_calc)\n",
    "# #bs_fire = af.bs_calc_new(fires_f.first()) #testing direct function call\n",
    "# #print(bs_fire.propertyNames().getInfo())\n",
    "# # while max() isn't really needed for 1 fire, keeping the pattern of the multi-fire runs\n",
    "# rdnbr_composite = ee.ImageCollection(rdnbr_coll).max().rename('RDNBR') \n",
    "\n",
    "# # To Asset\n",
    "# # get current date to put into asset name\n",
    "# cur_date = datetime.datetime.now().strftime(\"%Y%m%d\")\n",
    "# desc = f'mtbs_bs_rdnbr_{str(fire_name).replace(\" \", \"\")}_{str(cur_date)}'\n",
    "# print(desc)\n",
    "# #fails with Image.clipToBoundsAndScale: The geometry for image clipping must be bounded.\n",
    "# utils.exportImgtoAsset(rdnbr_composite,  \n",
    "#                     desc=desc,\n",
    "#                     region=rdnbr_coll.geometry(),\n",
    "#                     asset_folder='projects/pyregence-ee/assets/fires_bs_tool/dmn_testing', \n",
    "#                     export_type='single_fire', \n",
    "#                     export=False)\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.13"
  },
  "vscode": {
   "interpreter": {
    "hash": "5456adbb3bcda72c3f039b556ea69b5085f7c81f44c93a142b7acd907f0ad587"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}