The performance and functionality of integrated photonic devices can be enhanced by using complex structures controlled by a large number of design variables. However, the optimization of such high-dimensional structures is challenging, often limiting their realization. Global optimization algorithms and artificial neural networks are increasingly used to tackle these problems. Although these are exciting new developments, the outcome is a single optimized design meeting particular performance objectives selected upfront. The influences of the various design parameters remain hidden. Here we report on our strategy of using machine learning pattern recognition techniques to create a methodology for building the global performance map of a high-dimensional design space. As an example and demonstration, we study the design of a vertical grating coupler consisting of silicon and subwavelength metamaterial segments. We show how the relationship between designs with comparable primary performance can be clearly revealed by identifying the minimum number of characterizing parameters that defines the subspace of good designs, significantly scaling down the complexity of the problem. Moreover, the subspace can be identified using only a small number of good design solutions. We reveal design areas with comparable fiber coupling efficiency but with significant differences in other performance criteria, such as back-reflections, tolerance to fabrication uncertainty and minimum feature size. This novel approach provides the designer with a global perspective of the design space, enabling informed decisions based on the relative priorities of all relevant performance specifications and figures-of-merits for a particular application. Insights from the mapping exercise also inspired new design structures with enhanced characteristics.