Christoffel Transformation for a Matrix of Bi-variate Measures