mport pandas as pd
import numpy as np
import geopandas as gp
import matplotlib.pyplot as plt
import shapely
import fiona
from mpl_toolkits.axes_grid1 import make_axes_locatable
import scipy


#china_geod0 = gp.GeoDataFrame.from_file('D:\data\china\gadm36_CHN_0.shp', encoding = 'gb18030')
china_geod1 = gp.GeoDataFrame.from_file('D:\data\china city344')
#china_geod2 = gp.GeoDataFrame.from_file('D:\data\china\gadm36_CHN_1.shp', encoding = 'gb18030')
#china_geod3 = gp.GeoDataFrame.from_file('D:\data\china\gadm36_CHN_1.shp', encoding = 'gb18030')
china_geod1.head(6)
print(china_geod1.head(6))
fig, ax = plt.subplots(1, 1)
#divider = make_axes_locatable(ax)
#cax = divider.append_axes("bottom", size="5%", pad=0.1)

#china_geod1.boundary.plot()

#china_geod0.plot()
china_geod1.plot(column='NAME_1', ax=ax, legend=True)
#china_geod1.plot(column='NAME_1',
           #cmap='OrRd',
           #ax=ax,
           #legend=True,
           #legend_kwds={'label': "Population by Country",
                       # 'orientation': "horizontal"}
                #)
#china_geod2.plot()
#china_geod3.plot()

plt.show()