Without an example it is difficult to address your specific issue. Here is
an interpolation function:
data =
Table[
With[
{r = RandomReal[{0, 5}], t = RandomReal[{0, 2 Pi}]},
{r Cos[t], r Sin[t], Sin[r^2]/r^2}],
{2*10^4}];
f = Interpolation[{Most[#], Last[#]} & /@ data,
InterpolationOrder -> 1];
g[x_?NumericQ, y_?NumericQ] :=
f[x, y]*Boole[x^2 + y^2 <= 25]
Grid[{
{ListDensityPlot[data,
ImageSize -> 300],
Plot3D[g[x, y], {x, -5, 5}, {y, -5, 5},
PlotRange -> All,
RegionFunction ->
Function[{x, y, z}, x^2 + y^2 < 25],
ImageSize -> 300]},
{DensityPlot[g[x, y], {x, -5, 5}, {y, -5, 5},
RegionFunction ->
Function[{x, y, z}, x^2 + y^2 <= 25],
ImageSize -> 300],
ContourPlot[g[x, y] == 0, {x, -5, 5}, {y, -5, 5},
RegionFunction ->
Function[{x, y, z}, x^2 + y^2 <= 25],
ImageSize -> 300]}}] // Quiet
Bob Hanlon