For the problem of numerical differentiation of bivariate functions with mixed smoothness, we determine the exact orders for the minimal radius of Galerkin information and construct a truncation method, which is optimal in a sense of the indicated quantity.