{ "cells": [ { "cell_type": "markdown", "id": "731499a9", "metadata": {}, "source": [ "# BLS MSA stats\n", "\n", "By Andy Wheeler\n", "\n", "See blog post, [Use circles instead of choropleth for MSAs](https://andrewpwheeler.com/2022/03/13/use-circles-instead-of-choropleth-for-msas/) for overview.\n", "\n", "This scrapes the data from BLS website for particular employment codes." ] }, { "cell_type": "code", "execution_count": 1, "id": "above-transformation", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
OCodeDesc
0110000Management Occupations
1111011Chief Executives
2111021General and Operations Managers
3111031Legislators
4112011Advertising and Promotions Managers
.........
847537072Pump Operators, Except Wellhead Pumpers
848537073Wellhead Pumpers
849537081Refuse and Recyclable Material Collectors
850537121Tank Car, Truck, and Ship Loaders
851537199Material Moving Workers, All Other
\n", "

852 rows × 2 columns

\n", "
" ], "text/plain": [ " OCode Desc\n", "0 110000 Management Occupations\n", "1 111011 Chief Executives\n", "2 111021 General and Operations Managers\n", "3 111031 Legislators\n", "4 112011 Advertising and Promotions Managers\n", ".. ... ...\n", "847 537072 Pump Operators, Except Wellhead Pumpers\n", "848 537073 Wellhead Pumpers\n", "849 537081 Refuse and Recyclable Material Collectors\n", "850 537121 Tank Car, Truck, and Ship Loaders\n", "851 537199 Material Moving Workers, All Other\n", "\n", "[852 rows x 2 columns]" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Creating example bls maps\n", "from bls_geo import *\n", "\n", "# can check out https://www.bls.gov/oes/current/oes_stru.htm\n", "# or use this function\n", "os_cats = ocodes()\n", "os_cats" ] }, { "cell_type": "code", "execution_count": 2, "id": "04f48075", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVYAAADnCAYAAABFTn7nAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAABYi0lEQVR4nO2dd3yb1fWHnyuP2LHj7DgLskgICTNhhBX2KnuWDnYpq0DZ8INyORRaVoFCyyq7ZY8CZZcVRkIYhQBJCCRAyN7Die3Ylu7vj3MVy0K2ZVt25OQ++ehj69V9hyL5vPee8T3GOUcgEAgEMkdkXV9AIBAIrG8EwxoIBAIZJhjWQCAQyDDBsAYCgUCGCYY1EAgEMkzuur6AQCAQiLPf/lu6JYtXpTX200+/f805t38rX1KzCIY1EAhkDYsXl/HhR5LW2PycE3q08uU0m2BYA4FAFuFwrmZdX0SLCYY1EAhkDw6ci67rq2gxwbAGAoGsweGIhRlrIBAIZJIYsVjlur6IFhMMayAQyB6cw8XCjDUQCAQyS3AFBAKBQCYJWQGBQCCQYRzEqtf1RbSYYFgDgUDW4FyYsQYCgUBmcTGoCVkBgUAgkFnWgxlrULcKBAJZg8FhYjVpPRo9ljH3G2MWGmO+Stp+tjFmmjFmsjHmhoTtlxljpvvX9kvYPtoY86V/7TZjjGns3MGwBgKB7ME5iNWk92icB4E66lfGmD2AQ4EtnXMjgZv89hHAscBIv88dxpgcv9udwG+Bof7RqKJWMKyBQCCLyJxhdc69CyxN2nwGcJ1zbo0fs9BvPxR43Dm3xjn3PTAd2N4Y0wcocc5NcNp59WHgsMbOHQxrIBDIHpzDRNek9QB6GGM+SXj8No0zDAN2NcZMNMaMM8Zs57f3A2YljJvtt/Xzvydvb5AQvAoEAlmES8t/6lnsnNu2iSfIBboCY4DtgCeNMYOBVH5T18D2Rk8SCAQCWYKDWKvKBs4GnvXL+o+MMTGgh9++UcK4/sBcv71/iu0NElwBgUAge8hs8CoVzwF7AhhjhgH5wGLgBeBYY0wHY8wgNEj1kXNuHlBmjBnjswGOB55v7CRhxhoIBLIKk6EZqzHmMWB31Bc7G7DA/cD9PgWrCjjBz14nG2OeBKYANcBZrlZx+ww0w6AQeMU/GiQY1kAgkD24zLkCnHO/qOelX9cz/lrg2hTbPwE2b8q5g2ENBAJZhMPUVK3ri2gxwbAGAoGswTiXMVfAuiQY1kAgkF0EwxoIBAIZJMxYA4FAoBUIhjUQCAQyicPEYuv6IlpMMKyBQCB7cA5CVkAgEAhkFuPCjDUQCAQyRwYLBNYlwbAGAoHsIvhYA4FAIIM4FwxrIBAIZBaHqale1xfRYoJhDQQC2YMjzFgDgUAgswRXQCAQCGQWR8gKCAQCgUxi1pPKq9CaJRAIZBexWHqPRjDG3G+MWei7BSS/dqExxhljeiRsu8wYM90YM80Ys1/C9tHGmC/9a7f5Fi0NEgxrIBDIHpyDmpr0Ho3zILB/8kZjzEbAPsCPCdtGAMcCI/0+dxhjcvzLdwK/RftgDU11zGSCYQ0EAtmDA2IuvUdjh3LuXWBpipduAS6mbhvrQ4HHnXNrnHPfA9OB7Y0xfYAS59wE3xvrYeCwxs4dfKyBQCC7aEUfqzHmEGCOc25S0oq+H/BhwvPZflu1/z15e4MEwxoIBLKIJqVb9TDGfJLw/B7n3D31DTbGdAQuB/ZN9XLqi6l3e4Osd4ZVRPYDxgALgVestT+s2ysKBAJpE3cFpMdi59y2TTj6EGAQEJ+t9gf+Z4zZHp2JbpQwtj8w12/vn2J7g7Rbwyoi1wCvWmvfT9h2BHAFsAlwPLCviAwC1qDT/HetteXr4noDgUCatJJsoHPuS6BX/Lkx5gdgW+fcYmPMC8Cjxpibgb5okOoj51zUGFNmjBkDTETtyu2NnatdGlYRMcAOQKmIkGBctwS2Qe97V1trt/Tj89FZ7PkiUgCsBMYBn1pr0wovBgKBNsA5qMmMYTXGPAbsjroMZgPWOXdf6tO6ycaYJ4EpQA1wlnMuXqlwBpphUAi84h8Nn1sDXe0LEckBHgBOAK4GBqCRvq+B+Iz0CGvtv+vZvzOwGzAavbksBd4DPrPWtn8FiECgnbLtoIibaPPTGpt70ppPm+gKaDPaq2E1wDvW2t38826oYT0KKMgtr4lsNH7eUwUrq/8BvH/MpIcqGjleN2BXdLabixrnicBEa+2qVnwrgUAggdEDI27ilekZ1rxTgmHNOCJyu7X27KRtWwGfFyytZNA7a/3LFcDbwPXHTHro3TSPXQRsD5wNdAQ+A6YCk4AvrbXtv+YuEMhCRg+MuIlX5KU1Nu/Uqqw1rO3Sx+qpE4QSkVzgc4Du366IUVv8UIhWUxSme2Br7WrgbRE5Hp0Jf4DOaB8DnhKRCJrfNhn4FJhprW2fd6hAIJtwwHowbWnPhjU5v6wamJq3unqzkjmrkyvK7j1m0kOvNeMcmwDvW2udiGwH/MNaewusDYiNBPYCBooIQBT4BvgSmGatbf/tJgPrBBHZA9gc+AdwA5rZcoe19vt1emFtQfsXt2rXhvUiEbnXWvuNfz6tdNLizTou+ok7tQL4JHljuiTMRPcETknYXoW6CD6Lb/NBtWFodsLhfhYNOrue6h/fhUyEQEP4GMIRwCHAjsC9wATgZhG5Yb02rg5crFGNk6ynPRtWgN+gNb8AV3RYWf1Uwco6Qf15wJhjJj3040/2bAQRGQnMT9jUyVo7r6F9rLVRag1o4rGKgeFoitix3gAbdOEzF/gWrU2e448R2LAZBbwLXAnUWGvLAESXRT8H/roOr631Ca6AdYPPRQX4T3ybtfbpJ7c64WzqJu/mohVYzeF4IDFdq9m3UZ9Z8AlJM2c/M4knI+8L9Pf+2/gsuQKYharwzALmZpN7QUTyAEKKWsY5ELgpuZjFWjtfRHqvo2tqO8KMdZ1xMoC19r2k7fcAO6PyXwBnHjPpocpmnmMHdMYQZ6P6BjYX72aY4x/vJL8uIoX+vANQV0TfuDFLIAYsA6pQ/+4iYAmwxFq7JtPXnMQvgTNE5Bxr7UetfK4NAn9jLWigQnCRiPSy1jZ3wpDdBFfAOmUvYLmIPGutPSK+8ZhJD1UBv3hyqxMuA3qg6VHNJSfJMH0uIqdaa//RgmM2CWttBWosv6lvjP9D7Aq8BZwDbAp0B7qLSAc/LFlMwqEKPXPQWXG5/7km4VHtH4X+5wrUiLuE4/UHLkCrW4JhzQyl6OdSHy8ABwMpK4jaPybMWNch3wDvlJWtPuaII345ffORmxx09R+v/rogv383YA/gzcqq2T809+A+jzV5yb0dWqHVZoY1HXxO7RIRWWatHdeUfeOzIzRXtxDokPDI84+90aDdr1Fjmvit/xQ1qIe28G0EaqkhoZ49GWvtdBH5P9ZXw+rARdu/THR7NaydrbWXFeT3rwD+kROJPI2mphyM1vRGC/L7b1NZNfvLZh7/aDQKm8hpwD0i0t1au6S5F55NeKNcTlJOcCIiMgG42Fr7agNjvhGRe1GXxJtowszeaGDuEWvt4oxe+PrNQaiORUPsLiJHo6mADQZU2yWx9m9Y2907EJH+wPSC/P7FwIkAEyZMGlCQ3/9V4Ho/bCZJkfkmcjDwr8QN1tr/AH9BDXdW4X2x6yyAZK2911obz9AYDYwALkVns9NF5AgRSbtCRkRyReRnrXO1WU8RWlrdELujlYEvtfrVtDXOuwLSeWQx7XHGegFq4LZCE/RZvGhZMbBfwpg/VFbNbkmuaC9r7dfJG621N4nIuSLSz1rbkB+srRmGBq0yjrU26osf0hnrgGsTNk0QkW1QgZuPRaQKONWPORj9HB06QytDb4jbogb5b8DLGXob7YlPgN81NMBa+yOax/2WiNyim+zKNrm6NsC57Daa6dCuDGs8FenP1/5jM1TKq0vi67m5ObFDD99zcpcuJZeLyI7AdcAAa+14v38OajTrXT75MTn1vQ48i7oJNm7Rm8ksg9Cc3dairLk7+mT2/gAi8ju0RPgGIK57ORuVauuIfl6Xognx57TsktstJ6E3lXQ4HpXDfElExq43ZdXrgSugXRlWYL+amug4NPJdBDyDJvGfBVBTE40889R/f15ZNXuqiDwCPAJUi0hHNKKdA1T6aHkEDRSsRGdG/7LWrgDygeIGruF84AgRWWCtLW2Vd9l0NgFasxrHiMgAa+3MlhzEWrvWYPgc3rHAmcB5aDCmAPi7tXZ26iOs3/hKvaOstaelM95aO1tE3kRvVF1J3TivXeFCutU64ejc3JxTKqtmO/xdvSC//9F4w+qJJ63/Kp0Dikg/VNf1aV94YIDO9Y33S+OdgPHNewutwkBa19/2PHoTy4iSkPcJ/wn4d9LntNo/4r70G4BfANdaa6/IxLmzFX/zf4mfZqM0tE8H1JhOtda2e6OqmJAV0Jb4ipNpKZY7z6J/pMegPrrvmnJc7yv9k3/EU62+FJFO8VLCFLwJnNuU87Qy/YBprXj8vwD/y8SB/Ez1r8A13leYaswg4BJgBnC0tfbpTJw7W/Ez1SVoALKxwBUi0h24DS3QWGmtHdG6V9iGONYLV0B7egcnoD2961BZNTtaWTX78sqq2UOBIyurZrdImNpLBl5IQrlsCvpba59tyXkyTFda1xXweSYO4o3qTcCjwCwR2UFEBiSNGYhWvN0KbA386Pdr6Lj9fFAxPYXkLEJEegG/RWeqr6KfZX1jDxeRfwKLUaP6pbW23tVVe8XFTFqPxjDG3G+MWWiM+Sph243GmK+NMV8YY/5tjOmS8Nplxpjpxphpxpj9EraPNsZ86V+7zST1zU5Fu5ix+qBVz8Zy9iqrZq/I0Cm3pp7qF+8GyDqZiNYMXHjZxBYdwwcF70bdCqvR2WhfNAj5gM8eGIP6iy8A/ozmdE5G84fHJka+RWR/4Eg0aFeOCpfcJyJbokbnD/4mmVV4F8eeqD6EQY3kW2gO8J/R9sxH1bP734G4VsDieE+39QuTyayAB1GXYeKE7L/AZc65GmPM9cBlwCXGmBFoKfxI9Hv5hjFmmO97dSd68/sQjcfsTyN9r9qFYUVTqepNUM8k3ud6oLV2dD1DLgYWiMgs4CNr7ZFtcV3tGa9vMBXtT7YH2oY4ip+lichleANjrb3T73MjcKG1tkxE7gd+EJGLUMN7JCDAGXEJRhE5DLgRldl735/jizZ7k/Xg+6vtguadbu83nwf8M34z9O6nAdbaz0SkTwOHi78WY32tdsugK8A5964xZmDSttcTnn5I7U3sUOBx59wa4HtjzHRge9/JtcQ5NwHAGPMwcBjt3bD6ZeB+6JexLRiHznjq43A0cPVn4J9tckXtn+dRrYEVwCXW2moReQDN0pgLXJqswWCtnZ7w+zciMgoYjKbY3WWtfSRp/HPAc/77Mg04TURiaFPJNk2k9y6Jq4FOwA/o9+WPaNrZ4ylypPuh+gyg5cnDEnSG4yu2uI/7ITRz4uPWewfrDgfE0g9e9TDGJCrG3eOcu6cJpzsZeML/3g81tHFm+23V/vfk7Q2S9YYV9Xc+3hY5eiLyMPBVQ0pN/jp29OM3qFJNEdnIWjurGbvOA/6U0KYca+1J/pgXAdNEZHe0kODpRKOSMP4H1Ei9lXA9BWjQ8jJ0Bvw2P02u/yVtWKHkDfttwN+stV8lvfyDiHQTkQ5JAj/jUfcA6LLzQnTpiYhciS5RN0NXS3evT8UAP8GZpsxYFze355Ux5nI03TJ+g07lf0gWL0rc3iBZbVhFZDiqMvVho4Nbfq4/+3Md1oTdvhaRS62117XSZaXLIBE5Hc2I+A74vhUEsz9CW850aWjQXcNP2xQ49vSv7050yq6inhQ2a+0CoIuI9ESXZdeLyCbo7HQq8Gdr7TOw1midjIqGb4QWLkxC3Qg3owUGi1Gf63CvDtbWHA88lsKoxvkI2Am9CcTpGC9isda+IiKJqWXx/8ft19dZajKtncdqjDkB9d/v5Wq7qc6mrjRof3Q1Ndv/nry9QbLWsPo/otOB/2ujUx5P3f/AdPgCDbisM0SkE2pgugFboMGNrn752BDlaHHESr//Kv9zBRpcKgcqUTnBKNpI8cQ0LmkQcMFdw0+7+vSv73Yi8ivgOBpRvbfWLkJna3f6z340KqzzhIj81Vp7AVpQcK/f5Ttr7RARORn1Yf4eDSzchP7RPCoibwAvt1UrE/9ZDLPWPtTAsM74XN0EkkV9Kn1KlQPKrbVFGbzMrKc1S1qNMfujqXy7OecSxYdeAB41xsQrAocCHznnosaYMmPMGGAiaiduTz5uMllrWNHoXLQBwd9M8zkaMdy7Cfv8CZgiIr+31t6aaoCvjClHe2N9AkzwRiRTDAM+s9b+Kd0dvOHqhhrhHmiKT2f/6IWWl3ZEq9Dy0O/J5qSXK/s18MnpX9/trr7C9iTX/AFjNmpKhN67Wz4BPhGRx9Ak+AuAX6Hv9xlUKwKgp7/mUWiGwQXW2rtF5B60K8NTIvK4tfamdM/fHLxf9Ro0ct8Qr6ISjB/5/UrR7hCJvIXqKXzChubHb5oroEGMMY+hgjU9jDGzAYu6jToA//VZUx865053zk02xjyJllfXAGf5jADQ79WDqLTmKzQSuILsNqxLaNvrW4LO0NLGR6wvRFX0x6EG6mO//Va0RDMH9fvtgfrRzvBVNqD+m5Xokvcz1Hn+fRP9yYNJY2mSdN0O32UgnfF+mT4R2KehcXcNPy0f/dLuftfw09b0hNwlB/X68Iqb/tjstCdr7RoR6eBLlL9HOy1YnwKWCzyOBja3QeUkY34/B7wmIjOBqSLyXGJArBW4BpiYyj+cxF1ogOoCqy1tdkR9x2vZautOD0z6vOxzdNa+/U+OsJ6TKVeAc+4XKTbXq2PrnLuWuiJC8e2foBOLtMlmw3oqtTKAbcEuaICgqbyIOsDvQ5cQVT59phtqqB/zNfYPptrZp9fsgs64jkaX8VDbbHAlmvP5Iyrw/QUwO8H4DvavtSaL/XUcCdzfwLidqQ0u5Rug5KPlW4lIZ6/D0FwuRVOpPkL/j+N/HGeihulYNOiVqhXN936/1nYH/BJVnHqsoRujL4legi4nT0fFwi+O8kgHNItgxcGHlj68x17dP5gxvXzviROWjwTq89eudzjXpKyArCWbDesA6nZJbTW8Ify+nj/MBvFL3JT/jyJyH7o0bmj/ecBT/pHqGBuhM5fB6NL2t0DnhIT9zWg4PazF+NnhNqiqV0OG9adlvjFc6WNzf3fXY6dNPv3ru59r5iW8AXyMzt5uSjBcU4Hd/O8pyzr9jPdJNEWmNW9AEo1Gb7rmj9ftI1dd+0HMVf1klp6f13sv4Kt+/Uq3OvGkI+KZCrOBElRY6GNgnnOutLg4t6BXr/yV24wuyea/0VbA4FwwrK3JVNR/+R80daVFykqN8HN+2jEgE3ShbqfXJuPTm+pNcfIBmvfrez1TeOPamKtkJRroWiu7mL+4qhhdJr8IPNfM05+G+nzBazp4DkRns9eg7cPrYwENK5Y1i0ik0AADcK5/To45MxaLdQJeQxtRpkrZuxgYPWfOgt/h03j8/2sxMGXNmtgHHTpE9ohXTPbpW1CyenX0I9TXveEQ1K1alaFoxP1A4A5fwbICdeY/leF0ooPRSGGm6dHKfj2A3ObMtFuD07+++/i7hp82Fl1tJNOsyHZCdsiBwG7W2mUJL5dZa79N4zAXAWNFxDQ3HzonUpTjcJuC2wr1526DriS6YszSfv37Lvtx5qy4RVhQz2H+ha46Tn3llXeHgpzsn/f+xTFPznz9tU/6brFlH2646XC232EgALl5a9ZrAZpUBKHr1uVbtA1wXFc1LtDxO+As/wf3IXCLbbmaf/c0gg6Berhr+GkFaFL8UGpT1lYBxWt6dxjvcszWuWU1DYmH14uf0ZWi2QlDkl5ObgX+E0Tk58AHaCR+PxGZaK39SYCiIXIiRRsBj0SM6RVzbtMUQ7pVram+Aq2u6g68EDH5F8Vc1esA+Xm9DRrUi4uj7/nVl99UHnDA2IOBq6y1x+666+iHy8ujx038cDY33vBfjvn51nz+2Y8cd8LYo5Yuf/Dje++ZdWFFRew5a+1ZKc6//rCe6LFmszOjJ0ntRqy1P1hrL7TW7grsipb53Ssi74vIf0TkWC/20VTa/ye5bqlCA2vPoS6VueisdYsO89fssmJst5IlB/aqFJFm6T1YaxeiGRdrJwIish2NpH+J6pWeiKZmnW6tPQSo8TfotOnYseN5wK5gBqOujmrqVt/Mmb9g0RRql/9bAs9HTP5f/PNcNPJ/jX++fMgmG08G7r72mju/y8/rndehQ+69cc2k/772Def87lHmzlvO8OE984qKc27eelTJjsCZyWpg6xvO+1jTeWQzWTljFZExwJKGlvt+SfeofyAig9Eo8Vki0g2NpN8L/Keh5Z+IbI/quGYUH3Raf0sPEzj967tjaGI+dw0/bQmw5+lf372UWkX7KDpbfFBE7rHW/rYZpzkK1WeIMxetYGqIXDTlLQ/NZwQtb92dJjSFPP2Mk8bdcvMdZzvn8lCxjlfRKrRhfsjEWKxiXMTkf4vmQc9CU9m+BKiqnl+dn9f7C9TAXwY8dsQR+96OFgpcnJeXu9eHE+Zu2bVrUeXSpasLqmsq3erFFeZfD0/gzLP2JDenGxM+WP5ewv9D3GCvf4SsgFblAJrYN91a+x1aY42IvIMuy05AU2BAI8LPAi/4/ME4v6CFAaZ62AdNpdmgOP3ru/9FUofbONbaE0XkDyIyGdipiSlYnagr11iMuhvqxVq7WkT2A0qstUvycnua6ppFX4nIiU04L8XFRa9uv/2orydO/HRzQyQSja2uikQKE3OAD4pECgfGXNUPpNAM9uwH5FRVz68Zt/MRkXerazpO/OzzJwYM6J67YsXKLkuWLMgpL8+J5OZ2prCwsLKqqqIQ4MD9b+HqP54F5LyMph/u2pRrb48EH2srICJHAd/aetTl0yTX11Wvra32M9NfAuf4cs/VaBrPTmi0NtPsRK3AQ6vg5fgyrQnQqlhr/ygiC4Ar0KDSWooKhgxAZQE7A2+trpyxvLhgaDHQ+2cH7dTj1Zc/POTGP/+raFXltw+jubUNSezFz1cFLM7L7XkdcEhebk+76aaDt//nw//eevqMzz9P55qvlht6dOnSuZt/uhPq8liMfr53o26rBv29VdXzHVrRA3DLzuOm7HXdx+90qlizprB2VJRodDX5eV2iAJFIhP32HzN38aKcXmhbnPtQ99h6XY21PvhYs8qwikhfYAtrrW3BMYaTohW0V6z6KGFcKWpoo8BbXmJuKrpUfDEDalqb0PppUENIs3qqpUiSir9/HkFTgeJL7gK07K8QnWEWo0ayxP/s4n8OQctQk1mJahL0RJfr49C+V2e8/OJ40DrtNcUFQx9dVfntEhHJEZGeaZYIb4zm/D45bdp3nHjSUduSZmeEU35z3P7vvP3+tOXLV/yI1pSDugQisVjFT25sPuf352jebV/gVGvtewlDJubE3Dkb9+i58us5s6oNOcsh0g9iOc6tIb9DblHPnj2WdOvW9YxHH3npKRHZboMRYAl5rK3CQFpeZXIqWubYIF5V6Rb/iCvcH4imXl3kDccaVD3pMeDTxoxtJFJYCiyNxSqqgUiSy6E1GEYTy1lbwNnAdt7NEieG3phi6GysBg3srKFWwCVR3GUq6nddBDyTwihui1asQa2KVrL2ZQfUQH6JZiKcRW1QKCV5uT23RrUnag/SIT8t5SsR2a1fvz49vvn2qz0Tt8diFQ6IisgIVKN3V2pzbeeirX0EdQG8KdpDbc24nY84wm9jUEmX6Ndz5lRBZGNjDMbkMHbIsMpTu21kHstf8/YK4xYBbChGFVhvsgKyzbB+jEZxU1Yhpcl2NGNp7wNlL1A7I4lXZP0MrUUfdOSRRw/Nz8//sVevnh9369b1ZeC/cWm6SKTwQLTs8mY0aNYWPd6HA4eLyNb+eTmwEK1Ym4+2l/ke1UOd38Lc38VAjbV29xYcYy0iMg+flVJUMOQJNItgFLXFBV38z51T7D4a7fe0zEf+6yUvt+ckNEq/lgED+1FUVNigsRLtJHseWi67lb/x7oQqZ41CDTzo//cbwK+stalWD8/defrvH6/Kz537/EHHTe2ixQO5ANsUdenwCrGuYOjUoXD1FZtsGd2upHvJypJCju5WfMVxj9zdmg0is5bYetBMMKsMq1Vl+f+JyPAUKuvpYjJVPODLVZ9CCxJ++eKLr91fXV09ChiVn59/fPfu3SqfffY/Fb169Vjdq1fP7gsXLuo2evQ2PbyboUmCLs2kFDjXqnp+vHvnUNRI9UVngAejqUpFaUgJ5qNukFRKWY+hM8RMkRiIKkINTnzby9SWn14MVO+2xza7Ll1a9ur0SXM/HZ47IrF7boM3MOfcLD+m2BgzpEvXklXFxUVjr7/+ujp5y4X5G28LnFhR9ePvfCXUNcC76Cz0AzSl7DvUiF5j6+/g+xNGfPljF3QmXucm8WnFio79evVy52y02fIdOpR0BVjRueO0FV2LdjrukbvXk3bWTcSl1ygw28kqw+oZh84Sm2xYRWQrdHaWcRYuXHRYdXX12tlRVVVV4bx58wvnzZvfFSASiXwGdMvPzzsITUZ/r55DZZL+aCEFAH7GtIS6LSbSxrs/7haRL1D/6FLUIMT73WdyBvV3dIl8PuqPXIM20rsKuGZ15YwJIpJz0WW/fmfqQ8turfnc7VFaXvqL0vxNilHfdVrRcecq5gIHGpP3LuQNOeCA3W96+OH7PkscU5i/8b5oxghXXnnl/3Jycm4A9rLWTopvbyFnLupVsnuPhSuLTELu+Bkl/Skp2dhEjOkaM1Qv7N310crC/JNO+OedbbHayUocISugteiGLjubw2/wea2ZZt68+SMbej0WizmACRM+ennffffsRUILkVakB5qYnxG8D/m3XiZwZVyyD829vIIMKkRZax8Rkc+Aty++7PiZQPkPP8zrNOvHBTNGbj74ZhFZs3JmVadZr6/apKosVoyrU8xSDWv94o19h18F5joX++KAn+02tbi4Y50qvcL8jc9CZ+IRgM/+9/XKbbcbeTfqPpiUiff6zr5bzQEqd3990q3AH+Lbuxj1elQU5s9bVNr5jOMeufv5TJyvvROCV63DszRdyT/O1sA5mbuUWubOnd+YiMdqdLb4OprI/oeGh2eEnNYIkCUGlLwOwVUi8iUZvmlZa6eIyIFo+tDODzxwV533cljJ2buQeuYfX5UcgYq71Es0tvpZEfk36hYZCVwoImcA91x37X0XU1fUhTde/3DNttuNfBuVcMxIWpO11l17+R8cKtS9lqXdi7+rKOxwUb/ZS/593CN3b7Cz1GRiYcbaKrzTHGPhl7HNFtlo5Nh5W+VsmtuxJPf9mWvmdJtbtXBQtaspTBoWX5ruA+Q3xQfXHrDWPiMi74rI/tbajLUit9Z+IiJ/RMuTt4hvP6zk7Ajq40zFAj9b3RO4z7dnGYZKK8bTuDamrirYfNTd8Fe0b9aB3Xt0iS5ZvDz++nI0a6EGdUddAdC3026DgcjcsnHT+3bazcwtG9fk79eOHU8p7j84v3CHTvTIi/CjvzZqcnPn/vKJf2TC1bD+EHysrUaPZu63B5n1Aa7Fxdy+3ReUdOneaetdxnTaGuecixGbaTB5Ty15ddnMNXPiboJvSZFDux5xPdouuHsmD2qtfVREfp+47bmVt8cOKzl7ORp4q0Mkj/0Xfloxqtfowi6oL3gGmpP6IDDNC7e8Z63dLXlfERmLGtbl3bt3XpNgWHevqPpxUsI407fTbp3Rzg5z+3ba7S109txoKl8K7pj9XVWnP5d0X735ft2+Pfzbbxcbwyhw7aq4oy1wZC4rwBhzP5rFsdA5t7nf1g39Dg9E9RuOcc4t869dBpyCphCe45x7zW8fTW1rlpeBcxOaEKYkGw3rQhHpXk/qSkOcQhpNvpqDWxE9hmoXz1HEGGNyyBkA8LOuY+dGXWzip6smxz5d/dVzsVjFvSJygm/VAhoAmoQGlD5ozzNZa+0bIrJIRPp4ge6MICJHkFrU/G9obvGoxI2xajab8+7qCX9/4bpTGjhsLNVGa+27eNGd6669bzt01puP5tyuZVVZZR4ahCzxj+Fo94YmsWPHU7bAC5GXr4wVffTU4r2Wdyoed/Imq4jEXKPqXBsiGQxePYh+hxLLjC8F3nTOXWeMudQ/v8QYMwJNrRuJZtS8YYwZ5vte3YkKzH+IGtb9aaTvVTYa1qvQpl9rfaWbFR91NbpEuxedEW6BRqlvQdOL3t79mEEj77z/xlZpkx1bHB1V32udcoq3Bdiry46xXnnd9gXw6lsAiLZyHgscApzvE8qno3/c5WgKzyR0KTypDYoKWkoBWlXWYsPqK5ROBfYCjkl+/bmVt195WMnZUWoNa1kkzxTGql0usZSVW02iourHjwvzNx5cUfXjT2YfK5aW90f/KOOsBh7v22m3yNyycSmNdj3cRILwN8A3Zbljp5flTu3eLZbsTgo4QyxDwSvn3LvGmIFJmw+lVpDnIbSH2iV+++POuTXA98aY6cD2xpgfgBLn3AQAY8zDwGG0N8NqrZ0mIsVJs6LZaE321mgRwZ/QyHD8jn/i6uVVU1vrmmKLazo18PLzaKnsf7co2vQn7Ti80PV0fEsTERkXT7L3eafboYbjXKBfQq6pQZck89Bo/LeoD3Cqb1ZYSG3teVtyPRrE2qixgWnwNHAlcHYDucdz0Jzg059beftDp+978bfzP6zYBO0w2xBVIlJk6+kOm5fbc1fg/eqaRSmXdD1KOyW/vyL0s7gBvfmnRU6e6R+trnOKaWBunREtOKSgqENbVc21GxxNqrzqYYz5JOH5Pc65exrZp9Q5Nw/AOTfPGBO/QfejbpribL+t2v+evL1Bss6wAuTOqpyRs7j6j2j6FMCbwD3orOZwv63uMspQ2BKF+PoQkRyqXQ76x12QYsik62bdlnbr6US8u+NV/0h17jy0O+QINDizF9BbtNVyPtArqcR0NVoJtBBtid0cf2Bj13yniNyRocNVop1NG/I1PgI8+tzK28sBcosi8Rl9Y7PGmWgmwLjEjXm5PQ2qVfA2WmCR0uXUoSBvaX5+TkFVVbRjwuZ8GmgvIyIDKyvWHDt37qKD585dVFq+tKaQWGmiUMwdwB0Tyu+bLCLDpzdSiruh0gRXwGLn3LYZOm2qk7oGtjdIVhrWnMXVYztMLt/zL9ufM3nVgd3/M3XV09M3Kz7qdNS30QWdAb4NnITeUf5RunHxMrSc9LwMX87Y/LHFr1V8s+Z0tOHb1qiu5paoL+bldA/k28ukXZHl3QKf+UdjxzZoRHwQOpv8I80LtKTDtAzdxGagd/9629c8t/L2Or7PjqW5KzHcgaMx4z4FbZ+y1rDm5fYcgeYXz0cFXzbLy+25vLpmUSp9ih+6dC8auHDeyi/QKH5f4JW5ZeM+9//XO6CBkZOBZWZ5zSrz/tKBd01+taSypnrtDXinYWNxSzpWFvfIWzX6qF5bA8f61LX3rbXNzddef3EQbd2S1gXGmD5+ttoHnYSAzkQTVyn9Uc2H2dRN/4xvb5CsM6y3Dj2zpxlacH+sU852keU1Nxe+v+L3tw498/jOOTkVK6PRh5zOMl6auurpFzYrPuocwExd9XQFgIhMEZGR1trJGbykg4DXrpt1Ww0w2T+aKwc4lgwm9CfijdwC//hQRFqzhUc1Go1vdiDOC4H3t9aOa3RwAkV98mZuc36PO1J9xgM67bstMGNm2evL0EDTUQDG5EUAl5vTZRGqulXqd3kPIC+3Z1F1zaLypMP1OuCIraYWFOb94pPx3583+fPZpw4e1muMiLyL3gwG+3EvUxX7We43FS/mLWX7fTtt9u4Ly74YC1Bc3JHhY7t/+uH4SQVjD939YlQEewEaVGmxj3h9RNWtWjXd6gVUp/k6//P5hO2PGmNuRm+iQ4GPnHNRY0yZMWYMMBFVWGs0SJ51htXlmaqcZdHCmv4djuowtfzN/B8qNwbe+XnPbguMMYN+/+0da2cwU1c9nTz72wv4j5dZy5RLYGt8TmMG2A74pNFRCeTn9e6NppItAN6rqp6fDcGt+cC3IvI+uiSPL8vz0JSUxOVzvAw43vDQoCl1xWjgqqnMRb/0aw3rgE77dkaj+7eiPdH+g87yLzAm7yhU72FEdc2iqXm5PReiqTaJ7IIWdgBrXTCrCwrzfgDu2HanQR9vu9Og/VD/8r9Qf/8/UTGYvTuMW/muKas5EGCboo0Gz1yz5NNJ5XNGd+xYwGf/mzp6/oIl/zQR82v0s48Ch1lrGxTp3pDJlGE1xjyGBqp6GGNmo0Hx64AnjTGnoHoUR+s53WRjzJPoSqcGOMtnBACcQW261Ss0EriCLDSs5035+wrgwet/dkEdwQ9jzOxEo5oKa+08EXkCzRb4fYYuKTeuYJUBtqQJLUE8w6iteErl410XnIHmjV6KfofiUe8KVB5waSPtcP6CqkG9Xt+YBphFgmEc0Gnf91HVqfhfY39Q/7UP8MW3H2FMh5dzc0pS5Uk/RF3R7GeAG6y19yddN2imwNbW2kki8jrwJ7Oi5vZIuUZcckyk/5HdR/Xfqmijdx9e+OHYocMGUFMTPRV401r7SxG5GF+SG0hNpiqvnHO/qOelveoZfy3wk0aTzrlP0FhH2mSdYV2LYcEa58q/7RBxnVdEvyntYNISjbbW3ugrhFrsEohXc7XkGEl08S1kmsJ0VKSkR1X1/Ka0ue4rIkejkezZwLxMtcm21k4XkSnASdbay5txiEloulxz+A5f5Tag0747oj73n3xG+Xm9B5540hFF1M6Ur4HY0agP2qJB0K7o7PpGABHZFZ3dvp5sVD23AdXeqPYEdi18ZOFDJkrn5IGDOnQfesKo3ef22W1I3113G/Uiqt16JXoTGI3GCAIpCCIsrYiJuuFVHXNemTK9cl80CDG3CSKrPwdeFJFtW+gSGEXdssiW0uRvTFX1/LnAk03Zxwun5KFR8R6osE0nEcnk510OXCoizwKfN1Gq8RXg9Gae9xt8wj06+x+W9Ho88tF3wvjP+ufmdLsA3Eow71fXLDjQv3Y2cHZebs/RqJ94ht/+KKpqlfKGbK19CkBE9kd9c2eYKCWk6HacYyI5g5eVFFQuq8F1zY27S5ajieYNFTZs0DgXDGur4nLNRdHDe1zd+V8LL1kxr+pvQPlOHX9TOL783kaX5d4l8DgtdwkcCvy3BfuvK/qhxu6S1jyJL3Z4C8iLRAr7oX7TUcD8WKyi3iola+0in6vcnFXFNGrLnl9G/b1j/fNy1A8GMLmwsKDaGDPW38+G5+f13rSqev7asufqmkWf+vdhvL+4S5rXcznwW2vtQzc+fFZf4Hy0aEIF8EvzYq4w0qt6dPGVpjx2tYN9/X5HoD7g5U18zxsQprWzAtqErH0Hl7z4l/lEzC1jji8diYoLH4rmEaaFtfZGYIyItCTPbQdqo4Ytwvv72sq31o826IVlrZ2C5ovegy6nFwGv4bvlNsI+wPMiskWjI+syAo3aMrPs9fPQ78Vv0LSZ4pllr9/kx+38+edTE9NnBgOX5uf1ThWN/zXqUmiwc4WIfBDPCrDWPgRwyff371Tjovug5ZBflscqP/77rBfPqxlW+Loryjk/1jMPNKgGWt58NRpcC9SDcyatRzaTtYbVcyjqI3sWXdoOb+L+hwN3JTfCawKF1tplzdw3mT1phnh3M+lL8zVtm8Rbb72bc801N15I3aV9ozdA33NsD+BBETmhCac8JPHJzLLXl6NJ+48ljavs3adn8urmhAED+u4mIhuJyOUiEtfMHQ68ba09uZFzbwT8H5qmQyRSeBjw1OU/PLQd8F7MuUWPLHin09w5884q7+I+Q7+zV6Ipew/6Y1zkq/ECKXBo8CqdRzaT7Yb1bPSL/AAabX55p46/6dbwLrX4ktgnqJ0xrEt2RPPg2oJSahOfW5X33hu/OhqN5lLXrbRZJFJ4RmP7WmtnoX7g3zXhlDegAai1zCx7fcms1e/936zV7/WMRAr3Ovvsc+xFF58io7YZscIYE8+1je6yy+hXD/jZ2BuAf6N++6ki8hGam3hdQycVkdHAEmvt+wkdV+P6Budf8v0DW/3fDw/d/F3l/M2ATW64/q+TampqHkZnqJugQbf7aH7QbsPAhRlrW/AG0Hl8+b1VqI5mN5qYU+pdAqNEpMEOAMmIyDAyu5zenLZp1wKafL6gjc6VqkhgM+CAdHb2wcWIiMz0y+zGxlcDZSlWIfeg7/mNF1545fxbbrmj5D8v/ueEWCx2XMzVXAD0fuvtl3522223DrLWbmutPcpaexbqFz7HWjuHhjmPhMTwSKTwNHwOJDAGrcCKV4PlAI/+6dqbu1BXmSsKjG/sPW7YpDdbDTPWlnEOtcGnvwP/oIkJ9p5jgIeb6BI4AtUoyBRdUH9kW9CTtmuLfXPS89OArrFYxSGpBtfDGDQ39Sfaq8l4RaxYimyPtb2/Zs2aXVJZWbmlc+7SmujSztHoimOqa5bcb0zeNikOOYNataP6ztkZzRZITMF6ltrv4jcG8yZewDrOgIEb5aDFB6BuoN+iKmaBenD46qs0HtlMVhtWa+0KoKeImPHl9y4CThtffm+T24N4l8BjNM0lsAsZClx5WqW7QT10J7NpYg2RrCkajcUqUipK1YefhfalkRWCvzHGNSLqsMsuO2750z34EjWaO6DdandKMeYQ4GwRaehGsAVJSmKxWMUitOvqyWj2yUBqO8sCLDjwwH3XUNvZ9lzgFGttUyQHN0iisUhaj2wmu69OeRF1/jO+/N5mGyZr7U3ANk1wCTS371Y20IW2M6zJ1Ktd2wiXocv5htgH6GytfSm+QUS2FZHrttlmyyMGDNiookOHDnGN20fRG+kkNAB6PboUr4O/2W2OZii843OA48fuLiL/QvNu70reNxariMViFQ/EYhV3RWOrL0UNbfwcpXf8/b7JQGc0veol1CXVMfk4gUTS869mu481a/NYE3gTrc3ORIpKUwoHOgJ3+GUgqCrVZ2he67gsF6TOzVSVVRrkJD1vbtnt34AXRKTMWlvfZ30OUC0iVajxjCsSDevWrevIE0/8VbnIdYVASSxWMRHAmLzbgbecq/6yvhP7pob7onoB14nI/ajrqRdqFEvSXG3MQWeu56I3t5ustVER2RN1R52AzsyPSONYGyTOrR/NBE0jrVuyAhHZBxhmrf17Bo51IbCRtfbcBsZ0Bx6x1u6fsK0zGpDZA5UPjM/2y4BP0V7341IZNB8I+7O19siWXn86iMg7cTHt1iYSKdwMOH7nnccc88EHH94CfB+LVbzU2H6pEJGuqB/6a2DX5P9LEXkPnRXeA5yFCp5fXZ+YdRPO+yu0zHUfdGaZi852n29qCxoR2aSmJnpUTk5kF2NMd/TGMx2d9W6HZrqUtOcWPa3JsKJu7m8j90lr7H4fP/lpBvVYM0p7mLGCzgR+3tSdcnO6Rmqiy+r4tKy1N3ktgS2stfXNYn6GGsvE/Vag+qZ1NE5FpA+wH9qu4RKvjORQn9wMNMDRB/X3rXfEYhVTgctEZMf33nv7b43u0AA+Z7jEd139DC0GSKQKONBaG1cXanFlmYjsBgiwmV+FpCW24f29O1dWVB9UUJi3NbWKXktzc3PeAs6w1s7youTXohq+Z/sxK8msBsV6RWw9+K9pF4bVL9WqRKS4iXJrb+XmdO2AprhMACbURJfNQdNkXk50Cfho8zVomeRYNL0rnWubhyZ/P5i43fvqdkbzV08BykVkD2rVxxejM7NJaLuZmW0Y3MpqrLX3i8jvRaSzv6EhIm+gs79UTQebhYgMQHU4d6nPtePFVvZEm8yNRDMuANz/JszM+d/4H7Y75pTt+z9274TF6A10wdyycYm+3Ai1lWg1+L85EbnBWtsE+YsNA7eelLS2C8Pq6dgUo+oN6hhUD3RMwvYfgQ83GzF03tixY+6nNsJ8AZqycw5aejqhJRfrl7FvAW/5WdGR8eWfn+0MRaPV2wO/ArqLSDFaCAG1vXa+QzUiJwMzsty3m0lKgFwR+Rp974OAUXFD21K8kMo/0RSoUhE5Cp1V9kdnyh+hs8rV6GpjCvBwYtVU/057ng3s/MS9H3VGy2r/5PdfmxNrra0Ukd+i5b7v+59L0EyEK621aXeU2FDI9sBUOrQLwyoiRSS1J06D0dSKLCeyMbDx1CnfUrZy1bLu3a/c6uqrr56EVivdZ639fcuuNiUdEn1qfmb6DQ10E/DveQt0ljQGLx0otc0GE4n7JueiM7qNfPrQfDRpfn4bBrMyxSXorP554C/AH4AmiUN7X/lANGA0YMWKlZtHozWDu3XrNhQYgBrPU1Ef6EdohV8EmJ6mj7o/gCN2K5qe9+ncsnE/KTSw1v5DRHZAv8O3o376ZeiKJpO50u0eR+PNzNoD7cKwokuo3iJyCvCKtbZO8vvE3Q/aHf0jnLLDOy/GP5ddaITZs+d1/fjjSS+IyEBUMak5xQetgg/IfEjdzpEpEW0oeBY6Cx6IphrtiUamuwDFGZYMTMXQTB7MWvuEb99ylLV2PL6SS0T2Rr+3XYDh8+cv3HLp0uWDVq8u71FWtnpNTU1NZJNNBi0bPHjjVaghW4r6NGc8+MDDW6xcWdYnFotdjPpEV8Zc1b/j5/Sz2OvRIFOD9O+0Zy8g3+FwuNFoWlVp3067bTy3bNzafFYR+Tsq0JJYit0JDXo21P13wyTIBrYd1to1XrR5IHCmiPRCl2Vx8etjUFX75RN3P2gCMD4Ce6Vz53vrzfdLBw7c6MHevXuWobJz7RHj6+7XVe5q3Lhn8ngHou6ZN7wCVrxJ4rVog8ApwLQHH3zimKqq6s0S9/3ww0//HI2u+Ef8ecTkn4n6svf2x3kSzY8+2p9rX7QctRLNBKgTuKyHPdGqwA/9fiXoDPY3qPBKnHPRgGe8XHdj1Kieaq19Lo3zbHBkKt3KGHMe+nk41J1zEnpDfQK1JT8Axzjnlvnxl6HxkChwjnPuteaeu10YVh9pfxI1ntejfq8/ishya+1X1JYzdgEOWF5dvWNMf2+UaDTW4eWX3jzkpJN/XnPVVVc1Raw53WsvpYlL2CYevzOwRRqGrQKduZUBzyRE1jNFjp9hdgCK/KMjWovfCTU8ndDPpQSd4RX7R2IubHdr7Rbo512DSiBehorKzAGOtda+mpPTuS9qtOoYVU9yNsE2aF5q4r32rJirqvTugleBO4Hzm+Ay+QA42GBeMuQQI7oNulJ4NmlcKerOiKFuhh+BbdM03hskmShXNcb0Q+MlI5xzFb6f1bHod+NN59x1xphL0fZClxhjRvjXR6KuozeMMcMS+l41iXZhWK211SJyG9pmJM4vgXNeOOrXV5XqH+paPl+x4ge0CWBazJkzv8vMmbNbS7RkLDq7ai0Go+6RX9U3wAfLuqBR61J01pdpw9oHTfKvRiud1qAzuQr0RliGJtt/jwZvFqHL9PnxgJyI5ABLfCS+I9CtAdnGrYFUUfVp/FSo5yzgSGq1CBYCi0XkSDSb4yhrbbJBbJDZZW8lrxD+RwpZS2vtHN+SZR+0rPbtYFTrx2GoiWXMFZALFBpjqtHv01z0Jr27f/0h4B3Un38o8Lhzbg3wvTFmOhpYblYQu10YVgBr7dv+Dy8fnfn0BsZHYrE3qnNzNsqrqb2xvLZwUZMjx3PmzK9pfFSz2J7W9d0OopEUJB8sWwYsE5FF/vdM8y1wsrW2JYpg49AmhecCPzaihbtpA9t3RWehcaJod9UuqP902BFHHjYFncH3z1SmQSr8LN6iPnwHHNjwHoEmzFh7GGMS/7bucc7dA+Ccm2OMuQldIVQArzvnXjfGlDrn5vkx84wxceHzftSNZ8z225pFuzGsAF7AotI/FvnNIx454dTTNpk57y6AGatX//DukqVNrlefMvnbEhHp3kLDkIoR/FQBKpNsTEJ6Txr0pHVcE1+gPswnWnCMp9Gb5s7ozLohTqI2L3QJ2i5lGPAVukxfS8xVRdFlIVdcccVGZWVlE7t16zYysWLL37SfQ10UiS6Djqj+wa+ttY/7sXlocPR/DRllESlABVrOstY2qW/Zhkpc6DpNFtdXeWWM6YrOQgehK6WnjDG/buBYqU7a7LzydmVYE/EVLQOttd8Y+HZpl07zuy4v67VxYWH/XGPm1jjXpIjrokVL8mbPnvck9bTGbQHFTS2LbCL90FShdOlOag3VlvIB+n/XbMNqrb1VREx1dXXs8cef/Y3IdXmxWEV9ebt7oUb4DeCtaHTF3Y0dX0Q65eXlTe3WrZvEjao3qFcAV6F+3V8ml5v6MT+KSFy8O4qmt/1VRK5qwI3wR+A8H1gMpINTvYAMsDfwvXNuEYAx5llU3WyBMaaPn632oVYQPq47Eac/LZDebHeGVUSuQyO8TwHLHj3+VFtQVfWHA557LPf5o497uveiZUf+vF/fmY/MnrNxI4eqQ35+XmzZ0uWbtMpFty690SKCdOlK6xjWt2laJ4CfEIkU9v3Nb06Y9tJLrzJv3oJ7gU6RSOHTsVjF7OSx0eiKReiqZXBjx/UBvu1QY7wS2F9EDkZnKbugrppIfZVvVjvQ9ouY/L5Aj5ir+sIfNw94RkQOA05I3F9EdgHGB6PadDJU0vojMMYY0xF1BeyFfs6rUTGc6/zPuDToC8Cjxpib0eDVUJo2YalDuzOs1tpLE54ueuLXp1xvnDvlseN/0/8XT/3zqH8eeMK0yvLeBV0jb89YFls6JJ1jFhcXVZx8yrH333zzTWc3Pjrr6IkmuKdLN2qruzKGtbbML31bwjP33vtQccLzW1CXwGXJA/0scgBqWDf2vw9AZ/B56BK+BxpgqkSDZl+hhRYfxaP/IhJDS1pTGtWIyf872j5mJTqLmYIGwuI6sof4yqpPReRYa2286GM3tLAh0AQcmcljdc5NNMY8jQYVa1DtiXvQLJQnjTGnoMb3aD9+ss8cmOLHn9XcjABoh4Y1mZxo7PPDn3h47ex0YlXvPfJM5Yx9Ox1UEHOxuatiZTMX1Myrnl39Y/eF0flDSdHobvDgjaPJRrVzx81/gf5BfQG8sKL8q2WdO26ev6L8q7R7FolIJ1q/x1FBE4MvXWm99sstKpopKur42urV5WMStw0aNODE88+/YGznziVxl8AIdNm2EvWtLkCXc9+jylRfWmsrvFbDa8kVVEUFQ0pXV85ITKm6GVggIv+21v5EQBs1zGcmPD8veYC19h4ReQXVdP0GFd8ZGcpVm4MhmqE8VuecRQOHiayhHnefc+5aNGOmxawHhjU69Ilfn/yrn//r/kcA/vbf6+eePOz8mWsWRzeNmEjfkpzOfUtyOjO0w/BPnljxcAmwJbosHA2MjkQiw0duvmmqVKsjgKP87ys7d9z8fuDkzh03/xyNXt+8ovyr5Y1c3s40ULa6juhC01wHbcbee+9+5vPPv1xn2/ffz+z917/e+YtotPwdABE5FZV9vDLFIdbii0rWdjcoKhgSQZd9+xQVDPkX8NvVlTNi1toLReQi4D8ickTcXxox+Vuin19iu+xPY64qMdsgkaVoVPnv1trJTXjbgQSaGLzKWtq9jEzfhUtzB85eWKddi4utlXBLZFxNdNmamuiyj2uiy+6oiS47pSa6bOvL/u93lZtsMjBVV4GtE34vQatsStC81D+gs9nGGEP29TjqjBqBrEJE9vzuu5k9U7w0xzn3bcLzXNRPlg5TRWTfooIhI1F/7EFoAcMvUD8asDYd7VfATSIS97MfglZjnZNwvCERk1/fZORs4NpgVFtOLM1HNtPuDesO77zodnjnxST/mCtKMfSN5A0i0j0SiawgqZ9R546bFwON+WfP7txx88bat7RlZ9Z0WYlv0Jg4o8sCylGDfwN60xqC+k53jMUqEtPJCklfkOd2VGTFUTd1piMwo6hgyFNFBUMMrNXbPRx4TERuRJWqhqI30nvR4od5aGEKkUjhrpFIYWLmSWkanV4DabA+tGZp94Y1FV22Krio66iCF4qH5o8r7Jv7eV6XyHcDT+i8mYgcICIjve8Tn7O6GI0OJrI5jQsRH0LjLZ57ktA9NBvwS+gvRWQeUCUNN9FrE3zrkglHHHHw+bFYxSWxWMWnsVjFd7FYxaxYrCI5qp62YbXWTlq+rGwAWtKa/F3PR8XJT08cb63dDhhpr7r8KXvV5TNjrmp8zFWdGnNVZ6OrmJcjkcIhaO3/pEik8I99+gzY0Tm3shlvPZBEvDVLe29/3e59rKm45Zk/3w+sbVXsI8ilaNR4c+BAL8t3KBolfCbpEFuleaoxaEvu+nDZKF7tpRF/LyK7or7FhWiQbjHwjbW2NQsaUvEIsMBa+1AaYwtpgljOvHmL89Gqq1TkArcXFQx5f3XljLUdHqy1PxORj9CgWG1LbpMTQXuvxau+BgFXLFmy5Ip58+Y3Kb0vUD+ZCl6tS9ZLw5qMz0Oc6x8TAES7ZR4KbJ8iqp6uYd2vc8fNfwXkryj/6oEUr7fqN0RECklyYzQFa+17QBevHnU0mk50nYhsiqYqnY4Gb0YBLydWpYnI9mhd/kVoP69m9Z3yN73eqJFKh0LS97EycGCfR40xI51z9X3XL040qgnsAFSIyIfA0dbaWbFYRWUkUngC6k/9HUB+fn7V8OHDdr777rtCvmoGCHqs7RxrbbnX9rxARL5FZQhjAFU1yzrlRAq/j5gOA4wxDblL+qGzodOTX/D+y4yrZSUxgrrCNM3C9/76EkBE7kWrp96jtiplNVAkIk+jy+cqNGWrA9qV9XIRSauVTQoe9tfwQ5rjO6B5qWlR2LHgzj59e1wyd86iEup+38ejN9ovUu3nVxoFInILGoDsCRCLVXwTiRROATDGsNVWm8vEiR9kjY5v+yf7/afpsF76WNPFWrvYWns5mhT8VxHZxpi8SGX1wsNXr5k5qKzym+qVFd9MW1X5/fiKqnnjqmpWfBKLVc1L0dl2WorDb4eWPbYmwzN9DmvteGvtRtbaX1prDVqSWYz6JFcBB1hri6y1/ay1Pfxrv0GFU3YTESciJzbhlNvRtLzafJrQTcJaW7bTzlvG21L/iN4kPgSeW1054+jVlTN+EtRM2v88tHPD2nLH/fffe2Fpac8Fo0ZtZSdO/OBPTbj2QBqsD1kBG+yMNRFr7Uci8ilw8imnnHjSffc9GM8q6ABu05irIhatojq6Nj6xEswCY3KX5pj8yg553aemOOxOaOltazKEVshJ3avojDNRn3QFULbXDWf8Eq3Wuh6o3uuGM7Z5c/Wdn8XHW2vvA+7zZZwHoalLD6Z5utdRWb90aZJhBRg8pN/Miy87fuQNf374IWDh6soZTe0meyYw08sZlu6ww7aX7rDDtrtba79u4nECjZCpyqt1TTCsHu+H/cf//d/lPfbee4/e7703fvc1a9akyqsEKAFX4lw1Na56Rc2a1QtTjNka7U3fmgxEGxZmmtPQQor6+A+aFVEHa+37IjKNNI2q14l9D1VtT5cmG1ZUZOWZ1ZUzmqUFYa290zeEfA3Nh30uGNXWI5Z14d6mEwxrEn/607WLgWPOOefcHT74YPzNn302aRvnXGEDu0xzrjrVV6E/Wp/cmvQiw+lcexWdsRm1LZ7royFpxdXAmCa0ailF/bTpkoeWJaaNtXaSiMwTkZ9ba5urvnUCqjXQy1rbWJpdoJk4F7IC1mtuu+2vE0Vkly233OKUd9559/c//DAzVXUWQCo3AGgfqtYOXhVTG2DKFPuj3QDqYwEqh1cfFWi9/u7pnMzX9FeISKdkub56aM6MFTSl6+80U9bQl8jOBe5qzv6B9Ml2/2k6bNDBq8aw1roHHrj/3hNPPG67ww8/9NbOnTvPRCuXzkaX4XsDt7XkHCKS55fEzb7Glpw/Bamq1hKZ++bqO+v16zb1erzKlEEbB6ZDHs1r+rgZ0Khma334tuMjrbWPNfcYgfRYHyqvwow1Day1FcB5V1xxxU3Ll6+4pGfPHpVo2s+78X5NiXhDmajNuT9aGlmNGuaZaIQ6/hgvIvHh8X2XoMpK36FCLlPRdiWt7YH6voHXltM6yljjUP2Fq9IYG89OaCrb4Ut5m0lX6l+dBDJEyGPdALnmmmvmAOeIyGZoL6tjRST+f+jQFinfokY3UX18B+BKa22NF10egNbBj0bFdP9nrT0/Ptgb5kFoldgwNNJ+OprMn3hJqcRmWsrz6HJ/BpqGtgKdTT4BnP/m6jubrareALvRcAVbIrlN6KSayFx01tpc8eJCsk/3Yb0kQx0E1inBsDYDa+1UkmYvfqnYD9gELXlcICJXowa3v7W2xu+7Ak1Kjyem3ykiZ4vIsLhIsp+Vfsc6kPd7c/Wdq/YqOmPwm6vvLAfYq+gMAzz25uo7m9JZdJCIFKVTjeVFoqvRTIR0KBeRntbaRY0PrcNZaOuVdMpmU3EekKq6LpBhMtRBYJ0SDGuG8FVb8ZbIbzdx9/uBa0XkRWttgwnrbUHcqPrfHdDUds0voW2Ft2tokJ+93w1c2gQXxzfAjmgrjabQmWY0hxORwWgrmxp0JRJoRRwQzZAvwBjTBVUm29wf+mS0mOcJNEbyA3CMc26ZH38ZmvoXBc5xzr3W3HOH4FUW4Gd25wNHiMiwdX09GeAsYFQaQbnfokI11zfh2F+g2gVp49XMBgFbeW2CtPn042mTn3jkrTkv/+fDiyZ8MHmPpuwbaA6GWJqPNPgr8Kpzbjiq/zEVuBR40zk3FHjTP8cYMwI4FhiJZsbcYYxp0nclkTBjzRKstTEROQ94QER+b63NdBpVi4mY/BzfSrpBrLVORCqAC4EbGxh6Fk3Pw/0Q/eIDICJ9UB90XL+1F3W/1wbNe/0O7SR7s4hciHbkHOj36Yvm0/ZEu9h2KCsrz5v5/fxh+fl5y2bPWtQHYMjQfr8Cbmri9QaaQoa6tBpj4qL0JwI456qAKmPMocDufthD6MrqElSQ6XHn3Brge2PMdDSOMqE55w+GNYvwuZKC+mfHevWpdU7E5OcAJ6FfvoPT3G1P1JA1ZFirgAtFZCs0d7YUNYzxn8OBHZLcBFOp2y7l38Dn6B/IY6jsYco8Vz+DfgV4Fc3OWOAfc9Dy4x+BH6y1ZSWFI/oAswYP6bM2UDZvzuKN99v7aPPaG0+tB+GV7KSJWQE9jDGJAjj3OOfu8b8PRqvkHjDGbIW6s84FSp1z8wB8C+z4d6kfetOOM9tvaxYmhaBIYB3jfY/fAvtaaz9fx5dDxORvh0bTo0CvmKtKq7WLiDgv5FLf63egS6/ZwDI0xWwhqtg1G5iYan8RGYeqbO2K+rSfA0Y0V7owFSWFIwzwIppe9sv49m1GDz1t3PvP31PffoGW0Se/1J1Y+vO0xl43+/ZPnXMpWyQZY7ZFDeXOvmPrX/E56M65Lgnjljnnuhpj/g5McM79y2+/D3jZOZes1ZwWYcaahVhrV4jICPSO26Yh0s4dN98UDQy9B9y0ovyrr9G7/avoEnw/dGbYIH522KBWrLX2zIZeT0otS8Shs+fL0XS0CtR/Wol24Iy30P4GeMlau7ye4x+C5qc+Es/aiLOyYooDDiwpHFGI9sj6HNhs3rwl56FtlAOtgCNjJa2zgdnOuYn++dOoP3WBMaaPn632obZycTbqHorTn7opk00iBK+yFGvtYmClF7NuS65BjdUpwBudO24eQUWdd/KvFxcVDOlSVDCkMUGTx4H/tvBanhaRh0XkYhG5WUS+EZH30VYr/0CX8N/76z0S9bPeY621aLHB58C5InKRiNTRIxCRh9Gc3eXA7SKScpq0smJKBXDgyoopo4BLunUrmXTllVeOaOH7CjSAc+k9Gj6Gmw/MMsbEuz3sheZlv0BtK6YT0O8AfvuxxpgOxphBaL+z5uY8hxlrllOCfsApxZgzTeeOm5cCByZsKkNLXDv6a/nG/8wBphQVDJmHRlIPBopXV85I7Gi6G1pN1WystUeLiEP9oN+gPeIvQmv+R6J+1fHAR8npWv75ZGCyT5laKiKjgE6o/OFyYKi1djrwvIjcLSL9UrWlWVkx5RX/8zYRuRP4i4icm41td9YHMlh5dTbwiDEmHw1enoROJp80xpyC+tSPBnDOTTbGPIka3xrgLOdcs7U+gmHNbj5GZ4qtZlhzcjr3RyvBPirusNGJaIVRnOHol/PPwD0xV7W0qGDIAYCgNfsbo4YNYFVRwZCLV1fOiKv7L0Fnls3Gq/eXocGsHugMVay1zze4YxLW2u9EJF6UEUXVtCKJhtFae5qI3CYiv7XWplzqR0x+75irmi8ij6Mz+nub9cYC9aJ6rBk6lnOfk7pN/V71jL8WuDYT5w6GNbsZiyo/3Z1qduT9mL9Dy12XoxoEb6FlqE/XN6PKyencEzVWs4BjgL8Ad1VWL+3TIbcrxtTxceXFXJUrKhiyY1HBEEFTklI1zitGjWw8x9ShhQIp8bPIBY0EnF5CZ8gno4GkaTQzSd9aO8b/f21an5aqtfYcP3N9EG3Vk6x1e2PE5D8fc1VPi8i2IjLGWvthikMFWsD6oBUQsgKyHC+wfDmwf7wnV8Jrp6J+xetQxafB6B36JOBTa+3FqY6Zk9P5bFSVqwL4moSZZcTkjc/P7VySGyncHDW4F60o/8oVFQwpRGeu5zZwuTXAthdfdnxPNMB1p99u0Nniv71eQvxLFwOOtNY+18D7j4/NSX7/rYH3af8OFWzZOl46GzH5/wH2Qd0ST+bm5t5y+RWXnAPcZ63Nqhbn7ZnSvFL38x7HpjX29vm31ZsVsK4JhrUdICIHAqcC/7TWPuO3GeATa+3oFONz0Sj+X6y1ryS/7mesDwH1CTZ/WZTf7+SVFVPqNMkrKhiSgwaFrqjvWg88eOdHR24+5FPUAD3qiwUMKkRzBrq03xotM+wMzEe7tBb4PF6DzkpHULeU9mXgkDbQuAVARG4tKyufedffnv0XsLh8zczXUMMa540xO+5w6H777X2btfY3bXFNGwK98krd0d3TM6x3LAiGNdBCvMH5G7oUPwMNwMy21l7VwD7XAJOstU8lv5aT0/kZ4Ih6dv1ZNLriJwYZoKhgSAHwCRo8qsPY3bf5qKYmOva/b/w7LfUpX15qUBEWUP/uKjQ4BjqjfQn1rZ4O0FBebCbZbtQ+g6ZO+SEugvNc+ZqZ89D/90TG2KsuP9haW++NJtA0euWVuqO6pWdY71yYvYY1pFu1E6y1zlp7FhowmQx835BR9ftcAewvIsUpXl6OJuKvSPFaqvEA+ODU+YnbCgs7sOtu27zx7jufjUnXqPrri/r80Z5oWWEFtQIuzwB51tpDqM01TO8vLgPU1EQTxbQPzc3pdBcazEu84WyPqm2VttV1re/EK69Cl9ZAm2KtfR1dOqfLu2iA6v7EjdHoilMAcnI690YNLOgy/bhodEVjUfdxaI7f9h07Fqzed/8dDnvk0fubrcrlc3Zv8I9U3IxmA7TZ31O37iWPJzw1+bndOlXVLPkiYvKPRNWRDkZ92ha4WkTOC+lXmWF9aCYYZqzrP/8EftaA0lRewu/vp2FUufiy42t+e/rhj246fMD4YZtuvH1LjGo6WGtXtKVRBehU8pMONaMAYq6qAjgebds90Fq7Ep1lH9SW17c+49J8ZDPBsK7neIP0EfX/4c9Hhbm3oZ6ldrz6S0RyROQE4PouXTu9+7/P39r5gwmvTGmFy17nrFyx6soOHfSeE4mYhSSUFsdc1XLUP32M3/QGsHdCN4lAM3FO9VjTeWQzwbBuGNxLPUnR0eiK6mh0xTfR6IrPo9EVK5Nf9zPdcp/2VAN8Za290Frb2q291ylvvfP8e3vts13lb0475O6V5d+Urqr8tk7TyJirWh1zVdWwtsrrH/gAW6BlrA8+1mBYNwCstUuB90Tk/EYHJ+BFtxO/wxtba5vaTaDdMnKLwVt37VYSz15okDtufG3BF/+budNB+x7Xvy2ubX3FoT7WdB7ZTDCsGwg+/7W3b4TYKF4jNZ5e9JG11lhrZ7XaBWYh1tpp6Gz/VBEpTiWIU1q8042lxTstAha+/+bXvyhbWXlNm1/oekbwsQbaG7cCfxWRLdMY+zxafbTGWrtDa15UNmOt/QQVoinjp3msoDm3PeJPKiurGlP9CjRCmLEG2hXW2rmoeEg6rabjM9vBrXdF7YZbUJfIT4oigK8Sn1RXRUPDwRaSCdnAdU0wrBsYfjm/vKEItoj8iGoPlHtjvEHjMyv+CfxaRH6V9HIdw1q1puYneVqB9HFAjXNpPbKZYFg3TCagDfTq43b/8642uJZ2gbX2RGttB2BEUjDrKxICfDXVYcbaUoKPNdBeeRc4rb4XrbU3om0pZrTZFbUfXgH+JiLDARasGl+OiigDEMkxwXXSAlya/tXgYw1kHV5nNF9E8hsZumtbXE97wlr7PvAZcLkX4gbtHAtA2crKVaXFOzX2/xqoF5f2v2wmGNYNl1eAsxoodX0R7YQaSMJae4+19jjgWy/pOBXAGBMbMqy0w4JV46vW7RW2XzKdx2qMyTHGfGaMedE/72aM+a8x5lv/s2vC2MuMMdONMdOMMfu15H0Ew7rh8i4q5jJLRG4SkR5Jr/8VKGjA8G7wWGvvAE4avePgDj1LSzjgsK0je+y/eWjX0kIyXHl1LgkrCrRT65vOuaHAm/45xpgRaEn3SLQb8R3GmEYLQ+oj1DZvoHi5vktEe0z/AVgkIg+g7V2mou2CQQWpv1w3V9kueHmHXYaO3X7nTVYZY4rR7gOBZuKAaIYi/saY/mhzzGuplbo8FNjd//4QKlF5id/+uHNuDfC9MWY6Kgs5oTnnDjPWDRxrbbm19jIgH3gflRB8Am26B6oPEEiBiGwEPAJ8bIz5EOje1ipc6yPOubQeQA9jzCcJj98mHepW4GLqTnBLnXPz/HnmAb389n5oD7g4s/22ZhFmrAEArLXV1Gq23iMipwOPWGvL1uFlZSUikof2GTsf1abtDWwW9FgzQxPuTIvr6yBgjDkIWOic+9QYs3sax0rl8mr25xkMayAl1tqQwwrEfc9ejBsRuQktnvg9Knb9Slv14doQ0OBVRu5POwOHGGN+hrY7LzHG/AtYYIzp45ybZ4zpQ213itnARgn790dTDptFMKyBQMOsASoTnl8BRK21V66j61nvyUQqlXPuMuAyAD9jvdA592tjzI3ACeiK4wRUEwPgBeBRY8zNaPHMUFTHuFmEZoKBQCBr6JTT021bcGRaY98pvzutZoIJhvUgY0x34ElgY+BH4Gjn3FI/7nLgZDSu8HvnXMqGmukQDGsgEMgaOkV6ulEF9TUPrsu7FfdkbZfW4AoIBAJZRbZXVaVDMKyBQCBr0PbXwbAGAoFARgmGNRAIBDJK9guspEMwrIFAIGtwBmpM+08LDoY1EAhkFcEVEAgEAhlEHQHtX24hGNZAIJBVhBlrIBAIZJiYCTPWQCAQyBjqCAiGNRAIBDKII7oeSAAHwxoIBLIGR3AFBAKBQIYJroBAIBDIOI5QIBAIBAIZIwSvAoFAIOM4olSv64toMcGwBgKBrEFlA9v/jDW0vw4EAlmFI5bWoyGMMRsZY942xkw1xkw2xpzrt3czxvzXGPOt/9k1YZ/LjDHTjTHTjDH7teQ9BMMaCASyCEeMaFqPRqgBLnDObQaMAc4yxowALgXedM4NBd70z/GvHQuMBPYH7jDG5DT3XQTDGggEsgZHZmaszrl5zrn/+d/LgKlAP+BQ4CE/7CHgMP/7ocDjzrk1zrnvgenA9s19H8HHGggEsghHzKWdbtXDGPNJwvN7nHP3JA8yxgwEtgEmAqXOuXmgxtcY08sP6wd8mLDbbL+tWQTDGggEsghHLP2sgMWNdWk1xhQDz6DtrFcaY+odmvJimkkwrIFAIGtwQMxlJivAGJOHGtVHnHPP+s0LjDF9/Gy1D7DQb58NbJSwe39gbnPPHXysgUAgi3CZygowwH3AVOfczQkvvQCc4H8/AXg+YfuxxpgOxphBwFDgo+a+izBjDQQC2YMDl76PtSF2Bo4DvjTGfO63/R9wHfCkMeYU4EfgaADn3GRjzJPAFDSj4CzXggsxzrV/te5AILB+kBMpdMUFQ9Iau7Ji8qeN+VjXFWHGGggEsgqXIR/ruiQY1kAgkEU4Yi5oBQQCgUAGCV1aA4FAIKM4Mha8WqcEwxoIBLIIF3ysgUAgkGmCKyAQCAQyiXPEXOjSGggEAhkj+FgDgUAg42gPgfZOMKyBQCCrCMGrQCAQyCghjzUQCARagWBYA4FAIIM4XMgKCAQCgUwTZqyBQCCQWdYDKdNgWAOBQBbhcM1vNZU1BMMaCASyidegpkeaYxe36pW0gNBBIBAIBDJMaCYYCAQCGSYY1kAgEMgwwbAGAoFAhgmGNRAIBDJMMKyBQCCQYf4fl5/vF7NCCS0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# can check out https://www.bls.gov/oes/current/oes_stru.htm\n", "bio = '172031'\n", "bio_stats = oes_geo(bio)\n", "areas = get_areas() # this takes a few minutes\n", "state = state_albers()\n", "geo_bio = merge_occgeo(bio_stats,areas)\n", "\n", "ax = geo_bio.plot(column='Employment',cmap='inferno',legend=True,zorder=2)\n", "state.boundary.plot(ax=ax,color='grey',linewidth=0.5,zorder=1)\n", "ax.set_ylim(0.1*1e6,3.3*1e6)\n", "ax.set_xlim(-0.3*1e7,0.3*1e7) # lower 48 focus (for Albers proj)\n", "ax.set_axis_off()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 3, "id": "selective-mainstream", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
Make this Notebook Trusted to load map: File -> Trust Notebook
" ], "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Can export to excel (csv parsing is tough, maybe use tsv instead)\n", "bio_stats.to_excel('biojobs.xlsx')\n", "\n", "att = ['areaName','Employment','Location Quotient','Employment per 1,000 jobs','Annual mean wage']\n", "form = ['',',.0f','.2f','.2f',',.0f']\n", "\n", "map_bio = fol_map(geo_bio,'Employment',['lat', 'lon'],att,form)\n", "#map_bio.save('biomap.html')\n", "map_bio" ] }, { "cell_type": "code", "execution_count": 4, "id": "aquatic-video", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
areaNameEmploymentAnnual mean wageAnnual 75th percentile wage
120New York-Newark-Jersey City, NY-NJ-PA13860136040.0168040.0
156San Francisco-Oakland-Hayward, CA8890152080.0175580.0
102Los Angeles-Long Beach-Anaheim, CA7160120290.0147990.0
191Washington-Arlington-Alexandria, DC-VA-MD-WV6720126160.0164530.0
33Chicago-Naperville-Elgin, IL-IN-WI5100110800.0132900.0
44Dallas-Fort Worth-Arlington, TX5040108090.0135290.0
157San Jose-Sunnyvale-Santa Clara, CA4940233320.0NaN
10Atlanta-Sandy Springs-Roswell, GA4480110670.0134990.0
132Philadelphia-Camden-Wilmington, PA-NJ-DE-MD3680109460.0135590.0
165Seattle-Tacoma-Bellevue, WA3610142810.0174640.0
\n", "
" ], "text/plain": [ " areaName Employment \\\n", "120 New York-Newark-Jersey City, NY-NJ-PA 13860 \n", "156 San Francisco-Oakland-Hayward, CA 8890 \n", "102 Los Angeles-Long Beach-Anaheim, CA 7160 \n", "191 Washington-Arlington-Alexandria, DC-VA-MD-WV 6720 \n", "33 Chicago-Naperville-Elgin, IL-IN-WI 5100 \n", "44 Dallas-Fort Worth-Arlington, TX 5040 \n", "157 San Jose-Sunnyvale-Santa Clara, CA 4940 \n", "10 Atlanta-Sandy Springs-Roswell, GA 4480 \n", "132 Philadelphia-Camden-Wilmington, PA-NJ-DE-MD 3680 \n", "165 Seattle-Tacoma-Bellevue, WA 3610 \n", "\n", " Annual mean wage Annual 75th percentile wage \n", "120 136040.0 168040.0 \n", "156 152080.0 175580.0 \n", "102 120290.0 147990.0 \n", "191 126160.0 164530.0 \n", "33 110800.0 132900.0 \n", "44 108090.0 135290.0 \n", "157 233320.0 NaN \n", "10 110670.0 134990.0 \n", "132 109460.0 135590.0 \n", "165 142810.0 174640.0 " ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Data science jobs\n", "ds = '152051'\n", "ds_stats = oes_geo(ds)\n", "geo_ds = merge_occgeo(ds_stats,areas)\n", "geo_ds.sort_values(by='Employment',inplace=True,ascending=False)\n", "lim_fields = ['areaName','Employment','Annual mean wage','Annual 75th percentile wage']\n", "geo_ds[lim_fields].head(10)" ] }, { "cell_type": "code", "execution_count": null, "id": "b7aaa5f9", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.13" } }, "nbformat": 4, "nbformat_minor": 5 }