我是seaborn的新手,请原谅我这是一个愚蠢的问题......我打算用seaborn绘制我的空间数据的多个方面(作为栅格数据)。源数据来自netcdf文件...导入和一些清理后,我有一个以下列布局的pandas数据框:
import seaborn as sns
import cartopy
import cartopy.crs as ccrs
import xray
ds1 = xray.open_dataset("example1.nc")['value']
ds2 = xray.open_dataset("example2.nc")['value']
df1 = ds1.to_dataframe()
df2 = ds2.to_dataframe()
# make lat, lon regular cols
df1.reset_index(inplace=True)
df2.reset_index(inplace=True)
# add a year column
df1['year'] = 2010
df2['year'] = 2011
# produce the final dataframe
df = pd.concat( [df1, df2] )
g = sns.FacetGrid(df, col="year", hue="value",
subplot_kws=dict(projection=ccrs.PlateCarree()), size=4.5,
sharex=False, sharey=False, despine=False)
g.map(plt.scatter, "lon", "lat", s=5, transform=ccrs.PlateCarree())
for ax in g.axes.ravel():
ax.add_feature(cartopy.feature.COASTLINE)
如果我使用散点图(我使用cartopy进行地理变换绘图),我可以得到一些视觉输出:
{{1}}
然而,当我想要绘制栅格数据时我迷失了...我想我应该使用pandas.meshgrid,plot.contourf和其他一些魔法,但我不能为我的生活找出我如何从大熊猫面对光栅图......