The increasing number of satellite missions provides vast opportunities for continuous vegetation monitoring, crucial for precision agriculture and environmental sustainability. However, accurately estimating vegetation traits, such as nitrogen concentration (N%), from Landsat 7 (L7), Landsat 8 (L8), and Sentinel-2 (S2) satellite data is challenging due to the diverse sensor configurations and complex atmospheric interactions. To address these limitations, we developed a unified and physically based method that combines a soil–plant–atmosphere radiative transfer (SPART) model with the bottom-of-atmosphere (BOA) spectral bidirectional reflectance distribution function. This approach enables us to assess the effect of rugged terrain, viewing angles, and illumination geometry on the spectral reflectance of multiple sensors. Our methodology involves inverting radiative transfer model variables using numerical optimization to estimate N% and creating a hybrid model. We used Gaussian process regression (GPR) to incorporate the inverted variables into the hybrid model for N% prediction, resulting in a unified approach for N% estimation across different sensors. Our model shows a validation accuracy of 0.35 (RMSE %N), a mean prediction interval width (MPIW) of 0.35, and an R${}^2$ of 0.50, using independent data from multiple sensors collected between 2016 and 2019. Our unified method provides a promising solution for estimating N% in vegetation from L7, L8, and S2 satellite data, overcoming the limitations posed by diverse sensor configurations and complex atmospheric interactions.