# Imports
#
import matplotlib.pyplot as plt
from libpysal.weights.contiguity import Queen
from libpysal import examples
import geopandas as gpd
from esda.moran import Moran_Local
from splot.esda import plot_local_autocorrelation
#
# Data preparation and analysis
#
guerry = examples.load_example('Guerry')
link_to_data = guerry.get_path('guerry.shp')
gdf = gpd.read_file(link_to_data)
y = gdf['Donatns'].values
w = Queen.from_dataframe(gdf)
w.transform = 'r'
moran_loc = Moran_Local(y, w)
#
# Plotting with quadrant mask and region mask
#
fig = plot_local_autocorrelation(moran_loc, gdf, 'Donatns', p=0.05,
                                 region_column='Dprtmnt',
                                 mask=['Ain'], quadrant=1)
plt.show()