Skip to content

Instantly share code, notes, and snippets.

@danlamanna

danlamanna/gist:8c9a1f01bca28be1913c

Created September 16, 2015 16:21
Show Gist options
  • Save danlamanna/8c9a1f01bca28be1913c to your computer and use it in GitHub Desktop.
Save danlamanna/8c9a1f01bca28be1913c to your computer and use it in GitHub Desktop.
Download ZIP
Raw
gistfile1.txt
{
"cells": [
{
"cell_type": "code",
"execution_count": 59,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"b685670ffe9ff0eef2d26460523d3f04 /home/dan.lamanna/Downloads/tl_2012_us_cbsa/tl_2012_us_cbsa.shp\r\n"
]
}
],
"source": [
"data_file = '/home/dan.lamanna/Downloads/tl_2012_us_cbsa/tl_2012_us_cbsa.shp'\n",
"!md5sum $data_file"
]
},
{
"cell_type": "code",
"execution_count": 60,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import json\n",
"import pickle\n",
"import requests\n",
"import shapefile"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# What MSAs are we expected to have"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"with open('msanames.pickle', 'wb') as outfile:\n",
" pickle.dump(json.loads('{\"msaname\":[\"Fayetteville, NC MSA\",\"Worcester, MA MSA\",\"Killeen-Temple-Fort Hood, TX MSA\",\"Salinas, CA MSA\",\"Greeley, CO MSA\",\"Milwaukee-Waukesha-West Allis, WI MSA\",\"Athens-Clark County, GA MSA\",\"Tallahassee, FL MSA\",\"Waterloo-Cedar Falls, IA MSA\",\"Lexington-Fayette, KY MSA\",\"Stockton, CA MSA\",\"Topeka, KS MSA\",\"Louisville/Jefferson County, KY-IN MSA\",\"Allentown-Bethlehem-Easton, PA-NJ MSA\",\"Clarksville, TN-KY MSA\",\"Bellingham, WA MSA\",\"Reading, PA MSA\",\"Muncie, IN MSA\",\"Eugene-Springfield, OR MSA\",\"Memphis, TN-AR-MS MSA\",\"Gulfport-Biloxi, MS MSA\",\"Roanoke, VA MSA\",\"Santa Cruz-Watsonville, CA MSA\",\"Myrtle Beach-Conway-North Myrtle Beach, SC MSA\",\"Denver-Aurora, CO MSA\",\"Springfield, MO MSA\",\"Pittsburgh, PA MSA\",\"Binghamton, NY MSA\",\"Harrisburg-Carlisle, PA MSA\",\"Nashville-Davidson-Murfreesboro-Franklin, TN MSA\",\"Minneapolis-St. Paul-Bloomington, MN-WI MSA\",\"Anchorage, AK MSA\",\"Alexandria, LA MSA\",\"Sioux Falls, SD MSA\",\"Greenville-Mauldin-Easley, SC MSA\",\"Syracuse, NY MSA\",\"Redding, CA MSA\",\"Phoenix-Mesa-Scottsdale, AZ MSA\",\"Albany-Schenectady-Troy, NY MSA\",\"Philadelphia-Camden-Wilmington, PA-NJ-DE-MD MSA\",\"Oxnard-Thousand Oaks-Ventura, CA MSA\",\"Provo-Orem, UT MSA\",\"Houma-Bayou Cane-Thibodaux, LA MSA\",\"Atlantic City, NJ MSA\",\"Salt Lake City, UT MSA\",\"Merced, CA MSA\",\"South Bend-Mishawaka, IN-MI MSA\",\"Palm Bay-Melbourne-Titusville, FL MSA\",\"Miami-Fort Lauderdale-Pompano Beach, FL MSA\",\"Riverside-San Bernardino-Ontario, CA MSA\",\"Savannah, GA MSA\",\"Lubbock, TX MSA\",\"Portland-South Portland-Biddeford, ME MSA\",\"Orlando-Kissimmee, FL MSA\",\"Sacramento-Arden-Arcade-Roseville, CA MSA\",\"Amarillo, TX MSA\",\"Laredo, TX MSA\",\"Gadsden, AL MSA\",\"Lincoln, NE MSA\",\"Wausau, WI MSA\",\"Charlottesville, VA MSA\",\"Cape Coral-Fort Myers, FL MSA\",\"Johnstown, PA MSA\",\"Hagerstown-Martinsburg, MD-WV MSA\",\"Olympia, WA MSA\",\"Akron, OH MSA\",\"Las Vegas-Paradise, NV MSA\",\"Bridgeport-Stamford-Norwalk, CT MSA\",\"Seattle-Tacoma-Bellevue, WA MSA\",\"Williamsport, PA MSA\",\"Lansing-East Lansing, MI MSA\",\"Asheville, NC MSA\",\"Appleton, WI MSA\",\"New York-Northern New Jersey-Long Island, NY-NJ-PA MSA\",\"Davenport-Moline-Rock Island, IA-IL MSA\",\"Janesville, WI MSA\",\"Madison, WI MSA\",\"Cincinnati-Middletown, OH-KY-IN MSA\",\"Corpus Christi, TX MSA\",\"Detroit-Warren-Livonia, MI MSA\",\"San Antonio, TX MSA\",\"Salem, OR MSA\",\"Austin-Round Rock, TX MSA\",\"Lancaster, PA MSA\",\"Waco, TX MSA\",\"Chattanooga, TN-GA MSA\",\"Utica-Rome, NY MSA\",\"Pensacola-Ferry Pass-Brent, FL MSA\",\"Omaha-Council Bluffs, NE-IA MSA\",\"Baltimore-Towson, MD MSA\",\"New Orleans-Metairie-Kenner, LA MSA\",\"New Haven-Milford, CT MSA\",\"Medford, OR MSA\",\"Visalia-Porterville, CA MSA\",\"Little Rock-North Little Rock-Conway, AR MSA\",\"Macon, GA MSA\",\"Youngstown-Warren-Boardman, OH-PA MSA\",\"San Jose-Sunnyvale-Santa Clara, CA MSA\",\"Tulsa, OK MSA\",\"Indianapolis-Carmel, IN MSA\",\"Evansville, IN-KY MSA\",\"Fort Collins-Loveland, CO MSA\",\"Dayton, OH MSA\",\"Tampa-St. Petersburg-Clearwater, FL\",\"Houston-Sugar Land-Baytown, TX MSA\",\"Washington-Arlington-Alexandria, DC-VA-MD-WV MSA\",\"San Francisco-Oakland-Fremont, CA MSA\",\"Chicago-Naperville-Joliet, IL-IN-WI MSA\",\"Modesto, CA MSA\",\"St. Joseph, MO-KS MSA\",\"Kansas City, MO-KS MSA\",\"Abilene, TX MSA\",\"Boston-Cambridge-Quincy, MA-NH MSA\",\"Charlotte-Gastonia-Concord, NC-SC MSA\",\"Buffalo-Niagara Falls, NY MSA\",\"Wichita, KS MSA\",\"Providence-New Bedford-Fall River, RI-MA MSA\",\"Raleigh-Cary, NC MSA\",\"Flint, MI MSA\",\"Flagstaff, AZ MSA\",\"Joplin, MO MSA\",\"Knoxville, TN MSA\",\"Columbus, OH MSA\",\"Monroe, LA MSA\",\"Chico, CA MSA\",\"Cedar Rapids, IA MSA\",\"Tyler, TX MSA\",\"McAllen-Edinburg-Mission, TX MSA\",\"Cleveland-Elyria-Mentor, OH MSA\",\"Oklahoma City, OK MSA\",\"Toledo, OH MSA\",\"Albuquerque, NM MSA\",\"Hickory-Lenoir-Morganton, NC MSA\",\"Fargo, ND-MN MSA\",\"Pueblo, CO MSA\",\"Portland-Vancouver-Beaverton, OR-WA MSA\",\"Scranton--Wilkes-Barre, PA MSA\",\"Canton-Massillon, OH MSA\",\"Augusta-Richmond County, GA-SC MSA\",\"Richmond, VA MSA\",\"San Luis Obispo-Paso Robles, CA MSA\",\"Eau Claire, WI MSA\",\"Hattiesburg, MS MSA\",\"Yuma, AZ MSA\",\"Rocky Mount, NC MSA\",\"Greensboro-High Point, NC MSA\",\"Reno-Sparks, NV MSA\",\"Dallas-Fort Worth-Arlington, TX MSA\",\"Terre Haute, IN MSA\",\"San Diego-Carlsbad-San Marcos, CA MSA\",\"Beaumont-Port Arthur, TX MSA\",\"Boise City-Nampa, ID MSA\",\"Tucson, AZ MSA\",\"Birmingham-Hoover, AL MSA\",\"Jacksonville, FL MSA\",\"El Paso, TX MSA\",\"Shreveport-Bossier City, LA MSA\",\"Kokomo, IN MSA\",\"Altoona, PA MSA\",\"Brownsville-Harlingen, TX MSA\",\"Anniston-Oxford, AL MSA\",\"Bakersfield, CA MSA\",\"Virginia Beach-Norfolk-Newport News, VA-NC MSA\",\"Duluth, MN-WI MSA\",\"St. Louis, MO-IL MSA\"]}'), outfile)"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Reverting to pickle..\n"
]
}
],
"source": [
"r = requests.get('https://tempus-demo.ngrok.com/api/groups',\n",
" params={'table': 'escort_ads',\n",
" 'group_col': 'msaname'})\n",
"\n",
"if not r.ok:\n",
" print 'Reverting to pickle..'\n",
" \n",
" with open('msanames.pickle', 'rb') as infile:\n",
" msanames = pickle.load(infile)\n",
"else:\n",
" msanames = r.json()"
]
},
{
"cell_type": "code",
"execution_count": 22,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"required_msa = set([msa.replace(' MSA', '') for msa in msanames['msaname']])"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"165"
]
},
"execution_count": 25,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"len(required_msa)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# What MSAs do we actually have"
]
},
{
"cell_type": "code",
"execution_count": 26,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"shp_reader = shapefile.Reader(data_file)"
]
},
{
"cell_type": "code",
"execution_count": 31,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"obtained_msa = set([record[3] for record in shp_reader.records()])"
]
},
{
"cell_type": "code",
"execution_count": 32,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"(955, 165, 15)"
]
},
"execution_count": 32,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"len(obtained_msa), len(required_msa), len(required_msa - obtained_msa)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Generate GeoJSON from the obtained MSAs"
]
},
{
"cell_type": "code",
"execution_count": 37,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"!ogr2ogr -f GeoJSON -t_srs wgs84 tl_2012_us_cbsa.geojson $data_file"
]
},
{
"cell_type": "code",
"execution_count": 61,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"with open('tl_2012_us_cbsa.geojson', 'rb') as infile:\n",
" msas = json.load(infile)"
]
},
{
"cell_type": "code",
"execution_count": 62,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"{u'geometry': {u'coordinates': [[[-155.963465, 19.082072],\n",
" [-155.966489, 19.08848],\n",
" [-155.968892, 19.093265],\n",
" [-155.972538, 19.101382],\n",
" [-155.974041, 19.107789],\n",
" [-155.974474, 19.112743],\n",
" [-155.974277, 19.118636],\n",
" [-155.973452, 19.125],\n",
" [-155.974704, 19.125],\n",
" [-155.974738, 19.130161],\n",
" [-155.973882, 19.14257],\n",
" [-155.970554, 19.159255],\n",
" [-155.96945, 19.169695],\n",
" [-155.968079, 19.179963],\n",
" [-155.966067, 19.191258],\n",
" [-155.964162, 19.201012],\n",
" [-155.962087, 19.209824],\n",
" [-155.958526, 19.223],\n",
" [-155.95673, 19.230957],\n",
" [-155.956511, 19.234465],\n",
" [-155.956652, 19.23806],\n",
" [-155.956602, 19.242597],\n",
" [-155.956556, 19.24679],\n",
" [-155.956937, 19.25],\n",
" [-155.957202, 19.25201],\n",
" [-155.9562, 19.25995],\n",
" [-155.955298, 19.263451],\n",
" [-155.952859, 19.272932],\n",
" [-155.948894, 19.283696],\n",
" [-155.946461, 19.290445],\n",
" [-155.945104, 19.29668],\n",
" [-155.944644, 19.302785],\n",
" [-155.944275, 19.307692],\n",
" [-155.944256, 19.316418],\n",
" [-155.94425, 19.319471],\n",
" [-155.944328, 19.32519],\n",
" [-155.944229, 19.329629],\n",
" [-155.943755, 19.332093],\n",
" [-155.942326, 19.339535],\n",
" [-155.943213, 19.344227],\n",
" [-155.947229, 19.35224],\n",
" [-155.950711, 19.35906],\n",
" [-155.95437, 19.366904],\n",
" [-155.956065, 19.371338],\n",
" [-155.957535, 19.375],\n",
" [-155.957985, 19.375715],\n",
" [-155.961376, 19.381101],\n",
" [-155.963615, 19.386519],\n",
" [-155.964299, 19.389948],\n",
" [-155.965017, 19.393551],\n",
" [-155.965788, 19.40066],\n",
" [-155.967951, 19.407537],\n",
" [-155.968084, 19.40796],\n",
" [-155.97158, 19.418526],\n",
" [-155.972759, 19.428035],\n",
" [-155.974597, 19.435842],\n",
" [-155.978105, 19.445809],\n",
" [-155.983319, 19.450752],\n",
" [-155.989108, 19.454995],\n",
" [-155.990328, 19.455889],\n",
" [-155.994997, 19.461253],\n",
" [-155.999134, 19.466181],\n",
" [-156.0, 19.469029],\n",
" [-156.000111, 19.469308],\n",
" [-156.000447, 19.470145],\n",
" [-156.000559, 19.470425],\n",
" [-156.002939, 19.476373],\n",
" [-156.007959, 19.488914],\n",
" [-156.010271, 19.49414],\n",
" [-156.012864, 19.5],\n",
" [-156.012653, 19.5],\n",
" [-156.01202, 19.5],\n",
" [-156.011809, 19.5],\n",
" [-156.011808, 19.500376],\n",
" [-156.01366, 19.507812],\n",
" [-156.014375, 19.51068],\n",
" [-156.015774, 19.522086],\n",
" [-156.016806, 19.531701],\n",
" [-156.016907, 19.532641],\n",
" [-156.019396, 19.539285],\n",
" [-156.019858, 19.541947],\n",
" [-156.021246, 19.549927],\n",
" [-156.021246, 19.549934],\n",
" [-156.021589, 19.552615],\n",
" [-156.022699, 19.561306],\n",
" [-156.022997, 19.563632],\n",
" [-156.02583, 19.575382],\n",
" [-156.029921, 19.58651],\n",
" [-156.030277, 19.587478],\n",
" [-156.033929, 19.594293],\n",
" [-156.034738, 19.595803],\n",
" [-156.035895, 19.597961],\n",
" [-156.037815, 19.599847],\n",
" [-156.038935, 19.600948],\n",
" [-156.039075, 19.600979],\n",
" [-156.043237, 19.601896],\n",
" [-156.043279, 19.601914],\n",
" [-156.049329, 19.604467],\n",
" [-156.055054, 19.607234],\n",
" [-156.057121, 19.608233],\n",
" [-156.058855, 19.609256],\n",
" [-156.06321, 19.611823],\n",
" [-156.071087, 19.617461],\n",
" [-156.075293, 19.620962],\n",
" [-156.078592, 19.625001],\n",
" [-156.082882, 19.632659],\n",
" [-156.085427, 19.638842],\n",
" [-156.087628, 19.646083],\n",
" [-156.08795, 19.651952],\n",
" [-156.088287, 19.658435],\n",
" [-156.08882, 19.661408],\n",
" [-156.090902, 19.663653],\n",
" [-156.09332, 19.668346],\n",
" [-156.097823, 19.675374],\n",
" [-156.101682, 19.686358],\n",
" [-156.104427, 19.693151],\n",
" [-156.104742, 19.694812],\n",
" [-156.109583, 19.700432],\n",
" [-156.111219, 19.703297],\n",
" [-156.115179, 19.710685],\n",
" [-156.119075, 19.719301],\n",
" [-156.120001, 19.727176],\n",
" [-156.119298, 19.736482],\n",
" [-156.11792, 19.744571],\n",
" [-156.114727, 19.75],\n",
" [-156.115857, 19.75],\n",
" [-156.114858, 19.75459],\n",
" [-156.112599, 19.759765],\n",
" [-156.110521, 19.764601],\n",
" [-156.109163, 19.769515],\n",
" [-156.108432, 19.776788],\n",
" [-156.108061, 19.784142],\n",
" [-156.10724, 19.790402],\n",
" [-156.105518, 19.797262],\n",
" [-156.103078, 19.802777],\n",
" [-156.099464, 19.808978],\n",
" [-156.095672, 19.815013],\n",
" [-156.090616, 19.821143],\n",
" [-156.086104, 19.825831],\n",
" [-156.079517, 19.831131],\n",
" [-156.07248, 19.835251],\n",
" [-156.065627, 19.839115],\n",
" [-156.058774, 19.841881],\n",
" [-156.053544, 19.843534],\n",
" [-156.049208, 19.848643],\n",
" [-156.045052, 19.854425],\n",
" [-156.039363, 19.860643],\n",
" [-156.033943, 19.867028],\n",
" [-156.030059, 19.871457],\n",
" [-156.027348, 19.875],\n",
" [-156.024849, 19.875],\n",
" [-156.02306, 19.877898],\n",
" [-156.019654, 19.882247],\n",
" [-156.018069, 19.883959],\n",
" [-156.015709, 19.886511],\n",
" [-156.015002, 19.887201],\n",
" [-156.012391, 19.889753],\n",
" [-156.008081, 19.892741],\n",
" [-156.00408, 19.894793],\n",
" [-156.002419, 19.895646],\n",
" [-156.0, 19.896586],\n",
" [-155.998344, 19.896876],\n",
" [-155.986097, 19.902663],\n",
" [-155.985414, 19.902865],\n",
" [-155.97177, 19.906912],\n",
" [-155.962427, 19.909152],\n",
" [-155.959936, 19.915244],\n",
" [-155.95406, 19.929714],\n",
" [-155.953553, 19.931414],\n",
" [-155.950381, 19.942075],\n",
" [-155.94793, 19.94622],\n",
" [-155.94753, 19.946895],\n",
" [-155.946331, 19.948922],\n",
" [-155.945932, 19.949598],\n",
" [-155.940118, 19.959429],\n",
" [-155.938883, 19.961518],\n",
" [-155.924263, 19.978653],\n",
" [-155.918986, 19.983966],\n",
" [-155.917489, 19.985013],\n",
" [-155.91346, 19.987836],\n",
" [-155.908278, 19.991764],\n",
" [-155.908018, 19.991961],\n",
" [-155.907488, 19.993465],\n",
" [-155.907124, 19.9945],\n",
" [-155.905531, 19.99712],\n",
" [-155.904607, 19.998429],\n",
" [-155.9035, 20.0],\n",
" [-155.903157, 20.0],\n",
" [-155.90213, 20.0],\n",
" [-155.901788, 20.0],\n",
" [-155.899257, 20.004107],\n",
" [-155.897865, 20.005632],\n",
" [-155.897635, 20.005886],\n",
" [-155.8961, 20.008006],\n",
" [-155.894841, 20.009021],\n",
" [-155.893489, 20.010461],\n",
" [-155.891152, 20.012065],\n",
" [-155.889444, 20.013333],\n",
" [-155.887376, 20.014854],\n",
" [-155.885843, 20.016633],\n",
" [-155.884311, 20.018327],\n",
" [-155.885276, 20.020026],\n",
" [-155.886348, 20.021914],\n",
" [-155.887767, 20.024135],\n",
" [-155.888842, 20.025919],\n",
" [-155.889362, 20.026782],\n",
" [-155.89114, 20.02892],\n",
" [-155.893892, 20.032765],\n",
" [-155.894206, 20.033429],\n",
" [-155.895306, 20.035752],\n",
" [-155.897973, 20.038916],\n",
" [-155.900371, 20.042164],\n",
" [-155.903036, 20.045668],\n",
" [-155.906947, 20.050456],\n",
" [-155.911187, 20.055281],\n",
" [-155.911305, 20.055416],\n",
" [-155.91495, 20.059946],\n",
" [-155.916878, 20.062573],\n",
" [-155.918148, 20.064304],\n",
" [-155.92063, 20.068658],\n",
" [-155.923112, 20.072842],\n",
" [-155.925251, 20.074641],\n",
" [-155.927566, 20.076951],\n",
" [-155.929521, 20.0796],\n",
" [-155.930854, 20.081394],\n",
" [-155.932275, 20.083445],\n",
" [-155.9336, 20.086601],\n",
" [-155.934924, 20.090012],\n",
" [-155.935802, 20.093166],\n",
" [-155.936592, 20.095894],\n",
" [-155.936737, 20.096474],\n",
" [-155.93738, 20.099048],\n",
" [-155.937639, 20.100666],\n",
" [-155.93745, 20.102453],\n",
" [-155.937888, 20.104157],\n",
" [-155.938599, 20.105268],\n",
" [-155.940111, 20.106978],\n",
" [-155.941712, 20.108944],\n",
" [-155.94269, 20.110396],\n",
" [-155.943929, 20.112955],\n",
" [-155.944724, 20.114917],\n",
" [-155.945517, 20.117219],\n",
" [-155.946308, 20.119947],\n",
" [-155.946657, 20.121651],\n",
" [-155.947988, 20.123786],\n",
" [-155.947816, 20.125],\n",
" [-155.948782, 20.128744],\n",
" [-155.949482, 20.13213],\n",
" [-155.950001, 20.136276],\n",
" [-155.95079, 20.139495],\n",
" [-155.951398, 20.143557],\n",
" [-155.952802, 20.149063],\n",
" [-155.953768, 20.152536],\n",
" [-155.954996, 20.157533],\n",
" [-155.956312, 20.1627],\n",
" [-155.957719, 20.167614],\n",
" [-155.957882, 20.171757],\n",
" [-155.958838, 20.178442],\n",
" [-155.95962, 20.184027],\n",
" [-155.95996, 20.188509],\n",
" [-155.960217, 20.191554],\n",
" [-155.960295, 20.194513],\n",
" [-155.960818, 20.197729],\n",
" [-155.960095, 20.201019],\n",
" [-155.959817, 20.204481],\n",
" [-155.958916, 20.2076],\n",
" [-155.958377, 20.209117],\n",
" [-155.958821, 20.209544],\n",
" [-155.958971, 20.211379],\n",
" [-155.959077, 20.212673],\n",
" [-155.959247, 20.215041],\n",
" [-155.959507, 20.217071],\n",
" [-155.959492, 20.221466],\n",
" [-155.959215, 20.224422],\n",
" [-155.959029, 20.226871],\n",
" [-155.958125, 20.231004],\n",
" [-155.957667, 20.234972],\n",
" [-155.956408, 20.238765],\n",
" [-155.955508, 20.241546],\n",
" [-155.95443, 20.244579],\n",
" [-155.953621, 20.247023],\n",
" [-155.952636, 20.248789],\n",
" [-155.951919, 20.25],\n",
" [-155.951725, 20.25],\n",
" [-155.951143, 20.25],\n",
" [-155.95095, 20.25],\n",
" [-155.950539, 20.252328],\n",
" [-155.950285, 20.253772],\n",
" [-155.949994, 20.255978],\n",
" [-155.949255, 20.258095],\n",
" [-155.948844, 20.259178],\n",
" [-155.948517, 20.260043],\n",
" [-155.948047, 20.261404],\n",
" [-155.947962, 20.261652],\n",
" [-155.947228, 20.263175],\n",
" [-155.946133, 20.264864],\n",
" [-155.944853, 20.267062],\n",
" [-155.943664, 20.269176],\n",
" [-155.942841, 20.270528],\n",
" [-155.942381, 20.271628],\n",
" [-155.941013, 20.273655],\n",
" [-155.940285, 20.274499],\n",
" [-155.939375, 20.275595],\n",
" [-155.938734, 20.276779],\n",
" [-155.93673, 20.279395],\n",
" [-155.935637, 20.28083],\n",
" [-155.93409, 20.282686],\n",
" [-155.932365, 20.284285],\n",
" [-155.931178, 20.286144],\n",
" [-155.929538, 20.288338],\n",
" [-155.928222, 20.289883],\n",
" [-155.928173, 20.289941],\n",
" [-155.926446, 20.291795],\n",
" [-155.924178, 20.293729],\n",
" [-155.921817, 20.296003],\n",
" [-155.919366, 20.29819],\n",
" [-155.917642, 20.299704],\n",
" [-155.91592, 20.301049],\n",
" [-155.914291, 20.301969],\n",
" [-155.912479, 20.303228],\n",
" [-155.910853, 20.303893],\n",
" [-155.908501, 20.305147],\n",
" [-155.906873, 20.306067],\n",
" [-155.905247, 20.306733],\n",
" [-155.903168, 20.307649],\n",
" [-155.901541, 20.3084],\n",
" [-155.900006, 20.308896],\n",
" [-155.898651, 20.309394],\n",
" [-155.8963, 20.310648],\n",
" [-155.894673, 20.311313],\n",
" [-155.892413, 20.312398],\n",
" [-155.890607, 20.312977],\n",
" [-155.887989, 20.313888],\n",
" [-155.885912, 20.314549],\n",
" [-155.883656, 20.315209],\n",
" [-155.881674, 20.315362],\n",
" [-155.878878, 20.316017],\n",
" [-155.876622, 20.316591],\n",
" [-155.875, 20.316832],\n",
" [-155.870728, 20.318315],\n",
" [-155.866691, 20.319305],\n",
" [-155.860953, 20.320284],\n",
" [-155.856653, 20.320764],\n",
" [-155.851635, 20.321409],\n",
" [-155.84546, 20.321453],\n",
" [-155.840272, 20.321334],\n",
" [-155.836783, 20.321226],\n",
" [-155.83375, 20.320444],\n",
" [-155.831824, 20.320202],\n",
" [-155.829461, 20.319907],\n",
" [-155.825713, 20.318951],\n",
" [-155.821702, 20.317654],\n",
" [-155.8184, 20.316869],\n",
" [-155.814744, 20.315744],\n",
" [-155.812088, 20.31327],\n",
" [-155.809941, 20.313255],\n",
" [-155.80476, 20.312458],\n",
" [-155.80128, 20.311672],\n",
" [-155.798758, 20.310994],\n",
" [-155.798044, 20.310787],\n",
" [-155.795551, 20.309957],\n",
" [-155.79289, 20.30927],\n",
" [-155.787521, 20.306827],\n",
" [-155.783762, 20.305229],\n",
" [-155.780183, 20.30346],\n",
" [-155.778037, 20.301772],\n",
" [-155.776069, 20.300506],\n",
" [-155.774818, 20.298645],\n",
" [-155.771685, 20.29857],\n",
" [-155.768999, 20.298578],\n",
" [-155.766312, 20.298077],\n",
" [-155.764075, 20.297998],\n",
" [-155.761616, 20.297377],\n",
" [-155.758864, 20.296725],\n",
" [-155.755475, 20.295769],\n",
" [-155.753424, 20.295162],\n",
" [-155.751464, 20.29447],\n",
" [-155.75, 20.293527],\n",
" [-155.748148, 20.292489],\n",
" [-155.746839, 20.291286],\n",
" [-155.74318, 20.290134],\n",
" [-155.7404, 20.288925],\n",
" [-155.738333, 20.287493],\n",
" [-155.734791, 20.285413],\n",
" [-155.731427, 20.283255],\n",
" [-155.726865, 20.280122],\n",
" [-155.72474, 20.27793],\n",
" [-155.722541, 20.275316],\n",
" [-155.721296, 20.273161],\n",
" [-155.719571, 20.272081],\n",
" [-155.717968, 20.27008],\n",
" [-155.71475, 20.268769],\n",
" [-155.71266, 20.268008],\n",
" [-155.710324, 20.266733],\n",
" [-155.707991, 20.265376],\n",
" [-155.705145, 20.263575],\n",
" [-155.702217, 20.261602],\n",
" [-155.69991, 20.259572],\n",
" [-155.698205, 20.257902],\n",
" [-155.6966, 20.255985],\n",
" [-155.695008, 20.253731],\n",
" [-155.693154, 20.251384],\n",
" [-155.692566, 20.25],\n",
" [-155.690551, 20.247898],\n",
" [-155.690012, 20.246886],\n",
" [-155.689383, 20.24562],\n",
" [-155.689125, 20.245216],\n",
" [-155.688845, 20.244776],\n",
" [-155.688044, 20.244439],\n",
" [-155.687682, 20.244325],\n",
" [-155.68671, 20.244018],\n",
" [-155.684844, 20.243682],\n",
" [-155.683331, 20.243177],\n",
" [-155.682536, 20.243123],\n",
" [-155.682088, 20.243093],\n",
" [-155.680888, 20.242544],\n",
" [-155.680883, 20.242541],\n",
" [-155.680878, 20.242539],\n",
" [-155.680871, 20.242535],\n",
" [-155.680867, 20.242534],\n",
" [-155.680585, 20.242405],\n",
" [-155.679739, 20.242018],\n",
" [-155.679458, 20.241889],\n",
" [-155.679328, 20.241829],\n",
" [-155.677726, 20.241239],\n",
" [-155.675854, 20.23989],\n",
" [-155.674426, 20.238709],\n",
" [-155.672197, 20.237107],\n",
" [-155.670766, 20.235251],\n",
" [-155.668803, 20.233649],\n",
" [-155.666129, 20.231793],\n",
" [-155.664792, 20.230781],\n",
" [-155.663543, 20.229853],\n",
" [-155.662736, 20.228419],\n",
" [-155.660863, 20.226985],\n",
" [-155.659436, 20.225804],\n",
" [-155.658723, 20.225383],\n",
" [-155.656679, 20.224962],\n",
" [-155.655167, 20.224541],\n",
" [-155.653567, 20.224119],\n",
" [-155.65179, 20.223867],\n",
" [-155.649476, 20.222771],\n",
" [-155.648586, 20.222349],\n",
" [-155.647519, 20.222012],\n",
" [-155.646805, 20.221422],\n",
" [-155.645941, 20.221012],\n",
" [-155.645559, 20.220831],\n",
" [-155.644577, 20.219903],\n",
" [-155.643774, 20.219228],\n",
" [-155.642264, 20.218976],\n",
" [-155.641464, 20.218807],\n",
" [-155.640752, 20.21847],\n",
" [-155.639863, 20.218217],\n",
" [-155.639238, 20.217711],\n",
" [-155.63826, 20.217374],\n",
" [-155.637459, 20.216952],\n",
" [-155.636124, 20.216362],\n",
" [-155.63452, 20.21535],\n",
" [-155.633097, 20.214844],\n",
" [-155.631851, 20.214253],\n",
" [-155.630515, 20.213494],\n",
" [-155.628555, 20.212313],\n",
" [-155.626506, 20.210964],\n",
" [-155.625, 20.209951],\n",
" [-155.622555, 20.208451],\n",
" [-155.62085, 20.206858],\n",
" [-155.618158, 20.205352],\n",
" [-155.615734, 20.203509],\n",
" [-155.612862, 20.200994],\n",
" [-155.610259, 20.198731],\n",
" [-155.608734, 20.197474],\n",
" [-155.607387, 20.195879],\n",
" [-155.606399, 20.193948],\n",
" [-155.604066, 20.193029],\n",
" [-155.601823, 20.192026],\n",
" [-155.599312, 20.191445],\n",
" [-155.597787, 20.19086],\n",
" [-155.595095, 20.189018],\n",
" [-155.592673, 20.187426],\n",
" [-155.59043, 20.186339],\n",
" [-155.588008, 20.184748],\n",
" [-155.585765, 20.183829],\n",
" [-155.582087, 20.181316],\n",
" [-155.579844, 20.180397],\n",
" [-155.578768, 20.180147],\n",
" [-155.57599, 20.181333],\n",
" [-155.573749, 20.182264],\n",
" [-155.570611, 20.182946],\n",
" [-155.568549, 20.182951],\n",
" [-155.564872, 20.18313],\n",
" [-155.561734, 20.183643],\n",
" [-155.558237, 20.183063],\n",
" [-155.555816, 20.182733],\n",
" [-155.552946, 20.18232],\n",
" [-155.548642, 20.181827],\n",
" [-155.545772, 20.181246],\n",
" [-155.54075, 20.180249],\n",
" [-155.538262, 20.179512],\n",
" [-155.536804, 20.179081],\n",
" [-155.531244, 20.177413],\n",
" [-155.526222, 20.175659],\n",
" [-155.5212, 20.173989],\n",
" [-155.517881, 20.172483],\n",
" [-155.515819, 20.171899],\n",
" [-155.513846, 20.17081],\n",
" [-155.512501, 20.170309],\n",
" [-155.510348, 20.169304],\n",
" [-155.508376, 20.168468],\n",
" [-155.505417, 20.168054],\n",
" [-155.50443, 20.167804],\n",
" [-155.502995, 20.167471],\n",
" [-155.501471, 20.167222],\n",
" [-155.500485, 20.167225],\n",
" [-155.499997, 20.167056],\n",
" [-155.498362, 20.166308],\n",
" [-155.496042, 20.16596],\n",
" [-155.493187, 20.165356],\n",
" [-155.490064, 20.164836],\n",
" [-155.486137, 20.164397],\n",
" [-155.486026, 20.164357],\n",
" [-155.482037, 20.162944],\n",
" [-155.478824, 20.162507],\n",
" [-155.475884, 20.161228],\n",
" [-155.472228, 20.160198],\n",
" [-155.469379, 20.158582],\n",
" [-155.466705, 20.157472],\n",
" [-155.465016, 20.156114],\n",
" [-155.462877, 20.155176],\n",
" [-155.458772, 20.154735],\n",
" [-155.45458, 20.153872],\n",
" [-155.450477, 20.152924],\n",
" [-155.447178, 20.151896],\n",
" [-155.445581, 20.151439],\n",
" [-155.443879, 20.150952],\n",
" [-155.441384, 20.150096],\n",
" [-155.438177, 20.148815],\n",
" [-155.436042, 20.147201],\n",
" [-155.43346, 20.145669],\n",
" [-155.432862, 20.145207],\n",
" [-155.431593, 20.144226],\n",
" [-155.429451, 20.144046],\n",
" [-155.424455, 20.143263],\n",
" [-155.422804, 20.142711],\n",
" [-155.422405, 20.142578],\n",
" [-155.418127, 20.141122],\n",
" [-155.414207, 20.139499],\n",
" [-155.412073, 20.137885],\n",
" [-155.409755, 20.137198],\n",
" [-155.405565, 20.135996],\n",
" [-155.400219, 20.134028],\n",
" [-155.397547, 20.132833],\n",
" [-155.394256, 20.130538],\n",
" [-155.393924, 20.130349],\n",
" [-155.390963, 20.128664],\n",
" [-155.387315, 20.126536],\n",
" [-155.384829, 20.125003],\n",
" [-155.384608, 20.125003],\n",
" [-155.383947, 20.125003],\n",
" [-155.383727, 20.125003],\n",
" [-155.378513, 20.122865],\n",
" [-155.377078, 20.121918],\n",
" [-155.375, 20.120286],\n",
" [-155.373579, 20.119228],\n",
" [-155.37305, 20.118834],\n",
" [-155.37156, 20.117185],\n",
" [-155.369702, 20.115953],\n",
" [-155.367743, 20.115229],\n",
" [-155.364711, 20.11431],\n",
" [-155.362313, 20.113151],\n",
" [-155.359371, 20.112234],\n",
" [-155.357324, 20.111423],\n",
" [-155.3551, 20.110522],\n",
" [-155.35199, 20.109092],\n",
" [-155.348435, 20.107481],\n",
" [-155.345964, 20.105555],\n",
" [-155.343398, 20.103882],\n",
" [-155.342486, 20.103326],\n",
" [-155.34021, 20.101939],\n",
" [-155.337457, 20.100601],\n",
" [-155.335244, 20.099191],\n",
" [-155.332215, 20.098187],\n",
" [-155.332003, 20.09811],\n",
" [-155.330958, 20.09773],\n",
" [-155.329721, 20.09728],\n",
" [-155.327876, 20.096461],\n",
" [-155.327499, 20.096294],\n",
" [-155.326838, 20.096069],\n",
" [-155.325362, 20.095565],\n",
" [-155.322865, 20.094828],\n",
" [-155.321483, 20.094109],\n",
" [-155.320469, 20.093583],\n",
" [-155.317083, 20.092485],\n",
" [-155.314418, 20.091234],\n",
" [-155.31211, 20.090076],\n",
" [-155.309894, 20.088835],\n",
" [-155.307602, 20.086998],\n",
" [-155.306409, 20.086163],\n",
" [-155.305216, 20.085328],\n",
" [-155.301671, 20.083291],\n",
" [-155.301516, 20.083209],\n",
" [-155.298301, 20.081514],\n",
" [-155.294741, 20.080157],\n",
" [-155.291626, 20.07898],\n",
" [-155.291232, 20.078833],\n",
" [-155.286998, 20.077257],\n",
" [-155.284336, 20.075921],\n",
" [-155.279276, 20.073424],\n",
" [-155.276345, 20.072081],\n",
" [-155.271708, 20.070783],\n",
" [-155.269036, 20.069871],\n",
" [-155.265294, 20.068678],\n",
" [-155.261643, 20.067403],\n",
" [-155.259998, 20.066709],\n",
" [-155.258443, 20.066054],\n",
" [-155.254621, 20.064435],\n",
" [-155.252399, 20.063448],\n",
" [-155.25, 20.061777],\n",
" [-155.24961, 20.061341],\n",
" [-155.248505, 20.060104],\n",
" [-155.248422, 20.060063],\n",
" [-155.247897, 20.059808],\n",
" [-155.245501, 20.058643],\n",
" [-155.241259, 20.056494],\n",
" [-155.237552, 20.054006],\n",
" [-155.236732, 20.053454],\n",
" [-155.234109, 20.05169],\n",
" [-155.22837, 20.048085],\n",
" [-155.222371, 20.043454],\n",
" [-155.218132, 20.041135],\n",
" [-155.215039, 20.039758],\n",
" [-155.211768, 20.038466],\n",
" [-155.209475, 20.03658],\n",
" [-155.206123, 20.034092],\n",
" [-155.204331, 20.032785],\n",
" [-155.204006, 20.032548],\n",
" [-155.202948, 20.03152],\n",
" [-155.201982, 20.030236],\n",
" [-155.200578, 20.027926],\n",
" [-155.197662, 20.026635],\n",
" [-155.195364, 20.025688],\n",
" [-155.194378, 20.025473],\n",
" [-155.193771, 20.02534],\n",
" [-155.190328, 20.023194],\n",
" [-155.187858, 20.021563],\n",
" [-155.184503, 20.019501],\n",
" [-155.182387, 20.017872],\n",
" [-155.18125, 20.015477],\n",
" [-155.17798, 20.014099],\n",
" [-155.176846, 20.013106],\n",
" [-155.174759, 20.011892],\n",
" [-155.173736, 20.010735],\n",
" [-155.173143, 20.01014],\n",
" [-155.170818, 20.008945],\n",
" [-155.168753, 20.007852],\n",
" [-155.166803, 20.006424],\n",
" [-155.164306, 20.003925],\n",
" [-155.16269, 20.001574],\n",
" [-155.161184, 20.000178],\n",
" [-155.161058, 20.000001],\n",
" [-155.159393, 20.000001],\n",
" [-155.159166, 19.999744],\n",
" [-155.157486, 19.997608],\n",
" [-155.155805, 19.997671],\n",
" [-155.153701, 19.996108],\n",
" [-155.151058, 19.995015],\n",
" [-155.14622, 19.992666],\n",
" [-155.141911, 19.990402],\n",
" [-155.14033, 19.989481],\n",
" [-155.137543, 19.986315],\n",
" [-155.135281, 19.983567],\n",
" [-155.133977, 19.981984],\n",
" [-155.132837, 19.981231],\n",
" [-155.128796, 19.978886],\n",
" [-155.126867, 19.977549],\n",
" [-155.123905, 19.974298],\n",
" [-155.121683, 19.971052],\n",
" [-155.113585, 19.964665],\n",
" [-155.106264, 19.957824],\n",
" [-155.091533, 19.946651],\n",
" [-155.090018, 19.945281],\n",
" [-155.086049, 19.941743],\n",
" [-155.08317, 19.937459],\n",
" [-155.081796, 19.935415],\n",
" [-155.07992, 19.933167],\n",
" [-155.076109, 19.928598],\n",
" [-155.072665, 19.922765],\n",
" [-155.07133, 19.919586],\n",
" [-155.070843, 19.918425],\n",
" [-155.067451, 19.916068],\n",
" [-155.066487, 19.915398],\n",
" [-155.063597, 19.91339],\n",
" [-155.062752, 19.912804],\n",
" [-155.062642, 19.912711],\n",
" [-155.055552, 19.906672],\n",
" [-155.046811, 19.897704],\n",
" [-155.041831, 19.889371],\n",
" [-155.037401, 19.882876],\n",
" [-155.03162, 19.875],\n",
" [-155.032994, 19.875],\n",
" [-155.034618, 19.875],\n",
" [-155.031522, 19.866089],\n",
" [-155.030649, 19.857674],\n",
" [-155.030299, 19.848995],\n",
" [-155.030278, 19.839463],\n",
" [-155.032941, 19.828013],\n",
" [-155.034707, 19.822067],\n",
" [-155.034748, 19.820463],\n",
" [-155.034321, 19.818027],\n",
" [-155.033887, 19.8123],\n",
" [-155.033877, 19.812176],\n",
" [-155.033849, 19.811806],\n",
" [-155.03384, 19.811683],\n",
" [-155.033734, 19.810281],\n",
" [-155.033603, 19.807319],\n",
" [-155.033171, 19.797473],\n",
" [-155.033286, 19.794221],\n",
" [-155.033315, 19.793421],\n",
" [-155.029999, 19.792101],\n",
" [-155.029914, 19.792067],\n",
" [-155.02669, 19.790709],\n",
" [-155.025523, 19.790146],\n",
" [-155.024311, 19.790226],\n",
" [-155.013697, 19.79093],\n",
" [-155.007044, 19.790243],\n",
" [-155.006073, 19.790149],\n",
" [-155.000002, 19.789563],\n",
" [-154.996757, 19.78917],\n",
" [-154.99472, 19.788924],\n",
" [-154.987489, 19.786126],\n",
" [-154.985888, 19.785507],\n",
" [-154.984507, 19.784802],\n",
" [-154.977849, 19.781404],\n",
" [-154.977839, 19.781396],\n",
" [-154.971306, 19.776015],\n",
" [-154.963843, 19.769189],\n",
" [-154.961237, 19.766308],\n",
" [-154.958728, 19.763534],\n",
" [-154.956172, 19.760794],\n",
" [-154.954406, 19.758902],\n",
" [-154.95243, 19.756784],\n",
" [-154.950419, 19.752354],\n",
" [-154.94935, 19.749999],\n",
" [-154.949144, 19.749999],\n",
" [-154.948527, 19.749999],\n",
" [-154.948322, 19.749999],\n",
" [-154.942429, 19.739708],\n",
" [-154.940619, 19.735726],\n",
" [-154.936963, 19.727678],\n",
" [-154.931853, 19.715468],\n",
" [-154.928109, 19.704534],\n",
" [-154.926134, 19.691664],\n",
" [-154.926439, 19.689644],\n",
" [-154.927778, 19.680812],\n",
" [-154.928935, 19.673605],\n",
" [-154.925456, 19.671858],\n",
" [-154.920257, 19.669249],\n",
" [-154.915723, 19.665535],\n",
" [-154.912712, 19.663069],\n",
" [-154.911451, 19.662036],\n",
" [-154.911183, 19.661725],\n",
" [-154.907187, 19.657091],\n",
" [-154.905856, 19.655547],\n",
" [-154.90306, 19.652305],\n",
" [-154.897096, 19.642793],\n",
" [-154.895128, 19.636764],\n",
" [-154.894852, 19.636249],\n",
" [-154.892225, 19.634474],\n",
" [-154.884291, 19.625],\n",
" [-154.885722, 19.625],\n",
" [-154.886784, 19.625],\n",
" [-154.883694, 19.621899],\n",
" [-154.883208, 19.621481],\n",
" [-154.881428, 19.619951],\n",
" [-154.878265, 19.617844],\n",
" [-154.875, 19.616167],\n",
" [-154.870122, 19.61415],\n",
" [-154.86221, 19.610534],\n",
" [-154.854472, 19.606491],\n",
" [-154.852924, 19.605775],\n",
" [-154.848783, 19.60386],\n",
" [-154.844961, 19.602135],\n",
" [-154.839912, 19.601705],\n",
" [-154.835024, 19.599571],\n",
" [-154.827827, 19.596538],\n",
" [-154.822221, 19.59331],\n",
" [-154.814908, 19.587047],\n",
" [-154.811505, 19.582594],\n",
" [-154.803328, 19.578725],\n",
" [-154.797548, 19.575668],\n",
" [-154.792204, 19.571755],\n",
" [-154.788991, 19.568743],\n",
" [-154.787735, 19.566807],\n",
" [-154.778294, 19.559895],\n",
" [-154.771682, 19.552514],\n",
" [-154.769345, 19.549916],\n",
" [-154.764892, 19.544965],\n",
" [-154.759333, 19.536462],\n",
" [-154.757137, 19.527907],\n",
" [-154.755838, 19.520784],\n",
" [-154.755827, 19.519444],\n",
" [-154.755793, 19.515524],\n",
" [-154.755792, 19.515342],\n",
" [-154.756114, 19.511596],\n",
" [-154.756605, 19.506911],\n",
" [-154.757166, 19.503855],\n",
" [-154.757874, 19.5],\n",
" [-154.75828, 19.5],\n",
" [-154.759499, 19.5],\n",
" [-154.759906, 19.5],\n",
" [-154.762169, 19.492619],\n",
" [-154.76249, 19.491575],\n",
" [-154.767902, 19.476301],\n",
" [-154.770121, 19.470875],\n",
" [-154.773044, 19.46373],\n",
" [-154.775751, 19.457113],\n",
" [-154.781089, 19.447616],\n",
" [-154.785537, 19.441229],\n",
" [-154.786132, 19.440376],\n",
" [-154.79357, 19.432572],\n",
" [-154.810311, 19.416057],\n",
" [-154.828271, 19.400891],\n",
" [-154.842815, 19.389236],\n",
" [-154.847913, 19.385152],\n",
" [-154.859673, 19.377633],\n",
" [-154.864313, 19.374998],\n",
" [-154.86645, 19.373637],\n",
" [-154.872863, 19.369556],\n",
" [-154.875001, 19.368196],\n",
" [-154.882661, 19.363797],\n",
" [-154.889028, 19.360392],\n",
" [-154.895905, 19.35232],\n",
" [-154.90823, 19.34202],\n",
" [-154.918905, 19.333484],\n",
" [-154.927296, 19.326959],\n",
" [-154.931719, 19.319875],\n",
" [-154.939663, 19.312834],\n",
" [-154.949939, 19.307096],\n",
" [-154.962283, 19.302062],\n",
" [-154.971455, 19.297839],\n",
" [-154.980079, 19.294118],\n",
" [-154.991348, 19.28822],\n",
" [-154.997442, 19.284554],\n",
" [-155.0, 19.282037],\n",
" [-155.0, 19.28438],\n",
" [-155.005449, 19.282647],\n",
" [-155.010823, 19.280572],\n",
" [-155.01799, 19.277749],\n",
" [-155.018508, 19.277527],\n",
" [-155.026142, 19.274253],\n",
" [-155.031727, 19.271769],\n",
" [-155.032324, 19.271504],\n",
" [-155.040383, 19.269032],\n",
" [-155.046917, 19.267482],\n",
" [-155.049153, 19.267251],\n",
" [-155.052291, 19.265236],\n",
" [-155.056149, 19.262117],\n",
" [-155.061439, 19.258416],\n",
" [-155.067531, 19.255323],\n",
" [-155.070657, 19.254498],\n",
" [-155.072812, 19.25393],\n",
" [-155.07639, 19.253628],\n",
" [-155.079339, 19.254088],\n",
" [-155.083189, 19.25],\n",
" [-155.080098, 19.25],\n",
" [-155.079031, 19.248966],\n",
" [-155.085359, 19.245466],\n",
" [-155.090387, 19.242796],\n",
" [-155.10551, 19.23427],\n",
" [-155.112571, 19.232142],\n",
" [-155.116306, 19.230602],\n",
" [-155.12001, 19.228532],\n",
" [-155.123329, 19.227696],\n",
" [-155.125, 19.226887],\n",
" [-155.129596, 19.226479],\n",
" [-155.132272, 19.225443],\n",
" [-155.140876, 19.222117],\n",
" [-155.149854, 19.219609],\n",
" [-155.161301, 19.218375],\n",
" [-155.169958, 19.213358],\n",
" [-155.170286, 19.213168],\n",
" [-155.175617, 19.211386],\n",
" [-155.182098, 19.209954],\n",
" [-155.187864, 19.209466],\n",
" [-155.196552, 19.20968],\n",
" [-155.204441, 19.210561],\n",
" [-155.213474, 19.212171],\n",
" [-155.225031, 19.212226],\n",
" [-155.230051, 19.212851],\n",
" [-155.236866, 19.215207],\n",
" [-155.242446, 19.216352],\n",
" [-155.245808, 19.217659],\n",
" [-155.250001, 19.220088],\n",
" [-155.250001, 19.219991],\n",
" [-155.250001, 19.219172],\n",
" [-155.258015, 19.222171],\n",
" [-155.261333, 19.220645],\n",
" [-155.265096, 19.21964],\n",
" [-155.268405, 19.219495],\n",
" [-155.27318, 19.213752],\n",
" [-155.286462, 19.205747],\n",
" [-155.297856, 19.199453],\n",
" [-155.29885, 19.197649],\n",
" [-155.301562, 19.192493],\n",
" [-155.307406, 19.187016],\n",
" [-155.312707, 19.182398],\n",
" [-155.320261, 19.174778],\n",
" [-155.325377, 19.17128],\n",
" [-155.330756, 19.168647],\n",
" [-155.33578, 19.16558],\n",
" [-155.338831, 19.163533],\n",
" [-155.342151, 19.161315],\n",
" [-155.345555, 19.160047],\n",
" [-155.348155, 19.158773],\n",
" [-155.356931, 19.155993],\n",
" [-155.365346, 19.153814],\n",
" [-155.365729, 19.153736],\n",
" [-155.370267, 19.152817],\n",
" [-155.371521, 19.152481],\n",
" [-155.372952, 19.152233],\n",
" [-155.374298, 19.151294],\n",
" [-155.374999, 19.150913],\n",
" [-155.375538, 19.150599],\n",
" [-155.378633, 19.148798],\n",
" [-155.380511, 19.148385],\n",
" [-155.384004, 19.146956],\n",
" [-155.386153, 19.146031],\n",
" [-155.387406, 19.145612],\n",
" [-155.389376, 19.144943],\n",
" [-155.391434, 19.14436],\n",
" [-155.393505, 19.1412],\n",
" [-155.397914, 19.135056],\n",
" [-155.399903, 19.13304],\n",
" [-155.401865, 19.131054],\n",
" [-155.404294, 19.127726],\n",
" [-155.406541, 19.125002],\n",
" [-155.406515, 19.124582],\n",
" [-155.411831, 19.120069],\n",
" [-155.41883, 19.114193],\n",
" [-155.422438, 19.1117],\n",
" [-155.428934, 19.107215],\n",
" [-155.43922, 19.102885],\n",
" [-155.451195, 19.099414],\n",
" [-155.460243, 19.097045],\n",
" [-155.466097, 19.095352],\n",
" [-155.46769, 19.093561],\n",
" [-155.47247, 19.088789],\n",
" [-155.477454, 19.084407],\n",
" [-155.477606, 19.084274],\n",
" [-155.484075, 19.08053],\n",
" [-155.494275, 19.072085],\n",
" [-155.494878, 19.071586],\n",
" [-155.495471, 19.07034],\n",
" [-155.495936, 19.069367],\n",
" [-155.497611, 19.065441],\n",
" [-155.49822, 19.064478],\n",
" [-155.498583, 19.063905],\n",
" [-155.500002, 19.064338],\n",
" [-155.500203, 19.063825],\n",
" [-155.50081, 19.062289],\n",
" [-155.501012, 19.061777],\n",
" [-155.501383, 19.060837],\n",
" [-155.501369, 19.056819],\n",
" [-155.503307, 19.046483],\n",
" [-155.506596, 19.037605],\n",
" [-155.509352, 19.030178],\n",
" [-155.513265, 19.020361],\n",
" [-155.513664, 19.018922],\n",
" [-155.514594, 19.015579],\n",
" [-155.516554, 19.011483],\n",
" [-155.518785, 19.008071],\n",
" [-155.523256, 19.003813],\n",
" [-155.526053, 19.0],\n",
" [-155.532261, 18.99355],\n",
" [-155.537366, 18.989057],\n",
" [-155.541787, 18.981294],\n",
" [-155.543224, 18.9741],\n",
" [-155.545351, 18.96185],\n",
" [-155.551511, 18.951789],\n",
" [-155.55747, 18.943442],\n",
" [-155.5656, 18.934517],\n",
" [-155.574511, 18.928175],\n",
" [-155.583779, 18.922094],\n",
" [-155.591057, 18.91333],\n",
" [-155.591551, 18.912736],\n",
" [-155.603944, 18.90128],\n",
" [-155.61606, 18.890249],\n",
" [-155.625, 18.883908],\n",
" [-155.625, 18.884345],\n",
" [-155.625, 18.885658],\n",
" [-155.625, 18.886096],\n",
" [-155.630753, 18.879687],\n",
" [-155.639097, 18.875],\n",
" [-155.641812, 18.875],\n",
" [-155.649041, 18.870503],\n",
" [-155.662995, 18.866159],\n",
" [-155.67836, 18.86546],\n",
" [-155.694257, 18.866786],\n",
" [-155.695032, 18.867019],\n",
" [-155.697795, 18.867847],\n",
" [-155.705973, 18.870301],\n",
" ...]],\n",
" u'type': u'Polygon'},\n",
" u'properties': {u'ALAND': 10433603617.0,\n",
" u'AWATER': 2739477738.0,\n",
" u'CBSAFP': u'25900',\n",
" u'CSAFP': None,\n",
" u'GEOID': u'25900',\n",
" u'INTPTLAT': u'+19.5977643',\n",
" u'INTPTLON': u'-155.5024434',\n",
" u'LSAD': u'M2',\n",
" u'MEMI': u'2',\n",
" u'MTFCC': u'G3110',\n",
" u'NAME': u'Hilo, HI',\n",
" u'NAMELSAD': u'Hilo, HI Micro Area'},\n",
" u'type': u'Feature'}"
]
},
"execution_count": 62,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"msas['features'][0]"
]
},
{
"cell_type": "code",
"execution_count": 63,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# MSA Name -> GeoJSON\n",
"geoJSON = {}\n",
"\n",
"for feature in msas['features']:\n",
" if feature['properties']['NAME'] in required_msa:\n",
" geoJSON[feature['properties']['NAME']] = {'type': 'FeatureCollection',\n",
" 'features': [feature]}"
]
},
{
"cell_type": "code",
"execution_count": 64,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[u'Evansville, IN-KY',\n",
" u'Cleveland-Elyria-Mentor, OH',\n",
" u'Tyler, TX',\n",
" u'Tulsa, OK',\n",
" u'Duluth, MN-WI',\n",
" u'Cincinnati-Middletown, OH-KY-IN',\n",
" u'Allentown-Bethlehem-Easton, PA-NJ',\n",
" u'Little Rock-North Little Rock-Conway, AR',\n",
" u'Chattanooga, TN-GA',\n",
" u'Appleton, WI',\n",
" u'Sioux Falls, SD',\n",
" u'Louisville/Jefferson County, KY-IN',\n",
" u'Knoxville, TN',\n",
" u'Amarillo, TX',\n",
" u'Tallahassee, FL',\n",
" u'South Bend-Mishawaka, IN-MI',\n",
" u'Waco, TX',\n",
" u'St. Louis, MO-IL',\n",
" u'Cape Coral-Fort Myers, FL',\n",
" u'Macon, GA',\n",
" u'El Paso, TX',\n",
" u'St. Joseph, MO-KS',\n",
" u'Reading, PA',\n",
" u'Monroe, LA',\n",
" u'Portland-South Portland-Biddeford, ME',\n",
" u'Washington-Arlington-Alexandria, DC-VA-MD-WV',\n",
" u'Provo-Orem, UT',\n",
" u'Bellingham, WA',\n",
" u'Utica-Rome, NY',\n",
" u'Milwaukee-Waukesha-West Allis, WI',\n",
" u'Baltimore-Towson, MD',\n",
" u'Lansing-East Lansing, MI',\n",
" u'Fargo, ND-MN',\n",
" u'Laredo, TX',\n",
" u'Detroit-Warren-Livonia, MI',\n",
" u'Las Vegas-Paradise, NV',\n",
" u'Medford, OR',\n",
" u'Worcester, MA',\n",
" u'Santa Cruz-Watsonville, CA',\n",
" u'Visalia-Porterville, CA',\n",
" u'Kokomo, IN',\n",
" u'Virginia Beach-Norfolk-Newport News, VA-NC',\n",
" u'Scranton--Wilkes-Barre, PA',\n",
" u'Oklahoma City, OK',\n",
" u'Greeley, CO',\n",
" u'Johnstown, PA',\n",
" u'Syracuse, NY',\n",
" u'Muncie, IN',\n",
" u'Stockton, CA',\n",
" u'Buffalo-Niagara Falls, NY',\n",
" u'Gadsden, AL',\n",
" u'Fayetteville, NC',\n",
" u'Corpus Christi, TX',\n",
" u'Hattiesburg, MS',\n",
" u'Flint, MI',\n",
" u'Redding, CA',\n",
" u'Alexandria, LA',\n",
" u'Eau Claire, WI',\n",
" u'Savannah, GA',\n",
" u'Bridgeport-Stamford-Norwalk, CT',\n",
" u'Greenville-Mauldin-Easley, SC',\n",
" u'San Luis Obispo-Paso Robles, CA',\n",
" u'Gulfport-Biloxi, MS',\n",
" u'Abilene, TX',\n",
" u'Dallas-Fort Worth-Arlington, TX',\n",
" u'Boston-Cambridge-Quincy, MA-NH',\n",
" u'Miami-Fort Lauderdale-Pompano Beach, FL',\n",
" u'Janesville, WI',\n",
" u'Flagstaff, AZ',\n",
" u'Lexington-Fayette, KY',\n",
" u'Raleigh-Cary, NC',\n",
" u'Birmingham-Hoover, AL',\n",
" u'Terre Haute, IN',\n",
" u'Shreveport-Bossier City, LA',\n",
" u'Beaumont-Port Arthur, TX',\n",
" u'San Francisco-Oakland-Fremont, CA',\n",
" u'Waterloo-Cedar Falls, IA',\n",
" u'New Orleans-Metairie-Kenner, LA',\n",
" u'Toledo, OH',\n",
" u'Riverside-San Bernardino-Ontario, CA',\n",
" u'Minneapolis-St. Paul-Bloomington, MN-WI',\n",
" u'Pittsburgh, PA',\n",
" u'Harrisburg-Carlisle, PA',\n",
" u'Boise City-Nampa, ID',\n",
" u'Wausau, WI',\n",
" u'Wichita, KS',\n",
" u'Olympia, WA',\n",
" u'Altoona, PA',\n",
" u'Lincoln, NE',\n",
" u'Asheville, NC',\n",
" u'Canton-Massillon, OH',\n",
" u'Pueblo, CO',\n",
" u'Hagerstown-Martinsburg, MD-WV',\n",
" u'Philadelphia-Camden-Wilmington, PA-NJ-DE-MD',\n",
" u'Columbus, OH',\n",
" u'Seattle-Tacoma-Bellevue, WA',\n",
" u'Charlottesville, VA',\n",
" u'San Diego-Carlsbad-San Marcos, CA',\n",
" u'Clarksville, TN-KY',\n",
" u'Indianapolis-Carmel, IN',\n",
" u'Palm Bay-Melbourne-Titusville, FL',\n",
" u'Dayton, OH',\n",
" u'Omaha-Council Bluffs, NE-IA',\n",
" u'Tucson, AZ',\n",
" u'McAllen-Edinburg-Mission, TX',\n",
" u'Madison, WI',\n",
" u'Augusta-Richmond County, GA-SC',\n",
" u'Springfield, MO',\n",
" u'Youngstown-Warren-Boardman, OH-PA',\n",
" u'Lubbock, TX',\n",
" u'Anniston-Oxford, AL',\n",
" u'Rocky Mount, NC',\n",
" u'Eugene-Springfield, OR',\n",
" u'Richmond, VA',\n",
" u'Kansas City, MO-KS',\n",
" u'Reno-Sparks, NV',\n",
" u'Hickory-Lenoir-Morganton, NC',\n",
" u'Salem, OR',\n",
" u'Salinas, CA',\n",
" u'Roanoke, VA',\n",
" u'Modesto, CA',\n",
" u'Yuma, AZ',\n",
" u'Binghamton, NY',\n",
" u'Salt Lake City, UT',\n",
" u'Houma-Bayou Cane-Thibodaux, LA',\n",
" u'Anchorage, AK',\n",
" u'Providence-New Bedford-Fall River, RI-MA',\n",
" u'New Haven-Milford, CT',\n",
" u'Oxnard-Thousand Oaks-Ventura, CA',\n",
" u'Lancaster, PA',\n",
" u'Brownsville-Harlingen, TX',\n",
" u'Akron, OH',\n",
" u'Houston-Sugar Land-Baytown, TX',\n",
" u'Tampa-St. Petersburg-Clearwater, FL',\n",
" u'Topeka, KS',\n",
" u'Fort Collins-Loveland, CO',\n",
" u'New York-Northern New Jersey-Long Island, NY-NJ-PA',\n",
" u'Albuquerque, NM',\n",
" u'Greensboro-High Point, NC',\n",
" u'Killeen-Temple-Fort Hood, TX',\n",
" u'Williamsport, PA',\n",
" u'Davenport-Moline-Rock Island, IA-IL',\n",
" u'Jacksonville, FL',\n",
" u'San Jose-Sunnyvale-Santa Clara, CA',\n",
" u'Pensacola-Ferry Pass-Brent, FL',\n",
" u'Albany-Schenectady-Troy, NY',\n",
" u'Chico, CA',\n",
" u'Joplin, MO',\n",
" u'Cedar Rapids, IA',\n",
" u'Merced, CA']"
]
},
"execution_count": 64,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"geoJSON.keys()"
]
},
{
"cell_type": "code",
"execution_count": 65,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"with open('msa.geojson', 'wb') as outfile:\n",
" json.dump(geoJSON, outfile)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.10"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment