title: Scipy: how to obtain bivariate empirical density probability distribution (not by KDE)?
body: <p>My data are stored discretely, e.g. 1, 2, 3 ....</p>
  
  <p>For 1 dimension, the empirical density probability function would look like a histogram.</p>
  
  <p><a href="https://i.stack.imgur.com/cSb1D.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/cSb1D.png" alt="enter image description here"></a></p>
  
  <p>And for 2 dimensions, it will look like a bivairate histogram, and how can I obtain it? (I can obtain one by Kernel Density Estimation, but that result is often smoothed by the gaussian kernel. I want a un-smoothed pdf in this case.)</p>