如何找到包含所有給定框的框？

Question

我有一個boost幾何框類型的序列，並希望找到包含所有這些對象的框的尺寸。

我注意到boost幾何提供了encompass函數，這似乎是我要求的一般幾何概念，但我不知道如何對一系列方框執行此操作。

基本上我不得不推出我自己的，但我想知道是否有可能簡單地將一系列盒子變成一個「幾何」，以便我可以簡單地應用envelope函數。這是我目前寫的：

// elsewhere in my code: 
using Location = boost::geometry::model::d2::point_xy<double>; 
using Box = boost::geometry::model::box<Location>; 

// calculate the smallest box that fully encompasses all provided boxes    
template <template <typename...> class IterableContainer>       
Box Envelope(const IterableContainer<Box>& sequence)        
{                     
    Location minCorner(               
     std::numeric_limits<double>::max(),           
     std::numeric_limits<double>::max());          

    Location maxCorner(               
     std::numeric_limits<double>::lowest(),           
     std::numeric_limits<double>::lowest());          

    for (auto& box : sequence)              
    {                    
     if (box.min_corner().x() < minCorner.x())         
      minCorner.x() == box.min_corner().x();         

     if (box.min_corner().y() < minCorner.y())         
      minCorner.y() == box.min_corner().y();         

     if (box.max_corner().x() > maxCorner.x())         
      maxCorner.x() == box.max_corner().x();         

     if (box.max_corner().y() > maxCorner.y())         
      maxCorner.y() == box.max_corner().y();         
    }                    

    return Box(minCorner, maxCorner);            
}                     

Answer 1

你正在尋找的功能叫做std::accumulate。你需要爲它提供一個box-union函數。

using Location = boost::geometry::model::d2::point_xy<double>; 
using Box = boost::geometry::model::box<Location>; 
double pinfi = std::numeric_limits<double>::max(); 
double ninfi = std::numeric_limits<double>::lowest(); 

Box u = std::accumulate(container.begin(), container.end(), 
         Box(Location(pinfi,pinfi), Location(ninfi,ninfi)), 
         [](const Box& a, const Box& b) { 
          return Box(
            Location(
            std::min(a.min_corner().x(),b.min_corner().x()), 
            std::min(a.min_corner().y(),b.min_corner().y())), 
            Location(
            std::max(a.max_corner().x(),b.max_corner().x()), 
            std::min(a.max_corner().y(),b.max_corner().y()))); 

         }); 

更新：此函數的構建塊已經存在於boost :: geometry中。這裏是完整的測試代碼：

template <typename T> 
Box box_encompass (T beg, T end) 
{ 
    return std::accumulate(beg, end, boost::geometry::make_inverse<Box>(), 
       [](Box a, const Box& b) -> Box { 
        boost::geometry::expand(a,b); 
        return a; 
       }); 
}