import psycopg2
import shapely
import osgeo.ogr
import shapely.wkt
import geopandas as gpd
import matplotlib.pyplot as plt

try:
    conn = psycopg2.connect("dbname='strokes' user='postgres'    host='localhost' password='****'")
except:
print "I am unable to connect to the database"
cur = conn.cursor()
cur.execute("CREATE INDEX bassin_index ON bassin USING GIST(geom)")
connection.commit()
cur.execute("SELECT st_astext(geom) AS wkt, fid_limite, codebassin FROM bassin")
rows = cur.fetchall()

rows_list=[]
for geom,fid_limite,codebassin in cursor:
    data={'codebassin':codeb,'fidb':fidlim,'geom':shapely.wkt.loads(geom)}
    rows_list.append(data)
    gdf=gpd.GeoDataFrame(rows_list).set_index('codebassin')
    gdf.head()
    gdf.plot(column='rows_list', scheme='QUANTILES', k=5, colormap='gray')
    plt.show()

conn.commit()
conn.close()