from scipy.stats import gaussian_kde
import matplotlib.pyplot as plt
x,y=buoywinds, era5winds
xy = np.vstack([x,y])
z = gaussian_kde(xy)(xy)
idx = z.argsort()
plt.scatter(buoywinds, era5winds, c = z, s = 10, cmap = "Spectral")
plt.colorbar()
x_line = np.linspace(min(min(x), min(y)), max(max(x), max(y)), 100)
plt.plot(x_line, x_line, color='black', linestyle='--')
plt.xlabel('Buoy Wind Speed (m/s)')
plt.ylabel('DS Wind Speed (m/s)')
plt.title('Comparison of Wind Speeds')
plt.show()