We present the theory of inversion with two-dimensional regularization. We use this novel method to retrieve profiles of microphysical properties of atmospheric particles from profiles of optical properties acquired with multiwavelength Raman lidar. This technique is the first attempt to the best of our knowledge, toward an operational inversion algorithm, which is strongly needed in view of multiwavelength Raman lidar networks. The new algorithm has several advantages over the inversion with so-called classical one-dimensional regularization. Extensive data postprocessing procedures, which are needed to obtain a sensible physical solution space with the classical approach, are reduced. Data analysis, which strongly depends on the experience of the operator, is put on a more objective basis. Thus, we strongly increase unsupervised data analysis. First results from simulation studies show that the new methodology in many cases outperforms our old methodology regarding accuracy of retrieved particle effective radius, and number, surface-area, and volume concentration. The real and the imaginary parts of the complex refractive index can be estimated with at least as equal accuracy as with our old method of inversion with one-dimensional regularization. However, our results on retrieval accuracy still have to be verified in a much larger simulation study.