Understanding the determinants of urban forest diversity and structure is important for preserving biodiversity and sustaining ecosystem services in cities. However, comprehensive field assessments are resource-intensive, and landscape-level approaches may overlook heterogeneity within urban regions. To address this challenge, we combined remote sensing with field inventories to comprehensively map and analyze urban forest attributes in forest patches across the Minneapolis-St. Paul Metropolitan Area (MSPMA) in a multistep process. First, we developed predictive machine learning models of forest attributes by integrating data from forest inventories (from 40 12.5-m-radius plots) with Global Ecosystem Dynamics Investigation (GEDI) observations and Sentinel-2-derived land surface phenology (LSP). These models enabled accurate predictions of forest attributes, specifically nine metrics of plant diversity (tree species richness, tree abundance, and understory plant abundance), structure (average canopy height, dbh, and canopy density), and structural complexity (variability in canopy height, dbh, and canopy density) with relative errors ranging between 11% and 21%. Second, we applied these machine learning models to predict diversity metrics for 804 additional plots from GEDI and Sentinel-2. Finally, we applied Bayesian multilevel models to the predicted diversity metrics to assess the influence of multiple factors-patch dimensions, landscape attributes, plot position, and jurisdictional agency-on these forest attributes across the 804 predicted plots. The models showed all predictors have some degree of effect on forest attributes, presenting varying explanatory power with R