In this paper, we introduce two-dimensional hat functions and derive operational matrix of integration of these functions. Then, we utilize them to solve some classes of integral equations. The method is based upon expanding functions as their truncated hat functions. Also, an error analysis is provided under several mild conditions. Illustrative examples are included to demonstrate the validity, efficiency and applicability of the method.