Estimating Orientation Distribution Functions with Probability Density Constraints and Spatial Regularity