A high precision two-dimension numerical integral formula is given. The developed formula is constructed by means of an optimal binary interpolation formula given in a reproducing kernel space. An irregular node set is considered.