Leaf photosynthetic pigments play a crucial role in evaluating nutritional elements and physiological states. In facility agriculture, it is vital to rapidly and accurately obtain the pigment content and distribution of leaves to ensure precise water and fertilizer management. In our research, we utilized chlorophyll a (Chla), chlorophyll b (Chlb), total chlorophylls (Chls) and total carotenoids (Cars) as indicators to study the variations in the leaf positions of Lycopersicon esculentum Mill. Under 10 nitrogen concentration applications, a total of 2610 leaves (435 samples) were collected using visible-near infrared hyperspectral imaging (VNIR–HSI). In this study, a “coarse–fine” screening strategy was proposed using competitive adaptive reweighted sampling (CARS) and the iteratively retained informative variable (IRIV) algorithm to extract the characteristic wavelengths. Finally, simultaneous and quantitative models were established using partial least squares regression (PLSR). The CARS–IRIV–PLSR was used to create models to achieve a better prediction effect. The coefficient determination (R2), root mean square error (RMSE) and ratio performance deviation (RPD) were predicted to be 0.8240, 1.43 and 2.38 for Chla; 0.8391, 0.53 and 2.49 for Chlb; 0.7899, 2.24 and 2.18 for Chls; and 0.7577, 0.27 and 2.03 for Cars, respectively. The combination of these models with the pseudo-color image allowed for a visual inversion of the content and distribution of the pigment. These findings have important implications for guiding pigment distribution, nutrient diagnosis and fertilization decisions in plant growth management.