Increases in extreme weather events associated with climate change have the potential to put currently healthy forests at risk. One option to minimize this risk is the application of forest management measures aimed at generating species mixtures predicted to be more resilient to these threats. In order to apply such measures appropriately, forest managers need up-to-date, accurate and consistent forest maps at relatively fine spatial resolutions. Cost efficiency is a major factor when creating such maps. Taking European spruce (Picea abies) and Scots pine (Pinus sylvestris) as an example, this paper describes an innovative approach for mapping two tree species using a combination of commercial very high resolution WorldView-2 (WV2) images and Landsat time series data. As a first step, this study used a supervised object-based classification of WV2 images covering relatively small test sites distributed across the region of interest. Using these classification maps as training data, wall-to-wall mapping of fractional coverages of spruce and pine was achieved using multi-temporal Landsat data and Random Forests (RF) regression. The method was applied for the entire state of Bavaria (Germany), which comprises a total forested area of approximately 26,000 km2. As applied here, this two-step approach yields consistent and accurate maps of fractional tree cover estimates with a spatial resolution of 1 ha. Independent validation of the fractional cover estimates using 3780 reference samples collected through visual interpretation of orthophotos produced root-mean-square errors (RMSE) of 11% (for spruce) and 14% (for pine) with almost no bias, and R2 values of 0.74 and 0.79 for spruce and pine, respectively. The majority of the validation samples (75% (spruce) and 84% (pine)) were modeled within the assumed uncertainty of ± 15% of the reference sample. Accuracies were significantly better compared to those achieved using a single-step classification of Landsat time series data at the pixel level (30 m), because the two-step approach better captures regional variation in the spectral signatures of target classes. Moreover, the increased number of available reference cells mitigates the impact of occasional errors in the reference data set. This two-step approach has great potential for cost-effective operational mapping of dominant forest types over large areas.