pass constraints from delauneytri to triscatteredinterp

I have given delauneytri boundary constraints - but when I pass the results to triscatteredinterp (F=triscatteredinterp(dt,datavalues), evaluate the function (F(qx,qy)) on a meshgrid and then contourf it - the areas I have removed are still on the contour plot... I could probably remove the areas between evaluating the function and contouring it - however this seems to defeat the point of constraining the boundary before gridding... The problem is that the edges of the data are jagged - and I don't want the edges to influence each other across the boundaries. I am pretty sure the boundary constraints have worked ok within the delauneytri - because if I plot the triangulation - the triangles outside the boundaries have been removed... So the problem seems to be that triscatteredinterp interpolates across the boundaries even with the constraints in place.

Does anyone have an idea of how to work around this?? It is not something I'm very familiar with - so any advice would be appreciated....

I am using matlab 2012a

Thank you!!!