Downscaling of remote sensing precipitation products and the forecasting of circulation model are always the intense interests in hydrology and meteorology. The essence of downscaling is primarily to enhance resolution of observation or simulated rainfall field, and to appropriately increase its details or high frequency characteristics. Precipitation, as the main driving factors of the earth’s hydrologic cycle, not only affects the moisture and heat condition of a certain river basin, but also affects the global water and heat circulation. Based on the properties of rainfall self-similarity structure, the mathematically ill-posed precipitation problem solving method was used in low resolution downscaling precipitation for high resolution reconstruction. When solving the downscaling ill-posed problem, the greedy method of orthogonal matching pursuit was introduced so as to get the best high-resolution estimation in an optimal sense. It is hard to imagine that we might be able to find very similar (in mathematical norms) precipitation patterns over relatively large storm-scales. However, finding similar features over sufficiently small sub-storm scales seems more feasible. Based on the characteristics that small scale organized precipitation features tend to recur across different storm environments, the precipitation of both high and low resolution was obtained by training, which could be used to reconstruct the desired high-resolution precipitation field. Multi-source merged precipitation products were used in this experiment. Given the consideration of incompleteness of merged precipitation dataset, it was firstly interpolated based on the method of Fields of Experts (FoEs), which could solve the problem that common interpolation method could hardly work on the interpolation for dataset where consecutive missing data exists. Secondly, ideal experiments of precipitation products downscaling were carried out, where smooth coupling sampling and resampling operator were adopted respectively. Assessment based on the metrics of fidelity and spatial structural similarity demonstrates that the method used in this paper is feasible.