import matplotlib.pyplot as plt
import seaborn as sns
from scipy.stats import pearsonr
plt.figure(figsize=(10,6))
sns.scatterplot(x=data_merge_max_median['air_density_max'], y=data_merge_max_median['air_density_median'], s=6, color='royalblue')# , hue="royalblue"
corr, _ = pearsonr(data_merge_max_median['air_density_max'], data_merge_max_median['air_density_median'])
plt.title(f'Scatter plot of ws_max Air Density vs ws_median Air Density\nCorrelation: {corr:.2f}')
plt.xlabel('ws_Max Air Density(kg/m^3)')
plt.ylabel('ws_Median Air Density(kg/m^3)')
plt.tight_layout()# plt.legend()
plt.show()