We consider the coupling between massive and spinning particles and three dimensional gravity. This allows us to construct geometric operators (distances between particles) as Dirac observables. We quantize the system à la loop quantum gravity: we give a description of the kinematical Hilbert space and construct the associated spin-foam model. We quantize the physical distance operator and compute its spectrum.