Preferential concentration of inertial particles in turbulence is studied numerically by evaluating the Lagrangian compressibility of the particle velocity field using the "full Lagrangian method." This is compared with the "mesoscopic Eulerian particle velocity field" both in a direct numerical simulation of turbulence and in a synthetic flow field. We demonstrate that the Lagrangian method, in contrast to the Eulerian, accurately predicts the compressibility of the particle velocity field even when the latter is characterized by singularities. In particular we use the method to evaluate the growth rates of spatial moments of the particle number density which reflect the fractal structure of segregation and the occurrence of singularities.