An L-2-optimal identification method is extended to cope with MIMO errors-in-variables (EIV) model estimation based on a geometrical interpretation for the v-gap metric. The L-2-optimal approximate models are composed of system and noise models and characterised by a normalised right graph symbol (NRGS) and its complementary inner factor (CIF), respectively. This metric can be evaluated as the supreme of sine values of the maximal principal angles between NRGS frequency responses of two concerned models. In order to make full use of the angular cosine formula for complex vectors to reduce computational loads, a CIF of the NRGS of the perturbed model is introduced and thus, the system parameter optimisation can be efficiently solved by sequential quadratic programming methods. With the estimated system model, the associated noise model can be built by right multiplication of an inner matrix. Finally, a simulation example demonstrates the effectiveness of the proposed identification method.