#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Tue Nov 12 12:45:24 2019

@author: mitchellteal
"""
import geopandas as gpd
import matplotlib
import numpy as np
import matplotlib.pyplot as plt
##%matplotlib inline
import shapefile as shp
#
fp = "Cable.shx"
data=gpd.read_file(fp)
data.crs
data.head()
data.plot()
sf = shp.Reader(fp)

plt.figure()
for shape in sf.shapeRecords():
    x = [i[0] for i in shape.shape.points[:]]
    y = [i[1] for i in shape.shape.points[:]]
    plt.plot(x,y)
plt.show()

import pandas as pd
import seaborn as sns
from mpl_toolkits.basemap import Basemap

fig, ax = plt.subplots(figsize=(10,20))
m = Basemap(resolution='l', # c, l, i, h, f or None
            projection='merc',
            lat_0=52.5, lon_0=-1,
            llcrnrlon=-7.5,llcrnrlat=49.75,urcrnrlon=3,urcrnrlat=58, epsg=27700)

m.drawmapboundary(fill_color='#bee3ee')
m.fillcontinents(color='#ffffff',lake_color='aqua')
m.drawcoastlines()
for shape in sf.shapeRecords():
    x = [i[0] for i in shape.shape.points[:]]
    y = [i[1] for i in shape.shape.points[:]]
    m.plot(x,y)
plt.show()