In the present paper, a computational method for solving nonlinear Volterra-Fredholm Hammerestein integral equations is proposed by using compactly supported semiorthogonal cubic B-spline wavelets as basis functions. Dual functions and Operational matrices of B-spline wavelets via Galerkin method are utilized to reduce the computation of integral equations to some algebraic system, where in the Galerkin method dual of B-spline wavelets are applied as weighting functions. The method is computationally attractive, and applications are demonstrated through illustrative examples.