Gel electrophoresis serves as a basic analytical tool in the proteomic studies. However, processing of gel electrophoretic images is still the main bottleneck of data analysis, and there is an increasing need for the fully automated approaches. The proposed start-to-end strategy of analyzing the gel images consists of chemometric tools, which allow their effective preprocessing, automatic warping, and data modeling. The image preprocessing techniques: denoising in the wavelet domain and the penalized asymmetric least squares approach for the background estimation are proposed. Matching of images is based on fuzzy warping of features, extracted from the gel images. For the classification or calibration purpose, multivariate approaches such, as partial least squares (PLS) or kernel-PLS methods are used. Performance of the proposed strategy is demonstrated on the real set of the two-dimensional gel images.