boost::多边形的分割错误

我已经处理了多边形数据。现在我想看看处理后的数据与原始数据的匹配程度。对于这个任务，我想使用BOOST的多边形集操作符。下面的代码给了我一个segfault：

#include <iostream>
#include <boost/polygon/polygon.hpp>
using namespace boost::polygon::operators;
using namespace std;
typedef boost::polygon::polygon_data<double> BPolygon;
typedef boost::polygon::polygon_traits<BPolygon>::point_type BPoint;
typedef boost::polygon::polygon_set_data<double> BPolygonSet;
typedef std::vector<BPolygon> BPolygonVec;

double meassureError(BPolygonVec &polys1, BPolygonVec &polys2)
{
  BPolygonSet set1;
  BPolygonSet set2;
  assign(set1, polys1);
  assign(set2, polys2);
  return area(set1 ^ set2);
}
int main(int argc, char *argv[])
{
  BPolygonVec polys1;
  BPolygonVec polys2;
  loadPolysFromFile(polys1);
  loadPolysFromFile(polys2);
  cout << meassureError(polys1, polys2) << endl;
  return 0;
}

gdb输出：

Program received signal SIGSEGV, Segmentation fault.
0x08156ce7 in std::list<boost::polygon::point_data<double>, std::allocator<boost::polygon::point_data<double> > >::begin (this=0x0) at /usr/include/c++/4.8.2/bits/stl_list.h:759
759           { return iterator(this->_M_impl._M_node._M_next); }

我的数据由大约2000个多边形组成，每个多边形大约有10个顶点，我希望有足够的内存来处理它。我做错了什么？谢谢你的帮助！