A novel method is developed in this study for soil classification and zonation in a two-dimensional (2D) vertical cross section using cone penetration tests (CPTs). A CPT is usually performed vertically and the number of CPT soundings in a site is often limited in geotechnical engineering practice. It is, therefore, difficult to properly interpret CPT results along the horizontal direction or accurately estimate the horizontal correlation length of CPT data. The method proposed in this study bypasses the difficulty in estimating horizontal correlation length and provides proper identification of subsurface soil stratification (i.e., soil layer number is constant along horizontal direction) and zonation (i.e., soil layer number varies along horizontal direction) in a 2D vertical cross section directly from a limited number of CPT soundings. The proposed method consists of three key elements: 2D interpolation of CPT data using 2D Bayesian compressive sampling; determination of soil behavior type (SBT) using a SBT chart at every location in the 2D section, including locations with measurements and unsampled locations; and soil layer or zone delineation using an edge detection method. Both simulated and real data examples are used to illustrate the proposed method. Results show that the method performs well even when only five sets of CPT soundings are available.