Euclidean Random Matching in 2D for Non-constant Densities