Accurately estimating leaf chlorophyll content (Chl) at large spatial scales is crucial for monitoring agricultural production and plant photosynthesis. Sentinel-2 and Landsat-7/8 offer the potential to assess Chl with high spatial resolutions using various physically-based, empirical, and hybrid methods, but they still present challenges due to the confounding effects of canopy structure and soil background. The near-infrared reflectance of vegetation (NIRV) has been proposed to characterize canopy structural and soil background effects in sun-induced chlorophyll fluorescence. Therefore, in this study, we developed a novel strategy that incorporates NIRV to minimize canopy structural and soil background effects on Chl retrieval, using bidirectional reflectance factor (BRF) from Sentinel-2 and Landsat-7/8 imagery. We compared Chl estimation using empirical random forest (RF), the physically-based PROSAIL-5B inversion, and a hybrid method using RF to train simulated datasets from PROSAIL-5B. We validated Chl retrieval performance with measured datasets from the National Ecological Observatory Network, covering seven vegetation types. We also compared spatiotemporal Chl variations from Sentinel-2, Landsat-7/8, and MODIS. The results showed that NIRV effectively characterized the canopy structural effect on Chl estimation, associated with leaf area index, average leaf angle, and canopy fraction of senescent leaves. Compared to BRF, correcting BRF by NIRV improved Chl retrieval with RMSE values reduced by 12.2%–27.0% for Sentinel-2, 7.4%–16.4% for Landsat-7, and 6.9%–13.7% for Landsat-8, respectively. This also resulted in a significant increase in R2 values for Sentinel-2 (0.13–0.46), Landsat-7 (0.08–0.19), and Landsat-8 (0.12–0.2), respectively. Sentinel-2 outperformed Landsat-7/8 in Chl estimation, due to the contribution of red edge bands and a higher spatial resolution. Furthermore, the hybrid method performed best for Chl estimation from Sentinel-2 and Landsat-7/8 with relative RMSE values of 17.33%, 20.22%, and 19.23%, respectively. Our hybrid method demonstrated lower model uncertainty with smaller variations of RMSE values across varying biochemical and structural traits and solar and viewing angles, and exhibited more consistent spatiotemporal variations of Chl from Sentinel-2 with phenology compared to Landsat-7/8 and the published MODIS Chl data. Our hybrid method demonstrates significant potential in mitigating canopy structural effects, enhancing Chl retrieval, and facilitating the development of Chl data from various satellite imagery sources.