Existe uma maneira simples de descobrir se um ponto está dentro de uma célula voronoi?

Por exemplo, o código a seguir gera algo como o diagrama abaixo:

using namespace boost::polygon;

point_data<int> p1(0, 0);
point_data<int> p2(-10, 10);
point_data<int> p3(-10, -10);
point_data<int> p4(10, -10);
point_data<int> p5(10, 10);

std::vector<point_data<int>> pts = { p1, p2, p3, p4, p5 };
construct_voronoi(pts.begin(), pts.end(), vd);

Nesse caso, como posso descobrir se o ponto (5,5) está dentro da célula central?

Eu poderia criar um polígono a partir de cada célula e descobrir usando umponto no algoritmo de polígono, mas estou interessado em saber que a biblioteca oferece algo "de graça".