title: How can I get the region associated to a given point in voronoi_finite_polygons_2d function?
body: <p>I am trying to adapt <a href="https://stackoverflow.com/a/20678647/2443944">this code</a> I found on stackoverflow to create a voronoi cell with finite boundaries.</p>
  
  <p>However my problem is that I don't know how to get the region associated to a given point. This was done in normal Voronoi with the <code>point_region</code> method, but this doesn't work here because regions have changed.</p>
  
  <p>The data points I am using are:</p>
  
  <pre><code>points = array([[289255.176  , 921667.461  ],
         [289296.31699, 921687.13826],
         [289463.30305, 921770.12504],
         [289725.08002, 921905.75745],
         [289960.48198, 922099.46056],
         [290106.98928, 922361.79529],
         [289255.184  , 921646.244  ],
         [289307.48677, 921627.05485],
         [289500.80493, 921555.50067],
         [289825.14532, 921435.65147],
         [290141.79326, 921322.77935],
         [290454.91721, 921211.09355],
         [289255.187  , 921635.627  ],
         [289327.07776, 921558.85263],
         [289565.21795, 921298.17707],
         [289875.40176, 920978.013  ],
         [290192.86361, 920656.82017],
         [289255.185  , 921630.386  ],
         [289318.54181, 921453.18492],
         [289421.06861, 921167.57934],
         [289565.42462, 920770.1386 ],
         [289701.83141, 920376.28627],
         [289833.6501 , 919990.66467]])
  
  vor = Voronoi(points)
  
  min_x = vor.min_bound[0] - 100
  max_x = vor.max_bound[0] + 100
  min_y = vor.min_bound[1] - 100
  max_y = vor.max_bound[1] + 100
  
  regions, vertices, pts = voronoi_finite_polygons_2d(vor)
  </code></pre>