{
"cells": [
{
"cell_type": "raw",
"metadata": {},
"source": [
"%This notebook demonstrates the use of the workpackage template, replace with your own.\n",
"\n",
"\\documentclass[english]{workpackage}[1996/06/02]\n",
"\n",
"% input the common preamble content (required by the ipnb2latex converter)\n",
"\\input{header.tex}\n",
"\n",
"\n",
"% then follows the rest of the preamble to be placed before the begin document\n",
"% this preamble content is special to the documentclass you defined above.\n",
"\\WPproject{Computational Radiometry} % project name\n",
"\\WPequipment{} % equipment name\n",
"\\WPsubject{06 Diverse pyradi utilities} % main heading \n",
"\\WPconclusions{} \n",
"\\WPclassification{} \n",
"\\WPdocauthor{CJ Willers}\n",
"\\WPcurrentpackdate{\\today}\n",
"\\WPcurrentpacknumber{} % work package number\n",
"\\WPdocnumber{} % this doc number hosts all the work packages\n",
"\\WPprevpackdate{} % work package which this one supersedes\n",
"\\WPprevpacknumber{} % work package which this one supersedes\n",
"\\WPsuperpackdate{} % work package which comes after this one\n",
"\\WPsuperpacknumber{} % work package which comes after this one\n",
"\\WPdocontractdetails{false}\n",
"\\WPcontractname{} % contract name \n",
"\\WPorderno{} % contract order number\n",
"\\WPmilestonenumber{} % contract milestone number\n",
"\\WPmilestonetitle{} % contract milestone title\n",
"\\WPcontractline{} % contract milestone line number \n",
"\\WPdocECPnumber{} % ecp\\ecr number\n",
"\\WPdistribution{}\n",
"\n",
" \n",
"\n",
"% this is entered just before the end{document}\n",
"\\newcommand{\\atendofdoc}{\n",
"\\bibliographystyle{IEEEtran}\n",
"\\bibliography{references}\n",
"}\n",
"\n",
"%and finally the document begin.\n",
"\\begin{document}\n",
"\\WPlayout\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"# 6 Diverse pyradi utilities"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"This notebook forms part of a series on [computational optical radiometry](https://github.com/NelisW/ComputationalRadiometry#computational-optical-radiometry-with-pyradi). The notebooks can be downloaded from [Github](https://github.com/NelisW/ComputationalRadiometry#computational-optical-radiometry-with-pyradi). These notebooks are constantly revised and updated, please revisit from time to time. \n",
"\n",
"\n",
"[](http://dx.doi.org/10.5281/zenodo.9910)\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"The date of this document and module versions used in this document are given at the end of the file. \n",
"Feedback is appreciated: neliswillers at gmail dot com."
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"## Overview"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"The pyradi library has a reasonably complete collection of Planck radiator models, both spectral and wide band. A comprehensive collection of physical constants, pertinent to optical radiation is also included. This notebook introduces these functions in the [`pyradi.ryplanck`](http://nelisw.github.io/pyradi-docs/_build/html/ryplanck.html#pyradi.ryplanck) library."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [],
"source": [
"from IPython.display import display\n",
"from IPython.display import Image\n",
"from IPython.display import HTML"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"## File and directory services"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"### List files in a directory"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"The [`ryfiles.listFiles`](http://nelisw.github.io/pyradi-docs/_build/html/ryfiles.html#pyradi.ryfiles.listFiles) function \n",
"returns a list of file paths to files in a file system, searching a directory structure along the specified path, looking for files that matches the glob pattern. If specified, the search will continue into sub-directories. A list of matching names is returned.\n",
"The function supports a local or network reachable filesystem, but not URLs.\n",
"\n",
"The function signature is: \n",
"\n",
"`listFiles(root, patterns='*', recurse=1, return_folders=0, useRegex=False)`\n",
"\n",
"- `root (string)` directory root from where the search must take place.\n",
"- `patterns (string)` glob/regex pattern for filename matching.\n",
"- `recurse (unt)` flag to indicate if subdirectories must also be searched (optional).\n",
"- `return_folders (int)` flag to indicate if folder names must also be returned (optional).\n",
"- `useRegex (bool)` flag to indicate if patterns areregular expression strings (optional)."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"['05a-PlottingWithPyradi-GeneralAndCartesian.ipynb', '05b-PlottingWithPyradi-Polar-and-3D.ipynb', '.ipynb_checkpoints\\\\05a-PlottingWithPyradi-GeneralAndCartesian-checkpoint.ipynb', '.ipynb_checkpoints\\\\05b-PlottingWithPyradi-Polar-and-3D-checkpoint.ipynb']\n"
]
}
],
"source": [
"import pyradi.ryfiles as ryfiles\n",
"print(ryfiles.listFiles('.', '05*.ipynb'))"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"['05a-PlottingWithPyradi-GeneralAndCartesian.bbl', '05a-PlottingWithPyradi-GeneralAndCartesian.bib', '05a-PlottingWithPyradi-GeneralAndCartesian.blg', '05a-PlottingWithPyradi-GeneralAndCartesian.ipynb', '05a-PlottingWithPyradi-GeneralAndCartesian.log', '05a-PlottingWithPyradi-GeneralAndCartesian.pdf', '05a-PlottingWithPyradi-GeneralAndCartesian.tex', '05b-PlottingWithPyradi-Polar-and-3D.bbl', '05b-PlottingWithPyradi-Polar-and-3D.bib', '05b-PlottingWithPyradi-Polar-and-3D.blg', '05b-PlottingWithPyradi-Polar-and-3D.ipynb', '05b-PlottingWithPyradi-Polar-and-3D.log', '05b-PlottingWithPyradi-Polar-and-3D.pdf', '05b-PlottingWithPyradi-Polar-and-3D.tex', '05', '.git\\\\objects\\\\05', '.git\\\\objects\\\\06\\\\05114a0c65030f3ab63c80822e07a03700ef35', '.git\\\\objects\\\\51\\\\059c45e9c9d5cfd46021370481d725b6eecccd', '.git\\\\objects\\\\7d\\\\05a03633f8985bca427403a7ff84ba35e199e6', '.git\\\\objects\\\\99\\\\05b062b719b832d586eb36cbd3ec21fe74d328', '.git\\\\objects\\\\d8\\\\052cd1780d21a06306ec3090e7657f3ea3e17e', '.git\\\\objects\\\\d9\\\\05835262f64f8ca948425e72c5f4ee07305246', '.ipynb_checkpoints\\\\05a-PlottingWithPyradi-GeneralAndCartesian-checkpoint.ipynb', '.ipynb_checkpoints\\\\05a-PlottingWithPyradi-GeneralAndCartesian-checkpoint.pdf', '.ipynb_checkpoints\\\\05b-PlottingWithPyradi-Polar-and-3D-checkpoint.ipynb', '05\\\\05-testfile.txt', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_13_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_14_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_15_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_17_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_19_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_22_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_22_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_22_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_25_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_25_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_25_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_29_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_31_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_32_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_34_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_34_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_34_3.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_34_4.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_34_5.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_35_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_36_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_3.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_4.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_5.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_40_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_40_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_42_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_43_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_44_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_46_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_47_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_48_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_50_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_51_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_52_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_53_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_55_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_56_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_57_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_60_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_62_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_64_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_66_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_68_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_70_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_72_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_73_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_74_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_75_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_76_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_77_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_78_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_79_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_81_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_83_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_84_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_85_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_87_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_90_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_91_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_91_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_93_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_96_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_97_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_10_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_11_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_12_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_13_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_14_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_15_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_16_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_17_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_17_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_18_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_22_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_23_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_24_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_25_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_26_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_27_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_28_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_29_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_32_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_33_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_35_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_36_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_38_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_39_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_39_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_40_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_41_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_46_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_47_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_49_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_50_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_54_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_55_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_58_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_59_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_59_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_61_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_62_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_64_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_65_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_72_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_73_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_74_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_75_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_75_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_83_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_84_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_85_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_86_0.png']\n"
]
}
],
"source": [
"import pyradi.ryfiles as ryfiles\n",
"print(ryfiles.listFiles('.', '05*', return_folders=1))"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[]\n"
]
}
],
"source": [
"print(ryfiles.listFiles('.', '05*.ipynb', useRegex=1))"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"['05a-PlottingWithPyradi-GeneralAndCartesian.bbl', '05a-PlottingWithPyradi-GeneralAndCartesian.bib', '05a-PlottingWithPyradi-GeneralAndCartesian.blg', '05a-PlottingWithPyradi-GeneralAndCartesian.ipynb', '05a-PlottingWithPyradi-GeneralAndCartesian.log', '05a-PlottingWithPyradi-GeneralAndCartesian.pdf', '05a-PlottingWithPyradi-GeneralAndCartesian.tex', '05b-PlottingWithPyradi-Polar-and-3D.bbl', '05b-PlottingWithPyradi-Polar-and-3D.bib', '05b-PlottingWithPyradi-Polar-and-3D.blg', '05b-PlottingWithPyradi-Polar-and-3D.ipynb', '05b-PlottingWithPyradi-Polar-and-3D.log', '05b-PlottingWithPyradi-Polar-and-3D.pdf', '05b-PlottingWithPyradi-Polar-and-3D.tex', 'tape7-05b', '.git\\\\objects\\\\08\\\\2acc4bfa2505b8d6fe0f85446c972a05a42888', '.git\\\\objects\\\\0d\\\\7f4c205a5f765006ae18faf50ffcbe4a430767', '.git\\\\objects\\\\2f\\\\f2c19f11277ca6a4895b1fbb5fcaf805abafbc', '.git\\\\objects\\\\38\\\\4b27fe04aa105acd84ebdbb16611729d4342f4', '.git\\\\objects\\\\40\\\\dcbc05a06258a6fc781a76ef07a6ee877b65b4', '.git\\\\objects\\\\6b\\\\42dbd56bd6e9ba05a54d592cfba05733d3aa20', '.git\\\\objects\\\\78\\\\656e36efc86279d39fe23f5b7b73dac4f005b9', '.git\\\\objects\\\\7b\\\\dc9c6407a05b826590e66950ca58484a07c69b', '.git\\\\objects\\\\7d\\\\05a03633f8985bca427403a7ff84ba35e199e6', '.git\\\\objects\\\\99\\\\05b062b719b832d586eb36cbd3ec21fe74d328', '.git\\\\objects\\\\9a\\\\e05ae9a9ff4ed6b6d749d024585c86433296da', '.git\\\\objects\\\\b2\\\\bd7a8d1a1ba253d4b5cdc12a05a1b76c8a573c', '.git\\\\objects\\\\c0\\\\91f1f7728fb9da43f08705ad8ebe3968cefea6', '.git\\\\objects\\\\e8\\\\8705a4a9b8a5682b325f4fe715b42065c08d43', '.ipynb_checkpoints\\\\05a-PlottingWithPyradi-GeneralAndCartesian-checkpoint.ipynb', '.ipynb_checkpoints\\\\05a-PlottingWithPyradi-GeneralAndCartesian-checkpoint.pdf', '.ipynb_checkpoints\\\\05b-PlottingWithPyradi-Polar-and-3D-checkpoint.ipynb', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_13_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_14_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_15_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_17_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_19_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_22_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_22_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_22_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_25_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_25_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_25_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_29_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_31_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_32_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_34_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_34_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_34_3.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_34_4.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_34_5.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_35_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_36_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_3.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_4.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_38_5.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_40_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_40_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_42_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_43_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_44_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_46_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_47_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_48_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_50_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_51_2.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_52_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_53_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_55_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_56_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_57_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_60_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_62_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_64_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_66_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_68_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_70_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_72_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_73_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_74_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_75_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_76_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_77_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_78_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_79_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_81_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_83_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_84_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_85_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_87_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_90_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_91_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_91_1.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_93_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_96_0.png', 'pic\\\\05a-PlottingWithPyradi-GeneralAndCartesian_97_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_10_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_11_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_12_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_13_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_14_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_15_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_16_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_17_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_17_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_18_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_22_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_23_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_24_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_25_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_26_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_27_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_28_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_29_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_32_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_33_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_35_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_36_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_38_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_39_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_39_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_40_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_41_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_46_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_47_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_49_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_50_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_54_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_55_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_58_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_59_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_59_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_61_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_62_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_64_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_65_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_72_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_73_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_74_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_75_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_75_1.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_83_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_84_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_85_0.png', 'pic\\\\05b-PlottingWithPyradi-Polar-and-3D_86_0.png']\n"
]
}
],
"source": [
"print(ryfiles.listFiles('.', '05[ab].*', useRegex=1))"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"['05a-PlottingWithPyradi-GeneralAndCartesian.bbl', '05a-PlottingWithPyradi-GeneralAndCartesian.bib', '05a-PlottingWithPyradi-GeneralAndCartesian.blg', '05a-PlottingWithPyradi-GeneralAndCartesian.ipynb', '05a-PlottingWithPyradi-GeneralAndCartesian.log', '05a-PlottingWithPyradi-GeneralAndCartesian.pdf', '05a-PlottingWithPyradi-GeneralAndCartesian.tex', '05b-PlottingWithPyradi-Polar-and-3D.bbl', '05b-PlottingWithPyradi-Polar-and-3D.bib', '05b-PlottingWithPyradi-Polar-and-3D.blg', '05b-PlottingWithPyradi-Polar-and-3D.ipynb', '05b-PlottingWithPyradi-Polar-and-3D.log', '05b-PlottingWithPyradi-Polar-and-3D.pdf', '05b-PlottingWithPyradi-Polar-and-3D.tex']\n"
]
}
],
"source": [
"print(ryfiles.listFiles('.', '^05[ab].*\\.[^h].*', useRegex=1))"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"### Working with filenames and directories"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"The information in this section does not describe any pyradi functionality, but may be useful in this context.\n",
"\n",
"To split a qualified path into the folder and filename use [`os.path.split`](https://docs.python.org/2/library/os.path.html), yielding a tuple with the path as the first element and the filename as the second element."
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"('05', '05-testfile.txt')\n",
"05-testfile.txt\n",
"05\n"
]
}
],
"source": [
"import os.path\n",
"print(type(os.path.split('05\\\\05-testfile.txt')))\n",
"print(os.path.split('05\\\\05-testfile.txt'))\n",
"print(os.path.basename('05\\\\05-testfile.txt'))\n",
"print(os.path.dirname('05\\\\05-testfile.txt'))"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"To join a folder and a filename use [`os.path.join`](https://docs.python.org/2/library/os.path.html), but there is a catch here; `os.path.join` does not expect the tuple, but simply a number of arguments. Study the documentation carefully, because [`os.path.join`](https://docs.python.org/2/library/os.path.html) may discard some arguments.\n",
"\n",
"Alternatively, use `os.sep` or `os.path.sep` to get the pathname separator ('/' for POSIX and '\\\\' for Windows) and then use this character to join the list."
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"05\\05-testfile.txt\n",
"05\\05-testfile.txt\n",
"05\\05-testfile.txt\n"
]
}
],
"source": [
"print(os.path.join(* os.path.split('05\\\\05-testfile.txt')))\n",
"print(os.sep.join(os.path.split('05\\\\05-testfile.txt')))\n",
"print(os.path.sep.join(os.path.split('05\\\\05-testfile.txt')))"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"To test for the existance of a file use `os.path.lexists` or `os.path.exists`, see the documentation to confirm the difference between the two."
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"True\n",
"False\n"
]
}
],
"source": [
"print(os.path.lexists('05\\\\05-testfile.txt'))\n",
"print(os.path.lexists('06\\\\05-testfile.txt'))"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"### Downloading compressed and uncompressed files from the internet"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"The [`ryfiles.downloadFileUrl`](http://nelisw.github.io/pyradi-docs/_build/html/ryfiles.html#pyradi.ryfiles.downloadFileUrl) function downloads a file from a URL. The URL is used to download a file, to the saveFilename specified. If no saveFilename is given, the basename of the URL is used. The name of the downloaded file is returned, or None if the download file.\n",
"\n",
"The function signature is: \n",
"\n",
" `downloadFileUrl(url, saveFilename=None)`\n",
"\n",
"- `url (string)` the url to be accessed.\n",
"- `saveFilename (string)` path to where the file must be saved (optional).\n",
"\n",
"The [`ryfiles.downloadUntar`](http://nelisw.github.io/pyradi-docs/_build/html/ryfiles.html#pyradi.ryfiles.downloadUntar) downloads and untars a compressed tar archive, and saves all files to the specified directory. The tarfilename is used to open the tar file, extracting to the destinationDir specified. If no destinationDir is given, the local directory '.' is used. Returns a list of the intarred filenames, or None of not succesful. Before downloading, a check is done to determine if the file was already downloaded and exists in the local file system.\n",
"\n",
"The file signature is: \n",
"\n",
" `downloadUntar(tgzFilename, url, destinationDir=None, tarFilename=None)`\n",
"\n",
"- `tgzFilename (string)` the name of the tar archive file.\n",
"- `url (string)` url where to look for the file (not including the filename).\n",
"- `destinationDir (string)` to where the files must be extracted (optional).\n",
"- `tarFilename (string)` downloaded tar filename (optional).\n",
"\n",
"The [`ryfiles.untarTarfile`](http://nelisw.github.io/pyradi-docs/_build/html/ryfiles.html#pyradi.ryfiles.untarTarfile) untars a tar archive, and save all files to the specified directory. The tarfilename is used to open a file, extraxting to the saveDirname specified. If no saveDirname is given, the local directory '.' is used. Returns a list of filenames saved, or None if failed.\n",
"\n",
"The file signature is: \n",
"\n",
" `untarTarfile(tarfilename, saveDirname=None)`\n",
"\n",
"- `tarfilename (string)` the name of the tar archive.\n",
"- `saveDirname (string)` to where the files must be extracted\n",
"\n",
"\n",
"The [`ryfiles.unzipGZipfile`](http://nelisw.github.io/pyradi-docs/_build/html/ryfiles.html#pyradi.ryfiles.unzipGZipfile)\n",
"unzips a file that was compressed using the gzip format. The zipfilename is used to open a file, to the saveFilename specified. If no saveFilename is given, the basename of the zipfilename is used, but with the file extension removed. Returns the filename of the saved file, or None if failed.\n",
"\n",
"The file signature is: \n",
" \n",
" `unzipGZipfile(zipfilename, saveFilename=None)`\n",
"\n",
"- `zipfilename (string)` the zipfilename to be decompressed.\n",
"- `saveFilename (string)` to where the file must be saved (optional)."
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"The following example downloads a file, uncompress and untar the file, saving all the files in the tar archive to the destination directory (the current working directory in this case)."
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"filesAvailable are ['arrayplotdemo.txt']\n"
]
}
],
"source": [
"tgzFilename = 'arrayplotdemo.tgz'\n",
"destinationDir = '.'\n",
"tarFilename = 'arrayplotdemo.tar'\n",
"url = 'https://raw.githubusercontent.com/NelisW/pyradi/master/pyradi/data/'\n",
"dlNames = ryfiles.downloadUntar(tgzFilename, url, destinationDir, tarFilename)\n",
"\n",
"if dlNames:\n",
" print('filesAvailable are {}'.format(dlNames))"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"## Pulse detection: probability of detection and false alarm rate"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"The [`ryutils.detectSignalToNoise`](http://nelisw.github.io/pyradi-docs/_build/html/ryutils.html#pyradi.ryutils.detectSignalToNoise)\n",
"function solves for the signal to noise ratio, given the threshold to noise ratio and probability of detection. Using the theory of matched filter design, calculate the signal to noise ratio, to achieve a required probability of detection. The function returns the signal to noise ratio required to achieve the probability of detection.\n",
"Reference: \"Electro-optics handbook,\" Tech. Rep. EOH-11, RCA, 1974. RCA Technical Series Publication. \n",
"\n",
"When there is a signal present, the probability of detection (signal plus noise exceeds the threshold) is given by\n",
"\n",
"$$\n",
"P_d =\n",
"\\frac{1}{2}\\left[\n",
"1+{\\rm erf}\\left(\n",
"\\frac{i_s-i_t}{\\sqrt{2}i_n}\n",
"\\right)\n",
"\\right]\n",
"$$\n",
"\n",
"where ${\\rm erf}$ is the error function:\n",
"\n",
"$$\n",
"{\\rm erf}(z)=\\frac{2}{\\sqrt{\\pi}}\\int_0^z e^{-t^2}dt.\n",
"$$\n",
"\n",
"\n",
"\n",
"The function signature is:\n",
"\n",
" `detectSignalToNoise(ThresholdToNoise, pD)` \n",
"\n",
"- `ThresholdToNoise (float)` the threshold to noise ratio (unitless).\n",
"- `pD (float)` the probability of detection (unitless).\n",
"\n",
"The [`ryutils.detectThresholdToNoise`](http://nelisw.github.io/pyradi-docs/_build/html/ryutils.html#pyradi.ryutils.detectThresholdToNoise)\n",
"function solve for threshold to noise ratio, given pulse width and FAR, for matched filter. Using the theory of matched filter design, calculate the threshold to noise ratio, to achieve a required false alarm rate. The function returns the threshold to noise ratio. Reference R. D. Hippenstiel, Detection Theory: Applications and Digital Signal Pro-cessing, CRC Press, 2002.\n",
"\n",
"The average false alarm rate is given by\n",
"$$\n",
"{FAR}=\n",
"\\frac{1}{2 t_p \\sqrt{3}}\\exp^{-i_t^2/(2i_n^2)},\n",
"$$\n",
"where\n",
"$t_p$ is the pulse width,\n",
"$i_t$ is the threshold value, and\n",
"$i_n$ is the rms noise value at the input to the threshold detector.\n",
"\n",
"\n",
"The function signature is:\n",
" `detectThresholdToNoise(pulseWidth, FAR)`\n",
" \n",
"- `pulseWidth (float)` the signal pulse width in [s].\n",
"- `FAR (float)` the false alarm rate in [alarms/s].\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"For a laser pulse with width=1e-07, a FAR=15 and Pd=0.999,\n",
"the Threshold to Noise ratio must be 4.933070468777455\n",
"and the Signal to Noise ratio must be 8.023302774945268\n",
" \n"
]
}
],
"source": [
"import pyradi.ryutils as ryutils\n",
"\n",
"pulsewidth = 100e-9\n",
"FAR = 15\n",
"probDetection = 0.999\n",
"ThresholdToNoise = ryutils.detectThresholdToNoiseTpFAR(pulsewidth,FAR)\n",
"SignalToNoise = ryutils.detectSignalToNoiseThresholdToNoisePd(ThresholdToNoise, probDetection)\n",
"print('For a laser pulse with width={0}, a FAR={1} and Pd={2},'.format(pulsewidth,FAR,probDetection))\n",
"print('the Threshold to Noise ratio must be {0}'.format(ThresholdToNoise))\n",
"print('and the Signal to Noise ratio must be {0}'.format(SignalToNoise))\n",
"print(' ')"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"## Solving the range equation"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"It is frequently necessary to determine the operational detection distance of a source and sensor combination. The problem is usually stated as follows: \"What operating detection range can be achieved with a given source intensity, atmospheric attenuation, and sensor sensitivity?\" \n",
"The inband irradiance is given by (Sec. 7.5)\n",
"\n",
"$$\n",
"E_{\\cal S}=\\frac{v_{\\cal S} }{k \\,\\widehat{{\\cal R}}\\,Z_t\\,A_{1}}\n",
"=\n",
"\\int_{A_0} \\left(\n",
"\\frac{1}{R_{01}^2}\n",
"\\int_{0}^{\\infty}\n",
" \\epsilon_{0\\lambda} L_{0\\lambda}\n",
"\\,\\tau_{a\\lambda}(R) {\\cal S}_\\lambda\\,\n",
"d\\lambda\\right) dA_{0}\\,\\cos\\theta_0\n",
"$$\n",
"\n",
"The objective is to solve for $R$ where the irradiance $E_{{\\cal S}}=E_{{\\cal S}\\theta}$ is at the threshold at which the range must be determined. Normally, $E_{{\\cal S}\\theta}=SNR\\times NEE$ where the SNR is selected to yield a given probability of detection. \n",
"\n",
"Assuming a target smaller than the sensor field of view, and using the effective transmittance the above equation can be simplified to the form\n",
"\n",
"$$\n",
"E_{{\\cal S}\\theta}=\\frac{I_{{\\cal S}} \\tau_{\\rm eff}(R)}{R^2}\n",
"$$\n",
"\n",
"where $I_{{\\cal S}}$ is the source intensity. \n",
"This equation can be solved is $\\tau_{\\rm eff}(R)$ is available in lookup table form.\n",
"\n",
"The [`ryutils.rangeEquation`](http://nelisw.github.io/pyradi-docs/_build/html/ryutils.html#pyradi.ryutils.rangeEquation) function solves the range equation for arbitrary transmittance vs range, with the equation given by\n",
"\n",
"$$\n",
"E=\\frac{I_{{\\cal S}} \\tau_{\\rm eff}(R)}{R^n}\n",
"$$\n",
"\n",
"where $E$ is the threshold irradiance in [W/m2], and $I$ is the intensity in [W/sr]. This range equation holds for\n",
"the case where the target is smaller than the field of view.\n",
"\n",
"The range $R$ must be in [m], and $\\tau_a(R)$\n",
"is calculated from a lookup table of atmospheric transmittance vs. range.\n",
"The transmittance lookup table can be calculated from the simple Bouguer law,\n",
"or it can have any abritrary shape, provided it decreases with increasing range.\n",
"The user supplies the lookup table in the form of an array of range values and\n",
"an associated array of transmittance values. The range values need not be on\n",
"constant linear range increment.\n",
"\n",
"The parameter $n$\n",
"\n",
"* $n$=2 (default value) the general case of a radiating source smaller than the field of view.\n",
"\n",
"* $n$=4 the special case of a laser rangefinder illuminating a splot smaller than the sensor field of view, viewed against the sky. In this case there is an $R^2$ attenuation from the laser to the source and another $R^2$ attenuation from the source to the receiver, hence $R^4$ overall.\n",
"\n",
"The function signature is: \n",
" `rangeEquation(Intensity, Irradiance, rangeTab, tauTab, rangeGuess = 1, n = 2)`\n",
" \n",
"- `Intensity (float or np.array[N,] or [N,1])` in [W/sr].\n",
"- `Irradiance (float or np.array[N,] or [N,1])` in [W/m2].\n",
"- `rangeTab (np.array[N,] or [N,1])` range vector for tauTab lookup in [m].\n",
"- `tauTab (np.array[N,] or [N,1])` transmittance vector for lookup range in [m].\n",
"- `rangeGuess (float)` starting value range estimate in [m] (optional).\n",
"- `n (float)` range power (2 or 4) (optional).\n",
"\n",
"If the range solution is doubtful (e.g. not a trustworthy solution) the returned value is made negative. The following example attempts to solve the equation for three cases, only one of which provides a stable solution."
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Range equation solver with irradiance threshold of 1e-99 W/m2:\n",
" intensity of 200 W/sr at range -11381.128909552699 m the irradiance is 4.462603202968597e-07 W/m2, error is 4.46e+92. \n",
"Check maximum range in lookup table\n",
"\n",
"Range equation solver with irradiance threshold of 1e-05 W/m2:\n",
" intensity of 200 W/sr at range 3452.03409008682 m the irradiance is 1.000000000000012e-05 W/m2, error is 1.19e-14. \n",
"\n",
"\n",
"Range equation solver with irradiance threshold of 1.0 W/m2:\n",
" intensity of 200 W/sr at range -14.12716136943599 m the irradiance is 0.49750424774943763 W/m2, error is -5.02e-01. \n",
"Check range resolution in lookup table\n",
"\n"
]
}
],
"source": [
"import numpy as np\n",
"from scipy.interpolate import interp1d\n",
"\n",
"rangeTab = np.linspace(0, 10000, 1000)\n",
"tauTab = np.exp(- 0.00015 * rangeTab)\n",
"Intensity=200\n",
"Irradiancetab=[10e-100, 10e-6, 10e-1]\n",
"for Irradiance in Irradiancetab:\n",
" r = ryutils.rangeEquation(Intensity = Intensity, Irradiance = Irradiance, rangeTab = rangeTab,\n",
" tauTab = tauTab, rangeGuess = 1, n = 2)\n",
"\n",
" #test the solution by calculating the irradiance at this range.\n",
" tauTable = interp1d(rangeTab, tauTab, kind = 'linear')\n",
"\n",
" if np.abs(r[0]) < rangeTab[2]:\n",
" rr = rangeTab[2]\n",
" strError = \"Check range resolution in lookup table\"\n",
" elif np.abs(r[0]) > rangeTab[-1]:\n",
" rr = rangeTab[-1]\n",
" strError = \"Check maximum range in lookup table\"\n",
" else:\n",
" rr = r[0]\n",
" strError = \"\"\n",
"\n",
" irrad = Intensity * tauTable(rr) / rr ** 2\n",
"\n",
" print('Range equation solver with irradiance threshold of {} W/m2:'.format(Irradiance))\n",
" print(' intensity of {4} W/sr at range {0} m the irradiance is {1} W/m2, error is {2:.2e}. \\n{3}\\n'.format(\n",
" r[0],irrad, (irrad - Irradiance) / Irradiance, strError, Intensity))"
]
},