In the current study, a semi-analytical model for contact between a homogeneous, isotropic, linear elastic half-space with a geometrically anisotropic (wavelengths are different in the two principal directions) bisinusoidal surface on the boundary and a rigid base is developed. Two asymptotic loads to area relations for early and almost complete contact are derived. The Hertz elliptic contact theory is applied to approximate the load to area relation in the early contact. The noncontact regions occur in the almost complete contact are treated as mode-I cracks. Since those cracks are in compression, an approximate relation between the load and noncontact area can be obtained by setting the corresponding stress intensity factor (SIF) to zero. These two asymptotic solutions are validated by two different numerical models, namely, the fast Fourier transform (FFT) model and the finite element (FE) model. A piecewise equation is fit to the numerical solutions to bridge these two asymptotic solutions.