Пробовали ли вы предоставить большое скалярное число аргументу zorder в функции kdeplot? Это должно выглядеть так:
zorder
kdeplot
sns.kdeplot(df_toplot['longitude'], df_toplot['latitude'], \ shade=True, cmap='Purples', \ ax=ax, zorder=100);