Recently, I moved from R to python. Basically, I can quickly find a good alternative to R function in python except smoothScatter. I searched a lot and found these two answers helpful: 
- Generate a heatmap in MatPlotLib using a scatter data set
- how does 2d kernel density estimation in python (sklearn) work?
But both of them are relatively complicated and slow compared to R (I have about 2,000,000 data points and scipy.stats.gaussian_kde is really slow). So is there a simple python package can replace smoothScatter in R?

 
     
    
