How to improve "pointLocation" detection on the outer boundary of a "delaunayTriangulation" object?

I have a "delaunayTriangulation" object of a collection of points, and I am trying to use "pointLocation" to find the enclosing triangle of a query point. According to the "inpolygon" function, the query point is on the edge of the polygon defined by the convex hull of the collection of points. However, the "pointLocation" method of the "delaunayTriangulation" object indicates that the query point lies outside the convex hull. 

How can I correctly identify the enclosing triangle of a query point on the outer boundary of a "delaunayTriangulation" object?

Answer 1

[MathWorks Support Team]

To identify the enclosing triangle of a query point on the outer boundary of a "delaunayTriangulation" object, please convert your "delaunayTriangulation" object into a regular "triangulation" object. From there, you can use the "pointLocation" method of the "triangulation" object to obtain more accurate results.

For example, consider the following code, where "DT" is the "delaunayTrianguation" object and "P" is the query point:

T = triangulation(DT.ConnectivityList, DT.Points); % convert delaunayTriangulation into triangulation
[ID, B] = pointLocation(T, P); % find the enclosing triangle and barycentric coordinates of the query point