Spatialization of biomass and carbon stocks is essential for a good understanding of the forest stand and its characteristics, especially in degraded Mediterranean cork oak forests. Furthermore, the analysis of biomass and carbon stock changes and dynamics is essential for understanding the carbon cycle, in particular carbon emissions and stocks, in order to make projections, especially in the context of climate change. In this research, we use a multidimensional framework integrating forest survey data, LiDAR UAV data, and extracted vegetation indices from Landsat imagery (NDVI, ARVI, CIG, etc.) to model and spatialize cork oak biomass and carbon stocks on a large scale. For this purpose, we explore the use of univariate and multivariate regression modeling and examine several types of regression, namely, multiple linear regression, stepwise linear regression, random forest regression, simple linear regression, logarithmic regression, and quadratic and cubic regression. The results show that for multivariate regression, stepwise regression gives good results, with R2 equal to 80% and 65% and RMSE equal to 2.59 and 1.52 Mg/ha for biomass and carbon stock, respectively. Random forest regression, chosen as the ML algorithm, gives acceptable results, explaining 80% and 60% of the variation in biomass and carbon stock, respectively, and an RMSE of 2.74 and 1.72 Mg/ha for biomass and carbon stock, respectively. For the univariate regression, the simple linear regression is chosen because it gives satisfactory results, close to those of the quadratic and cubic regressions, but with a simpler equation. The vegetation index chosen is ARVI, which shows good performance indices, close to those of the NDVI and CIG. The assessment of biomass and carbon stock changes in the study area over 35 years (1985–2020) showed a slight increase of less than 10 Mg/ha and a decrease in biomass and carbon stock over a large area.