Leaf photosynthetic pigments play a crucial role in evaluating nutritional elements and physiological states. In facility agriculture, it is vital to obtain rapidly and accurately 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), chlorophyll (Chll), and carotenoid (Caro) as indicators to study the variations in leaf position 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 by using competitive adaptive reweighted sampling (CARS) and iteratively retained informative variable (IRIV) algorithm to extract 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, 2.38 for Chla, 0.8391, 0.53, 2.49 for Chlb, 0.7899, 2.24, 2.18 for Chll, and 0.7577, 0.27, 2.03 for Caro, respectively. The combination of these models with the pseudo-color image allowed for a visual inversion of the content and distribution of pigment. These findings have important implications for guiding pigment distribution, nutrient diagnosis, and fertilization decisions in plant growth management.