{ "cells": [ { "cell_type": "code", "execution_count": 70, "metadata": {}, "outputs": [], "source": [ "library(sp)\n", "library(maps)\n", "library(rgdal)\n", "library(dplyr)\n", "library(scales)\n", "library(rgeos)\n", "library(geosphere)\n", "library(rmapshaper)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Provinces et territoires\n", "\n", "Téléchargement du fichier SHP contenant les limites des provinces et des territoires\n", "\n", " wget http://www12.statcan.gc.ca/census-recensement/2011/geo/bound-limit/files-fichiers/gpr_000b11a_f.zip" ] }, { "cell_type": "code", "execution_count": 71, "metadata": { "scrolled": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "OGR data source with driver: ESRI Shapefile \n", "Source: \"donnees/provinces/gpr_000b11a_f.shp\", layer: \"gpr_000b11a_f\"\n", "with 13 features\n", "It has 6 fields\n" ] } ], "source": [ "gpr_000b11a_f <- readOGR(\"donnees/provinces/gpr_000b11a_f.shp\",encoding = 'ISO-8859-1') %>% ms_simplify()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Données tabulaires contenues dans le fichier SHP" ] }, { "cell_type": "code", "execution_count": 72, "metadata": { "scrolled": true }, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\n", "
PRIDUPRNOMPRFNOMPRFABREVPRANOMPRAABREVrmapshaperid
059 British Columbia / Colombie-Britannique Colombie-Britannique C.-B. British Columbia B.C. 0
124 Quebec / Québec Québec Qc Quebec Que. 1
262 Nunavut Nunavut Nt Nunavut Nvt. 2
311 Prince Edward Island / Île-du-Prince-Édouard Île-du-Prince-Édouard Î.-P.-É. Prince Edward Island P.E.I. 3
447 Saskatchewan Saskatchewan Sask. Saskatchewan Sask. 4
560 Yukon Yukon Yn Yukon Y.T. 5
646 Manitoba Manitoba Man. Manitoba Man. 6
735 Ontario Ontario Ont. Ontario Ont. 7
813 New Brunswick / Nouveau-Brunswick Nouveau-Brunswick N.-B. New Brunswick N.B. 8
961 Northwest Territories / Territoires du Nord-Ouest Territoires du Nord-Ouest T.N.-O. Northwest Territories N.W.T. 9
1048 Alberta Alberta Alb. Alberta Alta. 10
1110 Newfoundland and Labrador / Terre-Neuve-et-LabradorTerre-Neuve-et-Labrador T.-N.-L. Newfoundland and Labrador N.L. 11
1212 Nova Scotia / Nouvelle-Écosse Nouvelle-Écosse N.-É. Nova Scotia N.S. 12
\n" ], "text/latex": [ "\\begin{tabular}{r|lllllll}\n", " & PRIDU & PRNOM & PRFNOM & PRFABREV & PRANOM & PRAABREV & rmapshaperid\\\\\n", "\\hline\n", "\t0 & 59 & British Columbia / Colombie-Britannique & Colombie-Britannique & C.-B. & British Columbia & B.C. & 0 \\\\\n", "\t1 & 24 & Quebec / Québec & Québec & Qc & Quebec & Que. & 1 \\\\\n", "\t2 & 62 & Nunavut & Nunavut & Nt & Nunavut & Nvt. & 2 \\\\\n", "\t3 & 11 & Prince Edward Island / Île-du-Prince-Édouard & Île-du-Prince-Édouard & Î.-P.-É. & Prince Edward Island & P.E.I. & 3 \\\\\n", "\t4 & 47 & Saskatchewan & Saskatchewan & Sask. & Saskatchewan & Sask. & 4 \\\\\n", "\t5 & 60 & Yukon & Yukon & Yn & Yukon & Y.T. & 5 \\\\\n", "\t6 & 46 & Manitoba & Manitoba & Man. & Manitoba & Man. & 6 \\\\\n", "\t7 & 35 & Ontario & Ontario & Ont. & Ontario & Ont. & 7 \\\\\n", "\t8 & 13 & New Brunswick / Nouveau-Brunswick & Nouveau-Brunswick & N.-B. & New Brunswick & N.B. & 8 \\\\\n", "\t9 & 61 & Northwest Territories / Territoires du Nord-Ouest & Territoires du Nord-Ouest & T.N.-O. & Northwest Territories & N.W.T. & 9 \\\\\n", "\t10 & 48 & Alberta & Alberta & Alb. & Alberta & Alta. & 10 \\\\\n", "\t11 & 10 & Newfoundland and Labrador / Terre-Neuve-et-Labrador & Terre-Neuve-et-Labrador & T.-N.-L. & Newfoundland and Labrador & N.L. & 11 \\\\\n", "\t12 & 12 & Nova Scotia / Nouvelle-Écosse & Nouvelle-Écosse & N.-É. & Nova Scotia & N.S. & 12 \\\\\n", "\\end{tabular}\n" ], "text/markdown": [ "\n", "| | PRIDU | PRNOM | PRFNOM | PRFABREV | PRANOM | PRAABREV | rmapshaperid | \n", "|---|---|---|---|---|---|---|---|---|---|---|---|---|\n", "| 0 | 59 | British Columbia / Colombie-Britannique | Colombie-Britannique | C.-B. | British Columbia | B.C. | 0 | \n", "| 1 | 24 | Quebec / Québec | Québec | Qc | Quebec | Que. | 1 | \n", "| 2 | 62 | Nunavut | Nunavut | Nt | Nunavut | Nvt. | 2 | \n", "| 3 | 11 | Prince Edward Island / Île-du-Prince-Édouard | Île-du-Prince-Édouard | Î.-P.-É. | Prince Edward Island | P.E.I. | 3 | \n", "| 4 | 47 | Saskatchewan | Saskatchewan | Sask. | Saskatchewan | Sask. | 4 | \n", "| 5 | 60 | Yukon | Yukon | Yn | Yukon | Y.T. | 5 | \n", "| 6 | 46 | Manitoba | Manitoba | Man. | Manitoba | Man. | 6 | \n", "| 7 | 35 | Ontario | Ontario | Ont. | Ontario | Ont. | 7 | \n", "| 8 | 13 | New Brunswick / Nouveau-Brunswick | Nouveau-Brunswick | N.-B. | New Brunswick | N.B. | 8 | \n", "| 9 | 61 | Northwest Territories / Territoires du Nord-Ouest | Territoires du Nord-Ouest | T.N.-O. | Northwest Territories | N.W.T. | 9 | \n", "| 10 | 48 | Alberta | Alberta | Alb. | Alberta | Alta. | 10 | \n", "| 11 | 10 | Newfoundland and Labrador / Terre-Neuve-et-Labrador | Terre-Neuve-et-Labrador | T.-N.-L. | Newfoundland and Labrador | N.L. | 11 | \n", "| 12 | 12 | Nova Scotia / Nouvelle-Écosse | Nouvelle-Écosse | N.-É. | Nova Scotia | N.S. | 12 | \n", "\n", "\n" ], "text/plain": [ " PRIDU PRNOM \n", "0 59 British Columbia / Colombie-Britannique \n", "1 24 Quebec / Québec \n", "2 62 Nunavut \n", "3 11 Prince Edward Island / Île-du-Prince-Édouard \n", "4 47 Saskatchewan \n", "5 60 Yukon \n", "6 46 Manitoba \n", "7 35 Ontario \n", "8 13 New Brunswick / Nouveau-Brunswick \n", "9 61 Northwest Territories / Territoires du Nord-Ouest \n", "10 48 Alberta \n", "11 10 Newfoundland and Labrador / Terre-Neuve-et-Labrador\n", "12 12 Nova Scotia / Nouvelle-Écosse \n", " PRFNOM PRFABREV PRANOM PRAABREV\n", "0 Colombie-Britannique C.-B. British Columbia B.C. \n", "1 Québec Qc Quebec Que. \n", "2 Nunavut Nt Nunavut Nvt. \n", "3 Île-du-Prince-Édouard Î.-P.-É. Prince Edward Island P.E.I. \n", "4 Saskatchewan Sask. Saskatchewan Sask. \n", "5 Yukon Yn Yukon Y.T. \n", "6 Manitoba Man. Manitoba Man. \n", "7 Ontario Ont. Ontario Ont. \n", "8 Nouveau-Brunswick N.-B. New Brunswick N.B. \n", "9 Territoires du Nord-Ouest T.N.-O. Northwest Territories N.W.T. \n", "10 Alberta Alb. Alberta Alta. \n", "11 Terre-Neuve-et-Labrador T.-N.-L. Newfoundland and Labrador N.L. \n", "12 Nouvelle-Écosse N.-É. Nova Scotia N.S. \n", " rmapshaperid\n", "0 0 \n", "1 1 \n", "2 2 \n", "3 3 \n", "4 4 \n", "5 5 \n", "6 6 \n", "7 7 \n", "8 8 \n", "9 9 \n", "10 10 \n", "11 11 \n", "12 12 " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "gpr_000b11a_f@data" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Extraction des coordonnées limites de chacune des provinces" ] }, { "cell_type": "code", "execution_count": 73, "metadata": {}, "outputs": [], "source": [ "bbox_prov <- lapply(gpr_000b11a_f@polygons, bbox)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Sélection des limites de l'Alberta et affichage des polygones" ] }, { "cell_type": "code", "execution_count": 74, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAIAAAByhViMAAAACXBIWXMAABJ0AAASdAHeZh94\nAAAgAElEQVR4nO3dZ5hU9eH34TOzu4A0Rem9CUqx914Q7A01ttiisRuNRpM/xm5ssRs19ogN\nVAQxaKKigKIiCig2QATpvZdld+c8L/BRdxYj/cz89r5fwe6EfLkul/3saZOK4zgCACD/pZMe\nAADA+iHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAAC\nIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAA\nAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewA\nAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHs\nAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh\n7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAAC\nIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAA\nAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewA\nAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHs\nAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh\n7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAIh7AAAAiHsAAACIewAAAJRmPQAyCVL5yW9\nAFgr1eskvQBygiN2AACBEHYAAIEQdgAAgRB2AACBcPME/GDGjBmTx305ddr0BQsXJr0FWC2b\n1q7duFHDKIqabdmxfv36Sc+B5Ak7+EHfvn3PO++8pFcAa+ORRx4555xzkl4ByUvFcZz0BsgJ\ny5cvn/bd1607bP/yc//af9+9k54D/Ip3B7937Emnjf9yRBRFjVptVa1ataQXQfIcsYMfVKtW\nrVWLFlEU1axZs85mmyU9B/gVNWvWjKJo5ZdtpOogiiI3TwAABEPYAQAEQtgBAARC2AEABELY\nAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC\n2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAE\nQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEA\nBELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgB\nAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELY\nAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC\n2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAE\nQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEA\nBELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgB\nAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELY\nAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC\n2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAE\nQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEA\nBELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgB\nAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AEABELY\nAQAEQtgBAARC2AEABELYAQAEQthBtlQqlfQEAFgbwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAg\nEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOssVxnPQEAFgbwg4AIBDCDgAgEMIO\nACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDC\nDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQ\nwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAg\nEMIOACAQwg6AvJRKpZKeADlH2AEABELYAQAEQtgBAARC2AEABELYAQAEQtgBAARC2AGQl+I4\nTnoC5BxhBwAQCGEHABAIYQcAEAhhBwAQCGEHABAIYQcAEAhhBwAQCGEHABAIYQcAEAhhBwAQ\nCGEHABAIYQcAEAhhBwAQCGEHABAIYQcAEAhhBwAQCGEHABAIYQcAEAhhBwAQCGEHABAIYQcA\nEAhhBwAQCGEHQF7KZDJJT4CcI+wAyEtxHCc9AXKOsAMgLzliBxUJOwDykiN2UJGwAwIXx3Gb\nTjs0aNm+tLR05UdGfvb5TnsdEEXRMy/0PvWsc9fizyysXW99TmStZDKZVCqV9ArILcIOCNzw\nT0cWFxdvskm1Qe+9v+5/2tRp09f9D2G9cMQOKhJ2QOCef/GlU0884Tfdj3nhxT4VP7t4yZLu\nJ5++eZPWu+7b5bPRX6z84MBBg7fZZa8GLduf8NszJ0+ZGkXR6C+/2rfb4SedcfYBhxx5xHEn\nlZWVddppj+Li4nv/8XCTth1q1m+6675dvhk7bqP+xSo9YQcVCTsgZJlMpvfLfU875cSTTuje\np99rK1asyHpBv9cGnH7KiTMmfHNot4NOPO13ZWVls2bPPuHUsx646/bJY0e3atHitLPPX/nK\nD4cN36ZTx4/fG9j/pecLCgpGDx86e87cq/56/et9X5w1cdy2nTvddd8/Nvrfr1Jz8wRUJOyA\nkA15/4PGjRp22Kr9tp071au7xVvvDMp6wS477XDkYYcUFRX935/+OH3GzPHfTXj1328cdOB+\n++y1R1FR0Q1//cuwTz5ZGRC1ata86o9/qFWz5o//2y0233z8FyO26dQxjuOaNWssWLhwo/7d\nKr1MxhE7yFaY9ADIOc7vhOSFl/p8NvrLus3aRlG0cNGiXi+9cmi3g37+gubNmq78RVFRUbOm\nTWbMnDVp8pQ333630057rPx4y+bN586bF0VR/Xp10+lyPwwXFRU+9tTT/V57PZVK1ahRvVHD\nBhvjr8T/F8exWycgiyN2kE3YBaOkpOSlV14d+eHg2ZPGzZ40bvh7A1959bXly4t//pqJ309a\n+YvS0tJJk6c0a9qkQf16Rxx68OjhQ0cPHzrqoyEvPfevLTbfPIqirKqLoujlvv1f7NOvb69n\nhr838KzTTtk4fyl+FMdx5K5YKE/YQTYPUAjG2+8OrlNns/Zbtl35284dO9Sps9kbb77189d8\n/MmIvv3/XVZWduud97Tbsk3zZk0P7XZQ/wFvjBj1WUlJyW133Xv2BZes8j+J5cuLp02f3rhR\nw8aNGn4/afKTPZ9bvHjJxvhb8f+5xg4qEnZAsF54sc9Rhx3yY5alUqkjDzvkhZd+ujd209q1\n/3LFZY888a8GLdu/OfDd5558NJVKtWje7PGH7jvp9LPrt2j31sBBTz/6UMU/+YhDD27Tafvj\njz16eXFx4zYdjj/1jIvP+/2UqdP69Ou/kf5uOLgOq5LyhQE/WTovVWPzN/q91K3LAUlPAX7F\n871fPvV355Ytmh1FUVS9TtJzICc4YgfZnIqFvJCJvfMEZBN2AOQl19hBRcIOgLzkcSdQkbAD\nIC953AlUJOwAyEtOxUJFwg6AvBTHkeN1kEXYAZCXMpmMU7GQRdgBkJc8hxUqEnaQzZOxIC9k\nMp5jB9mEHQB5yRE7qEjYQTbHACAvZDIZX6yQRdgBkJfiOI7cFwvlCTsA8pJ3noCKhB0Aecnj\nTqAiYQfZfKeAvOCIHVQk7CCbO+0gL2Qy3isWsgk7yOYZCgDkKWEH2YQd5IU4jh2wgyzCDoC8\nFEcedwLZhB0AQCCEHQBAIIQdAEAghB0AQCCEHQBAIIQdAEAghB0AQCCEHQBAIIQdAEAghB0A\nQCCEHQBAIIQdAEAghB0AQCCEHQBAIIQdAEAghB0AQCCEHQBAIIQdAEAghB0AQCCEHQBAIIQd\nAHkpFaWiOE56BeQWYQdAXiooSGeEHZQn7ADIS0VFRXGcSXoF5BZhB0BeiuM4ilJJr4DcIuwA\nyEulpaXplLCDcoQdAHmprCyTEnZQnrADAAiEsAMACISwAwAIhLADAAiEsAMgL8VRHLl5AsoT\ndgAAgRB2AACBEHYA5C/vFQvlCDsA8pdr7KAcYQcAEAhhBwAQCGEHABAIYQcAEAhhBwAQCGEH\nABAIYQcAEAhhBwAQCGEHABAIYQcAEAhhBwAQCGEHABAIYQfZvpswMekJALA2hB2UU61q1fMu\n+eOkyVOSHgIAa0zYQTlD3nq9oCB98FHHddxx9zPPvSjpOQCwBoQdlLPTDtsN/u+ABg3qT5k2\n7fMvvkx6DgCsAWEH2fbYbZeBA/pVKSpq2KB+0lsAYA0IO1iFkZ+Nnj1nbrVqVZ965rmSkpJM\nJpP0IgD4dcIOVmHGzJlxHPfp99qZ515Uv0W7Aw49KulFAPDrhB2sQrcuB8yZPH7xzElbtmnd\nrm2bQUPe7/fagKRHAeVkMplUKukRkGMKkx4AOWrzOptFUTTms+FRFLXtvOMf/3x1JpM55sjD\nk94F/CCO41Sk7KAcR+zg173e98Wp06a98ebbSQ8BfhLHsa6DLMIOVmHFihV33HPff94cuPK3\nW7Zp3aRx4y+++jqO42SHAT/KZDKO2EEWp2JhFTrssNu3302IomjTTWu3adXq/LPP7Hrg/g89\n+sSw4Z/uuvOOSa8DoiiKCgoKMn7WgvKEHazCId0OeuDhR2+/+foZM2YOem/oORddGkXR7886\nffttOyc9DfhBQUGBg+iQJeWrAn6ydN6Pv2zStmNRUeFbr73SpnWrp597IZPJnH7KSem0qxcg\nV/z93geuu+nWxbMmR1EUVa+T9BzICb5LwaoNf+/tmTNnPfjo46lU6vRTTjrzt6eoOgBynG9U\nsGrHnXLmipKSTh06JD0E+GXunYDyhB2s2ohRn53Q/ehXXn3trvseTHoLsAqpVMrFRJDFzROw\napec//vb7rq3oKDgv2+/88VXXz18311FRUVJjwJ+kopSkbKD8oQdrNqtN17buFHDgoKC3i/3\nfeLpZ6/9v6uaN2ua9CjgJ+l0Ko6EHZQj7OAXXXLBuVEUNWncaPD7Q5cuW5b0HKAiF9lBOa6x\ng19x9BGHFRUV3XH3fUkPAYBfIezg153Q/ei+/QcsXrwk6SHAT1xfBxUJO/h1t1x/zaLFi6//\n221JDwHKcSIWsgg7+HXNmjbpfvQR9/zj4dYdt586bfp7Qz+85sZbFixcmPQuqNS8cxJU5OYJ\nWC3333l708ZN7vnHQ3sc0K2kpGTqtOn3P/RI61YtmzZpfPopJ86YOeuNN9/+x913NG3SOOml\nAFRewg5WS90tNr/jb9fvsP02z/d++fMvvjzskK5NGjUaM+7bz0Z/0f3k01OpVLVqVfc7+Ih3\n3+iv7WDjiKM4SjkZC+UIO1gDJx3f/aTju2d9cPny5dWqVZs9Z26rrbf981+vf+aJfyayDQBc\nYwfrqlq1alEU1d1i8yaNG02cNCnpOVCJOF4HWYQdrDeHHdztvaEfvjN4SNJDoFJw8wRUJOxg\nvTnx+GOjKFqwwN2yACRD2MF6s/OO2++4/Xa/Pfu8T0aMTHoLhC/lzgmoQNjB+jRs8FtlZWU3\n3vr3ip9avHjJlT2ubd1x+3bb7vzYUz03/jYIj3OxkEXYwfqUTqfvuuXm/gPeqHjQbtgnn9xx\nz/2LFi3atHatcy78w3/eGpjJZBIZCUCohB2sZ+edc2b9enUvu6rHy337FxcXDxv+6cqP77/P\n3scdc2RJaenHQwY2bFD/4KOOO/Pci1z9DcB6JOxg/evUscOQ9z847pTTazVovuu+XQYNeT+K\nolQqNX3GzC3q1ImiaNr4rx+6586nn3thm132WrhoUdJ7AQiEBxTD+rd1+3aTJk/+06WXzJg5\n86FHnzjulDMaN2p4/LFHvTf0w8cfun/la84758ztt9tm74MO3X2/rs//67FtOnVMdjPknTiO\n3T0BWVLOBMFPls5b73/kmLHjzr/0iokTv//2uwlRFLVs0fy7L3+6/G7st+O323XvY486oufj\nD6/3/2sI29/vfeC6m25dPGtyFEVR9TpJz4Gc4IgdbFjttmz79r/7lpWVFdauF0XRbTde+/PP\nbtmm9TFHHt7r5VcaNWzwt+v/WljoSxLWhEN2UJ7vIrAxFBQUXH7JhVWqVjmh+zFZn3rmiX8u\nWrz4jnvuP/WkE5yQhdXnjBNUJOxgI/n7LTf+0qe++mbMHrvtouoAWEfuioWEvfrv18d9O36f\nPfcoLi5OegvkkziOnYuFLMIOEjZn7rw4jm+9857Du5845P0PnF0CYK0JO0jY1GnToijabZed\n3npn0D5dD7v59juTXgT5IZPJpL1dLJQn7CBhPa68fMTQQUPeHHDphec3bdL4rXcGJb0I8kM6\nnc44wg3lCTtI3nbbdi4sLLz79pvbtW3z0cfDvYcsrI6qVapkysqSXgG5RdhBrhj95VcDBw35\n/VlnpNO+MOHXFRUVZWI/BUE5vn9ArmjXtk3VqlXmzJlbWlqa9BbIA+6KhYqEHeSKKlWq3Hzd\n1c/3funCy/6U9BbIA5lMxr0TkEXYQQ65/JKL9t9vnw8++jiKIo+1g//NETuoSNhBbqldq9bE\n7yddd/NttRu2aNK2wz5dD5sw8fukR0EucsQOKhJ2kFsaNqy/dNmy2+++t1nTJk0aN/7iy6/3\nOKDbkiVLk94FOSeOI0fsIIuwg9xy7llnNmvS5Nijjhgz6uNhg9/adecdPf0EVqmwsCD21QHl\nCTvILdt27jj+yxHPPP7PlQ892WO3XWbMnFWzftNd9+2S9DTILSknYqECYQc57eqrrhjy5oB6\n9eo2adw46S2QW+I4diYWsgg7yHV77r7rggULdthum6SHQG5xVyxUJOwg173/wUcrVpSccOwx\nSQ8BINcJO8h1Q4Z+EEXR4iWLkx4CuSWdTq+8Mxb4kbCDXHfRuee0bNF874MO7fl8r6S3AJDT\nhB3kulq1an47+tMu++932tnn9+nXP+k5kCvcPAEVCTvIA+l0ul/vZ6tWrTpt+oykt0CucPME\nVCTsIG9sWrvWw489OXPWrKSHAJCjhB3kjWFD3v5u4sQuhx1TUlKS9BZIXhzHjtdBFmEHeaNF\ns2YnndB90uQpK9+UAio519hBRYVJDwBWS2lpabejjhv47uA/XHheQUFB0nMAyEXCDvLDv//z\n5sB3Bz/5zwfOOPXkpLcAkKOc0IH80LRxo3Q6/fU3Y5MeAkDuEnaQH3bcfrvttun87pD3iouL\nk94COSGVSkXeeALKE3aQNw7t1mXY8E+bt9/mhRf7JL0FkpdOp3UdZBF2kDduvKbHlHFfzJ49\ne+KkSUlvgeSlUinvFQtZhB3kk+IVKzJx3KZVq6SHQPJSqZSsgyzCDvJJy+bNO3fscOXV1y5e\nvKTiZ5csWbps2fKNvwoSkUqlIhfZQXnCDvLMm/37TJo85XcXXDx7zpyff3zGzFnb7LrX2Rdc\nktQw2MjcPAEVCTvIMw0a1P/n/Xe/2Kdfm447LFy06MeP3/L3uxcsWPhS31cnT5ma4DzYaJyK\nhYqEHeSfBx95vKCgoMdVl9euVWvlR+6678GHH3uyQf16VatWOfCwoyZ+7+4KwpdKeUMxyCbs\nIM8sXbr0kxEjb73x2isv++ms65y5c4uLi4875sixo4bPmjVn74MOHTHqswRHwkbgGjuoSNhB\nnqlevXrTJo1HjPwsiqKysrKJ30/KZDKHH9ItiqKTju/eoEH98V+OLC0t22GP/eo0btVxx933\n6XpYr5deSXo1ABuDsIP8U1RU9P6HH038ftJpZ5/fcutta9RrcvypZxQWFm7Vvl0URZttWnvq\nt1/+t3+fc393RudOHaZPn3nuxZdNnzEz6dWwnrl5AioqTHoAsMb6PPf09nvsu0/Xw6ZMnXbT\ntT1Ky8o+Hv5ps6ZNfv6agw7Y76AD9ouiaMWKFZs1bnnx5Vf27vmka5IIjK6DLMIO8s+WbVtv\nskm1SZOnHH/s0T2uvPx/v7hKlSr33H7LuRdf1vP5XqedfOLGWQhAIoQd5J8aNWrMnTy+qKio\noKBgdV7f6+VX6my2Wd0tttjQw2BjcgQaKhJ2kJeqVau2mq+84v+ueXfwe/fc/rdDux20QSfB\nRibsoCI3T0DInuz53F33/aNVyxZn/vaUpLfAepZOp6LYVXZQjrCDYA39cNjFl1/VZf99x33+\nSc2aNZKeA+tZOp2WdZBF2EGw9ul62DadOjzx8ANJD4ENIp1Ox47YQXnCDoK1xeabb9WuXdMm\njZMeAhtEHMcusoMswg6CdeD++7wzeMiKFSuSHgIbRBzHkfsnoDxhB8E67ugjJ34/qf12uzzZ\n89mysrKk58B65jwsVCTsIFjHHnXExK8/a96s6VnnXbzNLnu1327nbkd2X7Zs+S+9/omnn33s\nqZ4bcyGsi0zGqVjIJuwgZM2aNhn0n9eOP/boRo0a7rzDDgMHDfndBRf/0ouvufFvTz3z3Mac\nB+vCqVioyAOKIXy9ez6x8he777rzJVf8+Tfdjxkx6vP3P/jwuGOOOvd3Z6z8VElJyYyZs9Jp\nP+yRNzKZTNITIOekXKMAP1k6L+kFG1Ymk6m2eaOSkpIaNarXrlVr2vQZ9evVq1+vbqOGDb4e\nM3bylKkFBQXL5kwtLPQjH3ngkSf+dfHlVxbPmxFFUVS9TtJzICf46RwqkXQ6PWzQW+NGf7J4\n5uTvv/n8oXvuPPmE49pt2XbpsmU7bb/dg3f/vbS0dMbMWUnPhNWSyWSiyKlYKMfP5VC5bLdt\n55W/KCwsPO+cM3/+qZkzZ51/6eVTp01v0rhREtNgzWQyGZfYQRZH7IAf1K27RSqVGjHqs6SH\nwGqJ49gRO8gi7IAffDN2bBzHX30zxjXp5AXvPAEVORUL/GDr9u3POfO0ex54qCCdvuziC5yQ\nJcd53AlU5Igd8JNHHrhnx+23u/O+f/zmtLOS3gKrw4MdoBxhB5Qz/L2Bm2yySc0aNZIeAr8i\nlUrpOsgi7IBsh3brMui99wcOGpz0EPhfUqmUroMswg7I1rvnkzWq1zju5DP26nLIJyNGJj0H\nVi2VSjkVC1mEHZAtnU6P+Wz4KSceP2bct3t1OeTFPv2SXgSr4FQsVCTsgFXYvM5m99952+cf\nvbfJJpv89Yab5y9YkPQiyOZULFQk7IBf1KBB/Q/ffXPC95MOO/Y3S5cuS3oOlJPyrBOoQNgB\n/0u7tm0+eOc/w4Z/etIZZy9ZsjTpOfAT19hBRcIO+BXbb7vN5Zdc+Oq/Xx/52edJb4GfpFKp\nWNdBecIO+HXvDB6y1x677bn7rkkPgXKci4Uswg74dY0aNPzq6zH3PPCQK+3IHXEcKzvIIuyA\nX/f8vx7dvM5ml13VY+hHw1b5gmXLlrtzliQoOyhH2AG/bpNNNnng7juiKFq8eHHWp26/+77t\ndttn570PaNtpx3eHvJfEOiqpsrIy98VCFmEHrJauB+7ftk3r2+669+cfnDN37k23/n3ylCkl\nJSXpgvT/XXtjUvOohEpLSz3xBLIUJj0AyBuZsrJyv81kzjrv4uIVxZPGjt60du39Dj5ixYoV\nSW2jEorj2KlYyOKIHbC6VpSUzJk7b+asWSt/+/RzL/Qf8Eb/F1/YtHbtKIqWL1v+3YSJ47+b\nkOREgMpN2AGra8hbA8Z/N+Gk089Z+dvHnupZp85mXbvsv/K3r/TqOWPmrGtvutVdFABJEXbA\n6mrZvPl1Pf48cNDgJ3s+e/X1N9erW3fu3HmlpaUrP9uoYcPjjj7y2V4v3nnvP5LdCVBppWLP\n7YYfLZ2X9IJct2LFis477zlu/HdVq1YtKSnp1uWA115+4ecvOPSYE6ZMnTbqoyFJLaTyeODh\nR//U45plc6ZFURRVr5P0HMgJbp4A1kCVKlW+GfXx/3jBt99NaNKo0UbbQ2WWTqcdm4AsTsUC\n69OsWbPnzJ27YOHCpIcQvoKCAl0HWYQdsD49+uC9o7/8qsMOuz3b60XvP8YGVVBQ4IgdZBF2\nwPrU/agjvh45rHbt2qeede6fr7k+6TmErKyszAOKIYuwA9azLdu0/urTD8/93Rn3P/TIG2++\nnfQcguVwHVQk7IAN4uH77kqn05MmT0l6CEAlIuyADWWTatV6vfTKmwPfTXoIYXIeFioSdsCG\n8tSjD47+8suuRxzb/eTTk94CUCkIO2BDOe7oI6d/981BB+zXf8AbS5YsTXoOoXHEDioSdsCG\ntWTp0t122alGjepJDyE0HlAMFQk7YMM6/JCDPxw2/JMRI5MeQmgKCtKRsIPyhN5Ns+sAABFg\nSURBVB2wYf3likvr1d3imBN/+8DDjw5+b2jScwhHOp2WdZBF2AEb3Bv9Xq5Vq9ZlV/XYt9vh\n474dn/QcApFKpRyxgyzCDtjgOnfc+ovhQ594+IEoitJp/+ywfmQymcj9E1Cef2GBjWTZsmVR\nFNWoUSPpIQDBEnbARvLhx8MLCwsXLVqU9BAC4ZZYqEjYARvJ871fPun47m3btE56CECwhB2w\nkRx0wH7PvND7lVdfS3oI4XCFHWQRdsBG8seLL4iiqFatWkkPIRDOxEJFwg7YSDps1T6dTr83\n9MOkhxAId8VCRYVJDwAqi/r163U9cP+bb7/zq2++SafTu+280yUXnOvtPllrZWVl/vuBLMIO\n2HgGvNL7j3+++tEn/pVOp194sc+Wbdsc2u2gpEeRr4QdVORULLBR3XXrTYtmThr7+SdRFD3f\n++UpU6dlMpmkR5GXysrKdB1kEXZAAurXq3vKice/2Kdv0y07tum0w6zZs5NeRP6J49h9sZBF\n2AHJeObxf86dMr53zyemTJ129gV/eLbXi0kvIs84DwsVCTsgMdWrVz/+2KPvuvWm/7z59mln\nn79s2fKkF5F3PPIEyhF2QMIuOu+cgw/qkslk5s2fn/QWgPwm7ICEfTb6y/4DXj//nLMaNWyQ\n9BbySSqVcsAOsgg7IGHVqlaJo2j6jJkumQJYR8IOSFi7LdvuvusuI0Z9Nmbst0lvAchvwg5I\n3r577zFh4vdb77Br/wFvJL2FvJFKpZyJhSzCDkje3677a9mi2Y0bNTzvkj9+OGz4O4OHLF9e\nXFxcnPQugDwj7ICckE6nX+/74tx583ffv+sBhxzVsFX7Bi3bP/6vZ7746uukpwHkDWEH5IpO\nHbYeNvitAa/0fuSBe7beqn1xcfHZF1zyt9vvSnoXQN4oTHoAwE86d+zQuWOHKIrOOfO0KIqK\nNq1fUFBQXFycTqeLioqSXgeQ6xyxA3LXkYce3PP5XlvvsNs+XQ9Legs5x/NxoCJhB+Sul59/\n+vzfnzV//oLhn4686urrSkpKkl4EkNOEHZDTHrz773OnjO/W5YA77rn/oUefSHoOQE4TdkAe\neO3lFw47pOt1N9+2YOHCpLcA5C5hB+SHR+6/e978+U8/+0LSQ8gpHlEM5Qg7ID80athws802\n7ffa60kPIVfEcRxF7p+AcoQdkDca1q8/+P2hr//3raSHkBPiOJZ1kEXYAXnjoAP3LykpubLH\ntUkPIScUFBTEsVOxUI4HFAN5oLi4uF7zLRctXnzYwV1P6H500nPICYWFhcIOsgg7IA/Mmz9/\n0eLFt9xwzdln/LbuFlskPYec4IgdVORULJAHGjZoULVKlR7X3VS/Rbubb78z6TnkhMLCAlkH\nWYQdkB9uvv7qzp067LX7bvc9+EjSW8gJqVTKETvIIuyA/HD5JReN/GDwHrvvsmDhwl4vvZL0\nHJLnrlioSNgB+eSW669p1aL5iaf/bpd9DlyyZGnSc0haStpBOcIOyCepVOqLTz747Um/Gf7p\nyPP/cPnb7w5KehGJyWQySU+AnCPsgDyTTqeffuyhLdu07vl8r0OOPmHsuG+TXkQyMpnYATvI\n4nEnQF76ZtTHmUxm00Ytbrjlju227fz+Bx917tjh+qv/nPQuNp5MJuMtxSCLsAPyVTqd3rxO\nnRde6vNsrxdr1qw5bPgnV1x6Ua2aNZPexUaSyWQcsYMsTsUCeWzg6/2OP/ao63r8eeCAflOm\nTnt38HtJL2LjSaVSkaedQHmO2AF5rE2rVs89+WgURZlMprCwcNTno4849OCkR7GRpFIpXQdZ\nHLEDAtGxw1ZvveMmWaBSE3ZACM69+LJRn41u1aJF0kPYeFKusIMKhB0Qgv332TuKon333iPp\nIQBJEnZACE7+zXG77LTDBZdeceOtd/z4/qGz58y5/6FH5sydO2Hi98nOY0NIpVKRuyegPDdP\nAIEYOvA/hx5zwjU33rLbLjvPmzf/rXfeHTHq809HjvrLtTcsX148/osRzZs1TXoj65uug/Ic\nsQMCUVBQcOuN10ZR9Fyvl35z2lkv9+0/YeL3Lz371N577B7HcZ9+/ZMeyHrmaSdQUerHcxZA\ntHRe0gtYV1XrNCwtLS0qLFw+b/qPHzy8+4kffzJi1EdDGjaon+A21q8nez77+4suK1kwM4qi\nqHqdpOdATnDEDgjKX//8p0wmc/aZp/38g717PjF33txnXuid1CqAjcMRO/gZR+yCsHDhotq1\na2V98KAjjv12/HffjBxWVFSUyCrWO0fsoCJH7IDQVKy6KIoevu/OCRO/v/amWzf+HoCNRtgB\nlUKbVq2OOfLwx57qOXee47JAsIQdUFk8+fADc+fNvezKHmPGfpv0FtaDTCb25hOQRdgBlUXt\n2rVuveHap5974bKr/i/pLawHS5YsKSjwXQzK8YBioBK54tKLvx4zttfLr2QymXRaE+S3ktLS\ndLog6RWQW/y7BlQu9evXW7p0Wcutt534/aSkt7BO4jh2IhayCDugcrnpmh49rrx88pSp1918\n2/eTJic9B2B9EnZA5ZJOp2/461/22HWXp555bt9uhyc9h7UXx3HkkB2UJ+yAyujFZ59Mp9Nb\nt2+X9BDWXhzHKWUH5Qk7oDIa/cXXmUzm5N8cl/QQ1l5paWnKHTBQni8JoDI66MD9tth88/P/\ncPnHn4xIegtrqbS0NJ12xA7KEXZAJfVGv5eWLVs+8rPPkx7CWiory6RTvotBOb4kgEpq5uxZ\nmUzGZXb5K51OxXGc9ArILcIOqKTeevvdwsLCHbbbNukhrKXCwsJMJpP0Csgtwg6opM4567TS\n0tJ7H3w46SGspUwm471iIYuwAyqprdu3b99uy14vvZL0ENaBroPyhB1Qed1647WjPh/9/gcf\nJT2EtZFKucYOsgk7oPKqVrVqFEW1atVMeghrI51O6zrIIuyAymvnHXeIomj+/AVJD2FtpNPp\nSNlBecIOqLyWLVsWRVGVKlWSHsLaSKfTTsVCFmEHVF5ffPV1FEULFi5MeghrI51Ox5Gwg3KE\nHVB5dT1w/3Zbtj3/D5eXlJQkvYU1lkqlHLCDLMIOqNTe6PfidxMm9u7TN+khrLF0yjV2kE3Y\nAZVaqxYt2rZufepZ5343YWLSW1gz3lIMKhJ2QGV38m+6R1G0cNGipIewZlxjBxUJO6Cy+8sV\nl6VTqeGfjkx6CGvGc+ygImEHVHbTZ87MxHHNGjWSHsKaSafTkSN2UJ6wAyq7ls2b16pV69OR\no5IewppxxA4qEnZAZdd55z0XLVpUt+4WSQ9hzaRSKUfsIIuwAyq7MWPHXX3VFZdfcmHSQ1gz\nqVRK10EWYQdUdiWlpdt27pRO+/cQyHv+IQMquypFRVOnTUt6BWsllfQAyDHCDqjsMnE8b/78\npFewxjydGCoSdkBld8ShB9/zwMPLli1PegjAuhJ2QGV36YXnzV+wYMB/3kx6CMC6EnZAZbf3\nnrs3bdL4saeeTnoIwLoSdgBR544dxn47PukVrAV3T0A5wg4gatG82ZIlS5JewZpx8wRUJOwA\nooKCgtlz5hYXFyc9BGCdCDuA6Imnnz34oAOrVq2a9BCAdSLsAKKuB+7/6chRjtgB+U7YAUS7\n7LTD1GnT58ydl/QQ1oBr7KAiYQdUdgP+++bV19987FFHNG7UMOktAOtE2AGV3eAh70dR1POx\nh5MeArCuhB1Q2e2+265xHH89ZkzSQwDWlbADKrV58xcMHjI0iqKvx4xNegvAuipMegBAYv77\n9jsnn3nOnDlzWzRvtv8+eyc9B2BdOWIHVF7nXPiHzTbd9J03Xv129KeNGjZIeg7AuhJ2QOXV\nrGmTKIr223uvgoKCpLcArAfCDqi8fn/m6eO/mzDk/Q9W/nboh8O++OrrZCcBrAthB1Rep51y\nYutWLbscfsweB3S7+PKr9jzw4L0OPMT7TwD5S9gBldqAPr1q1aw5dtz4Bx95vG3r1vMXLJg9\nZ27So1hdqaQHQK4RdkCl1m7LtrMnjZv1/diyRbM7ddiqVs2atWvXSnoUq8VbikFFwg7gB/8d\n+M6F555dq2bNpIcArCVhB/CDTCb20BMgrwk7gCiKokwms3z58vr16iU9BGDtCTuAKIqiSVOm\nRFHkiB2Q14QdQBRFUVlZWRRFceR6/LwRx7HbYiGLsAOIoihq3bJlOpWa+P2kpIcArL3CpAcA\n5IpWrVredue9ZWWZ4445snYtDz3JCw7ZQTmO2AH84O0Bfb/9bsLvzr/42ptuTXoLwNoQdgA/\naNGsWfG86dt27nTPAw9dc+MtSc8BWGPCDqCcER8MOqRrl38+/pQ3NgDyjrADKCeVSk2aPGXu\nvHljxn2b9Bb+F+UNFQk7gGx33nJTjerVTzr97KSHAKwZYQeQrWuX/c87+8yvx4yZNn1G0lsA\n1oCwA1iFnXbcftmy5Y3bbH3+Hy5PegvA6hJ2AKtw3NFHThn3ZfVNNnn136+7lgvIF8IOYNUa\nN2qYyWSOPuKwVMpTcIH8IOwAflHbtm0efOTx7Xff98BDj/5uwsQFCxd+9c2YpEcB/CJhB/CL\nRn04+KZre0yZNu3Dj4dvv/u+Ldpv02GH3YZ+OCzpXQCrJuwAflE6ne5x5eUzJ4z5ZuSw1q1a\ntmjerH69uvt0PezWO+9JehpRHMfOkUOWwqQHAOSBpk0afzr03ZW/Pun0s6++/ua9dt9trz12\nS3ITQAWO2AGsme227VxWVla9+iZJDwHIJuwA1szF552TTqdHfjY66SEA2YQdwJqpXr16YWHB\nJyNGJj0EIJuwA1hjN/y1x4OPPP6Xa26YN39+0lsAfuLmCYA1dtUfL5k3b95td907bPinb/37\nFU8wBnKEI3YAa+PWG6/t80LPgYMGj/rcxXbJiOM4UtRQniN28IPS0tJFTqslZMWKkiVLlqzO\nK2vWrFFUVLSh96ymzTbdNIqib8aMa9G8WdJbKqNly5fHcbTybHitKrUKC31Hgyjlza1hpaFD\nh+65555JrwDWxtChQ3ffffekV0DyhB38YPHixd98803SK9anOnXqrM7LiouLly5duqHH/G+F\nhYW1atVao//J/PnzE//nK5PJzJkzp169esnOqLRKS0sXLVq08r/zrbbaqkaNGkkvguQJOwCA\nQLh5AgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIO\nACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDC\nDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQ\nwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAg\nEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4A\nIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIO\nACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDC\nDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQ\nwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAg\nEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4A\nIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIOACAQwg4AIBDCDgAgEMIO\nACAQwg4AIBDCDgAgEMIOACAQ/w/I2d8WrSmRvQAAAABJRU5ErkJggg==", "text/plain": [ "plot without title" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "bbox_prov_alberta <- bbox_prov[gpr_000b11a_f@data$PRNOM==\"Alberta\"][[1]]\n", "\n", "map(gpr_000b11a_f,fill = TRUE, col = 'seashell',namefield = 'PRFNOM',\n", " xlim=c(bbox_prov_alberta[1,1],bbox_prov_alberta[1,2]),\n", " ylim=c(bbox_prov_alberta[2,1],bbox_prov_alberta[2,2]))\n", "invisible(text(coordinates(gpr_000b11a_f),\n", " labels=as.character(gpr_000b11a_f$PRFNOM),\n", " cex=0.6))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Parcs nationaux\n", "\n", "Source des données:\n", "\n", " mkdir donnees/limites_parcs_nationaux\n", " cd donnees/limites_parcs_nationaux \n", " wget http://ftp.cartes.canada.ca/pub/pc_pc/National-parks_Parc-national/national_parks_boundaries/limites_parcs_nationaux.shp.zip\n", " unzip limites_parcs_nationaux.shp.zip\n", " rm limites_parcs_nationaux.shp.zip" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Je transforme les données des parcs nationaux dans le même système de référence que les polygones des provinces" ] }, { "cell_type": "code", "execution_count": 75, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "OGR data source with driver: ESRI Shapefile \n", "Source: \"donnees/limites_parcs_nationaux/limites_parcs_nationaux.shp\", layer: \"limites_parcs_nationaux\"\n", "with 47 features\n", "It has 7 fields\n", "Integer64 fields read as strings: OBJECTID \n" ] } ], "source": [ "limites_parcs_nationaux <- spTransform(\n", " x = rgdal::readOGR(\"donnees/limites_parcs_nationaux/limites_parcs_nationaux.shp\"),\n", " CRSobj = CRS(proj4string(gpr_000b11a_f))) %>% ms_simplify()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Données tabulaires contenues dans le fichier SHP des parcs nationaux" ] }, { "cell_type": "code", "execution_count": 76, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\n", "
OBJECTIDShape_LengShape_Areaglobalidprecisionhnom_parc_enom_parc_frmapshaperid
01 399100.63 521622672 {F5BD4565-8591-41AA-99E4-7A6F41E47A09} 100 Pacific Rim National Park Reserve of Canada Réserve de parc national du Canada Pacific Rim 0
12 86148.58 187215384 {DCB2DCFE-17BC-4DFA-BE9C-FE9A4624FCA4} 100 Elk Island National Park of Canada Parc national du Canada Elk Island 1
23 239712.76 1331729306 {106EEED8-DDC6-447B-9226-006507154D40} 100 Glacier National Park of Canada Parc national du Canada des Glaciers 2
34 313783.71 764739692 {B2A0614F-A228-4244-9BCC-6B9B96433447} 100 Grasslands National Park of Canada Parc national du Canada des Prairies 3
45 814277.44 10942527257 {B8B6DB08-8E91-4216-9D19-982A7E857920} 100 Jasper National Park of Canada Parc national du Canada Jasper 4
56 282894.28 3827306859 {D365A5E5-1ED3-4BAF-9DDF-232A5E5EA4B9} 100 Prince Albert National Park of Canada Parc national du Canada de Prince Albert 5
67 83392.95 259236547 {571E0061-18B4-4485-9072-7138D567A683} 100 Mount Revelstoke National Park of Canada Parc national du Canada du Mont-Revelstoke 6
78 574328.45 10924352608 {C52E5D05-7345-42F8-8377-FEE3D4A7F522} 200 Wapusk National Park of Canada Parc national du Canada Wapusk 7
89 134090.44 498878286 {615C33BF-A6D6-4B69-8317-436C1D03FBD1} 100 Waterton Lakes National Park of Canada Parc national du Canada des Lacs-Waterton 8
910 74399.91 211648312 {EB67B4F3-0C53-45F9-974B-3BE819D84BCA} 100 Fundy National Park of Canada Parc national du Canada Fundy 9
1011 145056.93 546631833 {F4D2C860-F422-4A44-84C5-2A4BA233DFE6} 100 La Mauricie National Park of Canada Parc national du Canada de la Mauricie 10
1112 762246.23 95866783 {9290F5AD-5FD3-4A6C-A777-BDF080AD924A} 100 Mingan Archipelago National Park Reserve of Canada Réserve de parc national du Canada de l'Archipel-de-Mingan 11
1213 361707.33 1860082905 {AFD5D06C-3740-43CF-A3C2-3BE738AED300} 100 Pukaskwa National Park of Canada Parc national du Canada Pukaskwa 12
1314 103792.51 31184813 {CA4BBEAD-3BA3-4B9B-82C1-1CC1EFDCAFB4} 100 Sable Island National Park Reserve of Canada Réserve de parc national du Canada de l'Île-de-Sable 13
1415 524862.27 1257873345 {A0AC46C1-8A41-4718-AFA7-A0BC689B087B} 100 Saguenay-St. Lawrence Marine Park Parc marin du Saguenay-Saint-Laurent 14
1516 624231.91 11827725094 {E125CA8E-DCD5-4640-90A5-533BDFB19B8D} 200 Aulavik National Park of Canada Parc national du Canada Aulavik 15
1617 632131.31 17853403254 {FF705D51-F6F2-4CD6-84C6-E57AF223C39C} 200 Tuktut Nogait National Park of Canada Parc national du Canada Tuktut Nogait 16
1718 541881.60 4607855048 {60AB71D9-F267-4A78-B61F-844D5E4D0D17} 200 Nááts'ihch'oh National Park Reserve of Canada Réserve de parc national du Canada de Nááts'ihch'oh 17
1819 149762.52 161656326 {2EFC258D-7344-4EBE-8BF4-A8B200942C98} 100 Bruce Peninsula National Park of Canada Parc national du Canada de la Péninsule-Bruce 18
1920 70694.38 117124125 {3EB03934-6824-4097-B1A4-74E688A9C49F} 100 Fathom Five National Marine Park of Canada Parc marin national du Canada Fathom Five 19
2021 26326.70 16075792 {81607DCE-BD39-41E3-ADEE-C371D4667D7B} 100 Point Pelee National Park of Canada Parc national du Canada de la Pointe-Pelée 20
2122 422768.76 1799122547 {B577C7C8-F927-4C3A-A6B5-6DA9DCF1101A} 100 Gros Morne National Park of Canada Parc national du Canada du Gros-Morne 21
2223 2124790.50 9180689134 {B4CFB708-8ABD-4198-8340-231965E6901B} 200 Torngat Mountains National Park of Canada Parc national du Canada Monts-Torngat 22
2324 94220.18 14042598 {BC8D02DA-591B-4BA2-A34C-814B8C39E755} 100 Georgian Bay Islands National Park of Canada Parc national du Canada des Îles-de-la-Baie-Georgienne 23
2425 82455.80 397193668 {57653BC7-3C8F-47FA-8414-4E23474C6875} 100 Kejimkujik National Park and National Historic Site of Canada Parc national et lieu historique national Kejimkujik du Canada 24
2526 35755.34 25322482 {8BEB0C58-98D9-4946-AF21-3514ED879146} 100 Kejimkujik National Park Seaside Parc national Kejimkujik Bord de mer 25
2627 274964.28 1361036954 {092E9172-2FDC-453D-94F8-26BFF4132FB6} 200 Kootenay National Park of Canada Parc national du Canada Kootenay 26
2728 197069.12 1268469933 {D7164C4B-1609-4FBF-945B-2864EBBDAAD2} 200 Yoho National Park of Canada Parc national du Canada Yoho 27
2829 747138.11 6738767165 {914676D3-62D7-4ADD-A75C-ABF0607C6564} 200 Banff National Park of Canada Parc national du Canada Banff 28
2930 382133.80 3049951157 {691C8725-89CA-47D8-9C1E-69E4A99037CE} 100 Riding Mountain National Park of Canada Parc national du Canada du Mont-Riding 29
3031 1167897.57 18397952413 {590D0639-7974-47A1-9795-040E7D405338} 200 Auyuittuq National Park of Canada Parc national du Canada Auyuittuq 30
3132 872045.03 41167722479 {2D2C1EB6-0D7F-4A6D-A77F-31C99E499397} 200 Quttinirpaaq National Park of Canada Parc national du Canada Quttinirpaaq 31
3233 811615.27 19625619733 {09D47FA6-A6E7-4082-9813-1D4ED9AF06F3} 200 Ukkusiksalik National Park of Canada Parc national du Canada du Ukkusiksalik 32
3334 1355170.32 21424060989 {5FDB4A6D-8A01-4728-9BE8-EC005103F54B} 200 Sirmilik National Park of Canada Parc national du Canada Sirmilik 33
3435 1778463.99 1439600089 {978ABCB6-206D-47A4-85F7-D1A1E40C4710} 100 Gwaii Haanas National Park Reserve and Haida Heritage Site Réserve de parc national et site du patrimoine haïda Gwaii Haanas34
3536 175418.98 37052584 {F0104540-00C1-4A34-B9A7-AB5F3D0D23D7} 100 Gulf Islands National Park Reserve of Canada Réserve de parc national du Canada des Îles-Gulf 35
3637 118096.26 242115530 {4245A9AA-5234-4D5B-A419-D348FEB84526} 100 Kouchibouguac National Park of Canada Parc national du Canada Kouchibouguac 36
3738 1270941.90 43160601764 {323E3E60-B323-43C9-BE80-ABE621FBEE96} 200 Wood Buffalo National Park of Canada Parc national du Canada Wood Buffalo 37
3839 108077.13 246232998 {CE0D9EFB-F2AB-4F8B-B03E-44A57D1593C2} 100 Forillon National Park of Canada Parc national du Canada Forillon 38
3940 187503.01 35663421 {2C53DA06-A8E9-4FDE-9834-9FD67647D037} 10 Prince Edward Island National Park of Canada Parc national du Canada de l'Île-du-Prince-Édouard 39
4041 370575.29 5547513490 {F5CE7051-831B-4A3F-B820-F72A574C655D} 200 Kluane National Park of Canada Parc national du Canada Kluane 40
4142 553903.31 9260673207 {E02AF2ED-EC14-4FF3-B22F-9EE14C31CF5E} 200 Ivvavik National Park of Canada Parc national du Canada Ivvavik 41
4243 602842.41 4106240006 {339112BF-0C82-45B3-B52D-ED44737BF2AC} 200 Vuntut National Park of Canada Parc national du Canada Vuntut 42
4344 586086.66 15267021368 {2856F646-CF59-40F9-8D43-6E918F32564C} 200 Kluane National Park Reserve of Canada Réserve de parc national du Canada du Canada Kluane 43
4445 296692.13 406346468 {A8D35956-D43F-401A-BF53-797FA1063444} 100 Terra Nova National Park of Canada Parc national du Canada Terra-Nova 44
4546 188812.32 969640516 {CD356B3C-727E-46AC-96A8-F4A0887D9AF2} 100 Cape Breton Highlands National Park of Canada Parc national du Canada des Hautes-Terres-du-Cap-Breton 45
4647 1135818.33 28316219360 {1C3DDDEF-F7BB-44F7-BFDF-9BCED3036A70} 200 Nahanni National Park Reserve of Canada Réserve de parc national du Canada Nahanni 46
\n" ], "text/latex": [ "\\begin{tabular}{r|llllllll}\n", " & OBJECTID & Shape\\_Leng & Shape\\_Area & globalid & precisionh & nom\\_parc\\_e & nom\\_parc\\_f & rmapshaperid\\\\\n", "\\hline\n", "\t0 & 1 & 399100.63 & 521622672 & \\{F5BD4565-8591-41AA-99E4-7A6F41E47A09\\} & 100 & Pacific Rim National Park Reserve of Canada & Réserve de parc national du Canada Pacific Rim & 0 \\\\\n", "\t1 & 2 & 86148.58 & 187215384 & \\{DCB2DCFE-17BC-4DFA-BE9C-FE9A4624FCA4\\} & 100 & Elk Island National Park of Canada & Parc national du Canada Elk Island & 1 \\\\\n", "\t2 & 3 & 239712.76 & 1331729306 & \\{106EEED8-DDC6-447B-9226-006507154D40\\} & 100 & Glacier National Park of Canada & Parc national du Canada des Glaciers & 2 \\\\\n", "\t3 & 4 & 313783.71 & 764739692 & \\{B2A0614F-A228-4244-9BCC-6B9B96433447\\} & 100 & Grasslands National Park of Canada & Parc national du Canada des Prairies & 3 \\\\\n", "\t4 & 5 & 814277.44 & 10942527257 & \\{B8B6DB08-8E91-4216-9D19-982A7E857920\\} & 100 & Jasper National Park of Canada & Parc national du Canada Jasper & 4 \\\\\n", "\t5 & 6 & 282894.28 & 3827306859 & \\{D365A5E5-1ED3-4BAF-9DDF-232A5E5EA4B9\\} & 100 & Prince Albert National Park of Canada & Parc national du Canada de Prince Albert & 5 \\\\\n", "\t6 & 7 & 83392.95 & 259236547 & \\{571E0061-18B4-4485-9072-7138D567A683\\} & 100 & Mount Revelstoke National Park of Canada & Parc national du Canada du Mont-Revelstoke & 6 \\\\\n", "\t7 & 8 & 574328.45 & 10924352608 & \\{C52E5D05-7345-42F8-8377-FEE3D4A7F522\\} & 200 & Wapusk National Park of Canada & Parc national du Canada Wapusk & 7 \\\\\n", "\t8 & 9 & 134090.44 & 498878286 & \\{615C33BF-A6D6-4B69-8317-436C1D03FBD1\\} & 100 & Waterton Lakes National Park of Canada & Parc national du Canada des Lacs-Waterton & 8 \\\\\n", "\t9 & 10 & 74399.91 & 211648312 & \\{EB67B4F3-0C53-45F9-974B-3BE819D84BCA\\} & 100 & Fundy National Park of Canada & Parc national du Canada Fundy & 9 \\\\\n", "\t10 & 11 & 145056.93 & 546631833 & \\{F4D2C860-F422-4A44-84C5-2A4BA233DFE6\\} & 100 & La Mauricie National Park of Canada & Parc national du Canada de la Mauricie & 10 \\\\\n", "\t11 & 12 & 762246.23 & 95866783 & \\{9290F5AD-5FD3-4A6C-A777-BDF080AD924A\\} & 100 & Mingan Archipelago National Park Reserve of Canada & Réserve de parc national du Canada de l'Archipel-de-Mingan & 11 \\\\\n", "\t12 & 13 & 361707.33 & 1860082905 & \\{AFD5D06C-3740-43CF-A3C2-3BE738AED300\\} & 100 & Pukaskwa National Park of Canada & Parc national du Canada Pukaskwa & 12 \\\\\n", "\t13 & 14 & 103792.51 & 31184813 & \\{CA4BBEAD-3BA3-4B9B-82C1-1CC1EFDCAFB4\\} & 100 & Sable Island National Park Reserve of Canada & Réserve de parc national du Canada de l'Île-de-Sable & 13 \\\\\n", "\t14 & 15 & 524862.27 & 1257873345 & \\{A0AC46C1-8A41-4718-AFA7-A0BC689B087B\\} & 100 & Saguenay-St. Lawrence Marine Park & Parc marin du Saguenay-Saint-Laurent & 14 \\\\\n", "\t15 & 16 & 624231.91 & 11827725094 & \\{E125CA8E-DCD5-4640-90A5-533BDFB19B8D\\} & 200 & Aulavik National Park of Canada & Parc national du Canada Aulavik & 15 \\\\\n", "\t16 & 17 & 632131.31 & 17853403254 & \\{FF705D51-F6F2-4CD6-84C6-E57AF223C39C\\} & 200 & Tuktut Nogait National Park of Canada & Parc national du Canada Tuktut Nogait & 16 \\\\\n", "\t17 & 18 & 541881.60 & 4607855048 & \\{60AB71D9-F267-4A78-B61F-844D5E4D0D17\\} & 200 & Nááts'ihch'oh National Park Reserve of Canada & Réserve de parc national du Canada de Nááts'ihch'oh & 17 \\\\\n", "\t18 & 19 & 149762.52 & 161656326 & \\{2EFC258D-7344-4EBE-8BF4-A8B200942C98\\} & 100 & Bruce Peninsula National Park of Canada & Parc national du Canada de la Péninsule-Bruce & 18 \\\\\n", "\t19 & 20 & 70694.38 & 117124125 & \\{3EB03934-6824-4097-B1A4-74E688A9C49F\\} & 100 & Fathom Five National Marine Park of Canada & Parc marin national du Canada Fathom Five & 19 \\\\\n", "\t20 & 21 & 26326.70 & 16075792 & \\{81607DCE-BD39-41E3-ADEE-C371D4667D7B\\} & 100 & Point Pelee National Park of Canada & Parc national du Canada de la Pointe-Pelée & 20 \\\\\n", "\t21 & 22 & 422768.76 & 1799122547 & \\{B577C7C8-F927-4C3A-A6B5-6DA9DCF1101A\\} & 100 & Gros Morne National Park of Canada & Parc national du Canada du Gros-Morne & 21 \\\\\n", "\t22 & 23 & 2124790.50 & 9180689134 & \\{B4CFB708-8ABD-4198-8340-231965E6901B\\} & 200 & Torngat Mountains National Park of Canada & Parc national du Canada Monts-Torngat & 22 \\\\\n", "\t23 & 24 & 94220.18 & 14042598 & \\{BC8D02DA-591B-4BA2-A34C-814B8C39E755\\} & 100 & Georgian Bay Islands National Park of Canada & Parc national du Canada des Îles-de-la-Baie-Georgienne & 23 \\\\\n", "\t24 & 25 & 82455.80 & 397193668 & \\{57653BC7-3C8F-47FA-8414-4E23474C6875\\} & 100 & Kejimkujik National Park and National Historic Site of Canada & Parc national et lieu historique national Kejimkujik du Canada & 24 \\\\\n", "\t25 & 26 & 35755.34 & 25322482 & \\{8BEB0C58-98D9-4946-AF21-3514ED879146\\} & 100 & Kejimkujik National Park Seaside & Parc national Kejimkujik Bord de mer & 25 \\\\\n", "\t26 & 27 & 274964.28 & 1361036954 & \\{092E9172-2FDC-453D-94F8-26BFF4132FB6\\} & 200 & Kootenay National Park of Canada & Parc national du Canada Kootenay & 26 \\\\\n", "\t27 & 28 & 197069.12 & 1268469933 & \\{D7164C4B-1609-4FBF-945B-2864EBBDAAD2\\} & 200 & Yoho National Park of Canada & Parc national du Canada Yoho & 27 \\\\\n", "\t28 & 29 & 747138.11 & 6738767165 & \\{914676D3-62D7-4ADD-A75C-ABF0607C6564\\} & 200 & Banff National Park of Canada & Parc national du Canada Banff & 28 \\\\\n", "\t29 & 30 & 382133.80 & 3049951157 & \\{691C8725-89CA-47D8-9C1E-69E4A99037CE\\} & 100 & Riding Mountain National Park of Canada & Parc national du Canada du Mont-Riding & 29 \\\\\n", "\t30 & 31 & 1167897.57 & 18397952413 & \\{590D0639-7974-47A1-9795-040E7D405338\\} & 200 & Auyuittuq National Park of Canada & Parc national du Canada Auyuittuq & 30 \\\\\n", "\t31 & 32 & 872045.03 & 41167722479 & \\{2D2C1EB6-0D7F-4A6D-A77F-31C99E499397\\} & 200 & Quttinirpaaq National Park of Canada & Parc national du Canada Quttinirpaaq & 31 \\\\\n", "\t32 & 33 & 811615.27 & 19625619733 & \\{09D47FA6-A6E7-4082-9813-1D4ED9AF06F3\\} & 200 & Ukkusiksalik National Park of Canada & Parc national du Canada du Ukkusiksalik & 32 \\\\\n", "\t33 & 34 & 1355170.32 & 21424060989 & \\{5FDB4A6D-8A01-4728-9BE8-EC005103F54B\\} & 200 & Sirmilik National Park of Canada & Parc national du Canada Sirmilik & 33 \\\\\n", "\t34 & 35 & 1778463.99 & 1439600089 & \\{978ABCB6-206D-47A4-85F7-D1A1E40C4710\\} & 100 & Gwaii Haanas National Park Reserve and Haida Heritage Site & Réserve de parc national et site du patrimoine haïda Gwaii Haanas & 34 \\\\\n", "\t35 & 36 & 175418.98 & 37052584 & \\{F0104540-00C1-4A34-B9A7-AB5F3D0D23D7\\} & 100 & Gulf Islands National Park Reserve of Canada & Réserve de parc national du Canada des Îles-Gulf & 35 \\\\\n", "\t36 & 37 & 118096.26 & 242115530 & \\{4245A9AA-5234-4D5B-A419-D348FEB84526\\} & 100 & Kouchibouguac National Park of Canada & Parc national du Canada Kouchibouguac & 36 \\\\\n", "\t37 & 38 & 1270941.90 & 43160601764 & \\{323E3E60-B323-43C9-BE80-ABE621FBEE96\\} & 200 & Wood Buffalo National Park of Canada & Parc national du Canada Wood Buffalo & 37 \\\\\n", "\t38 & 39 & 108077.13 & 246232998 & \\{CE0D9EFB-F2AB-4F8B-B03E-44A57D1593C2\\} & 100 & Forillon National Park of Canada & Parc national du Canada Forillon & 38 \\\\\n", "\t39 & 40 & 187503.01 & 35663421 & \\{2C53DA06-A8E9-4FDE-9834-9FD67647D037\\} & 10 & Prince Edward Island National Park of Canada & Parc national du Canada de l'Île-du-Prince-Édouard & 39 \\\\\n", "\t40 & 41 & 370575.29 & 5547513490 & \\{F5CE7051-831B-4A3F-B820-F72A574C655D\\} & 200 & Kluane National Park of Canada & Parc national du Canada Kluane & 40 \\\\\n", "\t41 & 42 & 553903.31 & 9260673207 & \\{E02AF2ED-EC14-4FF3-B22F-9EE14C31CF5E\\} & 200 & Ivvavik National Park of Canada & Parc national du Canada Ivvavik & 41 \\\\\n", "\t42 & 43 & 602842.41 & 4106240006 & \\{339112BF-0C82-45B3-B52D-ED44737BF2AC\\} & 200 & Vuntut National Park of Canada & Parc national du Canada Vuntut & 42 \\\\\n", "\t43 & 44 & 586086.66 & 15267021368 & \\{2856F646-CF59-40F9-8D43-6E918F32564C\\} & 200 & Kluane National Park Reserve of Canada & Réserve de parc national du Canada du Canada Kluane & 43 \\\\\n", "\t44 & 45 & 296692.13 & 406346468 & \\{A8D35956-D43F-401A-BF53-797FA1063444\\} & 100 & Terra Nova National Park of Canada & Parc national du Canada Terra-Nova & 44 \\\\\n", "\t45 & 46 & 188812.32 & 969640516 & \\{CD356B3C-727E-46AC-96A8-F4A0887D9AF2\\} & 100 & Cape Breton Highlands National Park of Canada & Parc national du Canada des Hautes-Terres-du-Cap-Breton & 45 \\\\\n", "\t46 & 47 & 1135818.33 & 28316219360 & \\{1C3DDDEF-F7BB-44F7-BFDF-9BCED3036A70\\} & 200 & Nahanni National Park Reserve of Canada & Réserve de parc national du Canada Nahanni & 46 \\\\\n", "\\end{tabular}\n" ], "text/markdown": [ "\n", "| | OBJECTID | Shape_Leng | Shape_Area | globalid | precisionh | nom_parc_e | nom_parc_f | rmapshaperid | \n", "|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|\n", "| 0 | 1 | 399100.63 | 521622672 | {F5BD4565-8591-41AA-99E4-7A6F41E47A09} | 100 | Pacific Rim National Park Reserve of Canada | Réserve de parc national du Canada Pacific Rim | 0 | \n", "| 1 | 2 | 86148.58 | 187215384 | {DCB2DCFE-17BC-4DFA-BE9C-FE9A4624FCA4} | 100 | Elk Island National Park of Canada | Parc national du Canada Elk Island | 1 | \n", "| 2 | 3 | 239712.76 | 1331729306 | {106EEED8-DDC6-447B-9226-006507154D40} | 100 | Glacier National Park of Canada | Parc national du Canada des Glaciers | 2 | \n", "| 3 | 4 | 313783.71 | 764739692 | {B2A0614F-A228-4244-9BCC-6B9B96433447} | 100 | Grasslands National Park of Canada | Parc national du Canada des Prairies | 3 | \n", "| 4 | 5 | 814277.44 | 10942527257 | {B8B6DB08-8E91-4216-9D19-982A7E857920} | 100 | Jasper National Park of Canada | Parc national du Canada Jasper | 4 | \n", "| 5 | 6 | 282894.28 | 3827306859 | {D365A5E5-1ED3-4BAF-9DDF-232A5E5EA4B9} | 100 | Prince Albert National Park of Canada | Parc national du Canada de Prince Albert | 5 | \n", "| 6 | 7 | 83392.95 | 259236547 | {571E0061-18B4-4485-9072-7138D567A683} | 100 | Mount Revelstoke National Park of Canada | Parc national du Canada du Mont-Revelstoke | 6 | \n", "| 7 | 8 | 574328.45 | 10924352608 | {C52E5D05-7345-42F8-8377-FEE3D4A7F522} | 200 | Wapusk National Park of Canada | Parc national du Canada Wapusk | 7 | \n", "| 8 | 9 | 134090.44 | 498878286 | {615C33BF-A6D6-4B69-8317-436C1D03FBD1} | 100 | Waterton Lakes National Park of Canada | Parc national du Canada des Lacs-Waterton | 8 | \n", "| 9 | 10 | 74399.91 | 211648312 | {EB67B4F3-0C53-45F9-974B-3BE819D84BCA} | 100 | Fundy National Park of Canada | Parc national du Canada Fundy | 9 | \n", "| 10 | 11 | 145056.93 | 546631833 | {F4D2C860-F422-4A44-84C5-2A4BA233DFE6} | 100 | La Mauricie National Park of Canada | Parc national du Canada de la Mauricie | 10 | \n", "| 11 | 12 | 762246.23 | 95866783 | {9290F5AD-5FD3-4A6C-A777-BDF080AD924A} | 100 | Mingan Archipelago National Park Reserve of Canada | Réserve de parc national du Canada de l'Archipel-de-Mingan | 11 | \n", "| 12 | 13 | 361707.33 | 1860082905 | {AFD5D06C-3740-43CF-A3C2-3BE738AED300} | 100 | Pukaskwa National Park of Canada | Parc national du Canada Pukaskwa | 12 | \n", "| 13 | 14 | 103792.51 | 31184813 | {CA4BBEAD-3BA3-4B9B-82C1-1CC1EFDCAFB4} | 100 | Sable Island National Park Reserve of Canada | Réserve de parc national du Canada de l'Île-de-Sable | 13 | \n", "| 14 | 15 | 524862.27 | 1257873345 | {A0AC46C1-8A41-4718-AFA7-A0BC689B087B} | 100 | Saguenay-St. Lawrence Marine Park | Parc marin du Saguenay-Saint-Laurent | 14 | \n", "| 15 | 16 | 624231.91 | 11827725094 | {E125CA8E-DCD5-4640-90A5-533BDFB19B8D} | 200 | Aulavik National Park of Canada | Parc national du Canada Aulavik | 15 | \n", "| 16 | 17 | 632131.31 | 17853403254 | {FF705D51-F6F2-4CD6-84C6-E57AF223C39C} | 200 | Tuktut Nogait National Park of Canada | Parc national du Canada Tuktut Nogait | 16 | \n", "| 17 | 18 | 541881.60 | 4607855048 | {60AB71D9-F267-4A78-B61F-844D5E4D0D17} | 200 | Nááts'ihch'oh National Park Reserve of Canada | Réserve de parc national du Canada de Nááts'ihch'oh | 17 | \n", "| 18 | 19 | 149762.52 | 161656326 | {2EFC258D-7344-4EBE-8BF4-A8B200942C98} | 100 | Bruce Peninsula National Park of Canada | Parc national du Canada de la Péninsule-Bruce | 18 | \n", "| 19 | 20 | 70694.38 | 117124125 | {3EB03934-6824-4097-B1A4-74E688A9C49F} | 100 | Fathom Five National Marine Park of Canada | Parc marin national du Canada Fathom Five | 19 | \n", "| 20 | 21 | 26326.70 | 16075792 | {81607DCE-BD39-41E3-ADEE-C371D4667D7B} | 100 | Point Pelee National Park of Canada | Parc national du Canada de la Pointe-Pelée | 20 | \n", "| 21 | 22 | 422768.76 | 1799122547 | {B577C7C8-F927-4C3A-A6B5-6DA9DCF1101A} | 100 | Gros Morne National Park of Canada | Parc national du Canada du Gros-Morne | 21 | \n", "| 22 | 23 | 2124790.50 | 9180689134 | {B4CFB708-8ABD-4198-8340-231965E6901B} | 200 | Torngat Mountains National Park of Canada | Parc national du Canada Monts-Torngat | 22 | \n", "| 23 | 24 | 94220.18 | 14042598 | {BC8D02DA-591B-4BA2-A34C-814B8C39E755} | 100 | Georgian Bay Islands National Park of Canada | Parc national du Canada des Îles-de-la-Baie-Georgienne | 23 | \n", "| 24 | 25 | 82455.80 | 397193668 | {57653BC7-3C8F-47FA-8414-4E23474C6875} | 100 | Kejimkujik National Park and National Historic Site of Canada | Parc national et lieu historique national Kejimkujik du Canada | 24 | \n", "| 25 | 26 | 35755.34 | 25322482 | {8BEB0C58-98D9-4946-AF21-3514ED879146} | 100 | Kejimkujik National Park Seaside | Parc national Kejimkujik Bord de mer | 25 | \n", "| 26 | 27 | 274964.28 | 1361036954 | {092E9172-2FDC-453D-94F8-26BFF4132FB6} | 200 | Kootenay National Park of Canada | Parc national du Canada Kootenay | 26 | \n", "| 27 | 28 | 197069.12 | 1268469933 | {D7164C4B-1609-4FBF-945B-2864EBBDAAD2} | 200 | Yoho National Park of Canada | Parc national du Canada Yoho | 27 | \n", "| 28 | 29 | 747138.11 | 6738767165 | {914676D3-62D7-4ADD-A75C-ABF0607C6564} | 200 | Banff National Park of Canada | Parc national du Canada Banff | 28 | \n", "| 29 | 30 | 382133.80 | 3049951157 | {691C8725-89CA-47D8-9C1E-69E4A99037CE} | 100 | Riding Mountain National Park of Canada | Parc national du Canada du Mont-Riding | 29 | \n", "| 30 | 31 | 1167897.57 | 18397952413 | {590D0639-7974-47A1-9795-040E7D405338} | 200 | Auyuittuq National Park of Canada | Parc national du Canada Auyuittuq | 30 | \n", "| 31 | 32 | 872045.03 | 41167722479 | {2D2C1EB6-0D7F-4A6D-A77F-31C99E499397} | 200 | Quttinirpaaq National Park of Canada | Parc national du Canada Quttinirpaaq | 31 | \n", "| 32 | 33 | 811615.27 | 19625619733 | {09D47FA6-A6E7-4082-9813-1D4ED9AF06F3} | 200 | Ukkusiksalik National Park of Canada | Parc national du Canada du Ukkusiksalik | 32 | \n", "| 33 | 34 | 1355170.32 | 21424060989 | {5FDB4A6D-8A01-4728-9BE8-EC005103F54B} | 200 | Sirmilik National Park of Canada | Parc national du Canada Sirmilik | 33 | \n", "| 34 | 35 | 1778463.99 | 1439600089 | {978ABCB6-206D-47A4-85F7-D1A1E40C4710} | 100 | Gwaii Haanas National Park Reserve and Haida Heritage Site | Réserve de parc national et site du patrimoine haïda Gwaii Haanas | 34 | \n", "| 35 | 36 | 175418.98 | 37052584 | {F0104540-00C1-4A34-B9A7-AB5F3D0D23D7} | 100 | Gulf Islands National Park Reserve of Canada | Réserve de parc national du Canada des Îles-Gulf | 35 | \n", "| 36 | 37 | 118096.26 | 242115530 | {4245A9AA-5234-4D5B-A419-D348FEB84526} | 100 | Kouchibouguac National Park of Canada | Parc national du Canada Kouchibouguac | 36 | \n", "| 37 | 38 | 1270941.90 | 43160601764 | {323E3E60-B323-43C9-BE80-ABE621FBEE96} | 200 | Wood Buffalo National Park of Canada | Parc national du Canada Wood Buffalo | 37 | \n", "| 38 | 39 | 108077.13 | 246232998 | {CE0D9EFB-F2AB-4F8B-B03E-44A57D1593C2} | 100 | Forillon National Park of Canada | Parc national du Canada Forillon | 38 | \n", "| 39 | 40 | 187503.01 | 35663421 | {2C53DA06-A8E9-4FDE-9834-9FD67647D037} | 10 | Prince Edward Island National Park of Canada | Parc national du Canada de l'Île-du-Prince-Édouard | 39 | \n", "| 40 | 41 | 370575.29 | 5547513490 | {F5CE7051-831B-4A3F-B820-F72A574C655D} | 200 | Kluane National Park of Canada | Parc national du Canada Kluane | 40 | \n", "| 41 | 42 | 553903.31 | 9260673207 | {E02AF2ED-EC14-4FF3-B22F-9EE14C31CF5E} | 200 | Ivvavik National Park of Canada | Parc national du Canada Ivvavik | 41 | \n", "| 42 | 43 | 602842.41 | 4106240006 | {339112BF-0C82-45B3-B52D-ED44737BF2AC} | 200 | Vuntut National Park of Canada | Parc national du Canada Vuntut | 42 | \n", "| 43 | 44 | 586086.66 | 15267021368 | {2856F646-CF59-40F9-8D43-6E918F32564C} | 200 | Kluane National Park Reserve of Canada | Réserve de parc national du Canada du Canada Kluane | 43 | \n", "| 44 | 45 | 296692.13 | 406346468 | {A8D35956-D43F-401A-BF53-797FA1063444} | 100 | Terra Nova National Park of Canada | Parc national du Canada Terra-Nova | 44 | \n", "| 45 | 46 | 188812.32 | 969640516 | {CD356B3C-727E-46AC-96A8-F4A0887D9AF2} | 100 | Cape Breton Highlands National Park of Canada | Parc national du Canada des Hautes-Terres-du-Cap-Breton | 45 | \n", "| 46 | 47 | 1135818.33 | 28316219360 | {1C3DDDEF-F7BB-44F7-BFDF-9BCED3036A70} | 200 | Nahanni National Park Reserve of Canada | Réserve de parc national du Canada Nahanni | 46 | \n", "\n", "\n" ], "text/plain": [ " OBJECTID Shape_Leng Shape_Area globalid \n", "0 1 399100.63 521622672 {F5BD4565-8591-41AA-99E4-7A6F41E47A09}\n", "1 2 86148.58 187215384 {DCB2DCFE-17BC-4DFA-BE9C-FE9A4624FCA4}\n", "2 3 239712.76 1331729306 {106EEED8-DDC6-447B-9226-006507154D40}\n", "3 4 313783.71 764739692 {B2A0614F-A228-4244-9BCC-6B9B96433447}\n", "4 5 814277.44 10942527257 {B8B6DB08-8E91-4216-9D19-982A7E857920}\n", "5 6 282894.28 3827306859 {D365A5E5-1ED3-4BAF-9DDF-232A5E5EA4B9}\n", "6 7 83392.95 259236547 {571E0061-18B4-4485-9072-7138D567A683}\n", "7 8 574328.45 10924352608 {C52E5D05-7345-42F8-8377-FEE3D4A7F522}\n", "8 9 134090.44 498878286 {615C33BF-A6D6-4B69-8317-436C1D03FBD1}\n", "9 10 74399.91 211648312 {EB67B4F3-0C53-45F9-974B-3BE819D84BCA}\n", "10 11 145056.93 546631833 {F4D2C860-F422-4A44-84C5-2A4BA233DFE6}\n", "11 12 762246.23 95866783 {9290F5AD-5FD3-4A6C-A777-BDF080AD924A}\n", "12 13 361707.33 1860082905 {AFD5D06C-3740-43CF-A3C2-3BE738AED300}\n", "13 14 103792.51 31184813 {CA4BBEAD-3BA3-4B9B-82C1-1CC1EFDCAFB4}\n", "14 15 524862.27 1257873345 {A0AC46C1-8A41-4718-AFA7-A0BC689B087B}\n", "15 16 624231.91 11827725094 {E125CA8E-DCD5-4640-90A5-533BDFB19B8D}\n", "16 17 632131.31 17853403254 {FF705D51-F6F2-4CD6-84C6-E57AF223C39C}\n", "17 18 541881.60 4607855048 {60AB71D9-F267-4A78-B61F-844D5E4D0D17}\n", "18 19 149762.52 161656326 {2EFC258D-7344-4EBE-8BF4-A8B200942C98}\n", "19 20 70694.38 117124125 {3EB03934-6824-4097-B1A4-74E688A9C49F}\n", "20 21 26326.70 16075792 {81607DCE-BD39-41E3-ADEE-C371D4667D7B}\n", "21 22 422768.76 1799122547 {B577C7C8-F927-4C3A-A6B5-6DA9DCF1101A}\n", "22 23 2124790.50 9180689134 {B4CFB708-8ABD-4198-8340-231965E6901B}\n", "23 24 94220.18 14042598 {BC8D02DA-591B-4BA2-A34C-814B8C39E755}\n", "24 25 82455.80 397193668 {57653BC7-3C8F-47FA-8414-4E23474C6875}\n", "25 26 35755.34 25322482 {8BEB0C58-98D9-4946-AF21-3514ED879146}\n", "26 27 274964.28 1361036954 {092E9172-2FDC-453D-94F8-26BFF4132FB6}\n", "27 28 197069.12 1268469933 {D7164C4B-1609-4FBF-945B-2864EBBDAAD2}\n", "28 29 747138.11 6738767165 {914676D3-62D7-4ADD-A75C-ABF0607C6564}\n", "29 30 382133.80 3049951157 {691C8725-89CA-47D8-9C1E-69E4A99037CE}\n", "30 31 1167897.57 18397952413 {590D0639-7974-47A1-9795-040E7D405338}\n", "31 32 872045.03 41167722479 {2D2C1EB6-0D7F-4A6D-A77F-31C99E499397}\n", "32 33 811615.27 19625619733 {09D47FA6-A6E7-4082-9813-1D4ED9AF06F3}\n", "33 34 1355170.32 21424060989 {5FDB4A6D-8A01-4728-9BE8-EC005103F54B}\n", "34 35 1778463.99 1439600089 {978ABCB6-206D-47A4-85F7-D1A1E40C4710}\n", "35 36 175418.98 37052584 {F0104540-00C1-4A34-B9A7-AB5F3D0D23D7}\n", "36 37 118096.26 242115530 {4245A9AA-5234-4D5B-A419-D348FEB84526}\n", "37 38 1270941.90 43160601764 {323E3E60-B323-43C9-BE80-ABE621FBEE96}\n", "38 39 108077.13 246232998 {CE0D9EFB-F2AB-4F8B-B03E-44A57D1593C2}\n", "39 40 187503.01 35663421 {2C53DA06-A8E9-4FDE-9834-9FD67647D037}\n", "40 41 370575.29 5547513490 {F5CE7051-831B-4A3F-B820-F72A574C655D}\n", "41 42 553903.31 9260673207 {E02AF2ED-EC14-4FF3-B22F-9EE14C31CF5E}\n", "42 43 602842.41 4106240006 {339112BF-0C82-45B3-B52D-ED44737BF2AC}\n", "43 44 586086.66 15267021368 {2856F646-CF59-40F9-8D43-6E918F32564C}\n", "44 45 296692.13 406346468 {A8D35956-D43F-401A-BF53-797FA1063444}\n", "45 46 188812.32 969640516 {CD356B3C-727E-46AC-96A8-F4A0887D9AF2}\n", "46 47 1135818.33 28316219360 {1C3DDDEF-F7BB-44F7-BFDF-9BCED3036A70}\n", " precisionh nom_parc_e \n", "0 100 Pacific Rim National Park Reserve of Canada \n", "1 100 Elk Island National Park of Canada \n", "2 100 Glacier National Park of Canada \n", "3 100 Grasslands National Park of Canada \n", "4 100 Jasper National Park of Canada \n", "5 100 Prince Albert National Park of Canada \n", "6 100 Mount Revelstoke National Park of Canada \n", "7 200 Wapusk National Park of Canada \n", "8 100 Waterton Lakes National Park of Canada \n", "9 100 Fundy National Park of Canada \n", "10 100 La Mauricie National Park of Canada \n", "11 100 Mingan Archipelago National Park Reserve of Canada \n", "12 100 Pukaskwa National Park of Canada \n", "13 100 Sable Island National Park Reserve of Canada \n", "14 100 Saguenay-St. Lawrence Marine Park \n", "15 200 Aulavik National Park of Canada \n", "16 200 Tuktut Nogait National Park of Canada \n", "17 200 Nááts'ihch'oh National Park Reserve of Canada \n", "18 100 Bruce Peninsula National Park of Canada \n", "19 100 Fathom Five National Marine Park of Canada \n", "20 100 Point Pelee National Park of Canada \n", "21 100 Gros Morne National Park of Canada \n", "22 200 Torngat Mountains National Park of Canada \n", "23 100 Georgian Bay Islands National Park of Canada \n", "24 100 Kejimkujik National Park and National Historic Site of Canada\n", "25 100 Kejimkujik National Park Seaside \n", "26 200 Kootenay National Park of Canada \n", "27 200 Yoho National Park of Canada \n", "28 200 Banff National Park of Canada \n", "29 100 Riding Mountain National Park of Canada \n", "30 200 Auyuittuq National Park of Canada \n", "31 200 Quttinirpaaq National Park of Canada \n", "32 200 Ukkusiksalik National Park of Canada \n", "33 200 Sirmilik National Park of Canada \n", "34 100 Gwaii Haanas National Park Reserve and Haida Heritage Site \n", "35 100 Gulf Islands National Park Reserve of Canada \n", "36 100 Kouchibouguac National Park of Canada \n", "37 200 Wood Buffalo National Park of Canada \n", "38 100 Forillon National Park of Canada \n", "39 10 Prince Edward Island National Park of Canada \n", "40 200 Kluane National Park of Canada \n", "41 200 Ivvavik National Park of Canada \n", "42 200 Vuntut National Park of Canada \n", "43 200 Kluane National Park Reserve of Canada \n", "44 100 Terra Nova National Park of Canada \n", "45 100 Cape Breton Highlands National Park of Canada \n", "46 200 Nahanni National Park Reserve of Canada \n", " nom_parc_f \n", "0 Réserve de parc national du Canada Pacific Rim \n", "1 Parc national du Canada Elk Island \n", "2 Parc national du Canada des Glaciers \n", "3 Parc national du Canada des Prairies \n", "4 Parc national du Canada Jasper \n", "5 Parc national du Canada de Prince Albert \n", "6 Parc national du Canada du Mont-Revelstoke \n", "7 Parc national du Canada Wapusk \n", "8 Parc national du Canada des Lacs-Waterton \n", "9 Parc national du Canada Fundy \n", "10 Parc national du Canada de la Mauricie \n", "11 Réserve de parc national du Canada de l'Archipel-de-Mingan \n", "12 Parc national du Canada Pukaskwa \n", "13 Réserve de parc national du Canada de l'Île-de-Sable \n", "14 Parc marin du Saguenay-Saint-Laurent \n", "15 Parc national du Canada Aulavik \n", "16 Parc national du Canada Tuktut Nogait \n", "17 Réserve de parc national du Canada de Nááts'ihch'oh \n", "18 Parc national du Canada de la Péninsule-Bruce \n", "19 Parc marin national du Canada Fathom Five \n", "20 Parc national du Canada de la Pointe-Pelée \n", "21 Parc national du Canada du Gros-Morne \n", "22 Parc national du Canada Monts-Torngat \n", "23 Parc national du Canada des Îles-de-la-Baie-Georgienne \n", "24 Parc national et lieu historique national Kejimkujik du Canada \n", "25 Parc national Kejimkujik Bord de mer \n", "26 Parc national du Canada Kootenay \n", "27 Parc national du Canada Yoho \n", "28 Parc national du Canada Banff \n", "29 Parc national du Canada du Mont-Riding \n", "30 Parc national du Canada Auyuittuq \n", "31 Parc national du Canada Quttinirpaaq \n", "32 Parc national du Canada du Ukkusiksalik \n", "33 Parc national du Canada Sirmilik \n", "34 Réserve de parc national et site du patrimoine haïda Gwaii Haanas\n", "35 Réserve de parc national du Canada des Îles-Gulf \n", "36 Parc national du Canada Kouchibouguac \n", "37 Parc national du Canada Wood Buffalo \n", "38 Parc national du Canada Forillon \n", "39 Parc national du Canada de l'Île-du-Prince-Édouard \n", "40 Parc national du Canada Kluane \n", "41 Parc national du Canada Ivvavik \n", "42 Parc national du Canada Vuntut \n", "43 Réserve de parc national du Canada du Canada Kluane \n", "44 Parc national du Canada Terra-Nova \n", "45 Parc national du Canada des Hautes-Terres-du-Cap-Breton \n", "46 Réserve de parc national du Canada Nahanni \n", " rmapshaperid\n", "0 0 \n", "1 1 \n", "2 2 \n", "3 3 \n", "4 4 \n", "5 5 \n", "6 6 \n", "7 7 \n", "8 8 \n", "9 9 \n", "10 10 \n", "11 11 \n", "12 12 \n", "13 13 \n", "14 14 \n", "15 15 \n", "16 16 \n", "17 17 \n", "18 18 \n", "19 19 \n", "20 20 \n", "21 21 \n", "22 22 \n", "23 23 \n", "24 24 \n", "25 25 \n", "26 26 \n", "27 27 \n", "28 28 \n", "29 29 \n", "30 30 \n", "31 31 \n", "32 32 \n", "33 33 \n", "34 34 \n", "35 35 \n", "36 36 \n", "37 37 \n", "38 38 \n", "39 39 \n", "40 40 \n", "41 41 \n", "42 42 \n", "43 43 \n", "44 44 \n", "45 45 \n", "46 46 " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "limites_parcs_nationaux@data" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Affichage des polygones des parcs nationaux contenux dans les limites du polygone de l'Alberta" ] }, { "cell_type": "code", "execution_count": 77, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAACJVBMVEUAAAAAAQABAgECAwID\nBgMEBgQEBwQFCAUFCQUGCgYHCwcHDAcJDgkJDwkKEAoKEQoLEgsMEwwMFAwNDQ0NFg0OFw4P\nDw8PGA8PGQ8QGxARGxERHBESHRISHxITHxMUIRQUIhQVIxUWFhYXJhcXJxcYKBgZKRkZKhka\nKxobGxscLhwdMB0eMR4fMx8gNCAgNSAhNyEiNyIjOSMjOiMkJCQkOyQkPCQlPSUmQCYnJyco\nQigoQygqKiorRyssSCwtSy0vTi8wUDAxUDEyMjIyUzIzVDM1NTU2WTY2WjY4XTg5Xjk5Xzk6\nYDo7YTs7Yjs9ZT0+Zz5AQEBAakBBQUFBbEFCbUJDbkNDb0NEcERGdEZHdUdHdkdJeUlKekpK\ne0pMfkxOgU5Pg09RhVFSiFJTU1NTiVNUilRUi1RVjFVVjVVWjVZYkVhblltbl1tcmFxdml1e\nnF5fnV9fnl9hoGFhoWFiomJjo2NjpGNkpmRlpmVmqGZmqWZnqmdoq2horGhpaWlprWlprmlq\nrmpqr2pqsGprsWtrsmtssmxss2xutm5vuG9wuXBwunBxu3FyvHJyvXJzvnNzv3N0v3R0wHR2\nw3Z3d3d3xHd4xnh4x3h5x3l5yHl6yXp6ynp7y3t8zHx8zXyIiIiLi4uqqqqrq6uurq6vr6+9\nvb3ExMTMzMzj4+Pm5ubs7Ozt7e3u7u7v7+/z8/P29vb39/f5+fn6+vr8/Pz+/v7///99BYfo\nAAAACXBIWXMAABJ0AAASdAHeZh94AAASEElEQVR4nO3d+5tdV13H8Qnhog00sVBLIVBqrVSw\nCgKtUqGgXKWAcqulBQLYC6BFKkWhsVBAsBRK047NpZ3ZCmgVVLRK/j5nkpnJXM6ZnMvn7L3W\n3q/XL5NzJifPyvPs9/NdZ+3nzCydBea21PUCoA+EBAFCggAhQYCQIEBIECAkCBASBAgJAoQE\nAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQ\nIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQ\nEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQB\nQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAg\nQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBAS\nBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASMVpFqrr/11fCak4QqqR\nkIqjoxoJqTgyqpGQiiOiGgmpODqqkZCKs5iQZLVYQiqOkGokpOK0EZKW0oRUnJZCUlKUkIrT\nVkhSShJScdoLSUo5QipOmyEpKUVIxWk1JCWFCKk47YakpQwhFUdINRJScVoPSVQBQipOpyEp\naUZCKk63ISlpNkIqjpBqJKTidBySkmYipOLoqEZCKo6QaiSk4nQckpJmIqTidB2SkmYhpPJ0\nHZKSZiCk8nTdkZJmIKTydJ2RkGYgpPJ0nZGOZiCk8nTdkZBmIKTydN2RkGYgpPLoqEJCKk/X\nISlpBkIqj44qJKTyCKlCQqqfkAogpPoJqQBCqt/eFO64/NCNP2iae45esv5FSC0QUv32lPDQ\n87/52A3vaH74vM8/ftPbhNQKIdVvTwnHXtc0n31t8+UjTXPnq4XUCiHVb0QLTx+/+f3N8uFP\nP3zje6cPSUkzEFL9RqRw9wsue6Rpbl967qWPCqkVQqrfqBZOf/Ka5otXfeP0J64RUiuEVL89\nIdxxrGmeeM5TH3x306wcNJJaIaT67eng7qsfXfnIa5ovvPyrT378ZatCaoOQ6reng9UPXHb4\nDd9pmtteeuj6B1TUCiHVb4aRI6Q0IdVPRwUQUv2yIWlpJkKqXzwkJU1PSPUTUgGEVD8dFUBI\n9RNSAYRUPyEVQEj1E1IBhFQ/IRVASPWLh6Sk6QmpfjoqgJDqJ6QCCKl+QiqAkOoXD0lJ0xNS\n/fIhSWlqQqrfIkJS0pSEVL+FhKSk6QipfosJSUtTEVL9hFQAIdVPRwUQUv10VAAh1U9HBRBS\n/YRUACHVT0gFEFL9hFQAIdVvMSEpaSpCqp+QCiCk+gmpAEKq34JCUtI0hFQ/IRVASPVbVEhK\nmoKQ6rfQkNQ0GSHVb2EhbdXERQmpfkIqgJDqJ6MCCKl+QiqAkOqnowIIqX4yKoCQ6iejAgip\nfjIqgJDqJ6MCCKl+C4hISdMSUv3iGTE9IdUvGpKSZiOk+umoAEKqXzQkJc1GSPUTUgGEVD8h\nFUBI9RNSAYRUPyEVQEj1y4akpJkIqX5CKoCQ6hcOSUmzEFL90iEpaQZCqp+QCiCk+gmpAEKq\nXzwkKU1PSPUTUgGEVD8hFUBI9RNSAYRUPyEVQEj1E1IBhFQ/IRVASBAgJAgQEgQICQKEBAFC\nggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBA\nSBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIE\nCAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKC\nACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBI\nECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQI\nCQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIA\nIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQ\nICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJ\nAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAh\nQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAg\nJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkC\nhAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFB\ngJAgQEgQICQIENIQ/PiNh65/sOtF9JuQBuCnL35m5aOv+kXXy+g1IQ3ATz58tjlz4D+7Xkav\nCWkAfv6zprn3ShNpkYQ0BM3q547c1fUi+k1IA9CcuOmVX+p6ET0npAE4c907zzRdL6LnhDQA\nd1576tSp//EeaZGENAC3LK37r66X0WtCGoBmXdeL6DkhDYGQFk5IECAkCBDSaM2GrYebf+hq\nQZRNSCM1F2w+3np+67tdLpDCCGmUZqc9TwiJXYQ02shypMQ4QhpHSUxBSOOYSUxBSGNN2pGc\nENJ+JgrJYGKdkMabbHOnJM4KaT/ThKSkgRPSeFOFpKRhE9I+hMSkhLQfITEhIe1nWyRSYj9C\n2tf2RpTEeELa145ElMRYQrqoZtunkqTEaEKagJK4GCFNQ0iMIaTpKImRhDQVI4nRhDSVCUJS\n0iAJaRqTdKSkQRLSNCYLSUoDJKRpCYkRhDQtHTGCkKZkIDGKkKYlJEYQ0tQmmUlaGhohTW1v\nNb+1tLR0g5IGTUhT2xvSS/5+efmkkAZNSFPb84HZU7+0anM3dEKazfZkHjx07aHf+ZaSBk1I\ns9lezF+//ltn/vg3DKVBE9KsdjZz8sBjShoyIc3qQi733ds0pw8ujwhJSYMhpBltq+UvL/27\nlT95/aiOlDQYQprR9lpu/ZVLf+97Qho0Ic1qdDhKGighzU5JbBHS7CaeSWLqPyHNbsKQNseS\nnPpMSHOYKiSbvF4T0hymGEmbY4meEtLspuxISH0mpJlN25GQ+kxIs9MRW4Q0Ox2xRUhzEBKb\nhDQPIbFBSPPQERuENA8TiQ1CmouQOE9I8xES5whpPpuR6GjghDSXC5HoaNiENJftlQhpyIQ0\nl12V6GiwhDSfXYkIaaiEFKWjoRJSlJCGSkhRQhoqISXpaLCEFOSoYbiElKOjARNSjo4GTEg5\n4yaSlAZASEFCGi4hRSlpqISUpqRBElKakAZJSHFCGiIhLY6UBkRIC2QqDYeQFsj+bjiEtEDe\nKQ2HkBbJocNgCGmhlDQUQlooIQ2FkBbKHaWhENJiKWkghLRoQhoEIS2ckIZASIsnpAEQUgv2\nLUlZvSCkNuzXkSHVC0JqwT7vkuz2ekJILbhYSFqqn5BasbukjXKcQPSGkNqw+17SZk1C6g0h\nteJCL6OaUlL9hNSKC7XoqJ+E1IoLtYydSGKqmpDacWHsjCnJXKqbkNqxrRMh9ZGQWrfPSJJS\ntYTUuv1CUlKthNS6fUNSUqWE1D4l9ZCQWreVzJ1HD13/YNN85ujhm08qqXJCate2YB755ftW\nPvqq1Qde9LWTb32fkConpFZt38Pd89tNc+bA4x97e9N8+3Kbu8oJqVXbe1l+rGnuvXL1Q2sh\nfffA00qqm5BatfNcYfVzR+5q7n/hV068ZemEkOompFbt6OjETa/80tqXT131ilsPruiobkJq\n1faOzlz3zjPrO7zvN839VxtIlRNSi3Zu7O689tSa1eOXPPTkm27b+a2uF8rUhBQzweW/o5Zb\nltY90dx+5IpbVnRUOSGFbL/+x8ewcQ/2nqOX3PiDZqyWlkyQkEK2FTA2hs17sD983ucfv+lt\n40NSUn2ElDJBC5v3YL98ZO0d0quF1CdCipkgho17sMuHP/3wje8VUp8IKefiNWzcg21uX3ru\npY8KqU+ElHPxHjbuwX7xqm+c/sQ1QuoTIQWdr+COyw9tHcnt+gsb92A/+O6mWTk4fiR1snjm\nIqSk9Qoeev43H7vhHaOb+NH5e7BfePlXn/z4y1ZNpB4RUtRaBMde1zSffe3oIjbvwd720kPX\nP7DP1k5K1RFS1LkInj5+8/vH9LBvPHtKElQ9hBR1roG7X3DZI+dbGNHB5CmZTDURUtLG1X/6\nk9eMjUBI/SSkqLVL/45jTfPEc56KhNTy6pmdkKLWLv67r3505SOvGd+BjHpJSFlNs/qByw6/\n4Tv7lLCZiZD6REgRo39C/ti/3Oz6hX1Cqp6Q5vXP5+4N3bf5cPKQLlpSO+snQkhzas4sLy//\n7XVPbT2+eApb31JRfwhpTueu+jcfn2pvt+eXmquoekKa0/p1/xd/uHXxX0hhoiCU1BdCmtPa\nVb/y6/+wZyDt84HzXa8WUi8IaU5rV/1f/e6uS3/rKGGil19oR0n1EtKc1i769xwb/XHYSV8/\nIqToEmmBkOazdtGvXvHvs776x2/81T99sPmzcwfoz5zdWRVVEdJ81i76f/vNWV/80xc/83//\n9Aer5w/Q/3fjnxNSlYTUoZ98+OzZZw88vl7Om49vPCekOgmpQz//2dmz/3rl6vkD9I3nzKQ6\nCalTv/iXX7tr4wB98ykhVUlIXfqPP//99d+QtOMAXUhVElKHnn3X1589V817ju38CfxCqo6Q\nOvSjPzr307ma1Sse3vaskGokpA7948ZP5/qbl2zrRkdVElLHdu/k7OzqJKSu7SpHR3USUud2\np6OkGgmpc6M3dx0uiBkIqXtGUg8IqQB7DhiEVB0hFWDvZk5ItRFSAUbs5ZRUGSEVYMTxgo4q\nI6TuuQfbA0LqnE8g9YGQuiekHhBS94TUA0Lqnr1dDwipeyrqASEVQEf1ExIECAkChAQBQoIA\nIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQ\nICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJ\nAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAh\nQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChNQDTdcLQEi1a87rehmD\nJ6SqnYvoM0cP3/zfXa9k6IRUtfWOHnjR106+9X2GUreEVLe1kD729qb59uU66paQ6rYW0ofW\nQvrugae7XsnACalyTXP/C79y4i1LJ4ykTgmpdk3zqatecevBFYd3nRJS9Za/vzaVrnYK3i0h\nVe/4JQ89+abbhNQtIdWuaW4/csUtK0LqlpAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAg\nJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkC\nhAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFB\ngJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAk\nCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGAkCBASBAgJAgQEgQICQKE\nBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQIEBIECAkChAQBQoIAIUGA\nkCBASBAgJAgQEgQICQKEBAFCggAhQYCQIEBIECAkCBASBAgJAoQEAUKCACFBgJAgQEgQICQI\nEBIE/D9ik+kq8/P6OgAAAABJRU5ErkJggg==", "text/plain": [ "plot without title" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "map(limites_parcs_nationaux,\n", " fill= TRUE,\n", " col='palegreen3',\n", " namefield='nom_parc_f',\n", " border=FALSE,\n", " xlim=c(bbox_prov_alberta[1,1],bbox_prov_alberta[1,2]),\n", " ylim=c(bbox_prov_alberta[2,1],bbox_prov_alberta[2,2]))\n", "invisible(text(coordinates(limites_parcs_nationaux),\n", " labels=as.character(limites_parcs_nationaux$OBJECTID),\n", " cex=.6))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Combinaison" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "On combine les données du polygone de la province avec les données des parcs nationaux" ] }, { "cell_type": "code", "execution_count": 78, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAIAAAByhViMAAAACXBIWXMAABJ0AAASdAHeZh94\nAAAgAElEQVR4nOzdeXwN1//H8XNvVhEhEbHFTiwRu4Roi7aWohqKVqnY2ypavm3RBm1ja1FK\n1U6Rqgq1pvadWFvUnpAgSmwRhOx3fn/M9zu/6yYiiSQzd7yeD3/MnXvuzGcmyb1vZ+aca5Ak\nSQAAAMD6GdUuAAAAAHmDYAcAAKATBDsAAACdINgBAADoBMEOAABAJwh2AAAAOkGwAwAA0AmC\nHQAAgE4Q7AAAAHSCYAcAAKATBDsAAACdINgBAADoBMEOAABAJwh2AAAAOkGwAwAA0AmCHQAA\ngE4Q7AAAAHSCYAcAAKATBDsAAACdINgBAADoBMEOAABAJwh2AAAAOkGwAwAA0AmCHQAAgE4Q\n7AAAAHSCYAcAAKATBDsAAACdINgBAADoBMEOAABAJwh2AAAAOkGwAwAA0AmCHQAAgE4Q7AAA\nAHSCYAcAAKATBDsAAACdINgBAADoBMEOAABAJwh2AAAAOkGwAwAA0AmCHQAAgE4Q7AAAAHSC\nYAcAAKATBDsAAACdINgBAADoBMEOAABAJwh2AAAAOkGwAwAA0AmCHQAAgE4Q7AAAAHSCYAcA\nAKATBDsAAACdINgBAADoBMEOAABAJwh2AAAAOkGwAwAA0AmCHQAAgE4Q7AAAAHSCYAcAAKAT\nBDsAAACdINjhucTGxhr+JygoSO1y8FxCQkKUn2ZaWpra5WjF5s2bldNy586d/NjF7t27BwwY\n4OXlVaRIEVdXVz8/v169eu3fvz8/9pW3CuDkZOq3335T9lukSJHExMSMbebMmSM3cHZ2VlZ+\n+OGH8sqXXnqpwKoFChLBDgBU8+DBg4CAgJYtWy5YsCAyMjIhISE+Pv7IkSPLli17+eWXO3To\nEBMTo3aNWrRixQplOSEhYcuWLSoWA2gKwQ544aSlpSm9HeYfkChgjx8/btGixbp1657WICws\nrGXLlvHx8QVZlfbdu3dv06ZN5mtWrlypVjFZ428NBY9gB+C/mjZtuvh/jEbeHPLdiBEjjh8/\nLi87Ozt///33Bw4ciIqKWrduna+vr7z+0qVLX331lXo1atGaNWtSU1PN16xfv/7x48dq1QNo\niq3aBQDQiipVqlSpUkXtKl4U169fnzNnjrzs4eGxc+dOb29v+WGlSpVef/11f3//kydPCiGW\nLFny3Xffmd8o9oL77bff5IVmzZodOHBACPHo0aPNmzd37txZ1boATeA/5Sho27dvDwwMrFev\nXqFChUqVKtWpU6dt27ZZtDGZTGvWrGnTpk3lypUdHR0rV67cqlWrkJCQZ97RP2HChCZNmjRp\n0mTQoEEJCQmffvpp2bJlHR0dq1evPnr06IcPH1q0T01NXbx4cfPmzStUqODg4FCuXDl/f/8Z\nM2YkJCQobb7++mv5SkrVqlWFEKtWrapXr56Xl5fSIDIycujQoa+88oqbm1upUqVatGgxffp0\nix6FXB/RggUL5CPq0qVLenr6tGnTvL29CxUqVLFixS5dukREROT0iLp06WJnZ6e07969u8Fg\nOHTokHjKzeayO3fufPPNN+3atatQoYKLi4uvr2/v3r23bt0qSVK+ViuE6N+/v1xVixYtLPYl\nr7e1/e9/UA8dOmQ0GpVb45XaPvnkE+VymJKlnub06dOdOnUqXry4k5OTn5/fmjVrMm2W/aqe\n5tdff1V++iNGjFBSnczJyWnkyJHy8qNHj06cOCEvZ+eM5cdPIUcnJ/vbnDdv3tdff/3111/P\nmzcv69OluHnz5s6dO5XzVr9+fXk5F1djb9261b9//1KlSjk6Onp7e0+cODE5OTljs2e+ZT3t\nLSKLvzWRk9MO5IwEPIcbN24ov0tfffVV1o1NJtMnn3yS6e/h4MGDTSaT3Cw9PT0gICDTZg0a\nNHj48GEWu+jXr5/c0t/fv1mzZhYvr169+o0bN5TGSUlJTZs2zXRHPj4+CQkJcrOxY8fKK6tU\nqfLLL7/IyxUqVJCf/f333+3t7TNuoWbNmrdv337+I1L2Xrly5R49eli83MHB4dSpUzk6orff\nfjvjswcPHpQkafbs2fLDwoULm9ewZ8+e4sWLZ7rZnj17pqSk5F+15j/T5s2bm1c1f/58eb2N\njY2ycvjw4cpG1q1bJ0lSdHS08uHavHnz9PT0LM72qlWrChcubFFM165dlWXlZ5qjqjLVoUMH\n5bRk+juQlJR09H/u3r2b/TOWHz+FHJ2c7G+zYcOG8vqGDRtmfboUM2fOlF/i4uKSlJQ0YcIE\n+aGTk5P5lqWn/D5/8MEH8sratWtXq1bNojxfX9/79+8rjbP5lvW0t4gs/tayf4qAnCLY4bnk\nKNitXbvW/A20d+/eDRo0UNb8/vvvcrMZM2YoK6tUqRIQENC4cWNlTd++fbPYhfJxq6hQoYKN\njY3ysGfPnkrjL7/8Ulnv5eX1yiuvuLu7K2u++eYbuZnyru3h4eHq6qpsVpKkU6dOKanOaDTW\nrVu3bt26yhY6dOggv/U/zxEpe1eUKlXK/Ijatm2boyM6e/bsrl27lJVjx47dt2+f/GGW6Qfh\nlStXihUrZn4+Gzdu7ODgoKwZMWJE/lUr5TBCPXr0SLma7O3tnZaW1qtXL/lhoUKFIiMjszjV\nUVFR5hm9VKlSRYoUsTicPAx2SqqoXr161i1zesby46eQo5OTzW1KuQp2yn/Y3n//fUmSIiMj\nlS2Hhoaat8w62Cks3iIGDRqkNM7mW9bT3iKy+FvL/ikCcopgh+eS/WCXnJxcuXJlueXnn38u\nJx6TyfTFF1/IKytVqpSUlCRJUrt27eQ17733nvJy5b/pxYsXV/6jnJF5sKtXr96VK1ckSbp3\n756yTYPBcPr0abmxcvFrxIgRSj3vvvuuvPKNN96Qm5l/TJYpU2bRokUnT548e/asJEnvvfee\nvN7d3V3pBVH+yy6EuHDhwnMekfne27dvf/PmTUmSbt++3ahRI3mli4uL8vJsHpH5ZeLffvtN\n2VfWH4QGg2HJkiXyyhs3bjRp0kReb29vf/369fyrNqcRavfu3UoNn332mcFgkJenTp36tJMs\nU3bt7OwcFhZmMpnS0tImTZokzORhsFM6QVu3bp11S0UufmPz6qeQo5OTzW1KOQ92V65cUXa3\nfv16eaVyNbZr167mjZ8Z7GrXrh0VFSVJ0r1799544w15pb29/dWrV6WcvGVl8RbxtL+17J8i\nIKcYPIECcv78+aioKCGEwWAYOXKk/HFrMBhGjRo1derU9PT06OjogwcPtmjR4urVq/JLDh48\nGBoa2rZt2yJFigQGBnp6esrrU1NTM736aWHx4sXly5cXQhQrVmzJkiUVKlR4/PixJEkbNmzw\n9vY2mUxDhgyRJEkI8dZbb8n1SJKk3GRz69atjNtcuXKl+RXeffv2yQuDBw+uXbu2vNyrV69p\n06ZdvHhRCHHgwAEvL688OSKj0bhs2TK5P8Dd3X3o0KFyX9SDBw/u3r3r7u6euyN6JuV+pp49\neyq9X6VKlVq0aJG3t7ckSSkpKeHh4RZXndSqVgjRvHnzQYMG/fzzz0KIKVOmyCt9fX2fdk1N\nlpaWtnr1ann5s88+k7O4jY3NiBEjtm7dqpyE/GDeXZSFXJyxvPop5Ojk5KjOY8eOZfc0CSGE\n+P333+UFFxeX1q1by8tdu3aVBxdv3Ljx0aNHGa8XP83ixYsrVaokhChWrNjSpUsrVaqUkJCQ\nkpLy66+/jhw5MvtvWebbtHiLyFT+/fIDglGxKDAXLlyQFySzvgoLkZGRLVq0aNSo0enTp4UQ\n0dHR3bp1s7W19ff3b926datWrRo1apTNaTjKlStXr1495aG7u3vTpk137NghhJBvHjcajfL/\n3WNiYrZt23bixIkzZ86cOnXKvA/SgouLi7+/v/IwMTFRmTzWfBZ7g8Gg3Oouy6sjUq7yCCFK\nlSqlLMv34OfiiJ4pOTlZ/mwTQnTs2NH8qZo1a1atWlW+EHb+/HktVKuYNGnSxo0blTxtZ2e3\naNGirPPTpUuXlM6VTp06mT8VEBCQH8HO1dX17t27QohsTkGcizOWVz+FHJ2cfP3JKlPBdezY\nUbkfoGvXrvKVzcTExLCwsG7dumVnU2XKlFG6MIUQ7u7ur7zyyp9//in+92aV/bcs5aHFW8TT\n5OspAgh2KCDKp2wW5Mg1derUhISE1atXy/+jTUtL27t37969e4OCgry8vGbNmvX6668/c1Ol\nS5e2WFOmTBl54dq1a/LCpUuXBg4caPGx5OTk9LQJsdzd3ZXrehZHZP7xmVGeHJHFKEvzShQ5\nPaJn+vfff9PT0+XlsmXLWjzr6ekpBzsl/KlbraJIkSKDBg1SRpW2b9/eYsxpRpcvX1aWLX55\nMh54nvDy8pK7dS9fvpySkpKxyzY5Ofno0aPycuXKlcuUKZPTM5ZXP4Wcnpx8+slGRET8/fff\n8nJISEhISEjGNitXrsxmsDOPuTLlWOTf5+y/ZSks3iKykH+//ADBDgVEuezo7OxsMWu8Qv7M\ncHNzCw0NjYmJWbVq1caNG/fu3avMChEREdG+ffvz58/LF1CyEBsba7FG+d+wfH02JSUlICBA\n7kgrXrx4165dmzdv3qRJk4ULF44bNy7TbVq8ZZt/MGT9LZl5ckTPlIsjeqYyZcoYjUaTySSE\nuH79usWzyhrlh5uv1UpPTq2SxUwxd+7cmTp1qvJw/fr1hw8f9vPzy6Ie84ASGxtbokQJ84dZ\nvDD7VVl4+eWX5f6hhISE5cuX9+7d26LBnj172rRpIy+HhoZ27Ngxz3++2fwp5Ojk5MfvoSw7\n39wQFhaWkJCQnTn/Mlau/D7Lf9rZf8tSZDPV5d8pAgTBDgVGmfgtISGhZs2aT7u0ERcXd/bs\nWXn5k08+GTZs2IMHD7Zu3bp48WL5UzAlJWXz5s0fffRR1ru7evXqqVOnfHx85Id37tw5ePCg\nvFyzZk0hxKlTp+Q3ViHEgQMHqlevLi8re3+mokWLlixZ8ubNm0KI/fv3K3f8SJI0evRoOeq9\n8847devWzZMjeqbnP6KMHB0dK1WqdOnSJSHExo0bzW+kO3/+vDIgUbm/MD+qVS6hWtx4dOrU\nqadtfMiQIbdv31YemkymPn36/P33346Ojk97SZUqVQwGg5zSNm7cqPzmCCHCwsIyts9FVRZ6\n9OgRFBQkd4iOGzeuTZs25inBZDL99NNPysPGjRvnx883m9vM0cnJjzqFEJIkKfMSZyEpKWnj\nxo3KKIQsXL9+/fjx48rAi7i4OOWW2Vq1aolsv2XlQj6dIkDGBMUoIDVr1lQ+t77//nuln+Pi\nxYu1a9euUaNGjRo1jh07du/evZf/R55rwMXFpUuXLitWrFCuVWVzCv6+ffvKvXQPHjzo3bv3\no0ePhBA2Njbt27cXQph/8D948EBe2LZtm/kEB8/UsmVLeWHmzJnKRaJFixaNHz9+7ty5c+fO\ndXd3z8MjylrujujevXtZb1Y5xqVLly5fvlxevnnzZt++feWePDs7u+zcV5TrapW+zPPnz4eG\nhgohJElavXq1Mv7Uwtq1a5WunSlTpsiT2J07d+7bb7/Nop5ChQopE8tNmjRJvh1TkqQZM2bI\n+fs5q8qoXLlyytDaS5cuNW7c+Ndff42KioqLiwsPD+/cufOGDRvkZ9u0aVOhQoU8+Y21kM1t\n5ujk5KjO7E9Q/M8//yi3cq5cuTI1AyWHZX+m4n79+sm9dA8fPuzTp49SrdxRms23rGzuS5j9\nreXHjxL4f/k/8BZ6lp27fcuWLSs3Xrp0qbKyWbNmAwYMaNWqlTKErVmzZiaTyWQy1ahRQ15j\nY2NTt27dLl26tGzZUrn6YzQa5RkKMmUxj53RaKxatar55O8ffvih3NL8jnVbW9vGjRubT0En\nhGjQoIHc0nz2UYvdRUdHK51A8kaU6RuEEK1bt37+I3ra3s3nvpdnXc7+EZlMJicnJ3llvXr1\nZs6cKW8h0+khoqKiXFxclI1UrVrV399febl4yjx2eVit+fQlQoiSJUsWLVrUfI35xCJ3794t\nWbKkvL5p06bmE8za2NgcO3bsaedZkiSLIS8WIw9kyoweOarqaR4+fPjMzk5XV9eIiIgcnbH8\n+Clk/+Rkf5tSTqY7GTFihNzSzs7uwYMHGRsMGzZMbuDg4JDFvIwW89jJbxHmNzh2795daZyd\nt6wsTrj0lL+1HJ0iIKcIdnguOQp26enpPXv2zLRNkyZNlI/MCxcumN/EY+Hnn3/Ooh4l2NWr\nV69ixYoWr23dunVcXJzSOONdTZUqVVJm9ChWrNijR4+kLN+1JUl62nDLunXryt8W8JxHlP0P\n6ewfkSRJyoVjWdbfPLFz586njQ7p3bt3ampqflfbpUsXi5b29vbK5TbzCGX+C7Z//35Jkm7f\nvq0E0zp16iQnJ2dxti0mZpO1bdtWWVZ+S3NUVRbu3r3bqlWrTM+tEKJp06YnTpxQGj/nb+xz\n/hSyf3Kyv81sBjuTyaT8Obdq1SrTNuZH9+uvv0rPCnZVqlSRvwHMnL+/f2xsrNI4m29ZWb9F\nZPq3lv1TBOQUwQ7PJUfBTpIkk8m0atWqbt261axZ08HBoXz58q1atVq+fLnFx+39+/d//PHH\nl19+uVKlSo6Ojq6urj4+Pn369Dl+/HjW9ZhPGxsTE9O9e3c3NzdnZ+eWLVtOmTIlLS3NvHFi\nYuLEiRNr1Kjh5ORUv379YcOGxcfHHz16tOH/LFiwQHrWu7YkSSdPnvzggw+aNm1atGhR+bti\n58yZYx53nueIcvQhnc0jkiTpypUrnTt3dnd3t7GxcXFxOXr0qPT0YCdJUmxs7OjRo9u0aVOu\nXDlnZ+dGjRoFBgZu27bNYmrlfKo2NTX1+++/r1+/vpOTk6ura8eOHY8fP75mzRq5ma+vr9xM\nuXYphOjUqZOy9/Hjxyvrx44dm/UJ37p1a6tWrYoWLerm5tamTZs//vjjwIEDSkn37t1TWmaz\nqmcymUybNm16//33K1asWKhQIVdX10aNGr377rsbNmywOL3P+Rv7nD+F7J+c7G8zm8FOuUFW\nCPHjjz9m2iYpKUnpS+vYsaP0lN/n8ePHyzX07dv3zp07AwYM8Pb2LlKkyMsvvzxhwgTz78dT\nfjrPfMvK+i0i07+1HJ12IEf+ez8soA/9+/dfuHChEKJ58+YWF8sAANA9Bk8AAADoBMEOAABA\nJwh2AAAAOkGwAwAA0AkGTwAAAOgEPXYAAAA6QbADAADQCYIdAACAThDsAAAAdIJgBwAAoBME\nOwAAAJ0g2AEAAOgEwQ4AAEAnCHYAAAA6QbADAADQCYIdAACAThDsAAAAdIJgBwAAoBMEOwAA\nAJ0g2AEAAOgEwQ4AAEAnCHYAAAA6QbADAADQCYIdAACAThDsAAAAdIJgBwAAoBMEOwAAAJ0g\n2AEAAOgEwQ4AAEAnCHYAAAA6QbADAADQCYIdAACAThDsAAAAdIJgBwAAoBMEO2iXyWS6c+fO\n9evX09PT1a4FAAArQLCDtkiStHv37h49epQtW9bBwaFEiRJly5a1t7f39PR85513du/eLUmS\n2jUCAKBRBj4moR2JiYmdOnXasmWLEKJ06dLlypUrXry4ECIuLu7ff/+9du2aEKJ9+/arVq1y\ndHRUuVYAALTHVu0CgP83ceLELVu2tGnTZuLEifXq1TMYDMpTkiSdPXt23LhxK1asmDx58ujR\no1WsEwAAbaLHDhrSpEmT+Pj4U6dO2dnZZdrAZDK1aNEiLS0tPDw8z/ceExNz69atPN8sgPzm\n4eFRrlw5tasANIEeO2jI2bNnO3Xq9LRUJ4QwGo0vvfTSTz/9lB97b9eu3enTp/Njy8hb70x/\nR+0SNOr3T39XuwR11K5V89SZs2pXAWgCwQ4a4u3tfeTIkbS0NFvbzH8zJUk6dOiQt7d3fuw9\nJSVl7sxpA/sG5sfGkYdCrq1VuwSNWjFgttolqGDeoiVTZ+TLf/YAa8SoWGhImzZtzp8//9Zb\nb506dSrjTQIRERG9evXatWtX27ZtVSkP0LKengFqlwBAffTYQUNGjBhx5MiRP//8888//yxX\nrlz58uXd3NwMBsO9e/euXbsWHR0thGjbtu2IESPUrhRQn5zkQq6tJdIBUBDsoCGFChUKCwvb\nvXv3vHnzdu/effDgQZPJJIQwGo0eHh5du3b94IMPXn31VfPRssCLSQlzpDoA5gh20BaDwdCy\nZcuWLVsKIdLT02/fvi1JkoeHh42NzfNs9vr16zdu3Mi6zcOHD6Oiozf8ufn+gwfPsy/kt38r\n31G7BJWVjXIPESvVrkITirq4XImJSUpKvnXrloeHh9rlAOpjuhNYgZs3b8bExFSvXr1IkSK5\n20Lz5s337t2bt1VBLYyKzb4XZ5zsvHnzBgwYoHYVgPoIdtCWy5cvf/311/Xq1fv000+FEOHh\n4QMHDjxz5oz8bJs2bebMmVOxYsWcbjY1NTUhISHrNiVLlmzfttXaDX+uXr6kZfOXc147Csja\ne9vULsFqBLi2UruEfLR77/7O3XuN/zpowS/Lzp47zxfSAIJLsdCUyMjIJk2axMXFyROanD17\ntmXLlqmpqa1bt65cufK5c+e2bNni6+t79uxZd3f3HG3Zzs7O1dU16zYGg8G5sLMQwtnZ2bVY\nsVwfBfKbXcJTZzqEhbCE3fKCLm/Fc3Z2FkK4Fy9uZ2dLqgNkTHcCDRk1alR8fPyaNWs+++wz\nIcSXX36Zlpa2ffv2LVu2zJ49e/fu3aGhobdv3x4zZozalQJWhsn/gBcEPXbQkH379nXo0CEg\n4L9dC0ePHm3btu2rr76qNOjSpctrr722e/dudeoDrJmc7XTZdQdAQY8dNCQhIUG+tiJLTU0t\nXbq0RZvKlSvHxMQUbF2AfoRcW0vvHaBjBDtoiI+Pz86dO+Pj4+WHfn5+x44dMx/fk56eHh4e\n7uPjo1KBgE6Q7QC9IthBQ4YPHx4bG9u6devw8HBJkoKDgy9dujR27Nj09HQhRGJi4ieffHLm\nzBnzi7MAcodsB+gS99hBQ7p27RoZGTl27NhmzZp5enpWrFjRw8MjODh49uzZFSpUiIiIePjw\nob+//8iRI9WuFNADJdtx4x2gG/TYQUMMBsNXX311+fLloKAgZ2fnEydOREVFCSHu3LkTFRXl\n5+e3evXqPXv2mN+HBwAAFPTYQXM8PT2Dg4ODg4OFEA8ePIiPj/fw8GCSKiD/ZOeyLL16gFWg\nxw6a5uLiUr58eVIdoDruyQOsAsEOAJAtZDtA+wh2AAAAOkGwAwBkF512gMYR7AAA2cUQCkDj\nCHYAAAA6QbADAGQXl2IBjSPYAQBygGwHaBnBDoCV4TYv1ZHtAM0i2AEAcoxsB2gTwQ4AAEAn\nCHYAgBzjgjigTQQ7AAAAnSDYAQAA6ATBDgCQM1yHBTSLYAcAyBmGxAKaRbADAOQMPXaAZhHs\nAAAAdMJW7QIAIMe01mPEpUkAGkGPHQAAgE4Q7AAAAHSCS7EAoI6UxJRZw+fEXrlpNBr7fNvL\nq0G19HTTwq8WX/onysbGpvfX73s1rKZ2jQCsDD12AKCOA+sPFnUvOils3PtB7/0+ZZUQ4uSe\nfx7dfzQpbFy/8X2WT1qhdoEArA/BDgDU4d20ZqfBHZMeJd2Kue3qUUwIUaxE0QdxD+Nvxcdc\niHEt6ap2gU/FYBFAs7gUCwDq8CjvIYRYGPTL6fCz/YIDhRDlvDztHewm9Jocf+vep7OGqF0g\nAOtDjx0AqOPKuasJ9x/1G9d7Uljw4q+XSZIUtnBzDd/qk7dMmLZ7yrxRC1OTU9WuEYCVIdgB\ngDqObD62LWSHECL+1v20lFQhRGpyqindJISQTFJqcqokSSqX+HRcjQW0iUuxAKCOtr1b//Tp\n7FGbRhsMhoGT+hkMhvb92s4dufCrt8amp6X3/qaXvaO92jUCsDIEOwBQRxFX51FLPjdf4+Ti\nNOxnrd9ap7Wv/QBgjkuxAAAAOkGwAwBkF911gMZxKRYAkF3KmAkSHqBN9NgBAHKMUbGANhHs\nAAAAdIJgBwDIMS7FAtpEsAMAANAJgh0AAIBOEOwAAAB0gmAHAACgEwQ7AECOMd0JoE0EOwBA\njjEqFtAmgh0AAIBOEOwAAAB0gmAHAMgx7rEDtIlgBwDIjZBra4l3gNYQ7AAAuUe2AzSFYAcA\neC5kO0A7bNUuAABg9ZRsxzQogLrosQMAANAJgh0AIG/QXQeojmAHAMgDpDpACwh2AIDnRaoD\nNIJgBwAAoBMEOwAAAJ0g2AEAAOgEwQ4A8LyYoxjQCIIdAACAThDsAAAAdIJgBwDIA1yNBbSA\nYAcAAKATBDsAQN6g0w5QHcEOAJAvyHlAwbNVuwAAgH4Q5gB10WMHAACgE/TYAQDyXk/PALVL\nAF5E9NgBAADoBMEOAJDH6K4D1MKlWABAniHSAeoi2AEA8gCRDtACgh0A4LkQ6QDtINgBAHKJ\nSAdoDcEOAJAD5mEu5Npash2gKQQ7AEC2ZMxwpDpAa5juBACQLXxdGKB9BDsAQLbQPwdoH8EO\nAJAt9NgB2kewAwAA0AmCHQAAgE4Q7AAAAHSCYAcAAKATBDsAQHYxfgLQOIIdAACAThDsAAA5\nQKcdoGUEOwBAzpDtAM0i2AEAAOgEwQ4AAEAnCHawAqmpqSdPnoyOjpYkSe1aAAghRMi1tVyQ\nBTSIYAdtSUhImDhxYkBAQPfu3Q8ePCiEOHfunLe3d7169SpXrly1atXDhw+rXSMAABplq3YB\nwP+7d+9e06ZNL1y4ID8MCwvbuXPnO++8c+vWrffeey81NXX9+vUtWrQ4f/58hQoV1C0VAAAN\noscOGhIcHHzhwoUxY8ZcuXLl8OHDlStX9vf3//fffw8fPvzrr7+uXLly165dycnJ3377rdqV\nAgCgRQQ7aMiff/7p7+//zTfflC9f3tfXd+7cuampqW+//XatWrXkBk2bNm3VqlV4eLi6dQIA\noE0EO2jI1atXvb29lYe1a9cWQlSsWNG8TZUqVa5cuVLAhQEAYBUIdtCQ8uXLn0vPJnQAACAA\nSURBVD17Vnl45swZIYRFjLty5Qo32AEAkCkGT0BD2rVrN23atODg4AEDBty4cePDDz+0sbFZ\nvXr16NGjq1evLoQ4duzY5s2bAwMD1a4UeEJPzwC1SwAAIeixg6aMHj26WrVqY8aMKV26dIMG\nDSIjI/ft21eyZEk/P7/evXv36NHjlVdesbOzGz16tNqVAgCgRfTYQUNcXV3/+uuvadOmHTly\nxNHRcfjw4U2bNt2wYUPHjh2XLFkihKhQocLy5csrVaqkdqUAAGgRwQ7aUqRIkTFjxpiv8fHx\nuXDhwsmTJ4sVK1a1alWDwZCLzY4ePfrQoUNZt0lLS4u7dy8XGwcAQCMIdtC69PT06OhoJyen\nihUr5i7VCSGqV6+empqadZtdu3bZ29vnbvsAAGgBwQ4aEhQUVKlSpX79+skPU1JSvvvuu4kT\nJyYmJgoh7OzsevbsOWXKFDc3t5xuuWfPns9sM23aNOfChXO6ZQAAtIPBE9CQ8ePHL1u2THk4\nePDgMWPGODg4dO7ceeDAgQ0bNly8eHGTJk0eP36sYpEAAGgWPXbQqJMnT86fP79x48ZhYWEl\nSpSQVy5ZsqR3794TJkwYN26cuuUBAKBB9NhBow4ePCiEmDZtmpLqhBCBgYH+/v6bN29Wry4A\nALSLYAeNun79uhCiTp06Fut9fHzOnTunRkUAAGgdwQ4aVa1aNSFEdHS0xfrY2NhcDJ4AAOBF\nQLCDtkRERAQHB4eGhlavXt3d3f3bb781mUzKs+Hh4Rs3bvT19VWxQgAANIvBE9AQT0/Pa9eu\nmU9QvHr16s2bN7dr104IMWTIkHnz5tna2gYFBalXIwAA2kWwg4bExMQ8fvz44sWLkf8TERFh\nZ2cnPxsaGurl5bVgwYL69eurWycAANpEsIO2ODk51alTJ+OYCSHE/v37q1atWvAlAQBgLbjH\nDlaDVAcAQNYIdgAAADpBsAMAANAJgh0AAIBOEOwAAAB0gmAHAACgEwQ7AAAAnSDYAQAA6ATB\nDgAAQCcIdgAAADpBsAMAANAJgh0AAIBOEOwAAAB0gmAHAACgEwQ7AAAAnSDYAQAA6ATBDgAA\nQCcIdgAAADpBsAMAANAJgh0AAIBOEOwAAAB0gmAHAACgEwQ7AAAAnSDYAQAA6ATBDgAAQCcI\ndgAAADpBsAMAANAJgh0AAIBOEOwAAAB0gmAHAACgEwQ7AAAAnSDYAQAA6ATBDgAAQCcIdgAA\nADpBsAMAANAJgh0AAIBOEOwAAAB0gmAHAACgEwQ7AAAAnSDYAQAA6ATBDgAAQCcIdgAAADpB\nsAMAANAJgh0AAIBOEOwAAAB0gmAHAACgEwQ7AAAAnSDYAQAA6ATBDgAAQCcIdgAAADpBsAMA\nANAJgh0AAIBOEOwAAAB0gmAHAACgEwQ7wJLBYFC7BAAAcoNgBwAAoBMEOwAAAJ0g2AEAAOgE\nwQ4AAEAnCHYAAAA6QbADAADQCYIdAACAThDsAAAAdIJgBwAAoBMEOwAAAJ0g2AEAAOgEwQ6w\nJEmS2iUAAJAbBDsAAACdINgBAADoBMEOAABAJwh2AAAAOkGwAwAA0AmCHQAAgE4Q7AAAAHSC\nYAcAAKATBDsAAACdINgBAADoBMEOAABAJwh2AAAAOkGwAwAA0AmCHQAAgE4Q7AAAAHSCYAcA\nAKATBDsAAACdINgBAADoBMEOVuDmzZvHjh17+PCh2oUAAKBpBDtoy+XLl3v37j19+nT5YXh4\neO3atUuVKtW4cWMXF5e2bdtevnxZ1QIBANAuW7ULAP5fZGRkkyZN4uLivL29hRBnz55t2bJl\nampq69atK1eufO7cuS1btvj6+p49e9bd3V3tYgEA0ByCHTRk1KhR8fHxa9aseeutt4QQX375\nZVpa2vbt21999VW5wapVq7p27TpmzJiff/5Z1UoBANAigh00ZN++fR06dAgICJAfHj16tG3b\ntkqqE0J06dLltdde2717tzr1AQCgbdxjBw1JSEhwdnZWHqamppYuXdqiTeXKlWNiYgq2LgAA\nrAPBDhri4+Ozc+fO+Ph4+aGfn9+xY8ckSVIapKenh4eH+/j4qFQgAACaRrCDhgwfPjw2NrZ1\n69bh4eGSJAUHB1+6dGns2LHp6elCiMTExE8++eTMmTPmF2cBAICCe+ygIV27do2MjBw7dmyz\nZs08PT0rVqzo4eERHBw8e/bsChUqREREPHz40N/ff+TIkWpXCgCAFtFjBw0xGAxfffXV5cuX\ng4KCnJ2dT5w4ERUVJYS4c+dOVFSUn5/f6tWr9+zZY34fHgAAUNBjB83x9PQMDg4ODg4WQjx4\n8CA+Pt7Dw8PR0fF5thkdHX3p0qWs25hMpnST6Xn2AgCAugh20DQXF5fChQtfvHgxLS3Ny8vL\nzs4ud9vp37//zp07n9ksNvZm7rYPAIAWcCkWGhIUFLRw4ULlYUpKSnBwcJEiRWrUqFG7du3C\nhQv37ds3Li4uF1vesWOH9Cz29vZly1jOrgIAgBUh2EFDxo8fv2zZMuXh4MGDx4wZ4+Dg0Llz\n54EDBzZs2HDx4sVNmjR5/PixikUCAKBZXIqFRp08eXL+/PmNGzcOCwsrUaKEvHLJkiW9e/ee\nMGHCuHHj1C0PAAANoscOGnXw4EEhxLRp05RUJ4QIDAz09/ffvHmzenUBAKBdBDto1PXr14UQ\nderUsVjv4+Nz7tw5NSoCAEDrCHbQqGrVqgkhoqOjLdbHxsa6ubmpUREAAFpHsIO2REREBAcH\nh4aGVq9e3d3d/dtvvzWZzS0XHh6+ceNGX19fFSsEoBEGg0HtEgDNYfAENMTT0/PatWtjxoxR\n1qxevXrz5s3t2rUTQgwZMmTevHm2trZBQUHq1QgAgHYR7KAhMTExjx8/vnjxYuT/REREKJMS\nh4aGenl5LViwoH79+urWCQCANhHsoC1OTk516tTJOGZCCLF///6qVasWfEkAAFgL7rGD1SDV\nAQCQNYIdAACAThDsAAAAdIJgBwCwSpIkqV0CoDkEOwAAAJ0g2AEAAOgEwQ4AAEAnCHYAAAA6\nQbADAADQCYIdAACAThDsAAAAdIJgBwAAoBMEOwAAAJ0g2AEAAOgEwQ4AAEAnCHYAAAA6QbAD\nAADQCYIdAACAThDsAAAAdIJgBwAAoBMEOwAAAJ0g2AEAAOgEwQ4AAEAnCHYAAAA6QbADAADQ\nCYIdAACAThDsAABWyWQyqV0CoDkEOwCAVZIkSe0SAM0h2AEArBI9dkBGBDsAgFWixw7IiGAH\nQOckSapSu0HJitXT0tLkNSf+OdXopVeFECErVvbs+0EutmnrUiIvS0SumEwmg8GgdhWAthDs\nAOjcsb9PJCcnFyrkuGf/geff2vUbsc+/EeQJeuyAjAh2AHTut9BVPd/t9s7bnVaE/pHx2YRH\nj95+L9CtbGW/5q//c/qMvHLnnr11fF8qWbF6t/f7XPv3uhDi9Nlzzdt06N67/6tvdHyzS/f0\n9PTajfyTk5N/nDWnbNVazh6efs1fvxB5sUAP7IVHsAMyItgB0DOTybRy9dpePd7t3u3tP9Zt\nTElJsWiwbuOfgT3evXn5Qrs2rd7t1S89Pf32nTvdevb96Yfvr0WerlShQq/+H8ktDx05Vqe2\n99H9Ozes+s3Gxub0sfA7d+NGjP5m09rQ21cu1vWp/cOMWQV+fC80Bk8AGRHsAOjZvgMHy5Qu\nVatG9bo+tUu4F9++a49FA99GDTq2f8POzu7Lz4fH3rwVFX15fdjmVq+1eOUlfzs7u29Hjzry\n119ygCji7Dxi+CdFnJ2V1xZ3c4s6c7xObW9JkpydC99/8KBAj+2FZzLRYwdYslW7AEBzuL6j\nJytW/fHP6bPu5aoKIR48fPj7qjXt2rQyb1C+nKe8YGdnV86z7M1bt2Ou/bttx+7ajfzl9RXL\nl4+7d08I4VHC3Wh84j/Ddna2C35Zum7jJoPBULiwU+lSJQvikPA/kiQxdAKwQI8dYIlgpxup\nqamr1qw/cWjvnZiLd2IuHtu/c836jUlJyeZtrlyNkRfS0tJirv1bzrNsSY8Sb7Zre/pY+Olj\n4ScP71u1fElxNzchhEWqE0KsXrsh9I91a38PObZ/Z99ePQrmoKCQJEkwKhZ4EsEOsMQECrqx\nY/deV9di1atVlR/6eNdydS22edt28zZH/zq+dkNYenr6pKnTvapVKV/Os12bVhv+3Hz85D+p\nqanf/fBj/0FDM/2VSEpKvhEbW6Z0qTKlS12NubZ42fKEhEcFcVT4H+6xAzIi2AHQrRWhf7zV\n/g0llhkMho7t31ix6v/HxhZ1cRn12bB5i5aUrFh9287dyxfPNxgMFcqXWzh7RvfA/h4VvLbv\n3LN0/uyMW36zXdsqtet37RyQlJxcpkqtrj17D/lw4L/Xb/yxbkMBHRvoXAcyY+APA5A5ODh0\n6xwQsmLl5nWr2rz+qtrlAHiG31au7tnvg9k/Tp0646cLEcw1AwhBjx2QEZdiAatgkvjmCcAS\nwQ4AYJW4xw7IiGAHALBKTHcCZESwAwBYJaY7ATIi2AEArBKXYoGMCHYAAKskSYL+OsACwQ4A\nYJVMJhOXYgELBDsAgFViHlYgI4IdYImZsQCrYDIxjx1giWAHALBK9NgBGRHsAEv0AQBWwWQy\n8ccKWCDYAQCskiRJgnGxwJMIdgAAq8Q3TwAZEewAAFaJ6U6AjAh2gCU+KQCrQI8dkBHBDrDE\nSDvAKphMfFcsYIlgB1hiDgUAgJUi2AGWCHaAVZAkiQ47wALBDgBglSTBdCeAJYIdAACAThDs\nAAAAdIJgBwAAoBMEOwAAAJ0g2AEAAOgEwQ4AAEAnCHYAAAA6QbADAADQCYIdAACAThDsAAAA\ndIJgBwAAoBMEOwAAAJ0g2AEAAOgEwQ4AAEAnCHYAAAA6QbADAADQCYIdAACAThDsAAAAdIJg\nBwAAoBMEOwAAAJ0g2AEArJJBGIQkqV0FoC0EOwCAVbKxMZoIdsCTCHYAAKtkZ2cnSSa1qwC0\nhWAHALBKkiQJYVC7CkBbCHYAAKuUlpZmNBDsgCcQ7AAAVik93WQg2AFPItgBAADoBMEO1mHB\nggX79+9XuwoAADSNYAfrMGDAgJCQELWrAABA02zVLgD4fxs3bszi2StXrigNOnToUCAVAQBg\nTQh20JA333wzi2c3b968efNmeVliVlLghScJSTB4AngSwQ4a8ssvvwwaNOjx48e9evWqWbOm\n+VOjRo3y9fXt1KmTWrUBAKB9BDtoSGBgoJ+fX7du3UJDQ2fOnNm3b19lLoNRo0bVr19/5MiR\n6lYIAICWMXgC2lKjRo3Dhw8HBgb279+/e/fu9+/fV7siAFrGXRnAEwh20JxChQrNnj07NDR0\n06ZN9evXP3z4sNoVAdAs7rEDnkCwg0Z16dLlxIkTJUqUeOmllyZPnqx2OQAAWAGCHbSrUqVK\n+/bt+/TTT7/44gu1awEAwAoweAKaZm9vP3ny5LZt2x4/frxu3bq53s6IESP+/vvvrNukpqbe\njYvL9S4AAFAdwQ5al56e7unpWapUKS8vr1xvxM/Pz83NLes2u3btKlTIMde7AABAdQQ7aEhQ\nUFClSpX69esnP0xJSfnuu+8mTpyYmJgohLCzs+vZs+eUKVOeGdEy6ty58zPbjBkzxqmQU063\nDACAdnCPHTRk/Pjxy5YtUx4OHjx4zJgxDg4OnTt3HjhwYMOGDRcvXtykSZPHjx+rWCQAAJpF\njx006uTJk/Pnz2/cuHFYWFiJEiXklUuWLOndu/eECRPGjRunbnkAAGgQPXbQqIMHDwohpk2b\npqQ6IURgYKC/v7/yjbEAAMAcwQ4adf36dSFEnTp1LNb7+PicO3dOjYoAANA6gh00qlq1akKI\n6Ohoi/WxsbG5GDwBZF9aWpqhsFupStWVf+O+m2LeYNLU6dN/mv2ce7lz926Ner65e+3FS1GG\nwm4jR3+jPDQ6F5dL9Wv+et8Ph6SkpDxneQCsFMEO2hIREREcHBwaGlq9enV3d/dvv/3WZDIp\nz4aHh2/cuNHXN5cfh0D2xUZfUP4FjfhM7XKe8Ovvoc2a+i1fuUr56yhTupRc6uE92+Pv3w9Z\nEapuhQDUwuAJaIinp+e1a9fGjBmjrFm9evXmzZvbtWsnhBgyZMi8efNsbW2DgoLUqxEvqJSU\nlKGfjdy0dXtJjxJVq1T2bdhgf/ihydNnurm6Rl661KRxI4PBEHHxUkpKypoVIQ4O9l+P/y5k\nxUqDwfB+925jRn1x4ODh2fMXlSld6mJUtJ2d7aLZPw369LOr164N+vSzWdMmP7Oxs3NhpRJJ\nkkJWhP4yd1b33v3DDx15yb+JeZ0PHj68eeuW0cg3qAIvKIIdNCQmJubx48cXL16M/J+IiAg7\nOzv52dDQUC8vrwULFtSvXz9fy4i+fCVftw+r4FnNW1n+YdL4u3Fxl69cjTh5NP7+g8Yvv+rb\nsIEQ4sDBw9ciz6SnpxcvV+WP35ZOnvBt5+69duzeI4TYsn3HP4f3CyFavvGmX+NGzoUL/7F+\nw78Xz7q5unbs+t72Xbt/nj7ln1Nnfp4+JWzz1mc2DnizvVLMkWN/JycnN/Vr3PmtN5evXCUH\nuxuxNz2reUuSFHvz1ivN/Lt17lTApwuARhDsoC1OTk516tTJOGZCCLF///6qVavmdwGODg4f\nDh3erk2rcp5l83tf0LJrkWfMH74b2G/QwH4ODg4lPUp06xwgr2zq19jR0UEI4V7c7bUWrwgh\nypfzvBcff/Kf04E9uhcu7CSE6Plut91793d4o01TX183V1chROVKFZOSkpUt7967P/uNhRAh\nK1Z26xxgNBo7d+zQuXuvHydPFEKULlVSLjj68pU+Hw5eH7bp3a7PnpQbgP5wjx2sRgGkOiHE\nvu2bbGyMbd/q4t2waZ8PBhfAHmEVjEajwfDf65s2NjYWCxbLkiQpjQ0GQ3p6uhDCza1YplvO\nUePU1NTfV61ZvW5DVZ+GfT8akvDo0badu80bVKpY4T9DP17wy7JMXw5A9wh2wBMaNai3d+uf\nJUt6/HvjxqkzZ9UuB1rR/KVmcxf+kpqaejcubtXa9c9o/HKzkBUrExOTHj9ODFmxssUrL2Xa\nLC09PfuNZVt37Cpa1CXqzPGLp/66eOqvQQP7Lv99lXkDSZL27DvQsH7dHB4fAJ3gUixgyb+J\n784/13lUqFaqpIfatUA15vfYveTfZOn82cdP/lO9nm8J9+Lvd3/H1TXzHjWj0WBrY9vhjTaH\nj/5Vt8nLkiS9161LuzatDhw8/EQbW9tiRYsWK1q098CPF8/96ZmNlYchK1YO6NNL6eHrF/i+\nX/PXPx82RL7HTghhY2PzVod2349lgBHwgjJIkqR2DYAmODg4dOscsGzhHCHEiX9ON/Bv3vmt\nDh3eaNPjna42NjZGI93bgLZMnTFrTPDEad+NnzrjpwsRF9UuB9AEPquATNy8dUuSpD/Wbezz\nwWCPCl6vtntL7YqA/wq5tlb+p3YhALSIYAdkos3rr969FpVwK6ZalcpeVavs2Xdg3cY/1S4K\nL7S0tLTeAz+uXq/xF22/PLb1L7XL0QSTyWRgwj7gSdxjB2TOzbWYECLin2NCiKo+DYePDDKZ\nTJ06dlC7Lryg/li38c7du2f/OjTr+NKggLH1X8vf2RytgiRJBkGyA55Ajx3wbJvWhl6/cWPz\nth1qF4IXV8UK5b8dPeq3GxtcihdxcnES3B4thCRJ5DrAAsEOyERKSsrk6TO2bNspP6xWpXLZ\nMmXOnDvPYCOoxbdRg7Pu0bGXY6cMnN6+3xs2tjbPfo3emUwmeuwAC1yKBTJRq0GTS9GXhRBF\ni7pUqVTpo/59Wr/Wcvb8RUeO/e3XuKHa1cEqpaWlfT9txsrVa21sbPoF9hw0sF+OXh5ybe2B\n9Qf3rtrX7T9dqtSt3NMzIJ/qtCI2NjYm/q8FPIlgB2TijTatfpoz//vx39y8eWvP/vABgz8V\nQgzsG1i/ro/apcFaTZwy7cdZc0d9PiwxMfGLoLGOjo59e/XI/sv/3nH82La/v1j8mY0NV1r+\ny8bGhk50wALBDsjEzKnf/bFu46y5C7ZvXDN5wrdLl68wmUyBPbozmx1y7fdVa8aM+nzooA+E\nEAkJj1aErs5ZsNt5IuKvyC/afik/7HSstfz1si8ySZK4EAtYINgBmTu2f0cV7wY/z1/4w6Tx\ngT26q10OrJ7RaExLS5eX09LSzL9bNjv6j+9j/pBUByBTBDsgc1169ElJTa1dq5bahUAn3n/v\nnW8nfp/wKCExMWnWvAULfp6hdkW6QJcd8CSCHZC54yf/6fZ2wJr1G+Pj7w8fOkjtcmD1/jP0\nYwd7+xWr/rCzs5s7c1qPd7rm6OU9PQP4tgkLBoOBW+wACwQ7IHNDPxr43Q8/2tjYbN2x68y5\nc3Nm/GBnZ6d2UbBiRqNx6KAP5HvscoeRsBYMwsB8foAFgh2QuUnBY8uULmVjY7Ny9dpFS38d\n++WI8uU81S4KwP8zGg2SINgBTyDYAU8ld66ULVN674Hwx4mJT2v2tAtkmfavyI2f9hRdMkAO\ncZMd8ASCHfAMAW+2t7OzmzxtxsLZMzM+m8VtTznKcHJj85dkumWSHwAgCwQ74Nm6vR2wdsOf\nP06e5Oxc2Hx9dm5mV9pkmuSyfgmALHB/HZARs60CzzbxmzEPExK+mfCdxfqengHZ70ILubb2\n+RNbnmwE0A0uxAIWCHbAs5XzLPt2wJvTZ82p7F3/+o3Y/eGHxgRPvP/ggfxsjuJdniDbAUII\nvk8MyIhLsUC2zJz6vWeZstNnzfZ/tU1qaur1G7EzZ8+rXKmiZ9kygT3evXnr9nmPKz7eNR0d\nHQumHkZaAAAyItgB2eJe3G3yhG8a1K/z28rVp86cbf9G67KlS0dcvPTP6TNvvxdoMBgcHR3e\n/bG7v1/jgsl25qnuabfxAfomCUkYuBgLPIFgB+RA965vd+/6tsXKpKQkR0fHO3fjKtWs67q4\naP26dQqsHotrsnTjAcALjnvsgOcld9G5F3crW6Z0YmJSwez0aaMouP0OLxT66wALBDsgz7Rv\n2ybu3r27cXHqlkG2wwuCwRNARgQ7IM+827WzECI1NU3tQsh2APCCItgBeaZxw/oXl1w48c8/\n9+8/ULsWQP8MjJwAMmDwBJCXjuzd7uzh6b6oeKMG9S2eSktLj7x48cbNmwaDoUqlSuXLeapS\nIaAnXIsFLBDsgLxkNBp/mDh+z62j9+8/KFrUxfyp+PvxE976zr24W8UK5R/1elyokKN78eL5\n1+WQxTfVAgD0ikuxQB77cECf3RN2njl//kbsTZPJFH//vrzevXjxT1cMTU1LO7pv556Ju4Y3\n++LkqdP5XQw32wHAC4VgB+S92t614uLu/XX8xKat2/eHH1LGySYnpxR3dRVC3Ig6P3v61Kld\np+3ZfyAtLX8HW5DtAODFwaVYIO/VrO619eedn3869OatW7PnL7L/yt7RwaF06ZKzA+csnD1T\nbvPhgD7169V5uVU7MVvUr1fHpUiR/KuHy7LQJUmSGD0BWCDYAXlvxpRJynKXgI4fffrZlStX\nG33i9870d4InTe7bq4f8lF/jhmf+OljP72WXJS716/rkd1V8LwUA6B7BDshfXtWq7ghbm56e\nbutSQgjx+9KF5s9Wq1K5U8cO12/ccHSwr1Hdi+kbgJzhLwZ4EvfYAQXBxsbmP0M/HvX5sG5v\nd7J4KmTR3Mdb709467uHCQn5XQb320FP+OYJICN67IACMmVi8NOeOnch4sMlH+TrbXYAgBcB\nPXaAytaHbbp4Kaq4q6vJZCqA3dFpB92QJIlrsYAFgh2gsrtx97pN63YxKvrIsb/j4u4VwB7J\ndgCgVwQ7QGXXb9wQQrgWK3rn7t3ww0ciL10qgJ2GXFtLvIO1M5lMRsYbAU8i2AEq++qL/4z0\nHTSz07jSkW4HJu+7cydO7YoA62A0Gk2MnwCeRLAD1Fevro+tre2078d7Va1yLz6+YMb6Macd\nrJ2Dvb0pPV3tKgBtIdgBWnH67LkSncpUKF+O2eyA7LCzszNJBTHkCLAiBDtAK7yqVlk7Yk1K\nSgqzcwHZwahYICOCHaAV9vb2478Omv7Oj6fOnC2A3cnjJxhCAetlMpno3QYsEOwADfnP0MEt\nW7xyL/6+EKJgprUTzH4Cq0WPHZAR3zwBaItLkSKJiYkRkRcvRkXb29k5FXaqV8fHqVChrF9l\nSjd90Hiwc9HCQoiq9ap8PO3DAikWUBM9dkBGBDtAW0qV8riX/vhSdPSR6YfcXF0vRUU/mvD4\n1eYv29jYZPGquzfiavrVGD57aO52GnJtLYNkYXUkSdBjB1jgUiygLR/07XP4h4OGA2kRJ48e\n2bvdr3FDIZ49liL28s3Eh49nDZ8754sFsZdjc7FfLsjC6tja2kgFdccCYC0IdoC21PXxjjp7\nPGThXKPRKITwb+KbnJyyaev2/eGHsniVja2xVtNagWN7ejepOWXg9NyNq2UsBawLEwMBGRHs\nAE0LGvHZYO/AneN3ODo6ZtGsVpOanT7u6Fy08Eud/NNS0xPuJeR6j2Q7WAtJkrgSC1gg2AFa\n16yp3/3794sWLZJFmw1zw9bN3iiEuB1z25Ruci7m/Dx7JNvBKjAqFsiIYAdo3YGDhzt937lM\nqdJZtHntvVcj/or88s0x0wfPGjRloMH4vJ92ZDsAsEaMigW0bl/4QVFUpKWnZdHGqUihzxcM\nK7CSAC0wGo2Cr2kBnkSPHaB1gz8YcHjawfBDR679e70g90unHQBYHYIdoHVFijhfOv13wuZ7\no14Nio29WZC7JttByxg8AWREsAOsgNFoXLfyVwcHh6Tk5ALetZLtCHnQGgZPABkR7ACrUdSl\nyJWYmOSUlALerzK/nUW2I+oBgNYweAKwGkf27fBu2FTMMLzSrInBoPL/ykh1UJ0kSfTXARYI\ndoDVqFCuXPdubyclJqp4+Yk8B+3gHjsgIy7FAtYhLS3ttfYBD+umeHqW7R50CAAAIABJREFU\n5ZuUAACZoscOsA5hW7aVCChdt07tcmXLql0LAECj6LEDrINnmdKhw0MTEh6pXcj/6+kZoHYJ\nAIAnEOwA69Cwfr16dXzuxsWZTCa1axGCVAcNMBgMgi+eAJ5EsAOsRrs2r8/rM3/Hrj3Xb9xQ\nuxZAfUajkVwHWCDYAVYjeMxX/148E/JxSGJiktq1AOozGAx8VyxggWAHWJPklBSTJDk5Oalb\nBtdhoQUGg4FYB1gg2AHWpGL58j7etc5duJCWlp7x2fT09PT0grgDj9nsoAUGg0Fwkx3wJKY7\nAazMtg1/eHrVLvZr0dq1atrb2yvrk5NTDhw67FqsaP26dVQsDygwDJ4AMqLHDrAyJUt6zJ05\nbfq7M3bu2ZeWlqasvxgVlZaaeiP2ZlISd+DljZAVK9t3fufjYZ+//V7gX8dPZGxw+crVKT/+\nJC8v+GXZW916fDh0eN8Phzx+nJjNXRw6cqxMlZrySOcpP/7U7f0+vfp/9NEn/8mrQ9A3LsUC\nGRHsAOvz87yFNjY21apUtrX9b6d7VPTlK1dj/pp5dONXGw4eOZqYmN1gkTsvzj12Hw3oO2va\n5C+GDV23cdMzGw/5aOCcGT+80fr1pctXZNEsPT09JSVFXl4ftqlTxw5H/zouP/xi2NClC2Yn\nJDy6/+DB8xeve3wFC5ARl2IBK/P48eO/jp/4ZvOYypUqKitTUlN/G/rbmFGfDxrQr2aDJqZJ\npkYN6xd1ccmPAl6cVCeTJOnvEydrVvfaH35o3qIlRqMxaMR/Ii9FLVyyzGg0+jZqaNG+WVO/\nr8dPMm8ce/OWsvxLyG+Xr1wd++UX1apWEULEXPt31OfDflu52q9xQyHE1BmzCjs5JTx65FKk\niAqHam24xw7IiGAHWBknJyfPsmXkHh1JkpKSkhwdHUt6lBBCdO/6dsmSHlFnT9Rq0OSX2CU9\nf+rh6OhoZ29XsXz5MqVLqV24VZo9f9Gmrdv/OX1m2YI5XwSN9fGu9eDBwyW/roi9eXPlssW7\n9u47fvKUxUuMRqPBYJgxe67S+EJkpLIshBj52adyqouKvnz2/IUlIb/t3rc/eMyXQoj/DP24\ndq1ab3Z598rVmIoVyhf88QKwdgQ7wPrY2dnduxefmJh4PiLy3+s3jEajvb2dra1tjepeQohi\nRV2uXzq7befuHbv2XL569e/j//gOa1LczdXBweH5d/2iddd9NKBvhzfa7A8/tHnbjsJOhYcN\nHpSWnnbr1p0vv/5WCGE0ZnI3y4GDh+vV8Tly7G+l8cQp05TlpctXeJYtI7dcH7Zp7swfGjWo\nP2zEV1djrskrHR0dOrZ/o8AO0KoxeALIiGAHWJ8/li+t799cmistGxTyTdDItPT0o8f+7t/7\nffM2rV5t0erVFkKIlJSUYmUqllhcvGH9eqpUqwOeZcv8Frr6g369+3881LVYsaARn/V8t1v3\n3v2LODs3qFdXafbzvIVrN4Q9fpw4c+p39evWURqbv9B8swcOHR7y0UAhxOstm2/ctEUIMXXG\nLI8SJVq83Izuumwi1wEWDBLTdgNCCCEcHBy6dQ5YtnCO2oU826NHj0pUqJaUlNy1c8DvSxc+\ns/28RUt2ph+qV8dH6SvKnRetuw4aN3/x0o+Hff7TD99PnfHThYiLapcDaAI9doD1KVy4cNy1\nKDs7Oxsbm+y0/331mjJdy9vb2+V3YUBBYlQskBHTnQBWydHRMZup7rMvx+zeu796taoeJUrk\nd1VAQSLYARkR7AA9W7xs+Q8zZgXOCfQsW/Y5N8V12Oy7c/duUlKy2lXon9FoENxNBDyJYAfo\nVvihI0P+M+LjZR+92vxlW9tsde/h+Q39bGSJ8tVKVaq+Y/cetWvROaPRSKwDLBDsAN16pXX7\n92f1rOtT+/k3RXddNt2+c2fm7Hlrfw/p1LHD5Gkz1S5H54xGI+P/AAsEO0C3iru5OTsXdnR0\nfM7tkOqyz7mwc+HCTv/X3p3HRVE/fhz/7C6CAl4ECIoCKt63eWvq17xvK9PyKDPL0sooj6/3\nraVZambe5plXmbeppKl4G4qigCByiKjch7DH/P7YvvuDReVmlt3X8/H9w52dnXmvXwfefWbm\nMwePHLt563YlZ2e545g5SZK4yA4wQrEDzFaXzq89fRqrf8B8QWyN+H1rxO+FEsnslSlTes+2\nzcH3QuvU8lo8b5bcccycJEmC+yeArJjuBDBbbw7o9+a774nVwqtmDbcqlbmFsHj06NqlR9cu\ncqewCJyHBbJjxA6mKDo6+s6dO2q1Ovtbjx8/joiIKP5IJdGg/n3D7ty4tyXwv12mnzl73ufM\n3xcvX9FqXziAFx4RaXiwVXYM2sHU6HScigWMUexgWi5fvtygQQNXV9e6deu6ubmtX2/8WIWh\nQ4dWrVpVlmwlUVW3KqePHXxr0IDQbcHpJ5JXv7fmhr//i1a+GxQUERn5kq3R7WBSOBULZMep\nWJiQkJCQ1157LT09vUuXLtbW1idPnhw9enRSUtIXX3whd7SSbdeWDfo/tGnV4uzDa64uLomJ\nibFx8a6uldz/15IlSZeenuHk5RT5LLpKaRf5wgK5VfDrRwHzQ7GDCZk+ffqzZ88OHjzYu3dv\nIUR0dHTbtm0nTpzYtWvX+vXry53OHIz9cNQEh6nSt5JKpbKysnri/zQwMNjaxrq0jU1ycook\nSRH+kY2615M7JpArkiQxYAcYodjBhFy4cKFbt276VieEcHFxOXz4cJMmTSZOnHjo0CF5s5kH\npVJ56fSJsuXsa3h6ajSadRu33L0V/CAi4lFMjLOjY7cu//lLupSenv6SGVL0Z2OZAAWmQKfT\nCUGzA7Kg2MGEPHz4sFu3bpmX1KlT54svvli8ePFff/3VqVMnmXKZlSaNG+r/YGVl9fGH72d+\nKybm8VhP7/aXWuc49R31DqZAp9MxYgcY4eYJmBAPD4/r168bTWEwZcoUFxeXUaNGJSQkyBXM\nQjg6vqJQKBITk3K5PvdSQF6SJDFiBxih2MGEdO3a9eLFi5MnT05JSTEsLF++/E8//RQaGjpy\n5Mi4uDgZ45m9u0FBkiQlJSfnfnowuh1kxJMngOwodjAhs2bN8vDw+OabbxwcHLp3725YPmDA\ngMmTJ+/fv79atWrXr1+XMaF5q1u79ofvj1g04JuAu4HPnj3L5ad4LgXkwnQnQHYUO5iQihUr\n3rx5c9q0ac2aNYuOjs781sKFCzds2ODq6hobGytXPEuwZuX3zZs2Wdh/8bV//PL0QbodZMLD\nJ4AsuHkCpsXe3n7u3Llz587N/tb7778/cuTIsLCwe/fu5XWzd+/eDQ8Pf/k6Op1Oq9Xmdcvm\n58rZU7aOVVQqfjjA1CkUCnodYISf3TB1Wq02ODhYo9HUqlWrVKlSnp6enp6eed3IZ599dvz4\n8RxXexTzOF8ZzU2v7q8/jY198vSp4yuv5P5TWyN+5z5ZFCeFQkGvA4xwKhYmZNq0aZmfIZaR\nkTF37tyyZcvWqVOnQYMGdnZ2o0aNyt+p2GPHjkk5sba2ruzKExeEEGLXlo3HZx27et3v/IWL\nCQmJuf8gJ2RRnBQKBadiASMUO5iQ+fPnb9myxfBy3LhxM2bMsLGxGTRo0JgxY5o3b75x48bW\nrVunpqbKGNISKJXKwBtXKtwss+ervecuXHyY9XpHwERwKhbIjmIHE+Xn57d27doWLVoEBgbu\n3bv3559/9vX13bRpU1BQ0IIFC+ROZ/4cKlZYsXTxzYtnj848cjcwWK1Ry50IMMapWCA7ih1M\nlK+vrxBi2bJlTk5OhoUjR45s27bt0aNH5ctlWSpVcr7w159bx227dOUad5bA1CiY6wTIhmIH\nExUVFSWEaNSokdHyhg0bBgQEyJHIQtWqWcPX59i6D9Zf++dGbrpdXi+z47I85BvX2AHZUexg\nory8vIQQoaGhRsujo6MdHBzkSGS5mjZu5P3Zp8vfWZGQi6eN5enGWH2rY4pj5I9Cocj1Q1IA\nS0Gxg2kJDAycO3fu7t27a9eu7ejoOGfOHJ1OZ3j3/PnzBw8ebNmypYwJLZPPmb/Hbv7YoWKF\nHNfMfUujzKHgOBcLGGEeO5gQNze3iIiIGTNmGJbs3bv36NGjvXr1EkKMHz9+zZo1VlZW06ZN\nky+jhXKt5JKcnBx6P6xaVTeVSlUUuzD0PCbDQy5JkkSzA4xQ7GBCwsPDU1NTg4ODg/4nMDCw\nVKlS+nd3795dq1atdevWNW3aVN6cFmjH5rWNW3XYcm/rsvPfPnfWYq1Wp5O0paz+/T/LMFmx\nvq4ZdbWXjNXR6pBHNDsgC4odTIutrW2jRo2y3zMhhDh79mzNmjWLPxKEEGXKlFm57NsNIbs0\nGuP7J+6FhEZGPZQkKT09vXmzJq/87wrIzO3NqOcBhUKr1XJfLGCEa+xQYtDq5NWtS+erKy7f\nCwnJvDBDrQ66F7Lv671nFv11dNbRO3eDXvTxXF57R/ND7mk0GmY8AYwwYgcgt3RZpzuRJMnv\nhv/uL3fFhAWVL1euU4++Bd8Fp2KRe5IkcSoWMMKIHYDcylCrM9Tq9IwM/cuIqKgV7648sHtn\n+XLlhBDP0p6lpqUV8IFvhhE7hu4AIB8YsQOQW3+fOFyzQfMym0q3btlCCBEeHlGxYoVur3fW\nv/vbr1uq1Kzv+OsrDerXMdxFkQ9G3Y4xPADIPUbsAOSWR7Vqs6ZOfvI0Njwi8m5gkLW1dWxs\nnEaj0b/r6uLy5oB+y97+PiQkTN6cAGCxGLEDkAcTJ4xv2KLd7pDdNjY2arW6d49uVlb//2Nk\n19aNvQYOfhQTU7tWod3pwrgdXkSpVPJIMcAIxQ5AHlhbW9/1u/ySFe6F3vcs5VXo+zVMmAIY\nKJVKiWeKAVlxKhZAYXr8+IlarTacny1E3E4BIyqVil4HGKHYAShMa1f9sP7DDX+dORsZ9VCr\nNZ7NuIDodshMpVIxYgcYodgBKExv9O+76upyq1JW1/1uBNwNLPTt0+1goNVqmaAYMMI1dgAK\nh6Fy2dnadurQ/uat2/fDHlRycnJycjSso9Xq1k/dGHT9nkIhBn/5xqvdmmsyNOunb7rnF1La\ntvTH335YuYarTPFR8jBcB2THiB2AQpB9IK1h/XoKhSLt2bPMC68cv5oUl7T48LzJm77eMGOz\nVqvz2XVaoVQsOjTv7a/f3DJ/e/72BQDQo9gBKCpKpTLqYfTjJ08NSxyrOL7x2UClSlnulbK2\n5WyFJEUGRzVoW1+pUtZtWSf05v1cjsHQ7SCE4DwskB3FDkBRadKoYVJy0sXLV65c/0e/pEYj\nT4/67tH3o5eM+b73Bz1VVirPBh7nD1x4+jD2xLZTKYkp6anp8mYGgBKNa+wAFBVXl0quLpUu\nXr7y6FGMVqtVqVRCiHN/+J7Z8/dg7zdrNK4uhOgwsF18TPzaKRvqtKztWMXRxtYmN1tmTjsI\nRuyA56HYAShaGq22YsUK+lZ37eT1K39em7jxK5Xq39MFEUGRTlWd+n/S955fSMiNEH5VI/eY\noBjIjmIHoEByvNytkpPz3aCghITE8uXLXTv1T+DVoIk9/qt/a8Efs6vUrHx00/GTO3xsytiM\nmjMiTztl3M7CqVRKQbEDsqLYAWareB7DNcxtwMu7Xc0anvfDwq5cu16juucg7wGj579vtMKY\nRR/kb7/5+BTMiVKppNYBRih2gBnKzU2jmdcpYEnKsdu1avHqtX/8bgXckSSpc8cOdra2Bdkd\noKdQKBixA4xQ7ABz9qJTloU4XYh+U9lnHtYvuXcjRKVSvTdreMcO7SIio/65cbOwrqErnvFI\nmDKdTie4KBPIimIHmCGjIbRimPXNMPNw/OOEaQNmNu3S1O/0jZSElEWH5oXcvL9l7tZZu6fr\ndDohhErFjx0AKCrMYweYpzyNZhWk+el3lH3m4QpO5RNjk+Jj4sPvhlesVFEIERcfr1AoNFpN\nvvcFZMYtsUB2FDvAbBVbtxPPm3m4ai03a5tSC0Z8u23Bjtff+Y8QIjLqYZXKrlxgBwBFh3Mi\nAP5VwDlEjGYePrT+aJ2WtQeO65+ckDK134wlxxc5Ob4SERnlUsnZpVKlwswNC8YVdoARRuwA\nc5aPlrY14vd8jN4ZZh7WtzohhDpdrdPqhBCSTlKnqyVJqu7pIYSwsiqc/57kzglwJhbIjmIH\nmLnsBSgjLWPZ2OWTek2d0md64LWg534qT91umNsAw8zD3l0ne3ednJ6W3vuDHg/uRkztP3P+\nsEXvzR5hXdra3s5eoVDExsbl/8vkNyHMEnfFAtlxKhYwf0Y3yZ77w7e8Y/kvVo0PuHjn1yV7\npm+f8tyP5H77WyN+zz7tsCgjJqwan3mBjY21k+MrQfdCklNShBAVK1Tw9HDP/V6eu9+8RoU5\n0Wq1PIMOMMKIHWBx6repO3Bcv2cpz2LCH1d0rlDwDea+WrV8tbmnu3tMzOOYx09uBdyJefy4\n4Hs3GrrTn0pmPM8SUOyA7BixAyyOczVnIcT6aZv8z9/+YO7IzG/pK1o+WlHuP1ivbu16dWun\nZ2T8edInMiq6XNmyNjY2Bfz1XFg1Tq1WT501LyEx8WH0oy3rV5cvV87w1uTps6dN+sre3i6v\n27xw6cqgocMjgm4plcolP6y8dOVqaZvSdna2P/2wtFAyWzKtVkuvA4wwYgeYP6PeExbwIDkh\n5YN57y06NHfjrC3ZJwMrhpObNtbWVSq7PoyOPuFz2uf03xkZGUW9x9z4de9vdevU+nnFsqkT\nvfftP5Cnz2q12ud+iz8OHRnYr8/lq9f1LydO+OyXdT8lJ6ckJCYWQmLLJkkS98UCRhixA8yf\n0TV2l45esbK2Gvhpv/iYBE2GWq5UTRs3atRAG/P48bV/bvjdvFXZ1aVKZddC2XK+i+npv89N\nn/y1EKJVi+atWjS/ePnqpq3bk5KSly6ap1/h1OkzazZsVqs1nV9r36J5M8O7K35acz/swf0H\nD84cP/Tt9yt69+jWoF5d/UfCIyKnfD1hx669rVo0F0IsXf6jna1tckpKubJlC+O7WjTOwwLZ\nMWIHWJwe73W7c+nulD7Tfxi3csyiDzL/dizgOc28NiqVSuXq4lK/bp3HT578c+OmVqsryN71\nAQoy3KjRaB1feUUIIUlS8L2QSs5OH74/wsGh4uWr1/Qr7Nt/cPvGtSPfHSKEMHp38ldfvPv2\nW9f+8bsdcNfQ6kJC79++c3fz1h3HT57SL/H+7NOV332bnJwc9iC8IN8U/8OUJ0AWjNgBFiHz\noF3ZivZTNn9dDDvKzXIhhId7tSdPnkbHxKg1apXKpoiC5UabVi2Onzw1oG/v6343Nm3dYVum\nTJNGDevVqW1YQaPRCCGUSqUQYvW6jZnfdatSuV/vnlNnzatZw9Ow/h+Hjvy84rtXmzWdMGnq\ng/AI/cLSpW369e5ZfN8KgCWh2AFmxRRmAMne4V6eJzEp6VFMjHu1qqVt5Gx1QoiR7w7xnjL9\n+Emf8IjIH5d9e+jocf/bAYlJSeERkfoVBvbrM+yDj6xUVu3atHKvVtXwrn7Us0pl1zuBgVO+\nnmDY4LkLF8ePHSOEeL1zx4NHjgkhli7/0dnJqVOHdh7u1eT4imZFoVAwYAcYodgBZiLHs6gv\nGTMrHi8KoFIqJSHS0+W/f8LGxmbld98YXo79cJTRCgkJiQtnz9i5Z1+9OrU7tGtj9G5ScrJX\njRq1vWoaluzeukn/h949uhVJYgDIimIHlGBbI37f/9PBv3ad1r+Mf5ww77eZVbyqFHCbRTfg\n99xuZ2dnV7FihcTExJSUVDs72/xttjDS5ez1/3T8dc9v1T082rdtnf3dI8dOTBg/tniSAMBz\nUeyAkq3/2D79x/YRQjwKi9n57a7KNSuLIi5nuZF5WrvcJHnFwSH4Xshff599tVmTSs7OedpL\ncXKoWDH7MJ7B4Dd4BkaxUigUnIkFjHBXLGAOJEna8c2ud6cMeckEEC86D1vAO0lf4rlbfu7C\nOrW8evfoZmNjc9P/dlx8/NPYWJ1Op9PlcJMsj5cAACMUO8AcXDxyuUbj6o5VHPUvC1jUZClM\nCoWiVYvmGWp1aGJ4WHLk8ZM+x0/6PIiITEpOLv4wuXTO9+KUGXP2HzwsdxAA+BfFDjAHJ7ad\nem1Qu5evox8qM/zPsPC5axZ+xFzsoqy9fYe2rd3KuDoqKt77JXCv956pXaYH3wt5yXZkHLS7\nGxT8n179L125Nnj4qENHj8sVAwAyo9gBJd7TqKc6ra68Y3n9y1zWMqNTovo/F91p2ecGyL6w\nbNmyzk6O1aq6+focS4t9aGVlpVAodDqdJD3/tKyMlxL63fAvX67c8QN7W77a7NKVa3LFAIDM\nuHkCKKkMg1UXj1xu+p8mBa84xV+ScpyBpV+vHkvfWvb+mvdsrK3btXnOjagyeq19WyGEUzWv\n1LS0779ZIHccS8QjxYDsKHZAibd95mq5I+SNRqMZ/cnnvpcuKYRiwezpqc212SdtEW5CCLF3\nxy+fTPhq58R93Wf3CLgbWKdWTYXCVM4zuFRy9r9y7sxZ3+ZNG3t6uMsdBwCEoNgBJZe8E5oU\nxL79B588fXr76oWH0Y/qtWq94tz3z520RW/VsiWrli3p88aQhcd+nX10hklVKGcnpzcH9pM7\nBQD8P1P5b18AlsPDvdqc6VNUKpWzk6NtOVsh/TsZ2UsmbTm4d2fvnt0Cg+/pn9ZqsDXidyY9\nAQADRuwAFLeWrzYTQgQG3ft0wle9P+ipslLplxtN2mJkzYpllWvUq328Jk9ZRSZMUQxkwYgd\nABls3bnr48++bD+2Q+e3OxoWZp60Jfs4nKuLS4UK5aMfxRRfSpg2SZKE4P4JIAtG7AAUt/0H\nD//2x6HjB/bujD5oWJh50pYXXT7o4uz8NDY25vETZ6fnj+rBokiSRK0DjFDsABS3A4ePnjh7\nukrDBvqXC/6YbVPGRj9pi3jpTSFdu3Re8fka+/V2FDsIIVQqlSRxKhbIgmIHoLitW7V83arl\nRidbe33Q4yUfSU9Pd6rmlZScPH7buMquLpnfKrl3B6OArKysKHaAEa6xAyCPPBWyuPj4pOTk\nmUemN2nUwK1K5Zw/AAvAiB2QHcUOgGye2+2eO32JS6VKNtbWdwODjp/0Ccr09FiG6yyZlZWK\nWgcYodgBkFPum9n82dMC1vtH7Xpw/36YYSGT2FkyhULBiB1ghGIHwOQ8t655fzbuH98zbdu0\nVGs0UQ+jX74yLAF3xQLZUewAyCxPp1MXzp5xfeWVa//4nT1/QavVFl0qlAzZHlICWDiKHQD5\nZe92L3pWmEKhuHXV1/qiWDNq7c1bt588fVosAWGKdDqd3BEAk0OxA2CiXjSSp1Qqf1n3k1eN\n6kvfWnbx8tWUlNRiDgYTodNJDNgBRpjHDoApyvH87F2/yzqdrryre+UNrnNivv9n1bWG9evN\nnja5eOLBFOh0Oh4pBhhhxA6AKcrNLRFKpdKhYsWohw9n9Zx7wuf0+s1bkpKTiyEbTIROp2PE\nDjBCsQMgv3zf2XrqyH7FOfWsqZNPHd4fGfXwrzNnCzcYTJlCoRDMdgJkxalYADIryHwlNTw9\nt29cK4TQ6XRWVlZ+N/379nrZo8lgThQKBb0OMMKIHQCZFdbTI+rXq3PC53ShbAoASiiKHQBz\n8NH4CX43/D3d3eUOguKj4Ao7IBuKHQATladTtJ1f6yCE6NihbZHFAYASgGIHwOQMcxug/1/u\nP/LO22+2fLXZJ198NXfRt4bnhz55+nTFT2uexsbeD3tQNEkhJ4VCIbh7AsiKmycAyG+Y2wDD\n+Fy+L7k7f+pYr4GDZ8xd2Lpli7i4+BM+f133u3ntH78pM+c8e5Yecut6tapuhRcZpoFeB2TF\niB0Ak5CPUTojKpVq0dyZQojtv+55e8Sovb8fuB/2YM+2TR3atpEkad/+A4UXFiaB2U6A7Bix\nA2A+mjZuZG1t/cv2nTbW1k8j7ukXDuzXp88bQxYu+X7IW2+4VHKWNyEAFClG7ACYlemTv9bp\ndKPfH5F54a4tG2LjYrfu3CVXKgAoHhQ7AGZl2iTvhIdhK7/7JvNCW1vbTq91WLVmvVqtlisY\nABQDih0Ac1OuXNnsC1cvX3o/7MHMeYuKPw8AFBuKHQCLUMPTc2C/Pus2bYmNi5M7CwAUFYod\nAEuxcfXK2LjYCROnBgbdkzsLCoFOJ/HwCcAIxQ6ApShXruyiOTN/2b5zwqT/yp0FhSAlJUWl\n4rcYkAXTnQCwIF99Mf5OYNCve3/T6XRKJZ2gZFNrNEqlSu4UgGnh5xoAy+Ls7JSamuZRt3HY\ng3C5s6BAJEniRCxghGIHwLLMmzF16kTviMioWfMXPwiPkDsOABQmih0Ay6JUKudMn9K2VctN\nW7d37N5H7jjIP0mSBEN2QFYUOwCWaPe2jUqlsm7tWnIHQf5JkqSg2QFZUewAWCL/W3d0Ot07\nb78pdxDkn0ajUXAHDJAVhwQAS9S1S6dXHBzGfu59+ep1ubMgnzQajVLJiB2QBcUOgIU6un9P\nWtqzf27clDsI8kmr1SkV/BYDsuCQAGChYp481ul0XGZXcimVCkmS5E4BmBaKHQALdeLkX1ZW\nVs2aNJY7SEkyb/GSz7+eov9z935v7Pntj+zr3A97sOSHlfnb/srVay9cuiKEuBMYNGv+4pes\neeHSlaU//KjT6YQQcXHxgwcPzt8eATNDsQNgoT4cNUKj0fywarXcQUqYoOB7QojU1LToRzEv\nWe37lT9FRj18+ZI8Wb1uo//tAMPLPw4dqVunlr7YCSEmTpyY7y0D5oRHigGwUHVr165dy+vX\nPb9N+WqC3FlKkpo1qgcF3wu6F9LptfZCiLPnL6zZsFmpVE6b5B10L2T95i1KpbLlq82fPI21\ns7M9+dfp9Zu32tnarvp+iX7JtNnzHSpWDAsP/2DksMTEJMNna9aort/+spWrXF1cYmPjqnt6\nXLx8ddPW7UlJyUsXzTtz7nxiUlKDenX1q4VHRLZ6tbnfzVv6l0tM9TafAAAddUlEQVSXLt2x\nY4csfyGASaHYwXTpdLrY2NiMjIxKlSqpVDwREoVv0dyZAwa/e873Yrs2reTOUmL06t718LE/\nox/FdGjbRgix/KefG9avl5iYtHnbzuhHj3Zt2ehz5u/rfv/ekrJ24y/bNqwJCb0fn5Bg2EKf\nnt2fpT87f+HSqdNnDJ+dO+O/+ncnjPukdctX7wQG7dy9r5Kz04fvj9i0dcflq9dea9e2fdvW\n+nVCQu/fvnPXzraMRqPRL/H29i6+vwLAhFHsYFokSTp9+vTatWv/+uuvmJgY/U9tpVLp6ura\nrl27sWPHduzYUaFgggMUjtI2NkKIsmXt5Q5SktSs7nnwyDEbG+vSpW2ePUu3s7WbMO4TjVYT\nE/Pkv7PmCCGUWeeWUygUSqVSo9Ealjg7O0ZERgkhMn/2uftavW5jk0YN69WpbbT8j0NHfl7x\n3fGTPhevXIuNiyvkbwiUZFxjBxOSlpbWs2fPzp07b9++XZKkZs2a9ezZs2fPni1atFAoFLt2\n7ercuXPfvn2fPXsmd1KYiRbNmwkh4uMTclwTmdnYWNes/u+Z048+eG/0p59NmTHHzs522JDB\nQ98bvW3nbtsyZfTvjnhnyIjRY5cu/7FihQqZt6BQKMqVK5v5s8/dkXu1qv63A27fuet78XIl\nZ6c/Dh3RLz934WLTxo2USqVKqbzhf1sIsXTp0qL6tkCJwr3iMCEzZsyYO3du9+7dFy5c2KRJ\nk8wjc5Ik3b59e968eTt37pwzZ8706dMLfe82NjaDBw3Ysp5L6S1IZNRDN6/6vj7HW7d8Ve4s\nyLNvli2fu/DbpYvmLV2+8m5gsNxxAJPAiB1MyPHjx2vXrn3gwIGmTZsanW9VKBT169fftm1b\nhw4djhw5IldCmJlbAXeEEAmJiXIHQX4olUpJMDYBZEGxgwm5fft2q1atSpUq9aIVlEpl+/bt\n/f39iyrAnbtXrl3X6XTDP/j46vV/Mr81efrs5OSUfG95yQ8rBw9/f8TosWM/9xZCqNXqiVNn\nfjR+Qr+33qFVyKhbl861vGqO/dxbrVbLnQV5plAoOOcEGKHYwYTUr1//0qVLhtvcspMk6cKF\nC/Xr1y/SGLPmLx7Uv0/zpk0Kd7MTJ3z2y7qfkpNTEhITf937W906tX5esWzqRO99+w8U7o6Q\nJ0f37w69H7Zr3+9yB7E4qalpf/199mH0o3xvQalQCpodkBXFDiake/fud+7c6d+//82bN7Nf\n/RkYGDhixAgfH58ePXoUXYYdu/fevnNnYL8+Qoiz5y+MGD32vTGfBt8L0b/7199n3xr23qCh\nI+6HPRj54ScZGRmVPGprtdphoz6aNnv+d8tXff71lBv+ty5evjr2c+9hoz56FPPYsOWly38c\n/clnySkp5cqWPf33uS6dOgohWrVo/v7wd4vu6yBHnu7uNatXHzbqo9D7YXJnsSBJyclN2772\neu+BNRo0vXj5av42wiPFgOyY7gQmZNKkSZcuXTp8+PDhw4erVq1arVo1BwcHhUIRFxcXERER\nGhoqhOjRo8ekSZOKLoMkSeXKlgsKvudVs0bmCbr0767d8MvW9T/fCQzasuNXB4eKf5/3bd60\n8Zlz5yu7uohMs3P16NrFMPlWn57d9Z/1/uzTBvXq9X1zSNiDcI1G6/jKK/rd3QsJNUzNClm8\n8/YbcxZ+m5iUJHcQC+J78XJ4RGTCw7BBQ4fv3LOvVYvm+dgI19gB2TFiBxNSpkyZQ4cOnTp1\nasiQIWq12tfX98CBA3/88ce5c+fS0tLeeuutEydOHD58uHTp0kWX4Z3Bb86d8d9pc+aL/02y\nNXWS9/Chb+vf1f8WUSqVOp2uU4d2S3/4ceKXn327bEWXzh2FEM7OjvoZvFav2xgYdC/75Ful\nS9v0691TCNGmVYvjJ08JIa773Vj+05qi+zrIjSlfTVAqFFeu/ZPzqigkXjWqa7W6ER+OvXDp\nSsP69fK3EaVSyYAdYIQRO5gWhULRuXPnzp07CyG0Wu3jx48lSXJ2di7gkyfGjx/v6+v78nXU\nanVKaqoQokpl16aNGx08ckw/yVbFChWmTfpKv86H748Y8eFYtVqzdOFcB4eKk6bN6ti+3YjR\nYzu0bfP3OV/xv9m59JNvJSYlhUdEGkbsli7/0dnJqVOHdh7u1Ua+O8R7yvTjJ33CIyJ/XPZt\nQb4aCi46JkYnSfZ2dnIHsSCeHu6Hf/v11z2/LZo7871hQ/O3EaVSKRixA7LiAgWYOq1WGxwc\nrNFoatWq9ZIbZl/u4MGDt27devk6CxcuHP/xh4bnGsGilHNxHzv6/cXzZskdBHnw8/pNn301\necXSxcxjBxgwYgcTMm3aNE9Pzw8++ED/MiMjY/HixQsXLkxLSxNClCpVatiwYUuWLHFwcMjr\nlvv06dOnT5+Xr7Nhw4aqblXyERslXcMW7ZKSkhwdX5E7CPJGoVAwYgcY4Ro7mJD58+dv2bLF\n8HLcuHEzZsywsbEZNGjQmDFjmjdvvnHjxtatW6empsoYEuYnMCh42qSvvD/7VO4gyBuFQkGv\nA4xQ7GCi/Pz81q5d26JFi8DAwL179/7888++vr6bNm0KCgpasGCB3On+n1qt9rvpzyTDJZpa\no2ncsIHRo+sBoCTiBxlMlP5eh2XLljk5ORkWjhw5sm3btkePHpUvVxapqWmtOnZt0vo1jzqN\n/W4W1fMwUNSsS5WKeviw2Ha35IeVV65df9G702bPj09I0P953JcTX7Kd+2EPlvywMsfdXbh0\npXKNujqdTmR7Aoo5UOS8CmBRKHYwUVFRUUKIRo0aGS1v2LBhQECAHIme49yFiwF370aH3m3a\nuNEv23fKHQf5pJOkuPh4uVOI1es2+t9+2b9trVb7zbLlkVF566B/HDoysF+fy1f/rZKZn4CS\n/6ymgZv/gOy4eQImysvLSwgRGhpq1O2io6PzcfNEEXGvWlWj0f535pwb/rfeGtRf7jjIp769\neny/cvXECZ+XKVOEUyQ+19nzF9Zs2KxUKqdN8j5z7rx+kuQZcxc+efp06kRvtVo99nPv+2EP\nFsyevvf3A/fDHqg16jGjRh45fmL95i1KpbLlq81PnT6zZsNmtVrT+bX2TRo1NGzNMOt1eETk\nlK8n7Ni1Vz8J8NLlP9rZ2uqfgFLMXxZAMWDEDqYlMDBw7ty5u3fvrl27tqOj45w5c/SnkPTO\nnz9/8ODBli1bypgws1peNXZv3ZiUnOz9+adjRo2UOw7y6YtPP45PSDh87M/i3/Xyn372qlnd\nyfGVzdt2vtauba/uXYUQH49+/8vxn/hevJyRoZ41dfL6n1as3fiLEGLyV1941aghhNi3/8Cu\nLRs/+uA9IcS+/Qe3b1w78t0hRlvTbz8k9P7tO3c3b92hnxBbCOH92acrv/s2OTk57EF48X9f\nAEWNETuYEDc3t4iIiBkzZhiW7N279+jRo7169RJCjB8/fs2aNVZWVtOmTZMvo7EBfXsP6Ntb\n7hQokA7t2rhVqbxu0y9vDOhbzLvWP91Eo9XExDw5dfqMfmFlV5f7YQ8yr1bJ2Umt1rhVqax/\nqT9xrL/bQ6PRGP6ceWv6Nf84dOTnFd+92qzphElTH4RH6BcanoACwPxQ7GBCwsPDU1NTg4OD\ng/4nMDDQMCnx7t27a9WqtW7duqZNm8qbE+anYf16gcH3im13S5f/WMnZuU4tr8xPN6nk7PTH\noSOZV6tS2fW/M+c8eRr7zfzZW7b/alg+bMjgoe+NLmtv36xJ44H9+gz74CMrlVW7Nq2yPyvl\n3IWL48eOEUK83rnjwSPHRNYnoBTb9y1K3D0BZMGTJ1BiBAcH16xZs+i2X7t2be/xYzmjapnG\nfu79+4FDD0PuyB0kz3bt/b1Vi+Y79+xr26plh3Zt5I5TrHjyBJAd19ihxCjSVgcLp1KpnjyN\nTU9PlztInr3+n46Hj/1Z3cOjfdvWcmcBID9OxQKA2PDLth5du9jY2MgdJM8cKlYc++EouVMA\nMBWM2AGA6Nal87V//EriiB0AZEaxAwDR8tVmUQ+jn8bGyR0kZ8dOnPp+5U/hEZFyB5Ef14gD\n2VHsAFi6w8f/nDZ7/qD+fSu7usidJQebtm7vP/jdVWs3tHytS1JystxxAJgcih0AS3fm73NC\niC3rVssdJGc+p88OeXPQ9fOnox/F3A64K3ccACaHYgfA0rVp3UqSpDuBgXIHyVnXLp2279rT\nsGW7KpVd69etI3ccACaHu2IBWLS4+IQzf58XQtwJDGrWpLHccXIwbMhgtyqV79wNGtC3t729\nndxxAJgcih0Ay3X8pM8773/49Gmse7WqnV/rIHecXOnUoX2nDu3lTgHARHEqFoDl+vDTzyuU\nL+9z9I97/tdcXSrJHQcACopiB8ByVXWrIoTo1KG9SqWSOwsAFAKKHQDLNeb9kSGh9/8+56t/\nef7CpVsBJe9xsQBgQLEDYLlGvDukuqfH630Gtv1P9/Hek9p16dG+S0+ePwGg5KLYAbBoh/f9\nWtbePig4ZNWa9TWrV49PSHjyNFbuUMgthdwBAFNDsQNg0Wp51XwSHvz4QZA26UmDenXK2tuX\nK1dW7lDIFR4pBmRHsQOAfx0/5fPpR6PL2tvLHQQA8oliBwD/0ukkJj0BUKJR7ABACCF0Ot2z\nZ8+cnZzkDgIA+UexAwAhhAiPjBRCMGIHoESj2AGAEEJotVohhCS4Hr/EkCSJ22IBIxQ7ABBC\niOoeHkqFIuxBuNxBACD/rOQOAACmwtPTY/HSH7Ra3ZsD+5Ury6QnJQJDdkAWjNgBwL9OHv79\nXuj9D8aOnzlvkdxZACA/KHYA8C/3qlXT46IbN2zw/cqfZsxdKHccAMgzih0AZHHd93TPbq//\nvH4TDzYAUOJQ7AAgC4VCER4RGRsXFxh8T+4seBmaN5AdxQ4AjC1dOM/O1nboyNFyBwGAvKHY\nAYCxbq93/nj0+3cCAx9GP5I7CwDkAcUOAJ7j1eZN09KeVa5Rd+zn3nJnAYDcotgBwHO8OaBf\nZPBt2zJl/jh0hGu5AJQUFDsAeL7Kri46nW5A394KBbPgAigZKHYA8EI1a9ZYtWZ90zYdu/Qa\nEHo/LCExMeBuoNyhAOCFKHYA8EJ+F87Mmzk18uHDC5evNG3T0b12o3rNWp+/cCn7mlsjfi/+\neABghGIHAC+kVCqnTvSOuR94959L1T093KtVdXZyfK1b70VLvxdCbI343fA/w0u5I1sQSZI4\nRw4YsZI7AACUAG5VKl87/5f+z0NHjvYrH7j85kaHihXlzAQA2VDsACBvmjRueF26o1KpDEs0\nGZr10zfd8wspbVs66tuHlWu4CiGGuQ2QLyMAC8WpWADIm/Eff7j7y92JiUmGJT67TiuUikWH\n5r399Ztb5m8XtDoAMqHYAUDe2NraWlmpEhISDUsig6MatK2vVCnrtqwTevM+894BkAvFDgDy\nbM70qYsHfXvnbqBarRZCeDbwOH/gwtOHsSe2nUpJTElPTecuCgCy4Bo7AMizSV9+FhcXt6D/\n4k9++bh1yxYdBraLj4lfO2VDnZa1Has42tjaiKwToHBmFkDxoNgBQH4smjuzdasWA98eVu9y\nnfioeKeqTv0/6XvPLyTkRghPqigekiQJ/qaBrCh2wL8kSUpNTY2Lj5c7iCXKyFCnpKTkZk17\ne7tSpUoVdZ5cqlC+vBAiOSXF2d3p8IajJ7afsiljM2LGu/rzs5nx76oopD17JkkiNTVVq9Vp\nNBorK36jAULBRb6AXrVq1cLDw+VOgRLm7e/fzs1qv37xa1EnsXDnz59v06aN3CkA+VHsgH95\neXkNHjx40KBBcgcpNBVzN31uenp6ampqUYd5OSsrq7Jly+bpI/Hx8bL/+NLpdE+fPnVycpI3\nhsXSaDRJSUk+Pj5bt269ffu2nZ2d3IkA+TFwDfxLqVS6u7s3b95c7iAA8iAkJKR06dK0OkCP\n6U4AAADMBMUOAADATFDsAAAAzAQ3TwD/atiwob+/v9wpAORZgwYNbt68KXcKwCRQ7IB/hYeH\nx8TEtG7d+ssvv2zYsKHccfJMrVaPGjVq1qxZNWrUkDtLniUkJIwbN+6bb75xdXWVO0uePXz4\ncOLEiStWrKhQoYLcWfIsJCRk5syZ69evt7a2ljtLnvn7+y9ZsuTixYvOzs5Vq1aVOw5gEih2\nQBY2NjYHDhzo1q2b3EHyLC0tzdbW9sKFC61atZI7S55FR0e7uroGBATUqVNH7ix5dvfu3Tp1\n6kRFRZXEVnrp0qVWrVqlpKTY2trKnSXP/vzzz969e2dkZMgdBDAhXGMHAABgJih2AAAAZoJi\nBwAAYCYodgAAAGaCYgcAAGAmKHYAAABmgmIHAABgJih2AAAAZoJiBwAAYCYodkAW1tbWJfHZ\nSkIIlUqlUqlKaPhSpUopFIoSGt7a2lqhUJQqVUruIPlhbW2t/5cjd5D8KLlHK1B0eKQYkEVo\naKiHh4dCoZA7SH6EhIRUr15d7hT5RHi5lNzwkiTdv3/f09NT7iCACaHYAQAAmAlOxQIAAJgJ\nih0AAICZoNgBAACYCYodAACAmaDYAQAAmAmKHQAAgJmg2AEAAJgJih0AAICZoNgBAACYCYod\nAACAmaDYAQAAmAmKHQAAgJmg2AEAAJgJih0AAICZoNgBAACYCYodLNeqVasqVKiQfXl8fPyX\nX37ZsGFDOzs7Ly+v4cOHh4aGZl7h2bNnc+bMadu2bbly5dq2bTt79uxnz54VV+p/vSh8jivE\nxcV5e3vXq1fP1ta2Xr16X331VXx8fJHFzFu23Kywe/fudu3a2dvbu7q6Dh48ODg4uGgy5jlY\nLlcQQmzatEmhUBw8eLBQo+Us3+FzPCKKQb7Dm8LRChQrCbBIycnJdevWLV++fPblnp6eQohW\nrVp99NFHXbt2FUKULl368uXL+hV0Ol2vXr2EELVr1x4+fLiXl5cQokePHjqdTvbwOa6QmJhY\nq1YtIUSbNm1Gjx7dpk0b/RdJTk4u4sg5Z8vNCvPnzxdCuLi4DB06tG/fvkql0sHB4f79+0WZ\nN1fBcrmCJEkBAQG2trZCiAMHDhRBxhfKd/gcj4hikO/wpnC0AsWMYgeLc/To0UWLFtWuXVsI\nkf03wfTp04UQX3/9teFH/4EDBxQKRYMGDfQvT506JYTo1auXWq2WJCkjI6Nbt25CCB8fH9nD\n57jCnDlzhBCzZs3SfzudTqf/vgsXLjT98OHh4SqVqkWLFnFxcfol+/btE0K89957Jp7cIDU1\ntVGjRvr/qC62YlfA8DkeEUWqgOHlPVoBWVDsYHFsbGwMI9bZfxO0bt3a2traaASrS5cuQojo\n6GhJkoYMGSKEuHHjhuHdq1evCiHeffdd2cPnuELv3r2FEFFRUYYlYWFhQogBAwYUbe5cZMtx\nhWnTpgkhzp07Z1ii0+m++eabZcuWFW3uAic3+Pjjj8uUKTN8+PDiLHYFDJ/jEVGkChhe3qMV\nkAXX2MHixMfHp6WlpaWl1a1bN/u7qampnTp1srOzy7zQ2tpaCBEXFyeEOH36tJubW4MGDQzv\nNm3a1NXV1cfHp4iDC5FT+BxX0P8WfPTokWFJdHS0YXlRK2D4PXv2uLm56U8f6ykUiq+//vqL\nL74oqsS5C5abFYQQe/bsWb169YoVK+rVq1dkSfOTrYBHRJEqYHh5j1ZAFhQ7WJzS/6NUPuff\nv5+f37FjxzIviY6OPnnypJOTU40aNbRa7aNHj9zd3RUKhWEFhUJRtWrVmJgYnU4nb/gcV/jy\nyy/t7e1HjRp15cqV1NTUy5cvjxkzxt7e/vPPPy/i4Dlny3GFiIgIDw+PGzdu9OvXr1KlSm5u\nboMGDQoICCji1DkHy80KoaGho0ePfvvtt0eNGlWUSfOTrSBHRJEkznW2l68g+9EKyIJiB7zM\nnTt32rZtm5GRsWDBglKlSj1+/Fin0zk4OBit9sorr2g0mtjYWFlC5l67du2OHz9+48aNFi1a\n2NnZtWzZ8vbt26dOnco8DGaaEhMTk5OTo6Ki2rdvHxYW1qdPn4YNG+7fv79Jkybnz5+XO10O\nMjIyhgwZ4uDg8PPPP2cuGSWR0REhd5yXKelHK5A/VnIHAIpESkrK2rVrDS9r1qzZp0+fPG0h\nISFh0aJF33//vVar/e6770aPHi2EkCRJCPGi380ZGRkFiPz/Ch7+Rfz9/d955x0rK6vBgwd7\neHiEhobu27dv6NChBw8erFOnTqHsoojC6+dkCQkJmTZt2uzZs/XDM6dOnXr99dfHjBnj7+9f\n8F0U3V/71KlTr127dvbs2fLlyxfKBrMruvAGzz0iCkURhS+eoxUwOXJf5AcUiYcPH2b+d/7G\nG29kX6d+/fovus79999/d3V1FUL07NnT39/fsFyj0SiVyg4dOhit36pVK6VSqdFoTCH8i1ZI\nT0/39PQsV65cQECAYaG/v7+9vb2Xl5f+tkGTDZ+WliaEcHR0NPpL1t/kWChX8RdR8hMnTggh\nFi9ebFiycOFCUdg3TxRReIMXHRGFoojCF8/RCpgaRuxgnlxcXCRJyt9np06dumDBAk9PTx8f\nn06dOmV+S6VSOTs7R0REZF4oSVJUVJSrq6tKpcp34MwKEv4l/Pz8QkNDR44cmXlwrn79+n37\n9t2xY8etW7caN25c8L0UUfjSpUtXrFjR3d3d6C+5evXqQojw8PBKlSoVcBdFlPz69etCiEmT\nJk2aNCnz8r59+woh1q5dWyhDX0UUXu8lR0ShKKLwxXO0AqaGYgdksWnTpgULFvTv33/z5s3P\nPXHWoUOH3bt3BwUF6Sc7FUIEBASEh4cPHjy4eJPmmX5efrVabbRco9EIIcqVKydDprxo1qzZ\npUuXnj17Vrp0acNC/c0T+mnMTFPjxo0/+uijzEuuX79+6dKlHj16uLu7F9YZ8KKT4xFhykru\n0Qrkn3yDhYDMsp+70Wq1tWrVsre3j42NfdGn9GfWhg8fbpjjd+jQoaLYpzzNx5kpnU5XrVq1\nMmXKZH5mgK+vr42NTfXq1YtzLv78nRPcu3evEGLcuHGGk2h79uwRQnTr1q2oguYuWJ5WkIrm\nVGxu5CN8bo6I4pG/v3kTOVqB4sSIHfD/wsLCAgMDHRwc3nrrrezv7tixw8nJqXPnzt27d9+y\nZUtkZGTr1q3PnTt3+vTpXr16dezYsfgD54lCodi+fft//vOf1q1b9+zZs1q1aiEhIcePHy9V\nqtS2bdtM/27Nfv36tWnTZuXKladPn27Tpk1oaOiff/7p5OS0evVquaOZrdwcEcWfKvdK7tEK\n5BvFDvh/ISEhQojY2NiTJ09mfzc9PV0IoVQqf/vtt0WLFh07dkw/2eysWbMmTZpk+sVICNGu\nXbuAgIA5c+b4+vqeOHHC3d19xIgRM2fO9PDwkDtazqysrI4dO7Zo0aJTp05t27atatWqY8aM\nmTdvnol3ixItN0eEKSvRRyuQPwqpyK63BQAAQHFigmIAAAAzQbEDAAAwExQ7AAAAM0GxAwAA\nMBMUOwAAADNBsQMAADATFDsAAAAzQbEDAAAwExQ7AAAAM0GxAwAAMBMUOwAAADNBsQMAADAT\nFDsAAAAzQbEDAAAwExQ7AAAAM0GxAwAAMBMUOwAAADNBsQMAADATFDsAAAAzQbEDAAAwExQ7\nAAAAM0GxAwAAMBMUOwAAADNBsQMAADATFDsAAAAzQbEDAAAwExQ7AAAAM0GxAwAAMBMUOwAA\nADNBsQMAADATFDsAAAAzQbEDAAAwExQ7AAAAM0GxAwAAMBMUOwAAADNBsQMAADATFDsAAAAz\nQbEDAAAwExQ7AAAAM0GxAwAAMBMUOwAAADNBsQMAADATFDsAAAAzQbEDAAAwExQ7AAAAM0Gx\nAwAAMBMUOwAAADNBsQMAADATFDsAAAAzQbEDAAAwExQ7AAAAM0GxAwAAMBMUOwAAADNBsQMA\nADATFDsAAAAzQbEDAAAwExQ7AAAAM0GxAwAAMBMUOwAAADNBsQMAADATFDsAAAAzQbEDAAAw\nExQ7AAAAM0GxAwAAMBMUOwAAADNBsQMAADATFDsAAAAzQbEDAAAwExQ7AAAAM0GxAwAAMBMU\nOwAAADNBsQMAADATFDsAAAAzQbEDAAAwE/8HEFloKe+EH1cAAAAASUVORK5CYII=", "text/plain": [ "Plot with title “Les parcs nationaux du Canada: Alberta”" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "map(database = gpr_000b11a_f,\n", " fill = TRUE,\n", " col = 'seashell',\n", " namefield = 'PRFNOM',\n", " xlim=c(bbox_prov_alberta[1,1],bbox_prov_alberta[1,2]),\n", " ylim=c(bbox_prov_alberta[2,1],bbox_prov_alberta[2,2]))\n", "title(\"Les parcs nationaux du Canada: Alberta\")\n", "map.axes()\n", "map.cities(x = canada.cities,capitals=2,cex = 0.5)\n", "map.cities(x = canada.cities,capitals=0,cex = 0.5,minpop = 50000)\n", "plot(limites_parcs_nationaux,\n", " add = TRUE,\n", " col = alpha(\"palegreen3\",.7),\n", " border = FALSE,)\n", "invisible(text(coordinates(gpr_000b11a_f),\n", " labels=as.character(gpr_000b11a_f$PRFNOM),\n", " cex=.6))\n", "invisible(text(coordinates(limites_parcs_nationaux),\n", " labels=as.character(limites_parcs_nationaux$OBJECTID),\n", " cex=.5))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Un premier indicateur: distance d'un parc national" ] }, { "cell_type": "code", "execution_count": 79, "metadata": {}, "outputs": [], "source": [ "onePoint <- SpatialPoints(coordinates(readWKT(\"POINT(-114.0 50)\")),proj4string = CRS(projargs = proj4string(gpr_000b11a_f)))" ] }, { "cell_type": "code", "execution_count": 80, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\t\n", "\n", "
distancelonlatID
88554.93 -113.971849.20471 9
\n" ], "text/latex": [ "\\begin{tabular}{llll}\n", " distance & lon & lat & ID\\\\\n", "\\hline\n", "\t 88554.93 & -113.9718 & 49.20471 & 9 \\\\\n", "\\end{tabular}\n" ], "text/markdown": [ "\n", "distance | lon | lat | ID | \n", "|---|\n", "| 88554.93 | -113.9718 | 49.20471 | 9 | \n", "\n", "\n" ], "text/plain": [ " distance lon lat ID\n", "[1,] 88554.93 -113.9718 49.20471 9 " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "dist2Line(onePoint, limites_parcs_nationaux)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Ressources\n", "\n", "- [Making Maps with R](http://www.molecularecologist.com/2012/09/making-maps-with-r/)\n", "- [Topic 2(a): Spatial Data in R](https://science.nature.nps.gov/im/datamgmt/statistics/r/advanced/spatial.cfm)" ] } ], "metadata": { "kernelspec": { "display_name": "R", "language": "R", "name": "ir" }, "language_info": { "codemirror_mode": "r", "file_extension": ".r", "mimetype": "text/x-r-source", "name": "R", "pygments_lexer": "r", "version": "3.3.3" } }, "nbformat": 4, "nbformat_minor": 2 }