The role of intraspecific trait variation in functional ecology has gained traction in recent years as many papers have observed its importance in driving community diversity and ecology. Yet much of the work in this field relies on field-based trait surveys. Here, we used continuous canopy trait information derived from remote sensing data of a highly polymorphic tree species, Metrosideros polymorpha, to quantify environmental controls on intraspecific trait variation. M. polymorpha, an endemic, keystone tree species in Hawai’i, varies morphologically, chemically, and genetically across broad elevation and soil substrate age gradients, making it an ideal model organism to explore large-scale environmental drivers of intraspecific trait variation. M. polymorpha canopy reflectance (visible to shortwave infrared; 380–2510 nm) and light detection and ranging (LiDAR) data collected by the Global Airborne Observatory were modeled to canopy trait estimates of leaf mass per area, chlorophyll a and b, carotenoids, total carbon, nitrogen, phosphorus, phenols, cellulose, and top of canopy height using previously developed leaf chemometric equations. We explored how these derived traits varied across environmental gradients by extracting elevation, slope, aspect, precipitation, and soil substrate age data at canopy locations. We then obtained the feature importance values of the environmental factors in predicting each leaf trait by training random forest models to predict leaf traits individually. Of these environmental factors, elevation was the most important predictor for all canopy traits. Elevation not only affected canopy traits directly but also indirectly by influencing the relationships between soil substrate age and canopy traits as well as between nitrogen and other traits, as indicated by the change in slope between the variables at different elevation ranges. In conclusion, intraspecific variation in M. polymorpha traits derived from remote sensing adheres to known leaf economic spectrum (LES) patterns as well as interspecific LES traits previously mapped using imaging spectroscopy.