Method ClusterizeStreets

libs/search/geocoder.cpp:1262–1328

template <class FnT>
void Geocoder::CentersFilter::ClusterizeStreets(std::vector<uint32_t> & streets, Geocoder const & geocoder,
                                                FnT && fn) const
{
  std::vector<std::tuple<double, m2::PointD, uint32_t>> loadedStreets;
  loadedStreets.reserve(streets.size());

  // Calculate {distance, center, feature id}.
  for (uint32_t fid : streets)
  {
    m2::PointD ftCenter;
    if (geocoder.m_context->GetCenter(fid, ftCenter))
    {
      double minDist = std::numeric_limits<double>::max();
      for (auto const & c : m_centers)
        minDist = std::min(minDist, ftCenter.Length(c));
      loadedStreets.emplace_back(minDist, ftCenter, fid);
    }
    else
    {
      // In general, we don't have centers for newly created features, but editor doesn't support streets now.
      ASSERT(false, ("Street feature without table's center"));
    }
  }

  // Sort by distance.
  std::sort(loadedStreets.begin(), loadedStreets.end(),
            [](auto const & t1, auto const & t2) { return std::get<0>(t1) < std::get<0>(t2); });

  buffer_vector<m2::RectD, 2> rects(m_centers.size());
  for (size_t i = 0; i < rects.size(); ++i)
  {
    rects[i] = mercator::RectByCenterXYAndSizeInMeters(m_centers[i],
                                                       geocoder.m_params.m_filteringParams.m_streetSearchRadiusM);
  }

  // Find the first (after m_maxStreetsCount) street that is out of the rect's bounds
  size_t count = std::min(loadedStreets.size(), geocoder.m_params.m_filteringParams.m_maxStreetsCount);
  for (; count < loadedStreets.size(); ++count)
  {
    bool const outside =
        std::all_of(rects.begin(), rects.end(), [pt = std::get<1>(loadedStreets[count])](m2::RectD const & rect)
    { return !rect.IsPointInside(pt); });
    if (outside)
      break;
  }

  // Emit results.
  streets.clear();
  if (count == 0)
    return;

  double firstDist = std::get<0>(loadedStreets[0]);
  for (size_t i = 0; i < count; ++i)
  {
    double const currDist = std::get<0>(loadedStreets[i]);
    if (currDist - firstDist > geocoder.m_params.m_filteringParams.m_streetClusterRadiusMercator)
    {
      std::sort(streets.begin(), streets.end());