We introduce a new class of sine-Gordon models, for which interaction term is present in a region different from the domain over which quadratic part is defined. We develop a novel non-perturbative approach for calculating partition functions of such models, which relies on mapping them to statistical properties of random surfaces. As a specific application of our method, we consider the problem of calculating the amplitude of interference fringes in experiments with two independent low dimensional Bose gases. We calculate full distribution functions of interference amplitude for 1D and 2D gases with nonzero temperatures.