DemoNotebook.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "enclosed-desire",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "* Owlready2 * Warning: optimized Cython parser module 'owlready2_optimized' is not available, defaulting to slower Python implementation\n"
     ]
    }
   ],
   "source": [
    "import cromo\n",
    "from cromo.click_encapsulate.commands import find_models\n",
    "from cromo.catalogs.model_catalog import checkConfigViability, getModelRules"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "other-bracket",
   "metadata": {},
   "outputs": [],
   "source": [
    "RULES_DIR='/Users/varun/git/cromo/src/cromo/resources/rules' \n",
    "ONTOLOGY_DIR='/Users/varun/git/cromo/src/cromo/resources/ontology'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "painful-florence",
   "metadata": {},
   "outputs": [],
   "source": [
    "scenario = 'Fire'\n",
    "region_geojson = '/Users/varun/git/cromo/src/cromo/test/georgetown.geojson'\n",
    "start_date = '2020-10-01T00:00:00'\n",
    "end_date = '2020-10-31T00:00:00'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "covered-mobile",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'http://geojson.io/#data=data:application/json,%7B%22type%22%3A%20%22FeatureCollection%22%2C%20%22features%22%3A%20%5B%7B%22id%22%3A%20%220%22%2C%20%22type%22%3A%20%22Feature%22%2C%20%22properties%22%3A%20%7B%7D%2C%20%22geometry%22%3A%20%7B%22type%22%3A%20%22Polygon%22%2C%20%22coordinates%22%3A%20%5B%5B%5B-120.840449%2C%2038.920212%5D%2C%20%5B-120.835169%2C%2038.920212%5D%2C%20%5B-120.835169%2C%2038.923656%5D%2C%20%5B-120.840449%2C%2038.923656%5D%2C%20%5B-120.840449%2C%2038.920212%5D%5D%5D%7D%7D%5D%7D'"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import geojsonio\n",
    "import geopandas as gpd\n",
    "region = gpd.read_file(region_geojson)\n",
    "geojsonio.display(region.to_json())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "civil-religion",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "QUIC-Fire Flat Terrain High Wind \n",
      "=================================\n",
      "                                                                                                    \n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 782.3375756250456                                                  \n",
      "\t- mean_slope = 21.18185666441626                                                  \n",
      "\n",
      "Input: wind.csv\n",
      "- Dataset: Gusts-Georgia-Slide-Road-CA ( 1 files... )                                               \n",
      "\t- min_wind_10 = 0.0                                                  \n",
      "\t- mean_wind_10 = 2.0172165163081193                                                  \n",
      "\t- max_wind_10 = 13.88                                                  \n",
      "\t- min_wind_dir = 1.24                                                  \n",
      "\t- mean_wind_dir = 153.06256087981146                                                  \n",
      "\t- max_wind_dir = 359.4                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "- wind.csv : Gusts-Georgia-Slide-Road-CA\n",
      "\n",
      "✗ INVALID\n",
      "\t ✓ Mean wind speed >= 1 m/s\n",
      "\t ✗ Mean slope should be < 20%\n",
      "\n",
      "QUIC-Fire Max\n",
      "=============\n",
      "                                                                                                    \n",
      "Input: wind.csv\n",
      "- Dataset: Gusts-Georgia-Slide-Road-CA ( 1 files... )                                               \n",
      "\t- min_wind_10 = 0.0                                                  \n",
      "\t- mean_wind_10 = 2.0172165163081193                                                  \n",
      "\t- max_wind_10 = 13.88                                                  \n",
      "\t- min_wind_dir = 1.24                                                  \n",
      "\t- mean_wind_dir = 153.06256087981146                                                  \n",
      "\t- max_wind_dir = 359.4                                                  \n",
      "\n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 782.3375756250456                                                  \n",
      "\t- mean_slope = 21.18185666441626                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- wind.csv : Gusts-Georgia-Slide-Road-CA\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "\n",
      "✓ VALID\n",
      "\t ✓ Valid for any slope\n",
      "\t ✓ Valid for any wind speed\n",
      "\n",
      "QUIC-Fire Flat Terrain Low Wind\n",
      "===============================\n",
      "                                                                                                    \n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 782.3375756250456                                                  \n",
      "\t- mean_slope = 21.18185666441626                                                  \n",
      "\n",
      "Input: wind.csv\n",
      "- Dataset: Gusts-Georgia-Slide-Road-CA ( 1 files... )                                               \n",
      "\t- min_wind_10 = 0.0                                                  \n",
      "\t- mean_wind_10 = 2.0172165163081193                                                  \n",
      "\t- max_wind_10 = 13.88                                                  \n",
      "\t- min_wind_dir = 1.24                                                  \n",
      "\t- mean_wind_dir = 153.06256087981146                                                  \n",
      "\t- max_wind_dir = 359.4                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "- wind.csv : Gusts-Georgia-Slide-Road-CA\n",
      "\n",
      "✗ INVALID\n",
      "\t ✗ Mean slope should be < 20%\n",
      "\t ✗ Mean wind speed should be < 1 m/s\n",
      "\n",
      "QUIC-Fire Rough Terrain Low Wind\n",
      "================================\n",
      "                                                                                                    \n",
      "Input: wind.csv\n",
      "- Dataset: Gusts-Georgia-Slide-Road-CA ( 1 files... )                                               \n",
      "\t- min_wind_10 = 0.0                                                  \n",
      "\t- mean_wind_10 = 2.0172165163081193                                                  \n",
      "\t- max_wind_10 = 13.88                                                  \n",
      "\t- min_wind_dir = 1.24                                                  \n",
      "\t- mean_wind_dir = 153.06256087981146                                                  \n",
      "\t- max_wind_dir = 359.4                                                  \n",
      "\n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 782.3375756250456                                                  \n",
      "\t- mean_slope = 21.18185666441626                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- wind.csv : Gusts-Georgia-Slide-Road-CA\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "\n",
      "✗ INVALID\n",
      "\t ✓ Mean slope >= 20%\n",
      "\t ✗ Mean wind speed should be < 1 m/s\n",
      "\n",
      "DRAFT Hernan's Configuration: Example of problems with current UI\n",
      "=================================================================\n",
      "                                                                                                    \n",
      "Input: moisture.dat\n",
      "- No datasets found in data catalog matching input variables ['fuel_moisture_content_array'] for this region and time period.\n",
      "\n",
      "QUIC-Fire Flat Terrain\n",
      "======================\n",
      "                                                                                                    \n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 782.3375756250456                                                  \n",
      "\t- mean_slope = 21.18185666441626                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "\n",
      "✗ INVALID\n",
      "\t ✗ Mean slope should be < 20%\n",
      "\n",
      "QUIC-Fire Rough Terrain\n",
      "=======================\n",
      "                                                                                                    \n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 782.3375756250456                                                  \n",
      "\t- mean_slope = 21.18185666441626                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "\n",
      "✓ VALID\n",
      "\t ✓ Mean slope >= 20%\n",
      "\n",
      "QUIC-Fire Rough Terrain High Wind\n",
      "=================================\n",
      "                                                                                                    \n",
      "Input: wind.csv\n",
      "- Dataset: Gusts-Georgia-Slide-Road-CA ( 1 files... )                                               \n",
      "\t- min_wind_10 = 0.0                                                  \n",
      "\t- mean_wind_10 = 2.0172165163081193                                                  \n",
      "\t- max_wind_10 = 13.88                                                  \n",
      "\t- min_wind_dir = 1.24                                                  \n",
      "\t- mean_wind_dir = 153.06256087981146                                                  \n",
      "\t- max_wind_dir = 359.4                                                  \n",
      "\n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 782.3375756250456                                                  \n",
      "\t- mean_slope = 21.18185666441626                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- wind.csv : Gusts-Georgia-Slide-Road-CA\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "\n",
      "✓ VALID\n",
      "\t ✓ Mean slope >= 20%\n",
      "\t ✓ Mean wind speed >= 1 m/s\n"
     ]
    }
   ],
   "source": [
    "configs = find_models(scenario)\n",
    "for config in configs:\n",
    "        if config.has_input is not None:\n",
    "            rules = getModelRules(config.id, rulesdir=RULES_DIR)\n",
    "            if len(rules) > 0:\n",
    "                print(\"\\n{}\\n{}\".format(config.label[0], \"=\"*len(config.label[0])))\n",
    "                checkConfigViability(config, region_geojson, start_date, end_date, rulesdir=RULES_DIR, ontdir=ONTOLOGY_DIR)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "processed-amendment",
   "metadata": {},
   "outputs": [],
   "source": [
    "scenario = 'Fire'\n",
    "region_geojson = '/Users/varun/git/cromo/src/cromo/test/grassvalley.geojson'\n",
    "start_date = '2020-10-01T00:00:00'\n",
    "end_date = '2020-10-31T00:00:00'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "unauthorized-preparation",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'http://geojson.io/#data=data:application/json,%7B%22type%22%3A%20%22FeatureCollection%22%2C%20%22features%22%3A%20%5B%7B%22id%22%3A%20%220%22%2C%20%22type%22%3A%20%22Feature%22%2C%20%22properties%22%3A%20%7B%7D%2C%20%22geometry%22%3A%20%7B%22type%22%3A%20%22Polygon%22%2C%20%22coordinates%22%3A%20%5B%5B%5B-121.080751%2C%2039.21031%5D%2C%20%5B-121.019037%2C%2039.21031%5D%2C%20%5B-121.019037%2C%2039.233981%5D%2C%20%5B-121.080751%2C%2039.233981%5D%2C%20%5B-121.080751%2C%2039.21031%5D%5D%5D%7D%7D%5D%7D'"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import geojsonio\n",
    "import geopandas as gpd\n",
    "region = gpd.read_file(region_geojson)\n",
    "geojsonio.display(region.to_json())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "cleared-function",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "QUIC-Fire Flat Terrain High Wind \n",
      "=================================\n",
      "                                                                                                    \n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 779.2601915551985                                                  \n",
      "\t- mean_slope = 11.060850330723047                                                  \n",
      "\n",
      "Input: wind.csv\n",
      "- Dataset: Gusts-Maryland-Road-CA ( 1 files... )                                                    \n",
      "\t- min_wind_10 = 0.0                                                  \n",
      "\t- mean_wind_10 = 0.38171041666666666                                                  \n",
      "\t- max_wind_10 = 4.182                                                  \n",
      "\t- min_wind_dir = 0.75                                                  \n",
      "\t- mean_wind_dir = 169.23192700729928                                                  \n",
      "\t- max_wind_dir = 359.9                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "- wind.csv : Gusts-Maryland-Road-CA\n",
      "\n",
      "✗ INVALID\n",
      "\t ✓ Mean slope < 20%\n",
      "\t ✗ Mean wind speed should be >= 1 m/s\n",
      "\n",
      "QUIC-Fire Max\n",
      "=============\n",
      "                                                                                                    \n",
      "Input: wind.csv\n",
      "- Dataset: Gusts-Maryland-Road-CA ( 1 files... )                                                    \n",
      "\t- min_wind_10 = 0.0                                                  \n",
      "\t- mean_wind_10 = 0.38171041666666666                                                  \n",
      "\t- max_wind_10 = 4.182                                                  \n",
      "\t- min_wind_dir = 0.75                                                  \n",
      "\t- mean_wind_dir = 169.23192700729928                                                  \n",
      "\t- max_wind_dir = 359.9                                                  \n",
      "\n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 779.2601915551985                                                  \n",
      "\t- mean_slope = 11.060850330723047                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- wind.csv : Gusts-Maryland-Road-CA\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "\n",
      "✓ VALID\n",
      "\t ✓ Valid for any slope\n",
      "\t ✓ Valid for any wind speed\n",
      "\n",
      "QUIC-Fire Flat Terrain Low Wind\n",
      "===============================\n",
      "                                                                                                    \n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 779.2601915551985                                                  \n",
      "\t- mean_slope = 11.060850330723047                                                  \n",
      "\n",
      "Input: wind.csv\n",
      "- Dataset: Gusts-Maryland-Road-CA ( 1 files... )                                                    \n",
      "\t- min_wind_10 = 0.0                                                  \n",
      "\t- mean_wind_10 = 0.38171041666666666                                                  \n",
      "\t- max_wind_10 = 4.182                                                  \n",
      "\t- min_wind_dir = 0.75                                                  \n",
      "\t- mean_wind_dir = 169.23192700729928                                                  \n",
      "\t- max_wind_dir = 359.9                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "- wind.csv : Gusts-Maryland-Road-CA\n",
      "\n",
      "✓ VALID\n",
      "\t ✓ Mean slope < 20%\n",
      "\t ✓ Mean wind speed < 1 m/s\n",
      "\n",
      "QUIC-Fire Rough Terrain Low Wind\n",
      "================================\n",
      "                                                                                                    \n",
      "Input: wind.csv\n",
      "- Dataset: Gusts-Maryland-Road-CA ( 1 files... )                                                    \n",
      "\t- min_wind_10 = 0.0                                                  \n",
      "\t- mean_wind_10 = 0.38171041666666666                                                  \n",
      "\t- max_wind_10 = 4.182                                                  \n",
      "\t- min_wind_dir = 0.75                                                  \n",
      "\t- mean_wind_dir = 169.23192700729928                                                  \n",
      "\t- max_wind_dir = 359.9                                                  \n",
      "\n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 779.2601915551985                                                  \n",
      "\t- mean_slope = 11.060850330723047                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- wind.csv : Gusts-Maryland-Road-CA\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "\n",
      "✗ INVALID\n",
      "\t ✓ Mean wind speed < 1 m/s\n",
      "\t ✗ Mean slope should be >= 20%\n",
      "\n",
      "DRAFT Hernan's Configuration: Example of problems with current UI\n",
      "=================================================================\n",
      "                                                                                                    \n",
      "Input: moisture.dat\n",
      "- No datasets found in data catalog matching input variables ['fuel_moisture_content_array'] for this region and time period.\n",
      "\n",
      "QUIC-Fire Flat Terrain\n",
      "======================\n",
      "                                                                                                    \n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 779.2601915551985                                                  \n",
      "\t- mean_slope = 11.060850330723047                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "\n",
      "✓ VALID\n",
      "\t ✓ Mean slope < 20%\n",
      "\n",
      "QUIC-Fire Rough Terrain\n",
      "=======================\n",
      "                                                                                                    \n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 779.2601915551985                                                  \n",
      "\t- mean_slope = 11.060850330723047                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "\n",
      "✗ INVALID\n",
      "\t ✗ Mean slope should be >= 20%\n",
      "\n",
      "QUIC-Fire Rough Terrain High Wind\n",
      "=================================\n",
      "                                                                                                    \n",
      "Input: wind.csv\n",
      "- Dataset: Gusts-Maryland-Road-CA ( 1 files... )                                                    \n",
      "\t- min_wind_10 = 0.0                                                  \n",
      "\t- mean_wind_10 = 0.38171041666666666                                                  \n",
      "\t- max_wind_10 = 4.182                                                  \n",
      "\t- min_wind_dir = 0.75                                                  \n",
      "\t- mean_wind_dir = 169.23192700729928                                                  \n",
      "\t- max_wind_dir = 359.9                                                  \n",
      "\n",
      "Input: topo.dat\n",
      "- Dataset: CA-DEM-Placeholder ( 1 files... )                                                        \n",
      "\t- mean_elevation = 779.2601915551985                                                  \n",
      "\t- mean_slope = 11.060850330723047                                                  \n",
      "\n",
      "Constraint Reasoning Over MOdel:\n",
      "\n",
      "------ Data combination 1 -------\n",
      "- wind.csv : Gusts-Maryland-Road-CA\n",
      "- topo.dat : CA-DEM-Placeholder\n",
      "\n",
      "✗ INVALID\n",
      "\t ✗ Mean wind speed should be >= 1 m/s\n",
      "\t ✗ Mean slope should be >= 20%\n"
     ]
    }
   ],
   "source": [
    "configs = find_models(scenario)\n",
    "for config in configs:\n",
    "        if config.has_input is not None:\n",
    "            rules = getModelRules(config.id, rulesdir=RULES_DIR)\n",
    "            if len(rules) > 0:\n",
    "                print(\"\\n{}\\n{}\".format(config.label[0], \"=\"*len(config.label[0])))\n",
    "                checkConfigViability(config, region_geojson, start_date, end_date, rulesdir=RULES_DIR, ontdir=ONTOLOGY_DIR)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "anonymous-receiver",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.9.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}