
{
"cell_type": "markdown",
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"source": [
"The next example plots the detection range for a missile given the threshold irradiance, viewing an aircraft with signature of 200W/sr, through an atmosphere with an attenuation coefficient of 0.15 km-1. The graphs shows that when viewing the target through the atmosphere, an irradiance threshold (sensitivity) of 1 $\\mu$W/m$^2$ will provide a detection distance of 8 km. A missile with sensitivity of 0.1 $\\mu$W/m$^2$ provides a detection distance of 15 km. The example shows that you need ten times better sensitivity to double the detection distance. Increased detection distance comes at a tremendous cost!"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {
"run_control": {
"frozen": false,
"read_only": false
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAwsAAAFUCAYAAABr4aCxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABPcUlEQVR4nO3dd5wV1f3/8dcbUGyoiNgABUuwYIuKWIMN0ajYxVgw+vtiTIzfxBY1MZaIUb8aE7+WfFHsXWxorFHXkqjYYyWiqBBQxI4FFT+/P86sDte7u3f37u5seT8fj3lw98ycmc89O3u5n5lz5igiMDMzMzMzK9Wl6ADMzMzMzKxtcrJgZmZmZmZlOVkwMzMzM7OynCyYmZmZmVlZThbMzMzMzKwsJwtmZmZmZlaWkwUzK5Sk2ZJWbIH9DpT0jKRPJB3W3Pu3tkPSUEnTio4DQNIBkh4pOo6GtJc42ypJJ0q6sug4zFqDkwWzZiTpDUmfZ1+A35Z0qaRFio6rrZBUI+n/5csiYpGIeL0FDnc0UBMRPSLinGp3JmmUpKckfSxpmqQzJHXLrV9C0s2SPpX0pqSflNTfStIrkj6T9ICkFeo4zt6SXiopu7eOsmNyPy/XHF+Ys3P2lGr305LHlxSSVm6tmFpSQ+dNybaDJN0taZak702SlP19fZF9/syWNKllo68zzvkljc8+D0PS0Cr21WYSQbPOysmCWfPbMSIWAdYB1gWOLTacTmsF4MWmVMwnATkLAb8ClgQ2BLYCjsytPw/4Elga2Ae4QNIa2f6WBG4CjgeWAJ4Erqvj8A8Cq0nqnYtlbWChkrKNgIdy9bYH7mqG91mVlthnc5PUtegYcuo8b8r4CrgeOKie/R2aJeCLRMTA5g21UR4B9gXeLjCGdqE9/M1Y5+ZkwayFRMTbwN2kpAEAScdIei3rGvOSpF1y6w6Q9IikMyV9IGmKpO1y6wdIeiir+3dJ5+Vvg0saIumfkj6U9Fx9V/Oyq9A3Sno3O85huXULZld3P8hiPCp/Za/0qm7+SrCknpJuz/b7Qfa6b7ZuDLAZcG521fPc0v1JWkzS5Vn9NyX9TlKXStqn5P3dD2yRO9YPKtj3PySdLel94MQyv88LIuLhiPgyIv4DXAVsktVfGNgNOD4iZkfEI8AEYL+s+q7AixFxQ0R8ke1/bUmrljnOdOB1YPOs6IekpOfBkrIupKSj1vbAHVk8v5H0n+xcmSRpq6z8xOyK75WSPgYOKGm30aQvrEdn7XZbVt7QeTtP20nqJek2pbswT0g6RbkuL5JWVboz8n4W3571Hb8kxtoE6blsm71y646QNFPSDEk/zZVfKukCSXdI+hTYQtJqSlfiP5T0oqSdctvPcwdMJV12JA3L4v5I0vmSHlTJHbMKz9OGzpt5RMSkiBhHE5Pg+kj6n+zva7GS3+mHkl6XtHFWPjVr41F17Sv7G/lz9n7mVnDsn0p6OTu/Xpd0cFa+MHAnsJy+u1uyXHYe35Cdx59Iej77Gz82i22qpGG5/S8naUJ2vk2W9F+5dYMlPZmdq+9I+lNW3l/ps2m0pOnZOXVESejzK32mfJKdQ+uXHLOuz9jv/R1m7T4uO85/sr+ZtpTUWifmZMGshSh9Sd4OmJwrfo30hXkx4CTgSknL5tZvCEwiXb0+AxgnSdm6q4GJQC/Sl81vv1BI6gP8DTiFdOX6SOBGZVeiS+LqAtwGPAf0IV0h/5WkbbNNTgBWypZtgTq/FJTRBbiEdFV/eeBz4FyAiPgt8DDfXfk8tEz9/yW1zYrAj4D9gZ/m1tfXPt+KiC1LjvXvCvf9OrAUMKaC97o5331p+wEwNztOreeA2ivEa2Q/18b3KelcqOsK8kN8lxhsnr2XR0rKHouILwEkzZeV3StpIHAosEFE9CD9Dt/I7XsEMB5YnJTwfCsixmZlZ2TttmO2qpLzNt925wGfAsuQzp9vz6HsC+C9pPN5KWBv4HxJa9Rz/HyMtW2wdrZN7R2aZbL4+pCuvJ8nqWeu6k+y2HoAj5P+Bu7JYvglcFXWdvVSuks0nnTHsBfpfNy4ZLOKzlMaPm+a4o9K3ZT+oQq6/0jqIulCYC1gWER8lHsP/yK9x6uBa4ENgJVJdwzOVfN1sZwJ7AAsSvqbPFvSD7O/k+2A6bm7JdOzOjsCVwA9gWdIF2a6kH7/JwP/l9v/NcA0YDlgd+BUZQk08BfgLxGxKOkz7/qS2LYAVgGGAcdI2jq3bidSuyxOSvJqL4A09BkL3/87vAz4mtS+62bHmycBNSuKkwWz5neLpE+AqaT/BE+oXZFdWZ4eEd9kX3JeBQbn6r4ZERdGxFzSfx7LAktLWp70H/Xvs6t2tVcga+0L3BERd2T7vpd01Xn7MvFtAPSOiJOzfb0OXAiMzNbvCYyJiPcjYipQcX//iHgvIm6MiM8i4hPSl7MfVVI3u4q2F3BsRHwSEW8AZzHvVday7dNM+54eEf8bEV9HxOcN7O+nwPrAmVnRIsBHJZt9RPpiWsn6Uvm7CJuRkoWHS8oezG2/OfBc1uZzge7A6pLmi4g3IuK13LaPRsQt2XlS7/usVcF5+23bkbrU7AackJ0HL5F+V7V2AN6IiEuytn4auJH0Ja4aXwEnR8RXEXEHMBvIf/m/NSL+ERHfkO72LQKclv0N3A/cTkpcGrI96S7RTdn7PYfvd7Wp9Dxt7HnRkN+QkuE+wFjgNkkr1bP9fKQv0kuQuk9+lls3JfsdzSV1metHat85EXEP6ffcLONGIuJvEfFaJA+SkrjNGqj2cETcnf0ObgB6k36fX5G+wPeXtLikfsCmwG8i4ouIeBa4iO/+9r8CVpa0ZHZ357GS45wUEZ9GxPOkCyH5c+SR7DN3LilxWTsrb+gzFnJ/h6QkaTvgV9mxZgJnl2xvVhgnC2bNb+fsiu5QYFXS1UUAJO0v6dns1v6HwKD8enJfOnL/cS9CuiL2fsl/5lNzr1cA9qjdb7bvTUlfUkqtQLqtn9/2OL77MrNcyb7frOhdp/e3kKT/U+rm8zHpCvniFd5OXxKYv+R4b5K++NSqq32aY99TqYCknYHTgO0iYlZWPJv0H37eosAnFa4v9RCwVnZlfAjpi8UrwLJZ2aZ8f7zCHQARMZk0tuJEYKakayUtl9u2oveZV8F5m99nb6BbSVnpubphyfm3D+nOQDXey7441vqMec+NfAzLAVOzL2q1Ss+Huszz9xERQbpqnVfpedrY86JeEfF4lgzPiYjLgH9Q/oJBrZVJV7hPqr1LlfNO7vXn2f5LyxaRtHyui9DspsQtaTtJj2XdhD7MYl6ygWqlsczKvrR/Gy/zfnbm2zT/uz6IdIfnFaUuczuUHKf0szD/t5RPEj8DFlAaf9DQZ2zpflcgJW4zctv/H+mul1nhnCyYtZDsCtmlZFeflZ5+cyGpi0iviFgceAEo1z2h1AxgCUkL5cr65V5PBa6IiMVzy8IRcVqZfU0lXTXMb9sjImq/VMwo2ffyJfU/Iw32rZX/kncE6Wruhtlt/dor4bXv8XtPcMmZRbrKl39K0PLAf+qpU6lK9l1fbABIGk76He6YXWms9W+gm6RVcmVr8103pRf57qpjbVeclaij73l2JXI6MBp4KyJqv4Q9mpUtAuSvgG5P6oZWW//qiNiU9H4DOL0R73Oe9RWet/k675K6U/TNlZWeqw+WnH+LRMQhFcbXVPn9Tgf6Zd1FauXPh0+p+xyfQe69Zd2L8u+1MRo6b6oV1P/58jKp28+dlXTBKnuAiLfiuy5Cje6WJKk76c7SmcDS2fl1B5V9ZlRiOumzM3+35tvfdUS8GhF7k76Ynw6Mz/4+a5V+Fk6nYQ19xsK872sqMAdYMrf9ohFRTXc0s2bjZMGsZf0Z2EbSOsDCpP8g3oVvu7IMqmQnEfEmqVvRiUqPJdyI1Ge31pXAjpK2ldRV0gJKjxws9yVmIvCx0iDYBbPtB0naIFt/PXCs0mDlvqT+3HnPAj/J6g1n3m5GPUhX9T6UtAS5LliZd0jdJMq9x7nZscdI6pF9ST08e29VaY59S9qS1Ld4t4iYWLL/T0lPOzpZ0sKSNiFdsb0i2+RmYJCk3SQtAPwe+Fd2t6AuD2cxPpwreyQre7K2C5GkAUD32n0pzS+xZfYl7AvS76PBQaY5pb+jRp23WVvfRDpXF1IaxL1/bpPbgR9I2k/SfNmygaTV6jh+JTE21uOkhODo7PhDSX9P12brnwV2zeJfmXmfPvQ3YE1JO2dXkX9BE++KVHDezEPJAqS7ZGR/592z14tnf/8LSOomaR9Ssn53AzFcQ7rq/fcGuiw1iqTuWayQBgIvkCVWpeYndZt7F/haaTD4sNz6d4BekhZrShyRulL+kzSWYwFJa5F+n1dlce4rqXd2l+nDrFr+7+X47DxYg5RY1fUUs7yGPmNLY5xB6np1lqRFlcaRrCSpoi6cZi3NyYJZC4qId4HLSU87eYnUT/5R0n+Aa5K6CVRqH9LjMt8jDWS+jnQ1qvY/xBGk//TfJV2pOooyf+PZl7kdSf22p5Cuul9EGhwKaQDrm9m6e/j+F5f/zup/mMV0S27dn4EFs30+xvcf5fkXYHelp8SUGwvxS9KXuNdJX4yvBi4us11TVLvv40ltdEeu28WdufU/J733maR+4IdExIvw7XmwG2kMxwekwaMN9Ud+kHS1Mz9x1sNZWb4L0o/JuiBlupO6Sc0idZNYinReVGocabzDh5JuaeJ5eyiprd4mnT/X8N25+gnpy+BI0lXat0lXdLuXO34d+z8RuCzbZs9GvDeyGL4kDU7djtRO5wP755K3s0l98t8hjTm4Kld3FrAHaeDye8DqpER+TmPjyNR53uS6+NTe3VuBlPzV3nn4nDSQGlI3llNIf/+zSOf7zhHR4FwLWZelk4H7JfVv4vsoNSmLrw8pYfmcee/s1R77E+AwUjL/AWkg+oTc+ldI7fJ69vternQfFdgb6E86324mjae5N1s3HHhRqQvVX4CRkZ5YVutB0kMq7gPOzMZr1KuCz9hy9iclTi+R2mE85buRmrU6pe6WZtbeSLoOeCUiSq/eN/dxhgJXRkRTu1pYC5J0B3BupEG9bZKk04FlIqIxT9ZqF7KuTNOAfSLigaLjseaRJU1TgPlKxsKYdTq+s2DWTmRdNVbKblEPJ91JuKXgsKx4NUCb+pKqNI/CWlm3mcGkbh83Fx1Xc8m6+yyedQE6jtS/vvQpOmZmHYJnDTRrP5Yh9W/uRbqSeUhEPFNsSFa0iDij6BjK6EHqOrIcqXvNWcCthUbUvDYidWOr7Tayc1T4GFozs/bG3ZDMzMzMzKysVumGJOlipSnYX8iVraP0XOVnlaZaH5xbd6zSlOyTlJvxUNJ6StO6T5Z0Tu2TFbKnLlyXlT9e1wCtauubmZmZmXUmrTVm4VLSEwfyziBNBLMO6TGCZwBIWp30lIw1sjrn67sJnS4gPWN8lWyp3edBwAcRsTLpKRb5Z4rnVVvfzMzMzKzTaJUxCxHxUJmr9cF3M1cuxncTnYwAro2IOcAUSZOBwZLeABaNiEcBJF0O7AzcmdU5Mas/HjhXkiLXx0rSstXUL2fxxRePlVdultnurQk+/fRTFl544YY3tGbnti+W279Ybv/iuO2L5fYv1lNPPTUrInq39nGLHOD8K+BuSWeS7nBsnJX3Yd6nSkzLyr7KXpeW19aZChARX0v6iDQIdFZu+z5V1gdA0mjS3Ql69+7NmWeeWfEbtuY1e/ZsFlmk0ROGWjNw2xfL7V8st39x3PbFcvsXa4sttniziOMWmSwcAvw6Im7MJtUZB2xN+anp65qyvvbKf33rqGCbSuqnwoixwFiAgQMHxtChQ8ttZq2gpqYGt38x3PbFcvsXy+1fHLd9sdz+nVOR8yyMIj0GEuAGoHaA8zSgX267vqQuStOy16Xl89SR1I3Uren9kuNVW9/MzMzMrFMpMlmYDvwoe70l8Gr2egIwMntC0QDSQOSJETED+ETSkOwpRvvz3XO7J5CSD4DdgftLxxtUW9/MzMzMrLNplW5Ikq4BhgJLSpoGnAD8F/CX7Er+F2TjACLiRUnXkya6+Rr4RUTMzXZ1COnJSguSBibfmZWPA67IBkO/T3qaUu2xn82euNSk+mZmZmZmnVVrPQ1p7zpWrVfH9mOAMWXKnwQGlSn/Atijjn2tU019MzMzM7POqshuSGZmZmZm1oY5WTAzMzMzs7KcLJiZmZmZWVlOFszMzMzMrCwnC2bt2Ndfw1lnwezZRUdiZmZmHZGTBbN27Omn4eij4cADwTODmJmZWXNzsmDWjg0eDKeeCjfcAH/6U9HRmJmZWUfjZMGsnTv6aNhtt/Tv/fcXHY2ZmZl1JE4WzNo5CS65BH7wA9hrL5g6teiIzMzMrKNwsmDWAfToATffDHPmpLsMX3xRdERmZmbWEThZMOsgVl0VLrsMnngCDjus6GjMzMysI3CyYNaB7LILHHssXHhhWszMzMyq4WTBrIP5wx9gm23g0ENh4sSiozEzM7P2zMmCWQfTtStccw0su2wavzBzZtERmZmZWXvlZMGsA+rVC266CWbNgpEj00zPZmZmZo3lZMGsg/rhD+Gvf4UHHkjjGMzMzMway8mCWQc2ahT8/Odw5plplmczMzOzxnCyYNbBnX02bLQR/PSn8OKLRUdjZmZm7YmTBbMObv75012FRRZJj1b96KOiIzIzM7P2wsmCWSfQp09KGKZMgf33h2++KToiMzMzaw+cLJh1EpttBmedBRMmwB//WHQ0ZmZm1h60SrIg6WJJMyW9UFL+S0mTJL0o6Yxc+bGSJmfrts2Vryfp+WzdOZKUlXeXdF1W/rik/nXEUVV9s/bul7+En/wEjj8e7rqr6GjMzMysrWutOwuXAsPzBZK2AEYAa0XEGsCZWfnqwEhgjazO+ZK6ZtUuAEYDq2RL7T4PAj6IiJWBs4HT64ij2vpm7ZoEY8fCmmumpOH114uOyMzMzNqyVkkWIuIh4P2S4kOA0yJiTrZN7TyzI4BrI2JOREwBJgODJS0LLBoRj0ZEAJcDO+fqXJa9Hg9sVXvXoFa19c06ioUXThO2RaQZnj/7rOiIzMzMrK0qcszCD4DNsm4/D0raICvvA0zNbTctK+uTvS4tn6dORHwNfAT0KjletfXNOoyVVoKrr4bnnoOf/SwlDmZmZmaluhV87J7AEGAD4HpJKwLlruhHPeU0sI4KtqmkftpQGk3qykTv3r2pqakpt5m1gtmzZ7v9q7DggjBq1ApceukAevZ8lV12+U/Fdd32xXL7F8vtXxy3fbHc/p1TkcnCNOCmrEvQREnfAEtm5f1y2/UFpmflfcuUk6szTVI3YDG+3+2p2voARMRYYCzAwIEDY+jQoRW+XWtuNTU1uP2rs/nmMGsWnH/+Kuyxxypsumll9dz2xXL7F8vtXxy3fbHc/p1Tkd2QbgG2BJD0A2B+YBYwARiZPaFoAGkg8sSImAF8ImlINp5gf+DWbF8TgFHZ692B+7Mk5FvV1jfriLp0gSuugP79YY89YMaMoiMyMzOztqS1Hp16DfAoMFDSNEkHARcDK2aPU70WGBXJi8D1wEvAXcAvImJutqtDgItIg55fA+7MyscBvSRNBg4Hjskd+9lcKI2ub9bRLb443HwzfPxxShi+/LLoiMzMzKytaJVuSBGxdx2r9q1j+zHAmDLlTwKDypR/AexRx77Wqaa+WWcwaBCMGwd77w1HHAH/+79FR2RmZmZtQZFjFsysDRk5Ep54Av70Jxg8GPbbr+iIzMzMrGhFjlkwszbm9NNh6FAYPRqeeaboaMzMzKxoThbM7FvdusF110GvXrDrrvB+2WeCmZmZWWfhZMHM5rHUUnDjjTB9OvzkJzB3bsN1zMzMrGNysmBm37PhhmmQ8913w4knFh2NmZmZFcXJgpmV9V//BQceCKecArfe2vD2ZmZm1vE4WTCzsiQ47zxYf/30ZKRJk4qOyMzMzFqbkwUzq9MCC6TxC927pwHPs2cXHZGZmZm1JicLZlav5ZeHa6+FV15J3ZIiio7IzMzMWouTBTNr0FZbwR//CDfcAGedVXQ0ZmZm1lqcLJhZRY46CnbbDX7zG3j66cWLDsfMzMxagZMFM6uIBJdcAgMHwsknr87UqUVHZGZmZi3NyYKZVaxHD7j5Zvjqqy7stht88UXREZmZmVlLcrJgZo0ycCAce+wrPPEE/PKXRUdjZmZmLcnJgpk12qabzuK44+Cii+DCC4uOxszMzFqKkwUza5KTT4Zhw+DQQ2HixKKjMTMzs5bgZMHMmqRrV7j6alhuufSUpJkzi47IzMzMmpuTBTNrsl690gzPs2bBXnvB118XHZGZmZk1JycLZlaVH/4Q/vpXqKmBY48tOhozMzNrTt2KDsDM2r9Ro+CJJ+DMM2GDDWDPPYuOyMzMzJqD7yyYWbP4059g443hwAPhhReKjsbMzMyag5MFM2sW888PN9wAiywCu+4KH31UdERmZmZWrVZJFiRdLGmmpO9db5R0pKSQtGSu7FhJkyVNkrRtrnw9Sc9n686RpKy8u6TrsvLHJfWvI46q6ptZ/ZZbLiUMU6bA/vvDN98UHZGZmZlVo7XuLFwKDC8tlNQP2AZ4K1e2OjASWCOrc76krtnqC4DRwCrZUrvPg4APImJl4Gzg9DriqLa+mTVgs81Sl6QJE+DUU4uOxszMzKrRKslCRDwEvF9m1dnA0UDkykYA10bEnIiYAkwGBktaFlg0Ih6NiAAuB3bO1bksez0e2Kr2rkGtauubWeUOPRT22Qd+/3u4886iozEzM7OmKmzMgqSdgP9ExHMlq/oAU3M/T8vK+mSvS8vnqRMRXwMfAb3K7Lea+mZWIQnGjoW11oKf/ARef73oiMzMzKwpCnl0qqSFgN8Cw8qtLlMW9ZTXV6eS/VZaP20ojSZ1ZaJ3797U1NSU28xawezZs93+Bam07Y8+egEOPng9hg2bw7nnPs0CC3gQQ3PwuV8st39x3PbFcvt3TkXNs7ASMAB4Luvt0xd4WtJg0hX/frlt+wLTs/K+ZcrJ1ZkmqRuwGN/v9lRtfQAiYiwwFmDgwIExdOjQit6wNb+amhrc/sVoTNv37Ak//vF8XHnl5lxxRbrrYNXxuV8st39x3PbFcvt3ToV0Q4qI5yNiqYjoHxH9SV/WfxgRbwMTgJHZE4oGkAYiT4yIGcAnkoZk4wn2B27NdjkBGJW93h24PxuXkD9mVfXNrGm22w5OOgmuugrOPbfoaMzMzKwxWuvRqdcAjwIDJU2TdFBd20bEi8D1wEvAXcAvImJutvoQ4CLSoOfXgNqhk+OAXpImA4cDx+SO/Wxu942ub2bV++1vYccd4fDD4ZFHio7GzMzMKtUq3ZAiYu8G1vcv+XkMMKbMdk8Cg8qUfwHsUce+16mmvplVr0sXuOIK2GAD2GMPeOqpNCeDmZmZtW2ewdnMWsVii8FNN8HHH6eE4csvi47IzMzMGuJkwcxazaBBcPHF8M9/pi5JZmZm1rY5WTCzVrXXXnDEEXDeeXD55UVHY2ZmZvVxsmBmre6002DoUDj4YHjmmaKjMTMzs7o4WTCzVtetG1x3HSy5JOy6K7z3XtERmZmZWTlOFsysEEstBePHw/TpsM8+MHduw3XMzMysdTlZMLPCbLhhmqjt7rvhhBOKjsbMzMxKOVkws0L913/BQQfBmDFwyy1FR2NmZmZ5dU7KJunACvfxdUT4mSZm1mTnngvPPQf77w9PPAEDBxYdkZmZmUH9MziPBR6uYB8bAE4WzKzJFlgAbrwR1lsvDXh+7DHo0aPoqMzMzKy+ZOHziNiioR1I+qAZ4zGzTmr55dMTkrbZBg48EK6/HqSiozIzM+vc6huz8MMK97FBcwRiZrbllmkOhvHj4cwzi47GzMzM6kwWIuLVSnYQEZObLxwz6+yOPBJ23x2OOQbuv7/oaMzMzDq3+rohfUtSN2BvYF1gkfy6iBjdAnGZWSclwcUXw0svwV57wVNPpS5KZmZm1voqfXTqlcAxwDfAOyWLmVmz6tEDbroJvvwSdtsNvvii6IjMzMw6p4ruLADDgX4R8UlLBmNmVmvgQLj8cth5Zzj0ULjooqIjMjMz63wqvbPwErBESwZiZlZqxAj47W9h3Di48MKiozEzM+t8Kr2zsC9wkaR7KOl65AnZzKwlnXQSPPlkuruw1lqw4YZFR2RmZtZ5VJosHABsBvQEPs+VB56QzcxaUNeucPXVacK23XZLA56XXrroqMzMzDqHSpOF/wbWjYiXWzIYM7NyllgiDXjeeGMYORLuvRe6VfrpZWZmZk1W6ZiFd4C3WjIQM7P6rLsujB0LNTVpDgYzMzNreZVemzsbuErSacDM/IqIeL3ZozIzK2O//WDiRDjrLNhggzQPg5mZmbWcSu8snAfsBPwTmJxbKprlWdLFkmZKeiFX9j+SXpH0L0k3S1o8t+5YSZMlTZK0ba58PUnPZ+vOkaSsvLuk67LyxyX1ryOOquqbWfHOOit1RzrwQHjhhYa3NzMzs6arKFmIiC51LF0rPM6lpLka8u4FBkXEWsC/gWMBJK0OjATWyOqcL6n2OBcAo4FVsqV2nwcBH0TEyqS7IKfXEUe19c2sYPPPDzfcAIsuCrvuCh9+WHREZmZmHVeldxaqEhEPAe+XlN0TEV9nPz4G9M1ejwCujYg5ETGFdAdjsKRlgUUj4tGIqH0K0865Opdlr8cDW9XeNahVbX0zazuWWy4lDFOmwP77wzffFB2RmZlZx1RnsiDpsrrWlWx3STPEcSBwZ/a6DzA1t25aVtYne11aPk+dLAH5COhVcoxq65tZG7LppvCnP8Ftt8GYMUVHY2Zm1jHVN8B5d0mXAg1dYd8F+GlTA5D0W+Br4KraojKbRT3l9dWZ51BV1k8bSqNJXZno3bs3NTU15TazVjB79my3f0HaStsPGgTbbLMqJ5ywNPPP/zwbbvh+w5U6gLbS/p2V2784bvtiuf07p/qShXeBiyvYx9tNPbikUcAOwFZZ1yBIV/z75TbrC0zPyvuWKc/XmSapG7AYJd2emqE+ABExFhgLMHDgwBg6dGglb9VaQE1NDW7/YrSltt9wwzTg+bTT1uLJJ2GllYqOqOW1pfbvjNz+xXHbF8vt3znV2Q0pIvpHxIAKllWbcmBJw4HfADtFxGe5VROAkdkTigaQBiJPjIgZwCeShmTjCfYHbs3VGZW93h24P5d81L6fquqbWdu00EJpwjYpDXj+7LOG65iZmVllWmWAs6RrgEeBgZKmSToIOBfoAdwr6VlJfwWIiBeB64GXgLuAX0TE3GxXhwAXkQY9v8Z34xzGAb0kTQYOB76dsknSs7lQGl3fzNq+FVeEq6+G55+H0aPBqb6ZmVnzqHRStqpExN5lisfVs/0Y4HtDFiPiSWBQmfIvgD3q2Nc61dQ3s/Zh+HA4+WQ4/vjUNemXvyw6IjMzs/avVe4smJm1huOOg512gsMPh4cfLjoaMzOz9s/Jgpl1GF26wOWXw4ABsMceMH16w3XMzMysbhUnC5JWk3S8pPOyn1eVtFbLhWZm1niLLQY33wyzZ8Puu8OXXxYdkZmZWftVUbIgaQ/gQdLkZftlxYsAf2qhuMzMmmyNNeDii+HRR1OXJDMzM2uaSu8snAwMi4ifAbVPJnoOWLtFojIzq9Kee8KRR8J558FlFc1Hb2ZmZqUqTRaWIiUH8N3MxkEdsxybmbUFf/wjbLEF/Oxn8MwzRUdjZmbW/lSaLDzFd92Pao0EJjZvOGZmzadbN7j2WlhyyTRh23vvFR2RmZlZ+1JpsnAYcIqkB4GFJd0N/AH4dYtFZmbWDJZaCm68MT0Zaa+9YM6coiMyMzNrPypKFiLiFWBV4Dzgd8AlwJoR8WoLxmZm1iwGD4axY+G++2C//WDu3IbrmJmZWSNmcI6Iz4DrWzAWM7MWM2oUzJqVBj337Al//StIRUdlZmbWtlWULEh6mPKDmecA04CbIuK25gzMzKy5HXEEvP8+nHoqLLFEGgBtZmZmdat0zEIN0J8018KV2b8rAE8C7wAXSzq6BeIzM2tWp5wCBx8Mp50G//M/RUdjZmbWtlXaDWkYsG1EvFxbIOkq4LKI2FDSTcC1wBktEKOZWbOR0twLH34IRx+duiT9v/9XdFRmZmZtU6XJwqrA6yVlbwIDASJioqSlmjMwM7OW0rUrXH45fPRRusvQsyfstlvRUZmZmbU9lXZDegi4RNLKkhaQtDJwIfAIgKQ1gRktFKOZWbObf34YPx6GDIGf/ATuvbfoiMzMzNqeSpOFUdm2LwGfAi8CXYEDsvVfAns3d3BmZi1p4YXh9tth1VVhl13gsceKjsjMzKxtqXSehfcjYiSwALAcsGBE7B0Rs7L1kyLiyRaM08ysRfTsCXffDcssA9tvDy+8UHREZmZmbUeldxaQtBiwPrAGMFTSlpK2bLHIzMxayTLLpG5ICywAw4bBlClFR2RmZtY2VDrPwgGk2ZtnA5/lVgWwYvOHZWbWugYMgHvugc03h222gUceSUmEmZlZZ1bpnYUxwO4RsXREDMgtThTMrMMYNAjuuAPefhu23RY++KDoiMzMzIpVabLQDbinJQMxM2sLhgyBW26BV16BHXaATz8tOiIzM7PiVJosnA78TlLFYxzMzNqrrbeGq69OT0fafXf48suiIzIzMytGpV/+fw38DvhE0lv5pZLKki6WNFPSC7myJSTdK+nV7N+euXXHSposaZKkbXPl60l6Plt3jiRl5d0lXZeVPy6pfx1xVFXfzDqP3XaDsWPhrrtg//1h7tyiIzIzM2t9lSYL+wJbA9sD+5UslbgUGF5SdgxwX0SsAtyX/Yyk1YGRpKcuDQfOl9Q1q3MBMBpYJVtq93kQ8EFErAycTboTUk619c2sEznoIDjjDLjuOvjFLyCi6IjMzMxaV0VPQ4qIB6s5SEQ8VOZq/QhgaPb6MqAG+E1Wfm1EzAGmSJoMDJb0BrBoRDwKIOlyYGfgzqzOidm+xgPnSlLEd/+1S1q2mvpm1jkddRS8/z6cdhr06gVjxhQdkZmZWeupKFkAkLQOsBmwJKDa8oj4fROPvXREzMj2MUPSUll5HyA/j+q0rOyr7HVpeW2dqdm+vpb0EdALmJXbvk+V9c2skzr11JQwnHoqLLEEHHFE0RGZmZm1jkrnWRhN6p5zD7Ad6Wr8MODWFohJZcqinvL66lSy30rrpw1TW4wG6N27NzU1NeU2s1Ywe/Zst39BOmPb77kn/Pvfq3PkkUvxzjuvsP32bxcWS2ds/7bE7V8ct32x3P6dU6V3Fo4GhkfEw5I+iIhdJG1HGlvQVO9IWja7q7AsMDMrnwb0y23XF5ielfctU56vM01SN2Ax4P2S41VbH4CIGAuMBRg4cGAMHTq0ojdrza+mpga3fzE6a9tvthnsuCOcddaqDBmyKrvuWkwcnbX92wq3f3Hc9sVy+3dOlQ5wXioiHs5efyOpS0TcCexYxbEnAKOy16P47i7FBGBk9oSiAaSByBOzLkufSBqSPcVo/5I6tfvaHbi/dLxBtfXNzOafH266CTbcEPbeG/7+96IjMjMza1mVJgvTcgOU/w2MkLQZUNHTxyVdAzwKDJQ0TdJBwGnANpJeBbbJfiYiXgSuB14C7gJ+ERG1Dy08BLgImAy8RuoOBTAO6JUNhj6c7MlK2bGfzYXS6PpmZnkLLwx/+xsMHAg77wyPP150RGZmZi2n0m5IZwCrAW8AJ5OeGDQ/cFgllSNi7zpWbVXH9mOA7z1zJCKeBAaVKf8C2KOOfa1TTX0zs1I9e8Ldd8Omm8L228NDD8EaaxQdlZmZWfNr8M5C1mXnIeBegKz7UU+gZ0Rc0LLhmZm1TcsuC/feC927w7Bh8MYbRUdkZmbW/BpMFrK++88D3+TKvoyI2S0ZmJlZW7fiinDPPfD557DNNvB2cQ9IMjMzaxGVjll4BvhBSwZiZtYeDRqUxjBMnw7Dh8OHHxYdkZmZWfOpdMxCDXCXpEtJk5d9+6SgiLi4+cMyM2s/NtoIbr4ZdtghLXfdBYssUnRUZmZm1as0WdgEmAL8qKQ8ACcLZtbpDRsGV18NI0fCdtvBHXdAjx5FR2VmZladipKFiNiipQMxM2vvdt8drrkmzcEwfDjceScsumjRUZmZmTVdpWMWzMysAnvsAdddBxMnpoTh44+LjsjMzKzpnCyYmTWz3XZLCcMTT8C228JHHxUdkZmZWdM4WTAzawG77go33ABPPZXGM/gpSWZm1h45WTAzayE77wzjx8Mzz6R5GD74oOiIzMzMGqfSpyEhaTFgIDDPAwEj4v7mDsrMrKPYaSe46abUNWmbbdKszz17Fh2VmZlZZSpKFiQdAJwHzAY+y60KYMXmD8vMrOPYYYeUMOy6K2y9dUoYllii6KjMzMwaVmk3pDHA7hGxdEQMyC1OFMzMKvDjH8Mtt8CLL8JWW8F77xUdkZmZWcMqTRa6Afe0ZCBmZh3ddtvBrbfCyy+nhGHWrKIjMjMzq1+lycLpwO8keUC0mVkVtt0WJkyASZNgyy3h3XeLjsjMzKxulX75/zXwO+ATSW/llxaMzcysQxo2DG67DV59NSUMM2cWHZGZmVl5lT4Nad8WjcLMrJPZemv429/S4OcttoD774elly46KjMzs3lVlCxExIMtHYiZWWez5ZZwxx1p8HNtwrDMMkVHZWZm9p2KuiFJmk/SSZJel/RF9u9JkuZv6QDNzDqyoUNTwvDmmylhmDGj6IjMzMy+U+mYhTOArYGfAWtn/25JGvhsZmZV+NGP4M47YepU2HxzeOONoiMyMzNLKk0W9gB2ioh7ImJSRNwD7ALs2XKhmZl1HptvniZrmzULNt00PV7VzMysaJUmC2pkuZmZNdJGG8GDD8LXX8Nmm8GTTxYdkZmZdXaVJgs3ALdJ2lbSapKGA7cA11cbgKRfS3pR0guSrpG0gKQlJN0r6dXs35657Y+VNFnSJEnb5srXk/R8tu4cSWUTmWrrm5m1pLXWgkcegR490gDoB/14CTMzK1ClycLRwN+B84CngHOBB4DfVHNwSX2Aw4D1I2IQ0BUYCRwD3BcRqwD3ZT8jafVs/RrAcOB8SV2z3V0AjAZWyZbhZY5XVX0zs9aw8sopYejbF4YPh9tvLzoiMzPrrCpKFiLiy4j4fUSsHBELZf8eHxFzmiGGbsCCkroBCwHTgRHAZdn6y4Cds9cjgGsjYk5ETAEmA4MlLQssGhGPRkQAl+fq5FVb38ysVfTpAw89BGuuCbvsAldfXXREZmbWGdU5z4KkzSPioez1lnVtFxH3N/XgEfEfSWcCbwGfA/dExD2Slo6IGdk2MyQtlVXpAzyW28W0rOyr7HVpealq65uZtZoll4T77oMRI2DffeGjj+CQQ4qOyszMOpP6JmU7HxiUvR5XxzYBrNjUg2djEUYAA4APgRsk1TdbdLlxBFFPeXPXR9JoUnclevfuTU1NTdlAreXNnj3b7V8Qt33rOuaYLsyZszo///mSPPXU64wY4fYvks//4rjti+X275zqTBayMQS1rwe00PG3BqZExLsAkm4CNgbekbRsdldhWWBmtv00oF+ufl9St6Vp2evS8lLV1icixgJjAQYOHBhDhw5t+F1ai6ipqcHtXwy3fevbYgs48EAYN25FZs/uxjXXLI8fw1AMn//FcdsXy+3fOVU6g/OtdZTfVOXx3wKGSFooe/rQVsDLwARgVLbNKKD2+BOAkZK6SxpAGog8Meuy9ImkIdl+9s/Vyau2vplZIeabDy67DA49FK67bnlGj4a5c4uOyszMOrr6uiHlbVFH+dBqDh4Rj0saDzwNfA08Q7pqvwhwvaSDSAnFHtn2L0q6Hngp2/4XEVH73+UhwKXAgsCd2YKknUhPW/p9U+qbmbUVXbrAOefARx+9wUUX9efDD+GKK2CBBYqOzMzMOqp6kwVJJ2cv58+9rrUi8Ga1AUTECcAJJcVzSHcZym0/BhhTpvxJvhtjkS+fQLqj0KT6ZmZtiQQHHvgG667bn8MPh3ffhVtugcUXLzoyMzPriBrqhtQvW7rkXvcj9emfSnbF38zMWtevfw3XXAP//CdsuilMnVp0RGZm1hHVe2chIn4KIOmfEXFh64RkZmaVGDkSll4adt4ZNtoI7rwzzctgZmbWXCqdwXmOpLXyBZLWlrRfC8RkZmYV2mKLNNszwGabgZ9qaGZmzanSZOEPpG5HeVOBU5o3HDMza6w114RHH4W+fWHbbeG664qOyMzMOopKk4VFgY9Lyj4CFm/WaMzMrEn69YOHH4YhQ1L3pD/9qeiIzMysI6g0WXgJ2K2kbBfSnAhmZtYG9OwJd98Nu+8ORxwBhx8O33xTdFRmZtaeVTrPwm+AOyTtBbwGrEx6tOn2LRWYmZk13gILpG5Ihx8OZ58N06alydwWXLDoyMzMrD2q6M5CRDxCmoPgCWBhYCIwKCL+0YKxmZlZE3TpkhKFM8+E8ePTIOh33ik6KjMza48qvbNARLwl6Qxg6YiY0YIxmZlZlaTUFWnFFWGffdJYhttvhzXWKDoyMzNrTyq6syBpcUlXA18Ak7OynST5aUhmZm3YLrvAQw/BF1/AxhvDvfcWHZGZmbUnlQ5w/ivp6UcrAF9mZY8Ce7VEUGZm1nzWXx8mToT+/WG77eBCT7FpZmYVqjRZ2Ao4LOt+FAAR8S6wVEsFZmZmzadfvzR527BhMHo0HHWUn5RkZmYNqzRZ+AhYMl8gaXnAYxfMzNqJHj1gwgQ49NA0+Hm33WD27KKjMjOztqzSZOEi4EZJWwBdJG0EXEbqnmRmZu1Et27wv/8Lf/lLShw23hhef73oqMzMrK2qNFk4HbgeOA+YD7gYuBX4SwvFZWZmLeiww+Cuu9I8DBtsAPffX3REZmbWFlWaLCwdEX+OiNUjYuGIWC0i/gws3YKxmZlZC9pmmzTweZll0liGc86BiKKjMjOztqTSZOHfdZS/1FyBmJlZ61t5ZXjsMdhhB/jv/4aDDoI5c4qOyszM2opKkwV9r0BaFPCzNMzM2rkePeCmm+D3v4dLLoGhQ2GGH19hZmY0kCxImirpLWBBSW/lF9KTkG5pjSDNzKxldekCJ50EN94Izz+f5mZ4/PGiozIzs6J1a2D9vqS7CncA++XKA3gnIia1VGBmZtb6dt0VVlkFRoyAzTeHsWNh1KiiozIzs6LUmyxExIMAkpaMiM9aJyQzMyvSmmvCE0/AnnvCAQfAs8/C//xPeuyqmZl1LpWOWZgraYyk1yV9BCBpmKRDWzA2MzMrSK9ecPfd8KtfwZ//DMOHw3vvFR2VmZm1tkqThT8Dg4B9SF2QAF4EDqk2AEmLSxov6RVJL0vaSNISku6V9Gr2b8/c9sdKmixpkqRtc+XrSXo+W3eOpO8Nym6O+mZmnUW3bnD22WnQ88MPp3EMTz5ZdFRmZtaaKk0WdgZ+EhGPkj0BKSL+A/Rphhj+AtwVEasCawMvA8cA90XEKsB92c9IWh0YCawBDAfOl9Q1288FwGhglWwZXnqgauubmXVGBxyQkoVvvoFNNoG//tXzMZiZdRaVJgtfUjK+QVJvoKqb0tnjVzcHxgFExJcR8SEwArgs2+wyUrJCVn5tRMyJiCnAZGCwpGWBRSPi0YgI4PJcnbxq65uZdUqDB8PTT8NWW8Ehh8B++8GnnxYdlZmZtbRKk4UbgMskDQDIvlyfC1xb5fFXBN4FLpH0jKSLJC1MmjF6BkD271LZ9n2Aqbn607KyPtnr0vJS1dY3M+u0evWC22+HU06Ba65JCcQrrxQdlZmZtaRKn21xHHAG8DywEPAqcCFwUjMc/4fALyPicUl/IetyVIdy4wiinvLmro+k0aTuSvTu3ZuampqygVrLmz17ttu/IG77YhXd/ptsAmecsTinnLI6667blaOOmsSWW84sLJ7WVnT7d2Zu+2K5/TunipKFiPgS+BXwq6z70aysu061pgHTIqJ26p/xpGThHUnLRsSM7C7GzNz2/XL1+wLTs/K+ZcrLHa+a+kTEWGAswMCBA2Po0KENvEVrKTU1Nbj9i+G2L1ZbaP+hQ2HkSNhrL/jDH1bn/fdX56yzoHv3QsNqFW2h/Tsrt32x3P6dU0XdkCStLulgSccCuwKrNcfBI+JtYKqkgVnRVsBLwASgdhqgUcCt2esJwEhJ3bMuUasAE7OuSp9IGpI9xWj/XJ28auubmVmmTx944AE44gg47zzYbDN4882iozIzs+ZU752F7IvzONIX9mmkq+19gOUkXQEc2Ax3GH4JXCVpfuB14KekJOZ6SQcBbwF7AETEi5KuJyUUXwO/iIi52X4OAS4FFgTuzBYk7QSsHxG/b0p9MzOr23zzwZlnwsYbw09/CuuuC+PGwS67FB2ZmZk1h4a6IY0GhgJDIuKJ2kJJGwDXAAcDf60mgIh4Fli/zKqt6th+DDCmTPmTpLkgSssnkO4oNKm+mZk1bNddYe21U9ekXXdNT0w66yxYcMGiIzMzs2o01A1pP+CwfKIAkP38q2y9mZkZK60E//gHHHUUXHABbLghvPRS0VGZmVk1GkoWVgcerGPdg9l6MzMzAOafH844A+66C955J836PHasJ3EzM2uvGkoWukbEJ+VWZOWVztNgZmadyLbbwnPPwaabwsEHp6cmffhh0VGZmVljNTRmYT5JW1B+HoJK6puZWSe1zDLpDsOZZ8JvfwsTJ8IVV6SnJpmZWfvQ0Jf9mcDFDaw3MzMrq0sXOPpo+NGPYJ990r9HHw0nndQ55mQwM2vv6k0WIqJ/K8VhZmYd2IYbwrPPpjkZTj893XG48koY5GfQmZm1aR5zYGZmrWKRReD//g9uuw1mzID11kuPV/3mm6IjMzOzujhZMDOzVrXDDvDCC7D99nDkkbDllp752cysrXKyYGZmra53b7jpJrjkEnj6aVhrLbj8cj9i1cysrXGyYGZmhZDggAPgX/9Ksz+PGgW77w7vvlt0ZGZmVsvJgpmZFap/f3jggTSZ2+23w+qrw3XX+S6DmVlb4GTBzMwK17UrHHVU6pI0YACMHAm77QZvv110ZGZmnZuTBTMzazPWWAP++c90l+GOO9JdBo9lMDMrjpMFMzNrU7p1S3cZnnsuJQujRqUnKE2bVnRkZmadj5MFMzNrkwYOhAcfhL/8BWpq0l2HsWM9L4OZWWtysmBmZm1W165w2GHw/PNpEreDD4bNN4cXXyw6MjOzzsHJgpmZtXkrrgj33ZfmZXj5ZVhnHfjtb+Hzz4uOzMysY3OyYGZm7ULtvAyvvAL77AOnngprrgn33lt0ZGZmHZeTBTMza1d694ZLL4X770/dlIYNS8nDzJlFR2Zm1vE4WTAzs3Zpiy3SE5NOOAHGj4dVV4ULL/QAaDOz5uRkwczM2q0FFoATT0xJw1prwejRMGQITJxYdGRmZh2DkwUzM2v3Vl0VHngArrwyzcew4YZw4IHwzjtFR2Zm1r61iWRBUldJz0i6Pft5CUn3Sno1+7dnbttjJU2WNEnStrny9SQ9n607R5LqOFZV9c3MrG2S0tiFSZPSpG5XXgk/+AH8+c/w1VdFR2dm1j61iWQB+G/g5dzPxwD3RcQqwH3Zz0haHRgJrAEMB86X1DWrcwEwGlglW4aXHqTa+mZm1vb16AFnnJHmZthoI/j1r2HdddOAaDMza5zCkwVJfYEfAxflikcAl2WvLwN2zpVfGxFzImIKMBkYLGlZYNGIeDQiArg8Vyev2vpmZtZODBwId94Jt9wCn30GW20Fe+wBb75ZdGRmZu1H4ckC8GfgaCD//IqlI2IGQPbvUll5H2BqbrtpWVmf7HVpealq65uZWTsiwYgRacbnk06C229PScRvfgMffVR0dGZmbV+3Ig8uaQdgZkQ8JWloJVXKlEU95c1dH0mjSd2V6N27NzU1NWUDtZY3e/Zst39B3PbFcvs3zeabw6qrdueiiwZwxhnL8H//9yWjRr3JjjtOp1u3sh/5Zbn9i+O2L5bbv3MqNFkANgF2krQ9sACwqKQrgXckLRsRM7IuQrVT7UwD+uXq9wWmZ+V9y5SXqrY+ETEWGAswcODAGDp0aAVv01pCTU0Nbv9iuO2L5favzp57wtNPwxFHzM8556zCXXetwhlnwE47pTsRDXH7F8dtXyy3f+dUaDekiDg2IvpGRH/SwOP7I2JfYAIwKttsFHBr9noCMFJSd0kDSAORJ2ZdlT6RNCR7itH+uTp51dY3M7MO4Ic/TAOeJ0xICcLOO6dJ3p58sujIzMzalrYwZqGc04BtJL0KbJP9TES8CFwPvATcBfwiIuZmdQ4hDZKeDLwG3AkgaSdJJze1vpmZdUwS7LhjemrS+efDSy/BBhukx69Onlx0dGZmbUPR3ZC+FRE1QE32+j1gqzq2GwOMKVP+JDCoTPkE0h2FJtU3M7OObb754JBD4Cc/gdNPT/MyXH89HHQQHH889PHjLsysE2urdxbMzMxa1WKLwamnwmuvwcEHw8UXw8orw5FHwqxZRUdnZlYMJwtmZmY5yy4L556bZoLeay84+2xYcUU48UT4+OOiozMza11OFszMzMoYMAAuvTSNaRg2LM3TsOKKcM01/Zg9u+jozMxah5MFMzOzeqy+OowfD088AeuvD2PHrkT//vDHP/pOg5l1fE4WzMzMKrD++nDXXXDuuU8zeDAcdxz07w9/+AN8+GHR0ZmZtQwnC2ZmZo2wxhofc8cdMHEibLop/P73KWk44QR4//2iozMza15OFszMzJpggw3SpG5PPw1bbgknn5yShuOOg5kzi47OzKx5OFkwMzOrwrrrwk03wXPPwfDhcNppsMIKae4GT+5mZu2dkwUzM7NmsNZaaTK3l1+G/fZL8zT84Aewxx5pcLSZWXvkZMHMzKwZDRwIY8fCG2/AMcfAvffC4MGwxRZw550QUXSEZmaVc7JgZmbWApZdNs0IPXUqnHVW6pK0/faw9tpwySXwxRdFR2hm1jAnC2ZmZi2oRw84/HB47TW47LJ0Z+HAA6FfP/jd7+A//yk6QjOzujlZMDMzawXzzw/77w//+hfcdx9sskm689C/P+y9Nzz2mLsomVnb42TBzMysFUnpUau33JK6Jh12WBrLsNFGsOGGcOWV8OWXRUdpZpY4WTAzMyvIiium8QzTpsF558HHH6cnKfXtC7/5Teq6ZGZWJCcLZmZmBVtkEfj5z+Gll+Cuu1IXpbPOgpVXhm22gfHj4auvio7SzDojJwtmZmZtRJcusO22cPPN8NZbaVbof/87zdXQr1+aHXrKlKKjNLPOxMmCmZlZG7TccnD88fD66/C3v6XxDKefDiutlBKK667z41fNrOU5WTAzM2vDunZN8zPceiu8+SaccEKaJXrkyDSXwyGH+ElKZtZynCyYmZm1E337pmRhypQ0M/QOO6S5GzbaCFZbDf74xzRY2sysuThZMDMza2e6doWtt4YrroC334Zx42CppdKYhuWXh2HD4KqrYPbsoiM1s/bOyYKZmVk7tuiiaUbohx5K8zYcfzy8+irsu29KIPbaK83p4PENZtYUhSYLkvpJekDSy5JelPTfWfkSku6V9Gr2b89cnWMlTZY0SdK2ufL1JD2frTtHkuo4ZlX1zczM2qqVVoKTTkrzMzz0EPz0p/DAA7DLLrD00unnu++Gr78uOlIzay+KvrPwNXBERKwGDAF+IWl14BjgvohYBbgv+5ls3UhgDWA4cL6krtm+LgBGA6tky/DSg1Vb38zMrD3o0gU22yxN9DZ9ekoQdt0VbroJhg9PA6N//nOoqYG5c4uO1szaskKThYiYERFPZ68/AV4G+gAjgMuyzS4Dds5ejwCujYg5ETEFmAwMlrQssGhEPBoRAVyeq5NXbX0zM7N2pVu3NIbhkkvgnXfSHA5bbQWXXgpbbJESh//3/+COO2DOnKKjNbO2plvRAdSS1B9YF3gcWDoiZkBKKCQtlW3WB3gsV21aVvZV9rq0vFS19ZE0mnQHgt69e1NTU9Pwm7MWMXv2bLd/Qdz2xXL7F6u9t//ii8PPfgajRnVl4sQlePjhJbn22l6MG9eNhRb6miFD3mOzzWax4Ybvs+CCbeu2Q3tv+/bO7d85tYlkQdIiwI3AryLi43qGC5RbEfWUN3d9ImIsMBZg4MCBMXTo0LKBWsurqanB7V8Mt32x3P7F6kjtv9126d85c+D+++Hmm7txyy1Lc//9S9O9e7ojMWJEmudh2WWLjRU6Vtu3R27/zqnoMQtImo+UKFwVETdlxe9kXYPI/p2ZlU8D+uWq9wWmZ+V9y5SXqra+mZlZh9O9e0ocxo6FGTPgwQfT3YfnnktdlJZbDtZfP83xMHEifPNN0RGbWWsp+mlIAsYBL0fEn3KrJgCjstejgFtz5SMldZc0gDQQeWLWZekTSUOyfe6fq5NXbX0zM7MOrWtX2Hxz+POf4Y03UsJw6qkpoTjlFNhww3SX4YAD4IYb4KOPCg7YzFpU0d2QNgH2A56X9GxWdhxwGnC9pIOAt4A9ACLiRUnXAy+RnqT0i4io7VB5CHApsCBwZ7YgaSdg/Yj4fVPqm5mZdVYSrLVWWo49Ft57D+66C/72N5gwIc0e3a0bbLIJbLNNWtZbLyUcZtYxFJosRMQjlB8vALBVHXXGAGPKlD8JDCpTPoF0R6FJ9c3MzCzp1Qv22SctX38Njz0Gt9+eHs36u9+lZfHFYcstv0seVlqp6KjNrBpF31kwMzOzdqhbN9h007ScdhrMnAn33Qd//zvce2+a0wFgwICUNGy9NfzoR2lWaTNrP5wsmJmZWdWWWgr23jstEfDvf6ek4e9/h2uvTYOnAVZbLSUNtUtbeMqSmdXNyYKZmZk1KwkGDkzLoYemLktPPJGesvTgg3DVVfDXv6ZtV1553uRh+eWLjd3M5uVkwczMzFpUt26w0UZpOeaYlDw8++x3ycONN8K4cWnbFVaAjTf+bllrrVTfzIrhPz8zMzNrVd26pXkb1l8fjjgC5s6FF15IicPDD6d/r7kmbbvQQjB4cEocFllkCQYPTmVm1jqcLJiZmVmhunaFtddOy2GHpTEPU6fCP/+ZlkcfhdNPh7lz12LffZ0smLUmJwtmZmbWpkhp7MLyy8PIkans00/h4oufoV+/dYsNzqyTKXQGZzMzM7NKLLwwrLmmp4s2a21OFszMzMzMrCwnC2ZmZmZmVpaTBTMzMzMzK8vJgpmZmZmZleVkwczMzMzMynKyYGZmZmZmZTlZMDMzMzOzspwsmJmZmZlZWU4WzMzMzMysLCcLZmZmZmZWliKi6BjaLUmfAJNa+bCLAdXOd9/YfVSyfUPb1LW+XHmlZUsCsxqIq7lV2/4t0fYNbdfYdZWUtce2b8o+Wvvcr6vc7d+0bVri3IfWb/+22vYNbefPnqbvw589dcfQGvtoq589AyOiRwNxNb+I8NLEBXiygGOObe19VLJ9Q9vUtb5ceSPK2l37t0TbN7RdY9dVUtYe276l2r85z323f/O2f0uc+0W0f1tt+yLa3+d+5dv4s6e49u8onz21i7shtT+3FbCPSrZvaJu61pcrr7SsCNXG0RJt39B2jV3XVtu/M5z7dZW7/Zu2jc/9pu/Dnz3faavt78+eltuHP3ty3A2pCpKejIj1i46js3L7F8dtXyy3f7Hc/sVx2xfL7V+sotrfdxaqM7boADo5t39x3PbFcvsXy+1fHLd9sdz+xSqk/X1nwczMzMzMyvKdBTMzMzMzK8vJQh0kDZc0SdJkSceUWS9J52Tr/yXph5XWtfpV2fYXS5op6YXWjbrjqKD9V5X0qKQ5ko5sTF2rX0Pnrz93mle59pa0hKR7Jb2a/duzjrpl27vS+p1VY9tc0rFZG0+StG0d+6yqfkfWXO0taT1Jz2frzpGkOo5XVf2OoKXbXFJ3Sddl5Y9L6l9HHFXVn0cRj2Bq6wvQFXgNWBGYH3gOWL1km+2BOwEBQ4DHK63rpWXaPlu3OfBD4IWi30t7XCps/6WADYAxwJGNqeulwfav9/z1507LtzdwBnBM9voY4PQy9eps70rqd+alMW0OrJ61bXdgQNbmXcvss6r6HXlprvYGJgIbZZ89dwLblTlWVfU7ytLSbQ78HPhr9nokcF0dcVRVP7/4zkJ5g4HJEfF6RHwJXAuMKNlmBHB5JI8Bi0tatsK6Vrdq2p6IeAh4v1Uj7lgabP+ImBkRTwBfNbau1a+C89efO82ojvYeAVyWvb4M2LlM1frau5L6nVYj23wEcG1EzImIKcBkUtuXqrZ+h9Uc7Z19xiwaEY9G+oZ5OeXP62rrdwit0Ob5fY0Htiq9U1Nt/VJOFsrrA0zN/TwtK6tkm0rqWt2qaXurXjVt699Ly/PnTstbOiJmAGT/LlVmm/rau5L6Nq+62qzS87ra+p1NY9urT/a6tLxUtfU7suZs82/rRMTXpFmee5Ucr9r683CyUF65DKv0sVF1bVNJXatbNW1v1aumbf17aXn+3Gkb3N6to9p29u+pcar9fPHnU+M1pc2q/Z7U6N+Hk4XypgH9cj/3BaZXuE0lda1u1bS9Va+atvXvpeX5c6flvVPbrTH7d2aZbepr70rq27zqarNKz+tq63c2jW2vadnr0vJS1dbvyJqzzb+tI6kbsBjf7/ZUbf15OFko7wlgFUkDJM1PGgAyoWSbCcD+SoYAH2W3liqpa3Wrpu2tetWcvz73W54/d1reBGBU9noUcGuZbepr70rq27zqarMJwMjs6S0DgFVIgzabu35n06j2yj5jPpE0JOvbvj/lz+tq63dkzdnm+X3tDtyfjUv4VrX1v6elRoO394X01JF/k0am/zYr+xnws+y1gPOy9c8D69dX10urtf01wAzS4NtpwEFFv5/2tlTQ/stkbfsx8GH2etG66nppVNt/7/z1506rt3cv4D7g1ezfJbJtlwPuaKi966rvpfFtnm3/26yNJ5F7gg5wUe3535T6nWVpxvZeH3ghW3cu303quxNwclPrd8SlFdp8AeAG0mDoicCKuTrPVlO/rsUzOJuZmZmZWVnuhmRmZmZmZmU5WTAzMzMzs7KcLJiZmZmZWVlOFszMzMzMrCwnC2ZmZmZmVpaTBTMzMzMzK8vJgpmZmZmZleVkwcysAJIulXRK7ucXJQ0tLqKGSXpD0tZt+Tj11S1t8zLrQ9KnksY05dhFknS/pC8kPVJ0LGbWsThZMDMro7W+GNeKiDUioqa1jteQ1n7/bcjaEfHbhjaSdKykO0rKXq2jbGTu5+UkTWtMQJK6Sxon6U1Jn0h6RtJ2+W0iYkvSbN9mZs3KyYKZWSNI6lZJWWfWSdrjIWATSV0BJC0DzAf8sKRs5WzbWtsDdzXyWN2AqcCPgMWA44HrJfWv5g2YmVXCyYKZWQOyq+y/kfQv4FNJ3eooO0bSa9nV35ck7ZLbx7qSns7WXQcsUOYYW2ev69xPbtsjJf1L0keSrpO0QLaun6SbJL0r6T1J5+bqLSfpxmzdFEmH1fF+rwCWB26TNFvS0bnV69Rx3HLtUe/xsu3/k73PSZK2quA4q0mqkfRh1nVrp3p+b/W2eUMk/VbSBbmfe0r6KovlCVJysE62enPgAWBSSdlrETE9t9vtgTtybXZU9j4/ze4eLC3pzizmv0vqGRGfRsSJEfFGRHwTEbcDU4D1GvN+zMyawsmCmVll9gZ+DCweEV/XUfYasBnp6u9JwJWSlpU0P3ALcAWwBHADsFs9xyq7n5Jt9gSGAwOAtYADsivatwNvAv2BPsC1AJK6ALcBz2XlWwG/krRt6cEjYj/gLWDHiFgkIs6o77jl2gj4pr7jSRoIHApsEBE9gG2BNxp4f/Nl+7wHWAr4JXBVtq95NKHNy1kTeDb38zrApIj4IiK+BB4nJQRk/z4MPFJS9u1dhSz+zYF7c/vcDdgG+AGwI3AncBywJOn/6O8ldJKWzrZ/sZHvx8ys0ZwsmJlV5pyImBoRn9dVFhE3RMT07OrvdcCrwGBgCOkq9J8j4quIGE+6Ml1WPfspjWd6RLxP+gK9TrbNcsBR2dXoLyKidsDrBkDviDg5Ir6MiNeBC4GRNE6545Zrj4aONxfoDqwuab7sqvlrDRxnCLAIcFq2z/tJydHeZeJsVJvXoVyy8Fzu5wf5LjHYjJQsPFxS9mBu+82B5yLik1zZ/0bEOxHxn6zu4xHxTETMAW4G1s0HlCUcVwGXRcQrjXw/ZmaN5mTBzKwyUxsqk7S/pGezLjIfAoNIV4iXA/4TEZHb/M26DlTPfvLezr3+jPQluh/wZu7OR94KwHK1+8z2exywdF1x1KHccWvl26Pe40XEZOBXwInATEnXSlqugeMsB0yNiG9y694k3bko1ag2L5XdmVgJeD5XvDbzJg8PAZtK6klKjF4F/glsnJUN4vvjFeYZAA28k3v9eZmfv23f7O7QFcCXpLsyZmYtzsmCmVllor4ySSuQrpwfCvSKiMWBFwABM4A+kpSru3y5gzSwn4ZMBZZX+QHGU4EpEbF4bukREdvXsa9y77ch+ToNHi8iro6ITUmJRQCnN7D/6UC/7EtzreWB/5TZtuI2r8PqpGTjM4BsP0OZ987Co6SuYqOBfwBExMdZnKOB6RExJbf99sDfGhHDt7LjjyMlW7tFxFdN2Y+ZWWM5WTAzax4Lk77wvgsg6aekK8uQvlR+DRyWDfzdle93K6pkPw2ZSPqSfJqkhSUtIGmT3LqPs0HFC0rqKmmQpA3q2Nc7wIoVHreuWOo8nqSBkraU1B34gnQVfW4D+3wc+BQ4WtJ8SvNS7Eg2LqNEY9q8nDWBpSStJGlB4A+kpOaN2g2y7lZPAoeTuhDVeiQry49XGAB0r6Lr0AXAaqRxJJ83tLGZWXNxsmBm1gwi4iXgLNKX1HdIXzZrrzZ/CexKGgz8AbAXcFNj91NBDHNJX55XJg1QnpYdK79uHdKTdGYBF5GujJfzR+B3WReiIys5fh2x1HW87sBpWfnbpAHLxzWwzy+BnYDtsnrnA/uX+wLemDavw5rA3aQBx5NJv4vXgdI5GB7MYs9PhvZwVpbvgvRjvt8FqSLZ3aaDSW35ttITqmZL2qcp+zMzawzN253TzMysc5L0BTAHOAdYH7goIm5spn3fAZwbEU1KGCrY/72kQd0TI2KrhrY3M6tUZ5g4x8zMrEER8e08DEqzLL/cjLuvIc3D0CIiYpuW2reZdW6+s2BmZpaTPcnoHWBhDyQ2s87OyYKZmZmZmZXlAc5mZmZmZlaWkwUzMzMzMyvLyYKZmZmZmZXlZMHMzMzMzMpysmBmZmZmZmU5WTAzMzMzs7KcLJiZmZmZWVlOFszMzMzMrKz/D9a6RXfCZnueAAAAAElFTkSuQmCC\n",
"text/plain": [
"