1 Introduction  
2 Methodological Framework  
2.1 Gaussian Process Regression  
2.2 Horvitz-Thompson and Difference Estimator  
2.3 Bayesian Optimization  
3 Model-Assisted Sampling with Bayesian Optimization  
4 Case Study  
4.1 Response Variable: Field Plot Level Tree Volume  
4.2 Explanatory Variables: Laser Scanning Data  
4.3 Performance Metrics  
4.4 Results  
5 Discussion  
6 Conclusions