Issue with RegionFunction as an option in ListContourPlot

I’ve come across an issue using the RegionFunction option in the ListContourPlot function. The issue is illustrated in the attached notebook. I’ve defined a (concave) boundary mesh region for a set of points and I want to use the boundary mesh to restrict the display of contours such that they only appear within the mesh region. However, the functionality does not seem to be working as I expect it to. I’ve tried using a different region (i.e. a Disk) to confine the display of contours. And that does work as expected. I’d be most grateful if anyone can spot any issues with what I’ve done or suggest a workaround. It may also be that I’ve stumbled upon a bug.

Thanks in anticipation.

Update…

The attached notebook shows a workaround (scroll down to the notebooks penultimate section). The workaround is far from ideal. But it does at least work as I expect. And seems to suggest that there may be a bug in the way that RegionFunction is working in ListContourPlot.

However, I’m still not wholly convinced that I’m not doing something wrong and would welcome any comments. Especially from Wolfram insiders who might be able to confirm whether or not there’s a bug.

Update…

Wolfram Support has confirmed that this issue is a bug. Hopefully, it’ll be fixed soon.