Regression models based on spectral indices are typically empirical formulae enabling the mapping of biophysical parameters derived from Earth Observation (EO) data. Due to its empirical nature, it remains nevertheless uncertain to what extent a selected regression model is the most appropriate one, until all band combinations and curve fitting functions are assessed. This paper describes the application of a Spectral Index (SI) assessment toolbox in the Automated Radiative Transfer Models Operator (ARTMO) package. ARTMO enables semi-automatic retrieval and mapping of biophysical parameters from optical remote sensing observations. The SI toolbox facilitates the assessment of biophysical parameter retrieval accuracy of established as well as new and generic SIs. For instance, based on the SI formulation used, all possible band combinations of formulations with up to ten bands can be defined and evaluated. Several options are available in the SI assessment: calibration/validation data partitioning, the addition of noise and the definition of curve fitting models. To illustrate its functioning, all two-band combinations according to simple ratio (SR) and normalized difference (ND) formulations as well as various fitting functions (linear, exponential, power, logarithmic, polynomial) functions led to similar best correlations, though spatial differences emerged when applying the functions to HyMap imagery. We suggest that a systematic SI assessment is a strong requirement in the quality assurance approach for accurate biophysical parameter retrieval.