atlaas-dynamic1.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsYAAAE6CAYAAAAV0JRKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAD4FJREFUeJzt3X+o7/d90PHnK15Tlrr1jzqsJCZjlFIsDN0khU7dGRUT\nVlzwB5rqHBTsoGg6+odkjEFP/Mv/ROhQxCyxspphoXedFu2g+bZ0ZTObm51tQgrD/CA0c6Do6D9a\n3/5xT8fp3b2532/O+d7vPTePBxxyvp/7OZ/PGz4k95nPeX/en1lrBQAAb3R3HHoAAABwKxDGAACQ\nMAYAgEoYAwBAJYwBAKASxgAAUAljAACohDEAAFR7CuOZ+aGZ+cLM/LOZ+Yv7OAcAAJynfd0xXtX/\nrt5UvbyncwAAwLnZKoxn5vGZeXVmvnzV9gdn5rmZeX5mHv3W9rXWF9Za76t+qvpH5ztkAAA4f9ve\nMX6ieuD0hpm5o/rYyfZ3Ve+fmXde9XP/s7rzrIMEAIB9u7TNTmutL87MfVdtvr/62lrrhaqZeap6\nqHpuZv5qV4L5LV2JZwAAuKVtFcbXcXf10qnPL3cllltrfar61Gv98MysM5wbAAC2ttaaG+1zljC+\n1sF3it21tPFFdXx83PHx8aGHwevk+l1crt3F5vpdXK7dxTZzwyauzrYqxcvVvac+31O9cobjAQDA\nwewSxtO33yV+pnr7zNw3M3dWD1efPs/BAQDAzbLtcm2fqL5UvWNmXpyZD6y1vlk9Un22+kr11Frr\n2V1Ofnx83Gaz2XHI3AqOjo4OPQTOwPW7uFy7i831u7hcu4tps9nsNAVmDjXPd2aWOcYAAOzbzGz1\n8N2+3nwHAAAXijAGAICEMQAAVAcOYw/fAQCwLx6+AwCAUzx8BwAAOxDGAACQMAYAgMrDdwAA3KY8\nfAcAAKd4+A4AAHYgjAEAIGEMAACVMAYAgMqqFAAA3KasSgEAAKdYlQIAAHYgjAEAIGEMAACVMAYA\ngEoYAwBAJYwBAKCyjjEAALcp6xgDAMAp1jEGAIAdCGMAAEgYAwBAJYwBAKASxgAAUAljAACohDEA\nAFR16ZAn/+7v/lPdddd39eY3v+Xbtn/wg3+3j3zkQwcaFQAAt4PNZrPTy+QO+oKP+pVr/MkXOjr6\n1Z5++vJNHxMAALefbV/wcdA7xvWea2z73epXb/ZAAAB4gzPHGAAAEsYAAFAJYwAAqIQxAABUwhgA\nACphDAAAlTAGAIBKGAMAQHXwF3wcV0cnXwAAcH4u2Cuhr3Xuyx0dPemV0AAAnIttXwltKgUAACSM\nAQCgEsYAAFAJYwAAqIQxAABUwhgAACphDAAAlTAGAIBKGAMAQCWMAQCgEsYAAFAJYwAAqIQxAABU\nwhgAACphDAAAVV067OmPq6OTLwAAOD+bzabNZrP1/rPW2t9oXuvEM6uude7LHR092dNPX77pYwIA\n4PYzM6215kb7mUoBAAAJYwAAqIQxAABUwhgAACphDAAAlTAGAIBKGAMAQCWMAQCgEsYAAFAJYwAA\nqIQxAABUwhgAACphDAAAlTAGAIBKGAMAQCWMAQCgEsYAAFAJYwAAqIQxAABUwhgAACphDAAA1R7D\neGbumplfn5kf2dc5AADgvOzzjvGj1S/s8fgAAHButgrjmXl8Zl6dmS9ftf3BmXluZp6fmUdPbX9v\n9dXqd6s51xEDAMAebHvH+InqgdMbZuaO6mMn299VvX9m3nnyxz9cvbv629XfO5+hAgDA/lzaZqe1\n1hdn5r6rNt9ffW2t9ULVzDxVPVQ9t9b6mZNtP1793jmOFwAA9mKrML6Ou6uXTn1+uSux/AfWWh9/\n7UMcn/r+6OQLAABev81m02az2fnnzhLG15o7vHY7xPEZTg8AAH/Y0dFRR0dHf/D5scce2+rnzrIq\nxcvVvac+31O9cobjAQDAwewSxtO33yV+pnr7zNw3M3dWD1efPs/BAQDAzbLVVIqZ+URXJgC/dWZe\nrD661npiZh6pPtuVwH58rfXsbqc/ztxiAAD2Yde5xrPWjtOCz8nMrGtPSb7c0dGTPf305Zs+JgAA\nbj8z01rrhu/W2Oeb7wAA4MIQxgAAkDAGAIDqbOsYn4PjPHwHAMA+ePgOAABO8fAdAADsQBgDAEDC\nGAAAKg/fAQBwm/LwHQAAnOLhOwAA2IEwBgCAhDEAAFTCGAAAKqtSAABwm7IqBQAAnGJVCgAA2IEw\nBgCAhDEAAFTCGAAAKmEMAACVMAYAgMo6xgAA3KasYwwAAKdYxxgAAHYgjAEAIGEMAACVMAYAgEoY\nAwBAJYwBAKASxgAAUHnBBwAAtykv+AAAgFO84AMAAHYgjAEAIGEMAACVMAYAgEoYAwBAJYwBAKAS\nxgAAUAljAACohDEAAFReCQ0AwG3KK6EBAOAUr4QGAIAdCGMAAEgYAwBAJYwBAKASxgAAUAljAACo\nhDEAAFTCGAAAKmEMAACVMAYAgEoYAwBAJYwBAKASxgAAUAljAACohDEAAFR16bCnP66OTr4AAOD8\nbDabNpvN1vvPWmt/o3mtE8+suta5L3d09GRPP335po8JAIDbz8y01pob7WcqBQAAJIwBAKASxgAA\nUAljAACohDEAAFTCGAAAKmEMAACVMAYAgEoYAwBAJYwBAKASxgAAUAljAACohDEAAFTCGAAAKmEM\nAACVMAYAgEoYAwBAJYwBAKASxgAAUAljAACohDEAAFR1aR8HnZl3Vj9ZvbX63Frrn+/jPAAAcF72\ncsd4rfXcWutD1d+q3rOPcwAAwHnaKoxn5vGZeXVmvnzV9gdn5rmZeX5mHr3qz/5K9e+qz5zfcAEA\nYD+2vWP8RPXA6Q0zc0f1sZPt76refzKFoqq11i+ttd5X/dg5jRUAAPZmqznGa60vzsx9V22+v/ra\nWuuFqpl5qnqoem5mfqj6a9Wbqn9/juMFAIC9OMvDd3dXL536/HJXYrm11uerz9/4EMenvj86+QIA\ngNdvs9m02Wx2/rmzhPFcY9va7RDHZzg9AABc8ba3fU+vvvrCmY5xljB+ubr31Od7qlfONBoAAHgd\nrkTx9e7RXut+7h+2y3Jtc9VRn6nePjP3zcyd1cPVp3c4HgAA3DK2Xa7tE9WXqnfMzIsz84G11jer\nR6rPVl+pnlprPbvb6Y+rzW4/AgAAW9m0y9TdWWvHacHnZGbWtW93X+7o6MmefvryTR8TAAAX08z0\nWlMp1lo3nE+xlzffAQDARSOMAQAgYQwAANXBw/g4D98BALAfm3Z5+O4s6xifg+PDnh4AgNvY0cnX\nY1vtbSoFAAAkjAEAoBLGAABQHTyMj/PwHQAA+7HJw3cAAODhOwAA2J0wBgCAhDEAAFTCGAAAqoOH\n8XFWpQAAYD82WZUCAACsSgEAALsTxgAAkDAGAIBKGAMAQCWMAQCgEsYAAFAdPIyPs44xAAD7sck6\nxgAAYB1jAADYnTAGAICEMQAAVMIYAAAqYQwAAJUwBgCAShgDAEB18DA+zgs+AADYj01e8AEAAF7w\nAQAAuxPGAACQMAYAgEoYAwBAJYwBAKASxgAAUAljAACohDEAAFTCGAAAqoOH8XFeCQ0AwH5s8kpo\nAADwSmgAANidMAYAgIQxAABUwhgAACphDAAAlTAGAIBKGAMAQCWMAQCgEsYAAFAJYwAAqIQxAABU\nwhgAACphDAAAlTAGAIBKGAMAQHXwMD6uNocdAgAAt6lNV3pzO5f2NYztHB/29AAA3MaOTr4e22pv\nUykAACBhDAAAlTAGAIBKGAMAQCWMAQCgEsYAAFAJYwAAqIQxAABUwhgAACphDAAAlTAGAIBKGAMA\nQCWMAQCgEsYAAFAJYwAAqIQxAABUwhgAACphDAAAlTAGAIBKGAMAQCWMAQCgqkv7OvDMPFS9r/rO\n6ufWWr+8r3MBAMBZ7e2O8VrrF9daP1F9qPqb+zoPh7HZbA49BM7A9bu4XLuLzfW7uFy7N4atw3hm\nHp+ZV2fmy1dtf3BmnpuZ52fm0Wv86M9UP3vWgXJr8R+Ii831u7hcu4vN9bu4XLs3hl3uGD9RPXB6\nw8zcUX3sZPu7qvfPzDtP/fk/rj6z1vqtcxgrAADszdZhvNb6YvU/rtp8f/W1tdYLa63/Uz1VPVQ1\nM49U763+xsz8xDmNFwAA9mLWWtvvPHNf9Utrre87+fzXqwdO5hI3Mz9W3b/W+vAWx9r+xAAAcAZr\nrbnRPmddleJaJ9gqeLcZHAAA3CxnXZXi5ereU5/vqV454zEBAOCm2zWMp2+/S/xM9faZuW9m7qwe\nrj59XoMDAICbZZfl2j5Rfal6x8y8ODMfWGt9s3qk+mz1leqptdazWxzrRku8cYu63rJ93Ppm5p6Z\n+dzMfHVmfntmbvgsALeOmXnTzPzazPzmyfX76KHHxG5m5o6Z+c8z4wbSBTMz/21m/svJv3//6dDj\nYTcz85aZ+bcz8+zMfGVm3n3dfXd5+O48nCzx9nxXVqx4pSt3nR9eaz13UwfC6zIzf776/erj33oI\nk4thZt5WvW2t9Vsz88eq36ge8u/exTEzd621vjEzf6T6lerDay1/SV8QM/OR6geq71pr/eihx8P2\nZuZ3qh9Ya129OhcXwMw8WX1+rfXEzFyq7lpr/a9r7bu3N9+9husu8cat7zrL9nEBrLW+/q01xdda\nv189W9192FGxi7XWN06+fVNXHp62us8FMTP3VD9S/ctDj4XXZTpMM3FGM/Od1V9Yaz1Rtdb6v9eL\n4jrMRb67eunU55fzlzPcVDPzPdWfqX7tsCNhFye/iv/N6uvVL6+1njn0mNjaP6n+Yf5n5qJa1X+c\nmWdm5oOHHgw7+d7q92bmiZOpTP9iZr7jejsfIoxf9xJvwNmdTKP4ZPWTJ3eOuSDWWv9vrfVnu7IC\n0Ltn5k8fekzc2My8r3r15Dc2Vz/EzsXwnrXWn+vKXf+/fzKtkIvhUvX91c+utb6/+kb1U9fb+RBh\nbIk3OJCTuVWfrP71WusXDz0eXp+TXwNuqgcPPBS284PVj57MU/031Q/PzMcPPCZ2sNb6+sk//3v1\nqa5MC+VieLl6aa316yefP9mVUL6mQ4SxJd4uPnc8Lq6fq7661vqnhx4Iu5mZPz4zbzn5/juqv1R5\ncPICWGv99Frr3rXW93bl77zPrbV+/NDjYjszc9fJb9qamTdXf7n6r4cdFdtaa71avTQz7zjZ9N7q\nq9fb/6xvvtvZWuubM/MPurLE2x3V49ss8cat4WTZvqPqrTPzYvXRb01o59Y2Mz9Y/Z3qt0/mqa7q\np9da/+GwI2NLf7L6Vycr+9xR/cJa6zMHHhO8EfyJ6lMzs7rSTT+/1vrsgcfEbj5c/fzM/NHqd6oP\nXG/Hm75cGwAA3IosPQIAAAljAACohDEAAFTCGAAAKmEMAACVMAYAgEoYAwBAJYwBAKCq/w9AZJiD\nZ/rg2gAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fc140b0ead0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdAAAAK+CAYAAAASS0+3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X2UbHV95/v3BxCXDxFBBTIgoID4MGPQGQ2JmbF9iICT\nO3gzg9EkV1GTuGYkmpjJ+DB3FjA3d5Ym8fqwvA7LDDrojSLiRDDjyMFgT65GCFGOEnk6oAJH5CQK\nxPFhvHD6e//ofQ59muqq6n12V+2qer/WqnV27dq961fVp+tT32/99q5UFZIkaXMOmPYAJEmaRQao\nJEktGKCSJLVggEqS1IIBKklSCwaoJEktGKCSpLmW5KFJrk5ybZLrkpwzYJuDk1yUZEeSLyY5ZtR+\ntyRAk5yW5MYkNyd501bchyRJ46iqHwPPq6pnACcDpyd59rrNXgPcXVUnAu8Cfn/UfjsP0CQHAO8F\nTgWeBrw8yZO7vh9JksZVVT9sFh8KHASsP4vQGcCFzfIlwAtG7XMrKtBnAzuq6raqug+4qBmYJElT\nkeSAJNcCdwFXVNU16zY5CrgDoKp2A/cmOWzYPg/agnHuHURjJ6uhuo8knkNQklqoqkx7DIMc9/iH\n1G0775/2MHZV1ZHrV1bVCvCMJI8CPpnkqVV1/ZpN1j+n4cFV6j62IkAH/WIHDuKccx74HHdpaYml\npaUtGM7+Offcczn33HOnPYyRHGe3HGe3ZmGcfR3j8vIyy8vLe6+fd9550xvMCLftvJ/d3z5xqmM4\n8Cd3HDHs9qr6XpJl4DRgbYDeATweuDPJgcCjquqeYfvaigDdCaydvXQ0cOegDfv4n1WS+mR9cdHn\nAO2rJI8F7quqv0vyMOCFwNvWbfYp4JXA1cCZwJWj9rsVAXoNcEKSY4FvAy8DXr4F9yNJ6pkVVqY9\nhEF+EriwmeR6APCxqvp0kvOAa6rqT4ELgA8n2QF8l9XsGqrzAK2q3UnOBrY1A72gqm7o+n4mpY9t\n5UEcZ7ccZ7dmYZyzMEa1U1XXAc8csP6cNcs/Bl66mf1mWt8HmqT8LlJJ2pwkvZ1ElKTu+/bxUx3D\nQ37y1ok9P1vRwpUkLajd1csW7pbwVH6SJLVggEqS1IItXElSZ1aGn3tgrliBSpLUghWoJKkzPT0O\ndEtYgUqS1IIBKklSC7ZwJUmd2b1AJ8ixApUkqQUrUElSZzyMRZIkDWWASpLUgi1cSVJndtvClSRJ\nwxigkiS1YAtXktQZZ+FKkqShrEAlSZ3xTESSJGkoA1SSpBZs4UqSOrM43wZqBSpJUitWoJKkzngm\nIkmSNJQBKklSC7ZwJUmd2b04HVwrUEmS2jBAJUlqwRauJKkzHgcqSZKGsgKVJHVmN5n2ECbGClSS\npBYMUEmSWrCFK0nqzIrHgUqSpGEMUEmSWrCFK0nqjLNwJUnSUFagkqTOWIFKkqShDFBJklqwhStJ\n6sxK2cKVJElDWIFKkjrjJCJJkjSUASpJUgu2cCVJndm9QHXZ4jxSSZI6ZIBKktSCLVxJUmc8DlSS\nJA1lBSpJ6ozHgUqSpKEMUEmSWrCFK0nqzO5anLpscR6pJEkdsgKVJHVmpYd1WZKjgQ8BRwK7gT+q\nqves2+ZfA78CFPAQ4CnAY6vq3o32a4BKkubd/cAbq2p7kkcCX0qyrapu3LNBVf0h8IcASX4B+K1h\n4Qm2cCVJc66q7qqq7c3y94EbgKOG/MjLgY+O2q8VqCSpM30/DjTJccDJwNUb3P4w4DTgdaP2ZYBK\nkmbW9qt+yFeu+uFY2zbt20uANzSV6CD/C/D5Ue1bgFTV2APtUpKa1n1L0qxKQlU/TzibpP7sG0+a\n6hhe8ISbBz4/SQ4C/hT4b1X17o1+Psl/AS6uqotG3ZcVqCSpMz0+DvQDwPUjwvMQ4LmszsYdyQCV\nJM21JM9hNRSvS3Itq4eqvBU4Fqiqen+z6UuAy6vqR+Ps1wCVJHVmpYeTiKrqC8CBY2x3IXDhuPvt\nba0tSVKfGaCSJLVgC1eS1JndC1SXLc4jlSSpQwaoJEkt2MKVJHWmx8eBdm5xHqkkSR2yApUkdaaP\n3we6VUY+0iQXJNmV5Ktr1h2aZFuSm5Jc3pz+aM9t70myI8n2JCdv1cAlSZqmcd4qfBA4dd26NwOf\nraqTgCuBtwAkOR04vqpOBF4LnN/hWCVJ6o2RLdyq+nySY9etPoPVE+7C6mmPPsdqqJ4BfKj5uauT\nHJLkiKra1eGYJUk9tbufXxSzJdo2qw/fE4pVdRdweLP+KOCONdt9i+Hf+i1J0kzqehLRoLceG37p\n57nnnrt3eWlpiaWlpY6HI0mzbXl5meXl5WkPY2yLdCaisb5Qu2nhfqqqnt5cvwFYqqpdSY4EPldV\nT0lyfrP8sWa7G4HnDmrh+oXakrR5ff9C7YtveeZUx/DSE748sedn3LcKYd/q8jLgrGb5LODSNetf\nAZDkFOBeP/+UJM2jkS3cJB8BloDHJLkdOAd4G/DxJK8GbgfOBKiqTyd5cZJbgB8Ar9qqgUuS+mdl\ngc5ENM4s3F/e4KYXbrD92fs1IkmSZsDivFWQJKlDnspPktSZRZqFuziPVJKkDlmBSpI645mIJEnS\nUAaoJEkt2MKVJHXG7wOVJElDWYFKkjqze4HORLQ4j1SSpA4ZoJIktWALV5LUmZWBXws9n6xAJUlq\nwQCVJKkFW7iSpM44C1eSJA1lBSpJ6oxfZyZJkoYyQCVJasEWriSpMyt+H6gkSRrGAJUkqQVbuJKk\nzjgLV5IkDWUFKknqzIpnIpIkScMYoJIktWALV5LUmd1+H6gkSRrGClSS1BknEUmSpKEMUEmSWrCF\nK0nqjJOIJEnSUAaoJEkt2MKVJHXGWbiSJGkoK1BJUmd2W4FKkqRhDFBJklqwhStJ6syKx4FKkjQf\nkhyd5Mok1ye5Lsnrh2z7rCT3J/nFUfu1ApUkdaank4juB95YVduTPBL4UpJtVXXj2o2SHAC8DfjM\nODvt5SOVJKkrVXVXVW1vlr8P3AAcNWDT3wQuAf5mnP0aoJKkhZHkOOBk4Op16/8e8BLgfBjvg1xb\nuJKkzqzUZCcR3XbN33LbX31nrG2b9u0lwBuaSnStdwFvqqpKAmOEqAEqSZpZxz7rcRz7rMftvf75\n828cuF2Sg1gNzw9X1aUDNvlHwEVZTc/HAqcnua+qLtvovg1QSdIi+ABwfVW9e9CNVfXEPctJPgh8\nalh4ggEqSerQ7h5OrUnyHOBXgOuSXAsU8FbgWKCq6v3rfqTG2a8BKkmaa1X1BeDATWz/6nG2M0Al\nSZ2Z9CSiaepfrS1J0gwwQCVJasEWriSpMysLVJctziOVJKlDBqgkSS3YwpUkdWa3s3AlSdIwVqCS\npM54HKgkSRrKAJUkqQVbuJKkzqzU4tRli/NIJUnqkBWoJKkzu3ESkSRJGsIAlSSpBVu4kqTOeByo\nJEkaygCVJKkFW7iSpM54HKgkSRrKClSS1JkVjwOVJEnDGKCSJLVgC1eS1JndHgf6gCRHJ7kyyfVJ\nrkvy+mb9oUm2JbkpyeVJDlnzM+9JsiPJ9iQnb+UDkCRpGsZp4d4PvLGqngr8DPC6JE8G3gx8tqpO\nAq4E3gKQ5HTg+Ko6EXgtcP6WjFyS1DsrdcBUL5M08t6q6q6q2t4sfx+4ATgaOAO4sNnswuY6zb8f\nara/GjgkyREdj1uSpKnaVFwnOQ44GbgKOKKqdsFqyAKHN5sdBdyx5se+1ayTJGlujD2JKMkjgUuA\nN1TV95PURpsOWDdw2yfmqXuXD+VxHJbDB20mTVQStu2+eNrDkABYXl5meXl52sMY2yKdTH6sAE1y\nEKvh+eGqurRZvSvJEVW1K8mRwN8063cCj1/z40cDdw7a7/F5WrtRT8m23RfzogNfOu1haItVbfTe\nUJq8paUllpaW9l4/77zzpjcY7WPcFu4HgOur6t1r1l0GnNUsnwVcumb9KwCSnALcu6fVO+sMT0nS\nHiMr0CTPAX4FuC7Jtay2Y98KvB24OMmrgduBMwGq6tNJXpzkFuAHwKu2avBdsbKUpG4s0qn8RgZo\nVX0BOHCDm1+4wc+cvT+DmjTDc9WsvZGYtfFKmi+eyk97jRtGfZlgY3hK/bNSmeplkgxQbZrBJUkG\n6MzqSxU4C3yuJG0FA3RGWQWOb9BzZahKW8NT+UlTsjbYtjLkfAMiaX8ZoAuqjxXY+lm1hpykPjNA\n1+ljsGyFPoZT12NalN+l1CfOwl1gfQwWtZMszgHdkibPAJ1Bfays+jimnz/gzGkPQVo4K2Sql0ky\nQGdQH6vkPo5JkraSASpJUgsGqBZeH9vP0qxyEtGC8wV1sdh+ltSGATpA2xdUg1fSorMCVStWMpK0\nOAxQTdw0KnW7A5K6trAB6gvq9EyjUt/MfV6x8vEtHIk032zhTtkkws12qzbi/w1J4+hlgPoCpmmq\nqmkPQdIMOGjaA5AkzY9Jt1GnqZcVqNQ1P/OW1LVeBKgvbtpqfiwgTYYnk5+wrX5xM6AlSV3rRYBu\nNauP/TPuGxDfqEhaJAsRoBrfoBAc9w3IVrxRMZSl2dLH40CTHJ3kyiTXJ7kuyesHbHNSkr9I8j+T\nvHGcx2qAah99q9YnMZ71Ie2JFKS5cz/wxqp6KvAzwOuSPHndNt8FfhP4g3F3aoBq4a0P6Z8/4Mwp\njUSafX2sQKvqrqra3ix/H7gBOGrdNt+pqi+xGrZjMUAlSQsjyXHAycDV+7svT6TQU9t2X9y7dqok\n9c3d2+/gnu07x9o2ySOBS4A3NJXofjFAe2qc8JynkL1i5eNU1dw8HmlRTfpMRI/+qWN49E8ds/f6\nNy4cXFgmOYjV8PxwVV3axX0boDNsnsLGzx0lbbEPANdX1bvH2HasdwF+Btozkzxsw0NEJC2CJM8B\nfgV4fpJrk3w5yWlJXpvkN5ptjkhyB/DbwL9NcnvT8t2QFWjPTLKqnKcKVlI/9PFk8lX1BeDAEdvs\nAh6/mf1agar3rJQl9ZEVqHrPSlmaHdXDCnSrLFwFOo/VzDw+Jknqu4UL0LXVzLwEjxWaJE3ewgUo\nPBCcBo8kdcvvA51zBqckaX8tZIBKMD8tfEnTsfABOuxF1BfY+WYnQupeH7+NZassfIAOexH1BVaS\ntJGFD9B5MWvV8qyNV9J4qjLVyyQtdIDO04v4rFXLszZeSVpvoQPUF3FJUlsLHaCavHmq+iU9mJOI\npC1i1S9pXhigmggrT2kxOIlIvTEvwWPlKWneGKA9tyjBMy9vFCQtDgN0AfUxrBbljYI075xEpLlm\nWEnS/jNAZ1gfK8l5cMXKx6c9BEkzwACdYVtdSS5qQFuhS+1VTfcySQaoNrSoQbKobxwkbY4BKq2z\nqG8cpC6skKleJskA7ZjVy+yrSfeBJM0kA7RjVi+StBgOmvYApL5xFq7U3qRPpzdNVqALZhIt5llv\nY9tFkDQOA3TBTCIcph1A+xvgfgYqteeZiKQZNu0Al7QYDFCNNMst2Vkeu6R+M0C110ZhM8sV3SyP\nXZpFnolIM69N5WXYSNL4DNAWNhNO02oh9jEMbadKmie9CtBZeYHdTDj1McjGsRW/i1l9LiSNrypT\nvUxSrwLUF9j+GOd3MStveIYZ9BiSxTkQXFJ7nolIrc3DG55Bj8HjQKX2PBPRgpuHymp/LPrjl6Rx\nGKADzENltT8W/fFL0jgM0Ba6rNCs9sbj8yTNBk/lp6G6rNCs9sbj8ySpbxYyQK1m+mMrfxf+niVt\npYUMUKuZ/tjK38X6fRuo0tbzVH5rJHlokquTXJvkuiTnNOuPS3JVkpuSfDTJQc36g5NclGRHki8m\nOWarH4Q0Dt84SerSyACtqh8Dz6uqZwAnA6cn+Wng7cA7quok4F7gNc2PvAa4u6pOBN4F/P6WjFyS\n1DueiWidqvphs/hQVk++UMDzgE806y8EXtIsn9FcB7gEeEEnI5UkqUfGCtAkByS5FrgLuAK4Fbi3\nqlaaTXYCRzXLRwF3AFTVbuDeJId1OmpNjJ8bStJgY53KrwnKZyR5FPAnwFMGbdb8u76Gzprb9nFr\nfW3v8qE8jsNy+DjD2VLbdl/sZ2VrzNpz4e9P82Z5eZnl5eVpD2Nsi3Qqv02dC7eqvpfkvwOnAI9O\nckATrkcDdzab7QQeD9yZ5EDgUVV1z6D9HZ+ntR/5FvHFt99GBaS/P82bpaUllpaW9l4/77zzpjcY\n7WOcWbiPTXJIs/ww4IXA9cDngDObzV4JXNosX9Zcp7n9yi4HrP6ZZJvXgJT6raZ8maRxPgP9SeBz\nSbYDVwOXV9WngTcDb0xyM3AYcEGz/QXAY5PsAH6r2W7q/Cxv67QNNX8nkmbZyBZuVV0HPHPA+m8A\nPz1g/Y+B3pUJVi6jTeLzw7X30dffyRUrH5/2ECTNgIU8E5EG259AG7ea3HMffa4++xrs0izwOFBp\nkzYbOn0OKb9QW9I4DNAp63MlNot8PiVNigE6ZX2uxKZhfwPQ51OasgWahmuA9oBV0wMMQEmzYlMn\nUtDWMDQkzYtFOhORFWhPWZVKUr/1PkAXNUisSiWpO0kuSLIryVc3uP1RSS5Lsr357uuzRu2z9wE6\nT0Ey728G5v3xSRqtarqXIT4InDrk9tcBX6uqk1n9us53JBn6MWfvA3SezNObgUHm/fFJml1V9Xlg\n4Beb7NkE+Ilm+SeA71bV/cP2aYCOYFUlSeOb4TMRvRd4apI7ga8Abxj1A87CHcGqavL8Tk9J4/rR\n177Oj67/Rhe7OhW4tqqen+R44IokT6+q72/0A1agC2ptZd23KtvwlDSuhz3tiRx25gv2XvbDq4D/\nAlBVtwLfAJ487AcM0AW1NqTmMbD69qZAWhiV6V6GS3MZ5DZWv++aJEcATwK+PmxntnA1l+bxTYGk\n9pJ8BFgCHpPkduAc4GCgqur9wO8B/3nNYS7/pqruHrZPA1SSNPeq6pdH3P5thh/m8iAGqCSpM4v0\nbYB+BrqBK1Y+7udokqQNGaAb+PkDzuTUg35p4UJ00R6vpI75dWYCqKqFm4wyycdrWEuaZTMRoL7Q\nzqdFe3Miab7MRIBO64XW4JakzZnhU/lt2kwE6LRYIUmSNmKAaiircEkazACdQ12GnlW4pE1xFq5m\nmaEnSVvPMxFJkjoz6Yk802QFqoXl57uS9ocBqoVlq1vS/jBANZZxqzWrOmnBOYlI2te41ZpVnaRF\nYYBOgVVav12x8vFpD0GaYZnyZXJmMkD7HkCjxmeV1m/+fiSNYyYDtO8vcH0fn4arRfpGYEmtzWSA\nzpq+V8zaV7I4x7FJnXMS0fybZKhZkc4WK1BJ41jIAN22+2JDTZK0XxYyQCcZnrZvt57PsdQjtnBn\nUx9fSNeH9RUrH+/lOGeZ3QRJ0zBXJ5OfhRfSnz/gzGkPQZK2jieT11ayAp0+fweS9tdMBei8vOhN\nq1Kel+dvrbaPaRa6FZL6rdcBuv7Fcate9OYxWAYZ9/mbpefDIJT6pWq6l0nqdYBO6sXRF+F9+UZF\nkkbrdYCqvT6GlW9UpAXgYSwaVx+DCgwrSdpqBuh+MqgkaTEZoDOurxXwHn0fn6SOVaZ7mSADdMb1\nvQLu+/gkqS0DVJKkFgzQFmxLzje/D1RqLzXdyyQZoC3Ylpxvfh+opHEYoOqVzVT3dgKkHlqg40Cn\n+m0sn/nW9i3dfxWcfvTJW3ofamejLzXfTHVvJ0DSNE01QP2oaXEZfpJm3UK0cG31aTOcRCTtB48D\nnS9tqh1Dd3E5iUjSOBYiQNtYG7qGqSRpPQN0DF18XmcIS1oICzQL1wCdkP0NYQNYkvrFAO3YVgWd\ns1a74RsRaYtZgaqtSQddXwOhr+Pa6PfT1/FK6i8DdMb1tTIdNa6+BVZfn0dJ/bUQAdq3F2sZWNLc\nsoU7X3yxnh2+2ZE0KxYiQDU7tuLNjqEsTZBnIpLa6WNY2YGQtBUMUHVqVsKqj0EvabYYoHPMkNjY\nrAS9NGtS071MkgE6xwwJSdo6CxmgVmaSpP21kAE6L5WZbwQk9Y7HgWoWzNrZfiRpnhigc2weKu0u\n3wT4hkJSl8YO0CQHJPlyksua68cluSrJTUk+muSgZv3BSS5KsiPJF5Mcs1WD31/bdl/si2rPdfkm\nYB7eUEhqJ8kFSXYl+eoGtz83yb1Nzn05yf8+ap+bqUDfAFy/5vrbgXdU1UnAvcBrmvWvAe6uqhOB\ndwG/v4n7mKgXHfhSvz1FkhbDB4FTR2zz51X1zObye6N2OFaAJjkaeDHwn9asfj7wiWb5QuAlzfIZ\nzXWAS4AXjHMfi8IqqHu+KZH6o6/HgVbV54F7Rg1/M4913Ar0ncDv0sxxSvIY4J6qWmlu3wkc1Swf\nBdzRDHg3cG+SwzYzKGkzfFMiqSOnJLk2yX9N8tRRGx80aoMk/xTYVVXbkyztWc2Dk7rW3LbPLthg\ncvF5f/jdvcvP/dmHsfSzDx81nIWxbffFvOjAl+79t83PdrVdF/fb9X21vb9R49iKcUr7Y3l5meXl\n5WkPY3wTPqH7j3bcwv+85dYudvUl4Niq+mGS04FPAk8a9gOpGn7gTJL/APwqcD/wMOAnmh2/CDiy\nqlaSnAKcU1WnJ/lMs3x1kgOBb1fV4QP2W6PuW5K0ryRUTTilxpSknvCud0x1DN/4rd/Z8PlJcizw\nqap6+qj9JPkG8A+r6u6NthnZwq2qt1bVMVX1ROBlwJVV9avA54Azm81eCVzaLF/WXKe5/cpR9yFJ\n0gQM6p6u3pAcsWb52awWmBuGJ4zRwh3izcBFSf4P4Frggmb9BcCHk+wAvstq6EqSFkFPG4tJPgIs\nAY9JcjtwDnAwUFX1fuBfJPmXwH3Aj4BfGrnPabVRbeFK0ub1voX7zim3cH974xZu1zwTkSRJLexP\nC1eSpH0tUGPRClSSpBasQCVJnRl2NqB5YwUqSVILBqgkSS3YwpUkdccWriRJGsYAlSSpBVu4kqTu\n2MKVJEnDWIFKkjrjcaCSJGkoA1SSpBZs4UqSutPPb1rbElagkiS1YAUqSeqOk4gkSdIwBqgkSS3Y\nwpUkdcbjQCVJ0lAGqCRJLdjClSR1xxauJEkaxgpUktQZJxFJkqShDFBJklqwhStJ6o4tXEmSNIwB\nKklSC7ZwJUndsYUrSZKGsQKVJHXG40AlSdJQBqgkSS0YoJIktWCASpLUgpOIJEndcRKRJEkaxgCV\nJKkFW7iSpM54HKgkSRrKAJUkqQVbuJKk7tjClSRJw1iBSpK6YwUqSZKGMUAlSWrBFq4kqTMeBypJ\nkoayApUkdccKVJIkDWOASpLUgi1cSVJnnEQkSZKGMkAlSWrBFq4kqTu2cCVJ0jAGqCSpOzXlywaS\nXJBkV5KvbnD7Lyf5SpLtST6f5B+MeqgGqCRpEXwQOHXI7V8H/klVnQz8HvBHo3boZ6CSpLlXVZ9P\ncuyQ269ac/Uq4KhR+zRAJUmdmZPjQH8N+G+jNjJAJUkz6we33cIPbr+ls/0leR7wKuDnRm1rgEqS\nZtYjjj2BRxx7wt7r3/nCttb7SvJ04P3AaVV1z6jtDVBJUnf63cJNc3nwDckxwCeA/62qbh1nZwao\nJGnuJfkIsAQ8JsntwDnAwUBV1fuBfwccBrwvSYD7qurZw/ZpgEqSutPTCrSqfnnE7b8O/Ppm9ulx\noJIktWCASpLUgi1cSVJn5uQ40LFYgUqS1IIVqCSpO1agkiRpGANUkqQWbOFKkjrjJKJ1knyz+aLR\na5P8ZbPu0CTbktyU5PIkh6zZ/j1JdjRfTHryVg1ekqRpGbeFuwIsVdUz1pza6M3AZ6vqJOBK4C0A\nSU4Hjq+qE4HXAud3PGZJkqZu3ADNgG3PAC5sli9sru9Z/yGAqroaOCTJEfs5TknSLKgpXyZo3AAt\n4PIk1yT5tWbdEVW1C6Cq7gIOb9YfBdyx5me/xRjf7C1J0iwZdxLRz1bVXUkeB2xLchMbZ/2gr4oZ\nuO255567d3lpaYmlpaUxhyNJi2F5eZnl5eVpD2N8CzSJKFWbe7RJzgG+D/waq5+L7kpyJPC5qnpK\nkvOb5Y81298IPHdPtbpmP7XZ+5akRZeEqhr4nZbTlqT+/m//X1Mdw1+/840Te35GtnCTPDzJI5vl\nRwAvAq4DLgPOajY7C7i0Wb4MeEWz/SnAvevDU5KkWTdOC/cI4E+SVLP9H1fVtiR/BVyc5NXA7cCZ\nAFX16SQvTnIL8APgVVs0dklSz/SyNN4iIwO0qr4BPOhYzqq6G3jhBj9z9v4PTZKk/vJMRJKk7izQ\n1BbPhStJUgsGqCRJLdjClSR1xpPJS5KkoQxQSZJasIUrSeqOLVxJkjSMFagkqTtWoJIkaRgDVJKk\nFmzhSpI643GgkiRpKANUkqQWbOFKkrpjC1eSJA1jBSpJ6oyTiCRJ0lAGqCRJLdjClSR1xxauJEka\nxgpUktQZJxFJkqShDFBJklqwhStJ6o4tXEmSNIwBKklSC7ZwJUndsYUrSZKGsQKVJHXG40AlSdJQ\nBqgkSS3YwpUkdccWriRJGsYKVJLUmdTilKBWoJIktWCASpLUgi1cSVJ3FqeDawUqSVIbBqgkSS3Y\nwpUkdcZT+UmSNEeSnJbkxiQ3J3nTgNuPSfLZJF9JcmWSvzdqnwaoJKk7NeXLAEkOAN4LnAo8DXh5\nkiev2+wPgf9cVT8F/HvgbaMeqgEqSZp3zwZ2VNVtVXUfcBFwxrptngpcCVBVywNufxADVJI0744C\n7lhzfWezbq3twD8HSPKLwCOTHDpsp04ikiR1ZtKTiL531y38j7tuHbVZBqxbP9LfBd6b5Czgz4Fv\nAfcP26kBKkmaWY868gQedeQJe6/f+ZVtgzbbCRyz5vrRwJ1rN6iqb/NABfoI4J9X1f8Ydt+2cCVJ\n8+4a4IQkxyY5GHgZcNnaDZI8JsmeSvUtwAdG7dQAlSR1p4ezcKtqN3A2sA34GnBRVd2Q5Lwkv9Bs\ntgTclORG4HDg/xz1UG3hSpLmXlV9Bjhp3bpz1ix/AvjEZvZpgEqSOuOZiCRJ0lAGqCRJLdjClSR1\nxxauJEkaxgpUktQZJxFJkqShDFBJklqwhStJ6k4tTg/XClSSpBYMUEmSWrCFK0nqjLNwJUnSUFag\nkqTuWIHKFp1IAAAO0UlEQVRKkqRhDFBJklqwhStJ6kxWpj2CybEClSSpBStQSVJ3nEQkSZKGMUAl\nSWrBFq4kqTOeiUiSJA1lgEqS1IItXElSd/w+0H0lOSTJx5PckORrSX46yaFJtiW5KcnlSQ5Zs/17\nkuxIsj3JyVs3fEmSpmPcFu67gU9X1VOAnwJuBN4MfLaqTgKuBN4CkOR04PiqOhF4LXB+56OWJPVS\narqXSRoZoEl+AvjHVfVBgKq6v6r+DjgDuLDZ7MLmOs2/H2q2vRo4JMkRXQ9ckqRpGqcCfSLwnSQf\nTPLlJO9P8nDgiKraBVBVdwGHN9sfBdyx5ue/1ayTJGlujDOJ6CDgmcDrquqvkryT1fbtRsVyBqwb\nuO255567d3lpaYmlpaUxhiNJi2N5eZnl5eVpD2N8izOHiNSIGVNN+/WLVfXE5vrPsRqgxwNLVbUr\nyZHA56rqKUnOb5Y/1mx/I/DcPdXqmv3WqPuWJO0rCVU1qFCZuiT1nP/1D6Y6hi/8ye9O7PkZ2cJt\ngu+OJE9qVr0A+BpwGXBWs+4s4NJm+TLgFQBJTgHuXR+ekiTNunGPA3098MdJHgJ8HXgVcCBwcZJX\nA7cDZwJU1aeTvDjJLcAPmm0lSQtgkU7lN1aAVtVXgGcNuOmFG2x/9v4MSpKkvvNMRJKk7izQ3BbP\nhStJUgsGqCRJLdjClSR1ZpEmEVmBSpLUghWoJKk7VqCSJGkYA1SSpBZs4UqSOuMkIkmSNJQBKklS\nC7ZwJUndWVmcHq4VqCRJLViBSpK6szgFqBWoJEltGKCSJLVgC1eS1BmPA5UkSUNZgUqSulOLU4Ja\ngUqS1IIBKklSC7ZwJUmdcRKRJEkaygCVJKkFA1SS1J2a8mUDSU5LcmOSm5O8aYNtXprka0muS/L/\njHqofgYqSZprSQ4A3gu8ALgTuCbJpVV145ptTgDeBPxMVX0vyWNH7dcAlSR1Jv08DvTZwI6qug0g\nyUXAGcCNa7b5deD/rqrvAVTVd0bt1BauJGneHQXcseb6zmbdWk8CTkry+SR/keTUUTu1ApUkzax7\n7v46997z9VGbZcC69aXyQcAJwD8BjgH+3yRP21ORDmKASpK6szLZuzv00U/k0Ec/ce/1b37jzwZt\ntpPVUNzjaFY/C12/zReragX4ZpKbgBOBL21037ZwJUnz7hrghCTHJjkYeBlw2bptPgk8H6CZQHQi\nMLS0NUAlSXOtqnYDZwPbgK8BF1XVDUnOS/ILzTaXA99N8jXgz4B/XVX3DNuvLVxJUmd6OguXqvoM\ncNK6deesu/47wO+Mu08rUEmSWrAClSR1p58F6JawApUkqQUDVJKkFmzhSpK609NJRFvBClSSpBas\nQCVJncniFKBWoJIktWGASpLUgi1cSVJ3nEQkSZKGMUAlSWrBFq4kqTOZ8PeBTpMVqCRJLViBSpK6\n4yQiSZI0jAEqSVILtnAlSd1ZnA6uFagkSW1YgUqSOhMnEUmSpGEMUEmSWrCFK0nqji1cSZI0jAEq\nSVILtnAlSd3xZPKSJGkYK1BJUmc8DlSSJA1lgEqS1IItXElSd2zhSpKkYQxQSZJasIUrSeqOLVxJ\nkjSMFagkqTueiUiSJA1jgEqS1IItXElSZzyVnyRJGsoKVJLUHStQSZI0zMgATfKkJNcm+XLz798l\neX2SQ5NsS3JTksuTHLLmZ96TZEeS7UlO3tqHIEnS5I0M0Kq6uaqeUVXPBP4h8APgT4A3A5+tqpOA\nK4G3ACQ5HTi+qk4EXgucv1WDlyT1TNV0LxO02RbuC4Fbq+oO4Azgwmb9hc11mn8/BFBVVwOHJDmi\ng7FKktQbmw3QXwI+0iwfUVW7AKrqLuDwZv1RwB1rfuZbzTpJkubG2LNwkzwE+GfAm5pVG9XKGbBu\n4Lbnnnvu3uWlpSWWlpbGHY4kLYTl5WWWl5enPYzxLdAs3NSYDzbJPwP+VVWd1ly/AViqql1JjgQ+\nV1VPSXJ+s/yxZrsbgefuqVbX7K/GvW9J0qokVNWgQmXqktSpf/+tUx3D5X/9Hyb2/Gymhfty4KNr\nrl8GnNUsnwVcumb9KwCSnALcuz48JUlzamXKlwkaq4Wb5GGsTiD6jTWr3w5cnOTVwO3AmQBV9ekk\nL05yC6szdl/V7ZAlSZq+sQK0qn4EPG7durtZDdVB25+9/0OTJKm/PJWfJKkznkxekiQNZQUqSeqO\nFagkSRrGAJUkqQVbuJKk7qzYwpUkSUMYoJKkuZfktCQ3Jrk5yZsG3P7aJF9tvvf6z5M8edQ+DVBJ\nUnd6+H2gSQ4A3gucCjwNePmAgPzjqnp6VT0D+APgnaMeqgEqSZp3zwZ2VNVtVXUfcBEPfIc1AFX1\n/TVXH8kYZ9Z1EpEkqTv9PA50/fdU72Q1VPeR5F8BbwQeAjx/1E4NUEnSzPruD2/n7h/dMWqzsb6n\nuqreB7wvycuAf8cD3zg2kAEqSZpZj3n4MTzm4cfsvX7rPX8xaLOdwDFrrh8N3Dlktx8Dzh91334G\nKknqTg8nEQHXACckOTbJwcDLWP3u6r2SnLDm6i8AN496qFagkqS5VlW7k5wNbGO1cLygqm5Ich5w\nTVX9KXB2khcC/x9wD/DKUfs1QCVJc6+qPgOctG7dOWuWf2uz+zRAJUnd8VR+kiRpGCtQSVJ3auT5\nB+aGFagkSS0YoJIktWALV5LUnX6eym9LWIFKktSCFagkqTsexiJJkoYxQCVJasEWriSpO04ikiRJ\nwxigkiS1YAtXktQdW7iSJGkYK1BJUnesQCVJ0jAGqCRJLdjClSR1Z8XvA5UkSUNYgUqSuuMkIkmS\nNIwBKklSC7ZwJUndsYUrSZKGMUAlSWrBFq4kqTsrtnAlSdIQVqCSpM5UeSYiSZI0hAEqSVILtnAl\nSd1xEpEkSRrGAJUkqQVbuJKk7ngqP0mSNIwVqCSpOyseBypJkoYwQCVJasEWriSpO04ikiRJw1iB\nSpI6U04ikiRJwxigkiS1YAtXktQdJxFJkqRhDFBJklqwhStJ6o7fBypJkoaxApUkdac8DlSSJA1h\ngEqS1IItXElSZ8pJRJIkaRgrUElSd5xEJEmShjFAJUlqwRauJKkzTiLSXsvLy9MewlgcZ7ccZ7dm\nYZyzMEa1l+S0JDcmuTnJmwbcfnCSi5LsSPLFJMeM2qcBOsKs/FE5zm45zm7NwjhnYYxqJ8kBwHuB\nU4GnAS9P8uR1m70GuLuqTgTeBfz+qP0aoJKk7tTKdC+DPRvYUVW3VdV9wEXAGeu2OQO4sFm+BHjB\nqIdqgEqS5t1RwB1rru9s1g3cpqp2A/cmOWzYTlNT+vbwJIvzSbMkdaiqMu0xDJLkm8CxUx7Grqo6\ncu2KJP8CeFFV/UZz/VeBZ1XVG9Zs89fNNnc2129ptrlnozua2izcvv4HkCS1U1XHTXsMG9gJrJ0U\ndDRw57pt7gAeD9yZ5EDgUcPCE2zhSpLm3zXACUmOTXIw8DLgsnXbfAp4ZbN8JnDlqJ16HKgkaa5V\n1e4kZwPbWC0cL6iqG5KcB1xTVX8KXAB8OMkO4LushuxQU/sMVJKkWTaVFu6oA1qnKck3k3wlybVJ\n/rJZd2iSbUluSnJ5kkOmMK4LkuxK8tU16zYcV5L3NAcEb09y8pTHeU6SnUm+3FxOW3PbW5px3pDk\nRRMa49FJrkxyfZLrkry+Wd+r53PAOH+zWd+35/OhSa5u/mauS3JOs/64JFc1z+dHkxzUrN/0Aetb\nPM4PJvl6s/7LSZ6+5mem8nekGVFVE72wGtq3sDpT6yHAduDJkx7HkPF9HTh03bq3A/+mWX4T8LYp\njOvngJOBr44aF3A68F+b5Z8GrpryOM8B3jhg26cA17L6UcJxzf+LTGCMRwInN8uPBG4Cnty353PI\nOHv1fDb3/fDm3wOBq5rn6WPAmc36/wi8tln+l8D7muVfAi6a4P/PQeP8IPCLA7ad2t+Rl9m4TKMC\nHeeA1mkKD67M1x5geyHwkomOCKiqzwPrZ4StH9cZa9Z/qPm5q4FDkhwxxXHC6vO63hmsvnjeX1Xf\nBHaw+v9jS1XVXVW1vVn+PnADq7PyevV8bjDOPceu9eb5bMb3w2bxoawGeAHPAz7RrF/7d7PpA9a7\nMmCce4683+j5nMrfkWbDNAJ0nANap6mAy5Nck+TXmnVHVNUuWH1RAx43tdHt6/B14zq8Wb/+Of4W\n03+OX9e0wf7Tmtbo1MeZ5DhWK+arePDvuTfP55pxXt2s6tXzmeSAJNcCdwFXALcC91btPTXM2r/z\nTR+wvlXjrKprmpt+r3k+35HkIevH2ejD35F6ZBoBOuidXp9mMv1sVf0j4MWsvkj9Y/o1vnH07Tl+\nH3B8VZ3M6gvXO5r1Ux1nkkeyWgG9oanwNrrvvo2zd89nVa1U1TNYreSfzWo7eaOxrB9nmNI4kzwV\neHNVPQV4FvAYVtv3g8bJpMap2TCNAB3ngNapaSoPqupvgU+y+mKwa0/rJsmRwN9Mb4T72GhcO1k9\nIHiPqT7HVfW3VbXnheePeKCtOLVxNhNaLgE+XFWXNqt793wOGmcfn889qup7wH8HTgEendWTeK8f\ny95xZswD1rdwnKet6Trcx+rnob15PtVv0wjQcQ5onYokD2/e7ZPkEcCLgOtYHd9ZzWavBC4duIOt\nF/Z9V7x2XGfxwLguA14BkOQUVltpuyYzRGDdOJsw2uMXgb9uli8DXtbMynwCcALwlxMa4weA66vq\n3WvW9fH5fNA4+/Z8JnnsnjZykocBLwSuBz7H6gHpsO/fzWVs8oD1LRznjXuezyRh9XPatc/nNP+O\n1HfTmLkEnMbqjMIdrLZPpj6bqhnXE1idFXwtq8H55mb9YcBnmzFfATx6CmP7CKvvfn8M3A68Cjh0\no3Gx+tU9twBfAZ455XF+CPhq89x+ktXPGvds/5ZmnDeweh7KSYzxOcDuNb/rLzf/Jzf8PU/j+Rwy\nzr49n/+gGdv2Zlz/tln/BFY/s72Z1Rm5D2nWPxS4uPn7vwo4bsrj/LPm9/rV5rl9+DR/715m5+KJ\nFCRJasFz4UqS1IIBKklSCwaoJEktGKCSJLVggEqS1IIBKklSCwaoJEkt/P8zmzkABCh11wAAAABJ\nRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fc140bbf450>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import atlaas\n",
    "import numpy as np\n",
    "from atlaas.helpers.gdal2 import gdal2\n",
    "from atlaas.helpers.matplot import hist, show\n",
    "\n",
    "## clean up previous run\n",
    "import glob\n",
    "import subprocess\n",
    "subprocess.check_output(['rm'] + glob.glob('atlaas.*x*.tif'))\n",
    "####\n",
    "\n",
    "test = atlaas.Atlaas()\n",
    "test.init(120.0, 120.0, 0.1, 0, 0, 0, 31, True)\n",
    "\n",
    "transformation = np.identity(4, dtype=np.double)\n",
    "\n",
    "## cloud1 is a \"flat\" scan (all Z = 0)\n",
    "cloud1 = np.random.rand(200000, 4).astype('float32')\n",
    "cloud1[:,2] = 0 # all Z = 0\n",
    "cloud1 *= 20 # x,y = random [0..20]\n",
    "# merge it once\n",
    "test.merge(cloud1, transformation)\n",
    "####\n",
    "\n",
    "## cloud_obj is a \"dynamic\" object (all Z = 3m)\n",
    "cloud_obj = np.random.rand(20000, 4).astype('float32')\n",
    "cloud_obj[:,2] = 2\n",
    "cloud_obj *= 3\n",
    "# cloud1 + cloud_obj\n",
    "cloud2 = np.append(cloud1, cloud_obj, axis=0)\n",
    "# merge both\n",
    "test.merge(cloud2, transformation)\n",
    "####\n",
    "\n",
    "filepath = 'out2.tif'\n",
    "test.region('region.png')\n",
    "atlaas.merge('atlaas.*x*.tif', filepath)\n",
    "g = gdal2( filepath )\n",
    "band = g.dtm()\n",
    "hist(band)\n",
    "show(band, 0, 3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsYAAAE6CAYAAAAV0JRKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEedJREFUeJzt3WGoZGd5B/D/s00NprWWplBLlsRKkNBAsVoiaNURS5MY\nNNQKzaYiCNYirUo/SKwIO9tPgVpqMdZSXDe1mMYquGprWgV3lFTQtdXGqksCtnGX6LZSSytCsfbp\nhzums9d7s+fs3NnZufn9YMic986c88JJsv9973Oet7o7AADweHdg3RMAAIBLgWAMAAARjAEAIIlg\nDAAASQRjAABIIhgDAEASwRgAAJIIxgAAkGRFwbiqXlBVn6qqd1bV81dxDQAA2EurWjHuJP+V5PIk\nZ1Z0DQAA2DODgnFVHa2qs1X1wLbxm6rqVFU9WFV3fH+8uz/V3bckeVOS39vbKQMAwN4bumJ8LMmN\niwNVdSDJXfPx65Mcqqrrtn3vP5I8YdlJAgDAql025EPdfX9VXbNt+IYkD3X3w0lSVfcmuTXJqar6\nlWwF5idnKzwDAMAlbVAw3sVVSU4vHJ/JVlhOd38wyQcf68tV1UtcGwAABuvuOt9nlgnGO518VNjt\nlo031XQ6zXQ6Xfc0uEDu3+Zy7zab+7e53LvNVnXeTJxkua4UZ5JcvXB8MMkjY04wnU4zm82WmAIA\nAOxsNpuN+gvNmGBcOXeV+GSSa6vqmqp6QpLbknx4xPkynU4zmUzGfAUAAAaZTCZ7H4yr6p4kn07y\n9Kr6WlW9qru/l+R1ST6W5EtJ7u3ur4yfMpvIX2g2m/u3udy7zeb+bS737vGh1lXnW1V9+PDhTCYT\n/7IBALDnZrNZZrNZjhw5Mujhu7UGYw/fAQCwalU1KBivaktoAADYKGsNxrpSAACwKmO7UiilAABg\nX1NKAQAAIwjGAAAQNcYAAOxTaowBAGCBGmMAABhBMAYAgAjGAACQxMN3AADsUx6+AwCABR6+AwCA\nEQRjAABIctk6L/7GN/7ujuMvecmL8/znP+8izwYAgMeztQbjt771s0meNn9935fyuc/9QU6cEIwB\nALhws9lsVKOHtT58l+x07eOZTO7OiRPHL/qcAADYfzx8BwAAIwjGAAAQwRgAAJIIxgAAkEQwBgCA\nJGtu15ZMk0zmLwAA2DvatQEAwALt2gAAYATBGAAAIhgDAEASwRgAAJIIxgAAkEQwBgCAJIIxAAAk\nEYwBACCJne8AANin7HwHAAAL7HwHAAAjCMYAABDBGAAAkgjGAACQRDAGAIAkgjEAACQRjAEAIIlg\nDAAASQRjAABIssJgXFVXVNXnqurFq7oGAADslVWuGN+R5H0rPD8AAOyZQcG4qo5W1dmqemDb+E1V\ndaqqHqyqOxbGX5Tky0n+Ncl596UGAIB1G7pifCzJjYsDVXUgyV3z8euTHKqq6+Y/fmGSZye5Pcmr\n92aqAACwOpcN+VB3319V12wbviHJQ939cJJU1b1Jbk1yqrvfMh97ZZJv7uF8AQBgJQYF411cleT0\nwvGZbIXlR3X3e5Y4PwAAXDTLBOOdaod73CmmC+8n8xcAAFy42WyW2Ww2+nvLBOMzSa5eOD6Y5JFx\np5gucXkAAPhBk8kkk8nk0eMjR44M+t6YYFw5d5X4ZJJr57XHX09yW5JDI86XrWA8iZViAAD22tiV\n4+o+f/VDVd2TrfR6ZZKzSQ5397GqujnJ27LV3eJod985+MJVvXPlxfFMJnfnxInjQ08FAAC7qqp0\n93lbCA/tSnH7LuP3Jblv5NwAAOCSs0yN8R6YRikFAACrsJJSilVQSgEAwMUwtJRi6M53AACwryml\nAABgX1JKAQAAC5RSAADACIIxAABEjTEAAPuUGmMAAFigxhgAAEYQjAEAIGqMAQDYp9QYAwDAAjXG\nAAAwgmAMAAARjAEAIIlgDAAASXSlAABgn9KVAgAAFuhKAQAAIwjGAAAQwRgAAJIIxgAAkEQwBgCA\nJNq1AQCwT2nXBgAAC7RrAwCAEQRjAACIYAwAAEkEYwAASCIYAwBAEsEYAACSCMYAAJBEMAYAgCR2\nvgMAYJ+y8x0AACyw8x0AAIwgGAMAQARjAABIIhgDAEASwRgAAJIIxgAAkEQwBgCAJIIxAAAkEYwB\nACDJiraErqrrkrwhyZVJPtHdf7KK6wAAwF5ZyYpxd5/q7tcm+bUkz1nFNQAAYC8NCsZVdbSqzlbV\nA9vGb6qqU1X1YFXdse1nL0nyV0k+unfTBQCA1Ri6YnwsyY2LA1V1IMld8/Hrkxyal1AkSbr7I919\nS5JX7NFcAQBgZQbVGHf3/VV1zbbhG5I81N0PJ0lV3Zvk1iSnquoFSV6W5PIkf72H8wUAgJVY5uG7\nq5KcXjg+k62wnO7+ZJJPnv8U04X3k/kLAAAu3Gw2y2w2G/29ZYJx7TDW404xXeLyAADwgyaTSSaT\nyaPHR44cGfS9ZbpSnEly9cLxwSSPLHE+AABYmzErxpVzV4lPJrl2Xnv89SS3JTk07vLTKKEAAGAV\nxpZUVPf5qx+q6p5spdcrk5xNcri7j1XVzUnelq2V56PdfefgC1f1zpUXxzOZ3J0TJ44PPRUAAOyq\nqtLdO5UBn2NoV4rbdxm/L8l9I+cGAACXnJVsCT3cNEopAABYhZWUUqyCUgoAAC6GoaUUy3SlAACA\nfUMpBQAA+5JSCgAAWKCUAgAARhCMAQAgaowBANin1BgDAMACNcYAADCCYAwAAFFjDADAPqXGGAAA\nFqgxBgCAEQRjAACIYAwAAEkEYwAASKIrBQAA+5SuFAAAsEBXCgAAGEEwBgCACMYAAJBEMAYAgCSC\nMQAAJNGuDQCAfUq7NgAAWKBdGwAAjCAYAwBABGMAAEgiGAMAQBLBGAAAkgjGAACQRDAGAIAkNvgA\nAGCfssEHAAAssMEHAACMIBgDAEAEYwAASCIYAwBAEsEYAACSCMYAAJBEMAYAgCSCMQAAJBGMAQAg\niWAMAABJkstWdeKqujXJLUmelOTd3f3xVV0LAACWtbJg3N0fSvKhqvrxJL+fRDAGAOCSNbiUoqqO\nVtXZqnpg2/hNVXWqqh6sqjt2+Opbkrxj2YkCAMAqjakxPpbkxsWBqjqQ5K75+PVJDlXVdQs/vzPJ\nR7v7C3swVwAAWJnBwbi770/yrW3DNyR5qLsf7u7vJrk3ya1JUlWvS/KiJC+vqtfs0XwBAGAllq0x\nvirJ6YXjM9kKy+nutyd5+2N/fbrwfjJ/AQDAhZvNZpnNZqO/t2wwrh3GevjXp0teHgAAzjWZTDKZ\nTB49PnLkyKDvLdvH+EySqxeODyZ5ZMlzAgDARTd2xbhy7irxySTXVtU1Sb6e5LYkh4afbholFAAA\nrMLYkorqHlb5UFX3ZCvBXpnkbJLD3X2sqm5O8rZsrT4f7e47B56vd666OJ7J5O6cOHF80LwAAOCx\nVFW6e6cS4HMMXjHu7tt3Gb8vyX0j5gYAAJecle18N8w0SikAAFiFlZVS7DWlFAAAXAxDSymW7UoB\nAAD7glIKAAD2JaUUAACwQCkFAACMIBgDAEDUGAMAsE+pMQYAgAVqjAEAYATBGAAAosYYAIB9So0x\nAAAsUGMMAAAjCMYAABDBGAAAknj4DgCAfcrDdwAAsMDDdwAAMIJgDAAAEYwBACCJYAwAAEkEYwAA\nSKJdGwAA+5R2bQAAsEC7NgAAGEEwBgCACMYAAJBEMAYAgCSCMQAAJBGMAQAgiWAMAABJBGMAAEhi\n5zsAAPYpO98BAMACO98BAMAIgjEAAEQwBgCAJIIxAAAkEYwBACCJYAwAAEkEYwAASCIYAwBAEsEY\nAACSCMYAAJBkRcG4qn6mqt5VVX+5ivMDAMBeW0kw7u5/7u5Xr+LcAACwCoOCcVUdraqzVfXAtvGb\nqupUVT1YVXesZooAALB6Q1eMjyW5cXGgqg4kuWs+fn2SQ1V13bbv1dIzBACAi2BQMO7u+5N8a9vw\nDUke6u6Hu/u7Se5NcmuSVNVPVNU7kzzDSjIAAJvgsiW+e1WS0wvHZ7IVltPd/57ktec/xXTh/WT+\nAgCACzebzTKbzUZ/b5lgvFOZRI87xXSJywMAwA+aTCaZTCaPHh85cmTQ95bpSnEmydULxweTPLLE\n+QAAYG3GrBhXzl0lPpnk2qq6JsnXk9yW5NC4y0+jhAIAgFUYW1JR3eevfqiqe7KVXq9McjbJ4e4+\nVlU3J3lbtlaej3b3nYMvXNU7V14cz2Ryd06cOD70VAAAsKuqSneft1vaoBXj7r59l/H7ktw3cm4L\nprFiDADAKqxkxXgVrBgDAHAxDF0xXsmW0AAAsGkEYwAAyNqD8TTJbL1TAABg4z3lKU9NVe34GmqZ\nDT72wHS9lwcAYF84e/bh7L7X3LBwrJQCAACy9mA8jVIKAABWY5YxFQpKKQAA2Kcm89eRQZ9WSgEA\nABGMAQAgiWAMAABJ1h6Mp/HwHQAAqzGLh+8AAMDDdwAAMJ5gDAAAEYwBACDJ2oPxNB6+AwBgNWbx\n8B0AAHj4DgAAxhOMAQAggjEAACQRjAEAIIlgDAAASdYejKfRrg0AgNWYRbs2AADQrg0AAMYTjAEA\nIIIxAAAkEYwBACCJYAwAAEkEYwAASCIYAwBAEsEYAACSrD0YT2PnOwAAVmMWO98BAICd7wAAYDzB\nGAAAIhgDAEASwRgAAJIIxgAAkEQwBgCAJIIxAAAkEYwBACCJYAwAAEkEYwAASLKiLaGr6ookf5zk\nv5N8srvvWcV1AABgr6xqxfhlSd7f3b+Z5KUrugZrNJvN1j0FluD+bS73brO5f5vLvXt8GBSMq+po\nVZ2tqge2jd9UVaeq6sGqumPhRweTnJ6//94ezZVLiP9BbDb3b3O5d5vN/dtc7t3jw9AV42NJblwc\nqKoDSe6aj1+f5FBVXTf/8elsheMkqT2YJwAArNSgYNzd9yf51rbhG5I81N0Pd/d3k9yb5Nb5zz6Y\n5OVV9Y4kH9mryQIAwKpUdw/7YNU1ST7S3T83P/7VJDd292vmx69IckN3v37g+YZdGAAAltTd561i\nWKYrxU4nHxx2h0wOAAAulmW6UpxJcvXC8cEkjyw3HQAAWI8xwbhy7irxySTXVtU1VfWEJLcl+fBe\nTg4AAC6Woe3a7kny6SRPr6qvVdWruvt7SV6X5GNJvpTk3u7+ysDz7dbmjUvcbq37uPRV1cGq+kRV\nfbmqvlhVg54H4NJQVZdX1Weq6vPz+3d43XNinKo6UFX/UFUWkTZMVf1LVf3j/L+/z657PoxTVU+u\nqvdX1Veq6ktV9exdPzv04bu9Mm/z9mCSF2Wr9OJkktu6+9RFnQgXpKp+Mcm3k7zn+w9ishmq6ilJ\nntLdX6iqH03y90lu9d/e5qiqK7r7O1X1Q0n+Lsnru9sf0huiqn4nybOS/Fh32/xqg1TVV5M8q7u3\nd+hiA1TV3dnaiflYVV2W5Iru/s+dPruqne8ey2O1eeMSt0vrPjZAd3+ju78wf//tJF9JctV6Z8UY\n3f2d+dvLs/XwtO4+G6KqDiZ5cZJ3rXsuXJDKejITS6qqJyV5XncfS5Lu/p/dQnGynpt8Vf5/V7xk\n6yE+fzjDRVRVT03yjCSfWe9MGGP+q/jPJ/lGko9398l1z4nB/jDJG+MvM5uqk/xtVZ2sqt9Y92QY\n5WlJvllVx+alTH9aVU/c7cPrCMZLtXkDljMvo/hAkjfMV47ZEN39v93989nqAvTsqvrZdc+J86uq\nW5Kcnf/GZvuD7GyG53T3L2Rr1f+35mWFbIbLkjwzyTu6+5lJvpPkTbt9eB3BWJs3WJN5bdUHkvx5\nd39o3fPhwsx/DThLctOap8Iwz03y0nmd6l8keWFVvWfNc2KE7v7G/J//lq3dfW9Y74wY4UyS0939\nufnxB7IVlHe0jmCszdvms+Kxud6d5Mvd/UfrngjjVNVPVtWT5++fmOSXknhwcgN095u7++ruflq2\n/sz7RHe/ct3zYpiqumL+m7ZU1Y8k+eUk/7TeWTFUd59Ncrqqnj4felGSL+/2+WV2vrsg3f29qvrt\nbLV5O5Dk6NA2b6zfvHXfJMmVVfW1JIe/X9DOpa2qnpvk15N8cV6n2kne3N1/s96ZMdBPJ/mzeWef\nA0ne190fXfOc4PHgp5J8sKo6W7npvd39sTXPiXFen+S9VfXDSb6a5FW7ffCit2sDAIBLkdYjAAAQ\nwRgAAJIIxgAAkEQwBgCAJIIxAAAkEYwBACCJYAwAAEkEYwAASJL8H/VtZv6r+nMlAAAAAElFTkSu\nQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fc140c424d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdAAAAK+CAYAAAASS0+3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X20bHV95/n3BxCXDxFBBdIgoID40K1ot4a06fb4EAEn\n0zjpxmiSUdQkrm6JJqbTPvT04jKT6aVJHB+W47BMo41OFBE7gmkjF4OnMxohRLlKeLygAlfkJgrE\n9mEcuec7f5x9L4dDnao6++46tavq/Vqr1t21a59dv6pzT33q+63f3pWqQpIkbc4B0x6AJEmzyACV\nJKkFA1SSpBYMUEmSWjBAJUlqwQCVJKkFA1SSNNeSPDTJVUmuSXJtknMGbHNwkguT7EzypSTHjNrv\nRAI0yWlJbkxyc5I3T+I+JEkaR1X9GHh+VT0TOBk4Pclz1m32WuDuqjoReDfw+6P223mAJjkAeB9w\nKvA04BVJntz1/UiSNK6q+mGz+FDgIGD9WYTOAC5oli8GXjhqn5OoQJ8D7Kyq26rqJ8CFzcAkSZqK\nJAckuQa4C7i8qq5et8lRwB0AVbUHuDfJYcP2edAExrlvEI1drIbqAyTxHIKS1EJVZdpjGOS4xz+k\nbtt137SHsbuqjly/sqpWgGcmeRTwqSRPrarr12yy/jkND65SH2ASATroFztwEOecc//nuEtLSywt\nLU1gOPtn27ZtbNu2bdrDGMlxdstxdmsWxtnXMS4vL7O8vLzv+rnnnju9wYxw26772PPtE6c6hgN/\neucRw26vqu8lWQZOA9YG6B3A44E7kxwIPKqq7hm2r0kE6C5g7eylo4E7B23Yx/+sktQn64uLPgdo\nXyV5LPCTqvr7JA8DXgS8fd1mnwZeBVwFnAlcMWq/kwjQq4ETkhwLfBt4OfCKCdyPJKlnVliZ9hAG\n+WnggmaS6wHAx6vqM0nOBa6uqj8Fzgc+kmQn8F1Ws2uozgO0qvYkORvY3gz0/Kq6oev72Sp9bCsP\n4ji75Ti7NQvjnIUxqp2quhZ41oD156xZ/jHwss3sN9P6PtAk5XeRStLmJOntJKIk9ZNvHz/VMTzk\np2/dsudnEi1cSdKC2lO9bOFOhKfykySpBQNUkqQWbOFKkjqzMvzcA3PFClSSpBasQCVJnenpcaAT\nYQUqSVILBqgkSS3YwpUkdWbPAp0gxwpUkqQWrEAlSZ3xMBZJkjSUASpJUgu2cCVJndljC1eSJA1j\ngEqS1IItXElSZ5yFK0mShrIClSR1xjMRSZKkoQxQSZJasIUrSerM4nwbqBWoJEmtWIFKkjrjmYgk\nSdJQBqgkSS3YwpUkdWbP4nRwrUAlSWrDAJUkqQVbuJKkzngcqCRJGsoKVJLUmT1k2kPYMlagkiS1\nYIBKktSCLVxJUmdWPA5UkiQNY4BKktSCLVxJUmechStJkoayApUkdcYKVJIkDWWASpLUgi1cSVJn\nVsoWriRJGsIKVJLUGScRSZKkoQxQSZJasIUrSerMngWqyxbnkUqS1CEDVJKkFmzhSpI643GgkiRp\nKCtQSVJnPA5UkiQNZYBKktSCLVxJUmf21OLUZYvzSCVJ6pAVqCSpMys9rMuSHA18GDgS2AP8UVW9\nd902/xb4FaCAhwBPAR5bVfdutF8DVJI07+4D3lRVO5I8Evhyku1VdePeDarqD4E/BEjyC8BvDQtP\nsIUrSZpzVXVXVe1olr8P3AAcNeRHXgF8bNR+rUAlSZ3p+3GgSY4DTgau2uD2hwGnAa8ftS8DVJI0\ns3Zc+UO+euUPx9q2ad9eDLyxqUQH+R+BL4xq3wKkqsYeaJeS1LTuW5JmVRKq+nnC2ST159940lTH\n8MIn3Dzw+UlyEPCnwJ9V1Xs2+vkk/wW4qKouHHVfVqCSpM70+DjQDwLXjwjPQ4DnsTobdyQDVJI0\n15I8l9VQvDbJNaweqvI24FigquoDzaYvBS6rqh+Ns18DVJLUmZUeTiKqqi8CB46x3QXABePut7e1\ntiRJfWaASpLUgi1cSVJn9ixQXbY4j1SSpA4ZoJIktWALV5LUmR4fB9q5xXmkkiR1yApUktSZPn4f\n6KSMfKRJzk+yO8nX1qw7NMn2JDcluaw5/dHe296bZGeSHUlOntTAJUmapnHeKnwIOHXdurcAn6uq\nk4ArgLcCJDkdOL6qTgReB5zX4VglSeqNkS3cqvpCkmPXrT6D1RPuwuppjz7PaqieAXy4+bmrkhyS\n5Iiq2t3hmCVJPbWnn18UMxFtm9WH7w3FqroLOLxZfxRwx5rtvsXwb/2WJGkmdT2JaNBbjw2/9HPb\ntm37lpeWllhaWup4OJI025aXl1leXp72MMa2SGciGusLtZsW7qer6unN9RuAparaneRI4PNV9ZQk\n5zXLH2+2uxF43qAWrl+oLUmb1/cv1L7olmdNdQwvO+ErW/b8jPtWITywurwUOKtZPgu4ZM36VwIk\nOQW4188/JUnzaGQLN8lHgSXgMUluB84B3g58IslrgNuBMwGq6jNJXpLkFuAHwKsnNXBJUv+sLNCZ\niMaZhfvLG9z0og22P3u/RiRJ0gxYnLcKkiR1yFP5SZI6s0izcBfnkUqS1CErUElSZzwTkSRJGsoA\nlSSpBVu4kqTO+H2gkiRpKCtQSVJn9izQmYgW55FKktQhA1SSpBZs4UqSOrMy8Guh55MVqCRJLRig\nkiS1YAtXktQZZ+FKkqShrEAlSZ3x68wkSdJQBqgkSS3YwpUkdWbF7wOVJEnDGKCSJLVgC1eS1Bln\n4UqSpKGsQCVJnVnxTESSJGkYA1SSpBZs4UqSOrPH7wOVJEnDWIFKkjrjJCJJkjSUASpJUgu2cCVJ\nnXESkSRJGsoAlSSpBVu4kqTOOAtXkiQNZQUqSerMHitQSZI0jAEqSVILtnAlSZ1Z8ThQSZLmQ5Kj\nk1yR5Pok1yZ5w5Btn53kviS/OGq/VqCSpM70dBLRfcCbqmpHkkcCX06yvapuXLtRkgOAtwOfHWen\nvXykkiR1paruqqodzfL3gRuAowZs+pvAxcDfjrNfA1SStDCSHAecDFy1bv0/AF4KnAfjfZBrC1eS\n1JmV2tpJRLdd/Xfc9tffGWvbpn17MfDGphJd693Am6uqksAYIWqASpJm1rHPfhzHPvtx+65/4bwb\nB26X5CBWw/MjVXXJgE3+CXBhVtPzscDpSX5SVZdudN8GqCRpEXwQuL6q3jPoxqp64t7lJB8CPj0s\nPMEAlSR1aE8Pp9YkeS7wK8C1Sa4BCngbcCxQVfWBdT9S4+zXAJUkzbWq+iJw4Ca2f8042xmgkqTO\nbPUkomnqX60tSdIMMEAlSWrBFq4kqTMrC1SXLc4jlSSpQwaoJEkt2MKVJHVmj7NwJUnSMFagkqTO\neByoJEkaygCVJKkFW7iSpM6s1OLUZYvzSCVJ6pAVqCSpM3twEpEkSRrCAJUkqQVbuJKkzngcqCRJ\nGsoAlSSpBVu4kqTOeByoJEkaygpUktSZFY8DlSRJwxigkiS1YAtXktSZPR4Her8kRye5Isn1Sa5N\n8oZm/aFJtie5KcllSQ5Z8zPvTbIzyY4kJ0/yAUiSNA3jtHDvA95UVU8FfhZ4fZInA28BPldVJwFX\nAG8FSHI6cHxVnQi8DjhvIiOXJPXOSh0w1ctWGnlvVXVXVe1olr8P3AAcDZwBXNBsdkFznebfDzfb\nXwUckuSIjsctSdJUbSqukxwHnAxcCRxRVbthNWSBw5vNjgLuWPNj32rWSZI0N8aeRJTkkcDFwBur\n6vtJaqNNB6wbuO0T89R9y4fyOA7L4YM2k7ZUErbvuWjaw5AAWF5eZnl5edrDGNsinUx+rABNchCr\n4fmRqrqkWb07yRFVtTvJkcDfNut3AY9f8+NHA3cO2u/xeVq7UU/J9j0X8eIDXzbtYWjCqjZ6byht\nvaWlJZaWlvZdP/fcc6c3GD3AuC3cDwLXV9V71qy7FDirWT4LuGTN+lcCJDkFuHdvq3fWGZ6SpL1G\nVqBJngv8CnBtkmtYbce+DXgHcFGS1wC3A2cCVNVnkrwkyS3AD4BXT2rwXbGylKRuLNKp/EYGaFV9\nEThwg5tftMHPnL0/g9pqhueqWXsjMWvjlTRfPJWf9hk3jPoywcbwlPpnpTLVy1YyQLVpBpckGaAz\nqy9V4CzwuZI0CQbojLIKHN+g58pQlSbDU/lJU7I22CYZcr4BkbS/DNAF1ccKbP2sWkNOUp8ZoOv0\nMVgmoY/h1PWYFuV3KfWJs3AXWB+DRe0ki3NAt6StZ4DOoD5WVn0c088fcOa0hyAtnBUy1ctWMkBn\nUB+r5D6OSZImyQCVJKkFA1QLr4/tZ2lWOYlowfmCulhsP0tqwwAdoO0LqsEradFZgaoVKxlJWhwG\nqLbcNCp1uwOSurawAeoL6vRMo1LfzH1evvKJCY5Emm+2cKdsK8LNdqs24v8NSePoZYD6AqZpqqpp\nD0HSDDho2gOQJM2PrW6jTlMvK1Cpa37mLalrvQhQX9w0aX4sIG0NTya/xSb94mZAS5K61osAnTSr\nj/0z7hsQ36hIWiQLEaAa36AQHPcNyCTeqBjK0mzp43GgSY5OckWS65Ncm+QNA7Y5KclfJvl/k7xp\nnMdqgOoB+latb8V41oe0J1KQ5s59wJuq6qnAzwKvT/Lkddt8F/hN4A/G3akBqoW3PqR//oAzpzQS\nafb1sQKtqruqakez/H3gBuCoddt8p6q+zGrYjsUAlSQtjCTHAScDV+3vvjyRQk9t33NR79qpktQ3\nd++4g3t27Bpr2ySPBC4G3thUovvFAO2pccJznkL28pVPUFVz83ikRbXVZyJ69DOO4dHPOGbf9W9c\nMLiwTHIQq+H5kaq6pIv7NkBn2DyFjZ87SpqwDwLXV9V7xth2rHcBfgbaM1t52IaHiEhaBEmeC/wK\n8IIk1yT5SpLTkrwuyW802xyR5A7gt4F/n+T2puW7ISvQntnKqnKeKlhJ/dDHk8lX1ReBA0dssxt4\n/Gb2awWq3rNSltRHVqDqPStlaXZUDyvQSVm4CnQeq5l5fEyS1HcLF6Brq5l5CR4rNEnaegsXoHB/\ncBo8ktQtvw90zhmckqT9tZABKsH8tPAlTcfCB+iwF1FfYOebnQipe338NpZJWfgAHfYi6gusJGkj\nCx+g82LWquVZG6+k8VRlqpettNABOk8v4rNWLc/aeCVpvYUOUF/EJUltLXSAauvNU9Uv6cGcRCRN\niFW/pHlhgGpLWHlKi8FJROqNeQkeK09J88YA7blFCZ55eaMgaXEYoAuoj2G1KG8UpHnnJCLNNcNK\nkvafATrD+lhJzoPLVz4x7SFImgEG6AybdCW5qAFthS61VzXdy1YyQLWhRQ2SRX3jIGlzDFBpnUV9\n4yB1YYVM9bKVDNCOWb3MvtrqPpCkmWSAdszqRZIWw0HTHoDUN87Cldrb6tPpTZMV6ILZihbzrLex\n7SJIGocBumC2IhymHUD7G+B+Biq155mIpBk27QCXtBgMUI00yy3ZWR67pH4zQLXPRmEzyxXdLI9d\nmkWeiUgzr03lZdhI0vgM0BY2E07TaiH2MQxtp0qaJ70K0Fl5gd1MOPUxyMYxid/FrD4XksZXlale\ntlKvAtQX2P4Y53cxK294hhn0GJLFORBcUnueiUitzcMbnkGPweNApfY8E9GCm4fKan8s+uOXpHEY\noAPMQ2W1Pxb98UvSOAzQFrqs0Kz2xuPzJM0GT+Wnobqs0Kz2xuPzJKlvFjJArWb6Y5K/C3/PkiZp\nIQPUaqY/Jvm7WL9vA1WaPE/lt0aShya5Ksk1Sa5Nck6z/rgkVya5KcnHkhzUrD84yYVJdib5UpJj\nJv0gpHH4xklSl0YGaFX9GHh+VT0TOBk4PcnPAO8A3llVJwH3Aq9tfuS1wN1VdSLwbuD3JzJySVLv\neCaidarqh83iQ1k9+UIBzwc+2ay/AHhps3xGcx3gYuCFnYxUkqQeGStAkxyQ5BrgLuBy4Fbg3qpa\naTbZBRzVLB8F3AFQVXuAe5Mc1umotWX83FCSBhvrVH5NUD4zyaOAPwGeMmiz5t/1NXTW3PYAt9Z1\n+5YP5XEclsPHGc5Ebd9zkZ+VrTFrz4W/P82b5eVllpeXpz2MsS3Sqfw2dS7cqvpekv8GnAI8OskB\nTbgeDdzZbLYLeDxwZ5IDgUdV1T2D9nd8ntZ+5BPii2+/jQpIf3+aN0tLSywtLe27fu65505vMHqA\ncWbhPjbJIc3yw4AXAdcDnwfObDZ7FXBJs3xpc53m9iu6HLD6ZyvbvAak1G815ctWGucz0J8GPp9k\nB3AVcFlVfQZ4C/CmJDcDhwHnN9ufDzw2yU7gt5rtps7P8ianbaj5O5E0y0a2cKvqWuBZA9Z/A/iZ\nAet/DPSuTLByGW0rPj9cex99/Z1cvvKJaQ9B0gxYyDMRabD9CbRxq8m999Hn6rOvwS7NAo8DlTZp\ns6HT55DyC7UljcMAnbI+V2KzyOdT0lYxQKesz5XYNOxvAPp8SlO2QNNwDdAesGq6nwEoaVZs6kQK\nmgxDQ9K8WKQzEVmB9pRVqST1W+8DdFGDxKpUkrqT5Pwku5N8bYPbH5Xk0iQ7mu++PmvUPnsfoPMU\nJPP+ZmDeH5+k0aqmexniQ8CpQ25/PXBdVZ3M6td1vjPJ0I85ex+g82Se3gwMMu+PT9LsqqovAAO/\n2GTvJsBPNcs/BXy3qu4btk8DdASrKkka3wyfieh9wFOT3Al8FXjjqB9wFu4IVlVbz+/0lDSuH133\ndX50/Te62NWpwDVV9YIkxwOXJ3l6VX1/ox+wAl1QayvrvlXZhqekcT3saU/ksDNfuO+yH14N/BeA\nqroV+Abw5GE/YIAuqLUhNY+B1bc3BdLCqEz3MlyayyC3sfp91yQ5AngS8PVhO7OFq7k0j28KJLWX\n5KPAEvCYJLcD5wAHA1VVHwB+D/jPaw5z+XdVdfewfRqgkqS5V1W/POL2bzP8MJcHMUAlSZ1ZpG8D\n9DPQDVy+8gk/R5MkbcgA3cDPH3Ampx70SwsXoov2eCV1zK8zE0BVLdxklK18vIa1pFk2EwHqC+18\nWrQ3J5Lmy0wE6LReaA1uSdqcGT6V36bNRIBOixWSJGkjBqiGsgqXpMEM0DnUZehZhUvaFGfhapYZ\nepI0eZ6JSJLUma2eyDNNVqBaWH6+K2l/GKBaWLa6Je0PA1RjGbdas6qTFpyTiKQHGrdas6qTtCgM\n0CmwSuu3y1c+Me0hSDMsU75snZkM0L4H0KjxWaX1m78fSeOYyQDt+wtc38en4WqRvhFYUmszGaCz\npu8Vsx4oWZzj2KTOOYlo/m1lqFmRzhYrUEnjWMgA3b7nIkNNkrRfFjJAtzI8bd9Ons+x1CO2cGdT\nH19I14f15Suf6OU4Z5ndBEnTMFcnk5+FF9KfP+DMaQ9BkibHk8lrkqxAp8/fgaT9NVMBOi8vetOq\nlOfl+Vur7WOahW6FpH7rdYCuf3Gc1IvePAbLIOM+f7P0fBiEUr9UTfeylXodoFv14uiL8AP5RkWS\nRut1gKq9PoaVb1SkBeBhLBpXH4MKDCtJmjQDdD8ZVJK0mAzQGdfXCnivvo9PUscq071sIQN0xvW9\nAu77+CSpLQNUkqQWDNAWbEvON78PVGovNd3LVjJAW7AtOd/8PlBJ4zBA1Subqe7tBEg95HGg0mRt\nFH6bqe7tBEiaJgNUU2H4SZp1CxGgtvq0GU4ikvaDx4HOg+LPdu0A2lU7hu7ichKRpHHMcYCG049+\nRuufXhu6hqkkab05DlCAbsr5Lj6vM4QlLQRn4apr+xvCBrAk9YsB2rFJBZ2zVrvhGxFpwqxA1dZW\nB11fA6Gv49ro99PX8UrqLwN0xvW1Mh01rr4FVl+fR0n9tRAB2rcXaxlY0tyyhTtffLGeHb7ZkTQr\nehKgHriuVZN4s2MoS1vIMxFtNU+dNi/6GFZ2ICRNQk8CdH/cf8o+Td+shFUfg17SbJmDAN2/U/bN\nM0NiY7MS9NKsSU33spXmIEDBFvBghoQkTc6cBOjmWJlJkvbXQgbovFRmvhGQ1DseB6pZMGtn+5Gk\neTIDAeos27bmodLu8k2AbygkdWnsAE1yQJKvJLm0uX5ckiuT3JTkY0kOatYfnOTCJDuTfCnJMfs3\nxMnNst2+5yJfVHuuyzcB8/CGQlI7Sc5PsjvJ1za4/XlJ7m1y7itJ/pdR+9xMBfpG4Po1198BvLOq\nTgLuBV7brH8tcHdVnQi8G/j9TdzHBiYzy/bFB77Mb0+RpMXwIeDUEdv8RVU9q7n83qgdjhWgSY4G\nXgL8pzWrXwB8slm+AHhps3xGcx3gYuCF49zHorAK6p5vSqT+6OtxoFX1BeCeUcPfzGMdtwJ9F/C7\nNHOckjwGuKeqVprbdwFHNctHAXc0A94D3JvksM0MStoM35RI6sgpSa5J8l+TPHXUxiMDNMn/AOyu\nqh3cn87hwUlda257wC7YYHLxrXXdvsvd9bejhjKm+Zh0tLeqalNdjfszXVduw/Y3iSqxzf2NGofV\nrPpmeXmZbdu27bv03hafPP5HN9/KPZ/Zvu+yH74MHFtVzwTeB3xq1A+kaviBM0n+I/CrwH3Aw4Cf\nanb8YuDIqlpJcgpwTlWdnuSzzfJVSQ4Evl1Vhw/Yb426b0nSAyWhaou/dmRMSeoJ737nVMfwjd/6\nnQ2fnyTHAp+uqqeP2k+SbwD/uKru3mibkRVoVb2tqo6pqicCLweuqKpfBT4PnNls9irgkmb50uY6\nze1XjLoPSZK2wKDu6eoNyRFrlp/DaoG5YXgCHLQfA3kLcGGS/w24Bji/WX8+8JEkO4Hvshq6kqRF\n0NPGYpKPAkvAY5LcDpwDHAxUVX0A+FdJ/jXwE+BHwC+N3Oe02qi2cCVp83rfwn3XlFu4v71xC7dr\nM3AmIkmS+md/WriSJD3QAjUWrUAlSWrBClSS1JlhZwOaN1agkiS1YIBKktSCLVxJUnds4UqSpGEM\nUEmSWrCFK0nqji1cSZI0jBWoJKkzHgcqSZKGMkAlSWrBFq4kqTv9/Ka1ibAClSSpBStQSVJ3nEQk\nSZKGMUAlSWrBFq4kqTMeBypJkoYyQCVJasEWriSpO7ZwJUnSMFagkqTOOIlIkiQNZYBKktSCLVxJ\nUnds4UqSpGEMUEmSWrCFK0nqji1cSZI0jBWoJKkzHgcqSZKGMkAlSWrBAJUkqQUDVJKkFpxEJEnq\njpOIJEnSMAaoJEkt2MKVJHXG40AlSdJQBqgkSS3YwpUkdccWriRJGsYKVJLUHStQSZI0jAEqSVIL\ntnAlSZ3xOFBJkjSUFagkqTtWoJIkaRgDVJKkFmzhSpI64yQiSZI0lAEqSVILtnAlSd2xhStJkoYx\nQCVJ3akpXzaQ5Pwku5N8bYPbfznJV5PsSPKFJP9o1EM1QCVJi+BDwKlDbv868M+r6mTg94A/GrVD\nPwOVJM29qvpCkmOH3H7lmqtXAkeN2qcBKknqzJwcB/prwJ+N2sgAlSTNrB/cdgs/uP2WzvaX5PnA\nq4GfG7WtASpJmlmPOPYEHnHsCfuuf+eL21vvK8nTgQ8Ap1XVPaO2N0AlSd3pdws3zeXBNyTHAJ8E\n/uequnWcnRmgkqS5l+SjwBLwmCS3A+cABwNVVR8A/gNwGPD+JAF+UlXPGbZPA1SS1J2eVqBV9csj\nbv914Nc3s0+PA5UkqQUDVJKkFmzhSpI6MyfHgY7FClSSpBasQCVJ3bEClSRJwxigkiS1YAtXktQZ\nJxGtk+SbzReNXpPkr5p1hybZnuSmJJclOWTN9u9NsrP5YtKTJzV4SZKmZdwW7gqwVFXPXHNqo7cA\nn6uqk4ArgLcCJDkdOL6qTgReB5zX8ZglSZq6cQM0A7Y9A7igWb6gub53/YcBquoq4JAkR+znOCVJ\ns6CmfNlC4wZoAZcluTrJrzXrjqiq3QBVdRdweLP+KOCONT/7Lcb4Zm9JkmbJuJOI/mlV3ZXkccD2\nJDexcdYP+qqYgdtu27Zt3/LS0hJLS0tjDkeSFsPy8jLLy8vTHsb4FmgSUao292iTnAN8H/g1Vj8X\n3Z3kSODzVfWUJOc1yx9vtr8ReN7eanXNfmqz9y1Jiy4JVTXwOy2nLUn9w9/+P6Y6hr9515u27PkZ\n2cJN8vAkj2yWHwG8GLgWuBQ4q9nsLOCSZvlS4JXN9qcA964PT0mSZt04LdwjgD9JUs32f1xV25P8\nNXBRktcAtwNnAlTVZ5K8JMktwA+AV09o7JKknullaTwhIwO0qr4BPOhYzqq6G3jRBj9z9v4PTZKk\n/vJMRJKk7izQ1BbPhStJUgsGqCRJLdjClSR1xpPJS5KkoQxQSZJasIUrSeqOLVxJkjSMFagkqTtW\noJIkaRgDVJKkFmzhSpI643GgkiRpKANUkqQWbOFKkrpjC1eSJA1jBSpJ6oyTiCRJ0lAGqCRJLdjC\nlSR1xxauJEkaxgpUktQZJxFJkqShDFBJklqwhStJ6o4tXEmSNIwBKklSC7ZwJUndsYUrSZKGsQKV\nJHXG40AlSdJQBqgkSS3YwpUkdccWriRJGsYKVJLUmdTilKBWoJIktWCASpLUgi1cSVJ3FqeDawUq\nSVIbBqgkSS3YwpUkdcZT+UmSNEeSnJbkxiQ3J3nzgNuPSfK5JF9NckWSfzBqnwaoJKk7NeXLAEkO\nAN4HnAo8DXhFkiev2+wPgf9cVc8A/lfg7aMeqgEqSZp3zwF2VtVtVfUT4ELgjHXbPBW4AqCqlgfc\n/iAGqCRp3h0F3LHm+q5m3Vo7gH8JkOQXgUcmOXTYTp1EJEnqzFZPIvreXbfw3++6ddRmGbBu/Uh/\nF3hfkrOAvwC+Bdw3bKcGqCRpZj3qyBN41JEn7Lt+51e3D9psF3DMmutHA3eu3aCqvs39FegjgH9Z\nVf992H3bwpUkzburgROSHJvkYODlwKVrN0jymCR7K9W3Ah8ctVMDVJLUnR7Owq2qPcDZwHbgOuDC\nqrohyblJfqHZbAm4KcmNwOHA/z7qodrClSTNvar6LHDSunXnrFn+JPDJzezTAJUkdcYzEUmSpKEM\nUEmSWrCFK0nqji1cSZI0jBWoJKkzTiKSJElDGaCSJLVgC1eS1J1anB6uFagkSS0YoJIktWALV5LU\nGWfhSpILWzUkAAAO20lEQVSkoaxAJUndsQKVJEnDGKCSJLVgC1eS1JmsTHsEW8cKVJKkFqxAJUnd\ncRKRJEkaxgCVJKkFW7iSpM54JiJJkjSUASpJUgu2cCVJ3fH7QB8oySFJPpHkhiTXJfmZJIcm2Z7k\npiSXJTlkzfbvTbIzyY4kJ09u+JIkTce4Ldz3AJ+pqqcAzwBuBN4CfK6qTgKuAN4KkOR04PiqOhF4\nHXBe56OWJPVSarqXrTQyQJP8FPDPqupDAFV1X1X9PXAGcEGz2QXNdZp/P9xsexVwSJIjuh64JEnT\nNE4F+kTgO0k+lOQrST6Q5OHAEVW1G6Cq7gIOb7Y/Crhjzc9/q1knSdLcGGcS0UHAs4DXV9VfJ3kX\nq+3bjYrlDFg3cNtt27btW15aWmJpaWmM4UjS4lheXmZ5eXnawxjf4swhIjVixlTTfv1SVT2xuf5z\nrAbo8cBSVe1OciTw+ap6SpLzmuWPN9vfCDxvb7W6Zr816r4lSQ+UhKoaVKhMXZJ67v/0B1Mdwxf/\n5He37PkZ2cJtgu+OJE9qVr0QuA64FDirWXcWcEmzfCnwSoAkpwD3rg9PSZJm3bjHgb4B+OMkDwG+\nDrwaOBC4KMlrgNuBMwGq6jNJXpLkFuAHzbaSpAWwSKfyGytAq+qrwLMH3PSiDbY/e38GJUlS33km\nIklSdxZobovnwpUkqQUDVJKkFmzhSpI6s0iTiKxAJUlqwQpUktQdK1BJkjSMASpJUgu2cCVJnXES\nkSRJGsoAlSSpBVu4kqTurCxOD9cKVJKkFqxAJUndWZwC1ApUkqQ2DFBJklqwhStJ6ozHgUqSpKGs\nQCVJ3anFKUGtQCVJasEAlSSpBVu4kqTOOIlIkiQNZYBKktSCASpJ6k5N+bKBJKcluTHJzUnevME2\nL0tyXZJrk/zfox6qn4FKkuZakgOA9wEvBO4Erk5ySVXduGabE4A3Az9bVd9L8thR+zVAJUmdST+P\nA30OsLOqbgNIciFwBnDjmm1+Hfg/q+p7AFX1nVE7tYUrSZp3RwF3rLm+q1m31pOAk5J8IclfJjl1\n1E6tQCVJM+ueu7/Ovfd8fdRmGbBufal8EHAC8M+BY4D/J8nT9lakgxigkqTurGzt3R366Cdy6KOf\nuO/6N7/x54M228VqKO51NKufha7f5ktVtQJ8M8lNwInAlze6b1u4kqR5dzVwQpJjkxwMvBy4dN02\nnwJeANBMIDoRGFraGqCSpLlWVXuAs4HtwHXAhVV1Q5Jzk/xCs81lwHeTXAf8OfBvq+qeYfu1hStJ\n6kxPZ+FSVZ8FTlq37px1138H+J1x92kFKklSC1agkqTu9LMAnQgrUEmSWjBAJUlqwRauJKk7PZ1E\nNAlWoJIktWAFKknqTBanALUClSSpDQNUkqQWbOFKkrrjJCJJkjSMASpJUgu2cCVJnckWfx/oNFmB\nSpLUghWoJKk7TiKSJEnDGKCSJLVgC1eS1J3F6eBagUqS1IYVqCSpM3ESkSRJGsYAlSSpBVu4kqTu\n2MKVJEnDGKCSJLVgC1eS1B1PJi9JkoaxApUkdcbjQCVJ0lAGqCRJLdjClSR1xxauJEkaxgCVJKkF\nW7iSpO7YwpUkScNYgUqSuuOZiCRJ0jAGqCRJLdjClSR1xlP5SZKkoaxAJUndsQKVJEnDjAzQJE9K\nck2SrzT//n2SNyQ5NMn2JDcluSzJIWt+5r1JdibZkeTkyT4ESZK23sgAraqbq+qZVfUs4B8DPwD+\nBHgL8LmqOgm4AngrQJLTgeOr6kTgdcB5kxq8JKlnqqZ72UKbbeG+CLi1qu4AzgAuaNZf0Fyn+ffD\nAFV1FXBIkiM6GKskSb2x2QD9JeCjzfIRVbUboKruAg5v1h8F3LHmZ77VrJMkaW6MPQs3yUOAfwG8\nuVm1Ua2cAesGbrtt27Z9y0tLSywtLY07HElaCMvLyywvL097GONboFm4qTEfbJJ/AfybqjqtuX4D\nsFRVu5McCXy+qp6S5Lxm+ePNdjcCz9tbra7ZX41735KkVUmoqkGFytQlqVP/4dumOobL/uY/btnz\ns5kW7iuAj625filwVrN8FnDJmvWvBEhyCnDv+vCUJM2plSlfttBYLdwkD2N1AtFvrFn9DuCiJK8B\nbgfOBKiqzyR5SZJbWJ2x++puhyxJ0vSNFaBV9SPgcevW3c1qqA7a/uz9H5okSf3lqfwkSZ3xZPKS\nJGkoK1BJUnesQCVJ0jAGqCRJLdjClSR1Z8UWriRJGsIAlSTNvSSnJbkxyc1J3jzg9tcl+Vrzvdd/\nkeTJo/ZpgEqSutPD7wNNcgDwPuBU4GnAKwYE5B9X1dOr6pnAHwDvGvVQDVBJ0rx7DrCzqm6rqp8A\nF3L/d1gDUFXfX3P1kYxxZl0nEUmSutPP40DXf0/1LlZD9QGS/BvgTcBDgBeM2qkBKkmaWd/94e3c\n/aM7Rm021vdUV9X7gfcneTnwH7j/G8cGMkAlSTPrMQ8/hsc8/Jh912+95y8HbbYLOGbN9aOBO4fs\n9uPAeaPu289AJUnd6eEkIuBq4IQkxyY5GHg5q99dvU+SE9Zc/QXg5lEP1QpUkjTXqmpPkrOB7awW\njudX1Q1JzgWurqo/Bc5O8iLg/wPuAV41ar8GqCRp7lXVZ4GT1q07Z83yb212nwaoJKk7nspPkiQN\nYwUqSepOjTz/wNywApUkqQUDVJKkFmzhSpK6089T+U2EFagkSS1YgUqSuuNhLJIkaRgDVJKkFmzh\nSpK64yQiSZI0jAEqSVILtnAlSd2xhStJkoaxApUkdccKVJIkDWOASpLUgi1cSVJ3Vvw+UEmSNIQV\nqCSpO04ikiRJwxigkiS1YAtXktQdW7iSJGkYA1SSpBZs4UqSurNiC1eSJA1hBSpJ6kyVZyKSJElD\nGKCSJLVgC1eS1B0nEUmSpGEMUEmSWrCFK0nqjqfykyRJw1iBSpK6s+JxoJIkaQgDVJKkFmzhSpK6\n4yQiSZI0jBWoJKkz5SQiSZI0jAEqSVILtnAlSd1xEpEkSRrGAJUkqQVbuJKk7vh9oJIkaRgrUElS\nd8rjQCVJ0hAGqCRJLdjClSR1ppxEJEmShrEClSR1x0lEkiRpGANUkqQWbOFKkjrjJCLts7y8PO0h\njMVxdstxdmsWxjkLY1R7SU5LcmOSm5O8ecDtBye5MMnOJF9KcsyofRqgI8zKH5Xj7Jbj7NYsjHMW\nxqh2khwAvA84FXga8IokT1632WuBu6vqRODdwO+P2q8BKknqTq1M9zLYc4CdVXVbVf0EuBA4Y902\nZwAXNMsXAy8c9VANUEnSvDsKuGPN9V3NuoHbVNUe4N4khw3baWpK3x6eZHE+aZakDlVVpj2GQZJ8\nEzh2ysPYXVVHrl2R5F8BL66q32iu/yrw7Kp645pt/qbZ5s7m+i3NNvdsdEdTm4Xb1/8AkqR2quq4\naY9hA7uAtZOCjgbuXLfNHcDjgTuTHAg8alh4gi1cSdL8uxo4IcmxSQ4GXg5cum6bTwOvapbPBK4Y\ntVOPA5UkzbWq2pPkbGA7q4Xj+VV1Q5Jzgaur6k+B84GPJNkJfJfVkB1qap+BSpI0y6bSwh11QOs0\nJflmkq8muSbJXzXrDk2yPclNSS5LcsgUxnV+kt1JvrZm3YbjSvLe5oDgHUlOnvI4z0myK8lXmstp\na257azPOG5K8eIvGeHSSK5Jcn+TaJG9o1vfq+Rwwzt9s1vft+Xxokquav5lrk5zTrD8uyZXN8/mx\nJAc16zd9wPqEx/mhJF9v1n8lydPX/MxU/o40I6pqSy+shvYtrM7UegiwA3jyVo9jyPi+Dhy6bt07\ngH/XLL8ZePsUxvVzwMnA10aNCzgd+K/N8s8AV055nOcAbxqw7VOAa1j9KOG45v9FtmCMRwInN8uP\nBG4Cnty353PIOHv1fDb3/fDm3wOBK5vn6ePAmc36/wt4XbP8r4H3N8u/BFy4hf8/B43zQ8AvDth2\nan9HXmbjMo0KdJwDWqcpPLgyX3uA7QXAS7d0REBVfQFYPyNs/bjOWLP+w83PXQUckuSIKY4TVp/X\n9c5g9cXzvqr6JrCT1f8fE1VVd1XVjmb5+8ANrM7K69XzucE49x671pvnsxnfD5vFh7Ia4AU8H/hk\ns37t382mD1jvyoBx7j3yfqPncyp/R5oN0wjQcQ5onaYCLktydZJfa9YdUVW7YfVFDXjc1Eb3QIev\nG9fhzfr1z/G3mP5z/PqmDfaf1rRGpz7OJMexWjFfyYN/z715PteM86pmVa+ezyQHJLkGuAu4HLgV\nuLdq36lh1v6db/qA9UmNs6qubm76veb5fGeSh6wfZ6MPf0fqkWkE6KB3en2ayfRPq+qfAC9h9UXq\nn9Gv8Y2jb8/x+4Hjq+pkVl+43tmsn+o4kzyS1QrojU2Ft9F9922cvXs+q2qlqp7JaiX/HFbbyRuN\nZf04w5TGmeSpwFuq6inAs4HHsNq+HzROtmqcmg3TCNBxDmidmqbyoKr+DvgUqy8Gu/e2bpIcCfzt\n9Eb4ABuNaxerBwTvNdXnuKr+rqr2vvD8Efe3Fac2zmZCy8XAR6rqkmZ1757PQePs4/O5V1V9D/hv\nwCnAo7N6Eu/1Y9k3zox5wPoEx3namq7DT1j9PLQ3z6f6bRoBOs4BrVOR5OHNu32SPAJ4MXAtq+M7\nq9nsVcAlA3cweeGB74rXjuss7h/XpcArAZKcwmorbffWDBFYN84mjPb6ReBvmuVLgZc3szKfAJwA\n/NUWjfGDwPVV9Z416/r4fD5onH17PpM8dm8bOcnDgBcB1wOfZ/WAdHjg382lbPKA9QmO88a9z2eS\nsPo57drnc5p/R+q7acxcAk5jdUbhTlbbJ1OfTdWM6wmszgq+htXgfEuz/jDgc82YLwcePYWxfZTV\nd78/Bm4HXg0cutG4WP3qnluArwLPmvI4Pwx8rXluP8XqZ417t39rM84bWD0P5VaM8bnAnjW/6680\n/yc3/D1P4/kcMs6+PZ//qBnbjmZc/75Z/wRWP7O9mdUZuQ9p1j8UuKj5+78SOG7K4/zz5vf6tea5\nffg0f+9eZufiiRQkSWrBc+FKktSCASpJUgsGqCRJLRigkiS1YIBKktSCASpJUgsGqCRJLfz/IHFm\nPnSxnWwAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fc1409c44d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# merge cloud1 back again\n",
    "test.merge(cloud1, transformation)\n",
    "####\n",
    "\n",
    "filepath = 'out3.tif'\n",
    "test.region('region.png')\n",
    "atlaas.merge('atlaas.*x*.tif', filepath)\n",
    "g = gdal2( filepath )\n",
    "band = g.dtm()\n",
    "hist(band)\n",
    "show(band, 0, 3)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}