The potential of visible-near-infrared spectroscopy to determine selected individual and total glucosinolates in broccoli has been evaluated. Modified partial least-squares regression was used to develop quantitative models to predict glucosinolate contents. Both the whole spectrum and different spectral regions were separately evaluated to develop the quantitative models; in all cases the best results were obtained using the near-infrared zone between 2000 and 2498 nm. These models have been externally validated for the screening of glucoraphanin, glucobrassicin, 4-methoxyglucobrassicin, neoglucobrassicin, and total glucosinolates contents. In addition, discriminant partial least-squares was used to distinguish between two possible broccoli cultivars and showed a high degree of accuracy. In the case of the qualitative analysis, best results were obtained using the whole spectrum (i.e., 400-2498 nm) with a correct classification rate of 100% in external validation being obtained.