# coding: utf-8
import geopandas as gpd
import matplotlib.pyplot as plt

#le fichier du cadastre pour example
# a télécharger à l'avance et à mettre qq part
#https://cadastre.data.gouv.fr/data/etalab-cadastre/2018-10-01/geojson/communes/92/92023/cadastre-92023-parcelles.json.gz


dfgeo = gpd.read_file('fichier.json')  #le fchier précédement telecharger

# trier par surface
dfgeo.sort_values(by='contenance',inplace=True)

surf_min = 1100  #surface minimum à chercher
surf_max = 1100 #surface max à chercher

f, ax = plt.subplots(1, figsize=(20, 20))
ax.set_aspect('equal')
dfgeo.geometry.plot(ax=ax,alpha=.7)

selection = dfgeo[dfgeo.contenance.between(surf_min,surf_max)]['geometry']
selection.centroid.plot(ax=ax,color='red')

cpt=0
for index,row in dfgeo[dfgeo.contenance.between(surf_min,surf_max)].iterrows():
    cpt+=1
    lat = (row['geometry'].centroid.x)
    lon = (row['geometry'].centroid.y)
    ax.annotate(cpt,(lat,lon),size=20)
    surf = row['contenance']
    url = 'https://www.google.fr/maps/place/'+str(lon)+','+str(lat)
    print (cpt, url , surf )



plt.show()