Abstract Combining forest growth models with remotely sensed data is possible under a generalized hierarchical model-based (GHMB) inferential framework. This implies the existence of two submodels: the growth model itself ($\mathcal{M}_{1}$) and a second submodel that links the growth predictions to some remotely sensed variables ($\mathcal{M}_{2}$). Analytical GHMB estimators are available to fit submodel $\mathcal{M}_{2}$ and account for the uncertainty stemming from submodel $\mathcal{M}_{1}$, i.e. the growth model. However, when the growth model is individual based, it is usually too complex to be differentiated with respect to its parameters. As a result, the analytical GHMB estimators cannot be used. In this study, we developed a bootstrap approach for the GHMB inferential framework in order to combine individual-based forest growth models with remotely sensed data. Through simulation studies, we showed that the bootstrap estimators were nearly unbiased when both submodels were linear. The estimator of the parameter estimates remained nearly unbiased when submodel $\mathcal{M}_{1}$ became complex, i.e. non-differentiable, and submodel $\mathcal{M}_{2}$ was nonlinear with heterogeneous variances and correlated error terms. The variance estimator showed some biases but these were relatively small. We further demonstrated through a real-world case study that the predictions of a complex individual-based model could be linked to a Landsat-8 near-infrared spectral band in the boreal forest zone of Quebec, Canada.