The electronic density of states of graphene is equivalent to that of relativistic electrons [1-3]. In the absence of disorder or external doping the Fermi energy lies at the Dirac point where the density of states vanishes. Although transport measurements at high carrier densities indicate rather high mobilities [4-6], many questions pertaining to disorder [7-14] remain unanswered. In particular, it has been argued theoretically, that when the average carrier density is zero, the inescapable presence of disorder will lead to electron and hole puddles with equal probability. In this work, we use a scanning single electron transistor to image the carrier density landscape of graphene in the vicinity of the neutrality point. Our results clearly show the electron-hole puddles expected theoretically [13]. In addition, our measurement technique enables to determine locally the density of states in graphene. In contrast to previously studied massive two dimensional electron systems, the kinetic contribution to the density of states accounts quantitatively for the measured signal. Our results suggests that exchange and correlation effects are either weak or have canceling contributions. The kinetic energy of Dirac particles in graphene increases linearly with momentum [1,15,16]. The total energy per particle however, also referred to as the chemical potential, µ , contains additional exchange and correlations contributions that arise from the Coulomb interactionHere K E , ex E and c E are the kinetic, exchange and correlation terms respectively. Therefore, the chemical potential and its derivative with respect to density, known as the inverse compressibility or density of states, provide direct insight into the properties of the Coulomb interaction in such a system. Compressibility measurements of conventional, massive, two dimensional electron systems made for example of Si or GaAs have been carried out by several groups [19,[23][24][25][26]. It has been shown that at zero magnetic field the chemical potential may be described rather accurately within the Hartree-Fock approximation. Quantitatively it has been found that Coulomb interactions add a substantial contribution to the compressibility and become dominant at low carrier densities. In a perfectly uniform and clean graphene sample, the inverse compressibility is expected to diverge at the Dirac point in view of the vanishing density of states. This divergence, however, is expected to be rounded off by disorder on length scales smaller than our tip size. Long range disorder on the other hand will cause local shifts in the Dirac point indicative of a non-zero local density. In this work we measure the spatial dependence of the local compressibility versus carrier density across the sample. Fluctuations in the Dirac point across the sample are translated into