The epidermal leaf patterns of plants exhibit remarkable diversity in cell shapes, sizes, and arrangements, driven by environmental interactions that lead to significant adaptive changes even among closely related species. The Solanaceae family, known for its high diversity of adaptive epidermal structures, has traditionally been studied using qualitative phenotypic descriptions. To advance this, we developed a workflow combining multi-scale computer vision, image processing, and data analysis to extract digital descriptors for leaf epidermal cell morphology. Applied to nine wild potato species, this workflow quantified key morphological parameters, identifying descriptors for trichomes, stomata, and pavement cells, and revealing interdependencies among these traits. Principal component analysis (PCA) highlighted two main axes, accounting for 45% and 21% of variance, corresponding to features such as guard cell shape, trichome length, stomatal density, and trichome density. These axes aligned well with the historical and geographical origins of the species, separating southern from Central American species, and forming distinct clusters for monophyletic groups. This workflow thus establishes a quantitative foundation for investigating leaf epidermal cell morphology within phylogenetic and geographic contexts.