#include <boost/geometry.hpp>
#include <boost/geometry/geometries/geometries.hpp>
#include <boost/geometry/multi/geometries/multi_point.hpp>
#include <boost/geometry/algorithms/perimeter.hpp> 
#include <vector>
using namespace std;

#define ADD(X,Y) append(all_points_in_radius,make<point_xy<float> >(X,Y))

int main()
{
	using boost::geometry::model::d2::point_xy;
	using boost::geometry::append;
	using boost::geometry::make;

	boost::geometry::model::polygon<point_xy<float> > all_points_in_radius;

	ADD(0,0); 
	ADD(0,2);

	cout << "num inner rings in poly: " << all_points_in_radius.inners().size() << endl;

	
	typedef vector<point_xy<float> > pointvector;
	pointvector &r = all_points_in_radius.outer();
	cout << "points in poly:" << endl;
	for (pointvector::iterator pi = r.begin(); pi!=r.end(); ++pi)
		cout << pi->x() << "," << pi->y() << endl;

	boost::geometry::model::polygon<point_xy<float> > hull;
	boost::geometry::convex_hull(all_points_in_radius,hull);
	pointvector &p = hull.outer();
	cout << "points in hull:" << endl;
	for (pointvector::iterator pi = p.begin(); pi!=p.end(); ++pi)
		cout << pi->x() << "," << pi->y() << endl;
	cout << "hull perimiter:" << endl;
	cout << boost::geometry::perimeter(hull) << endl;
	cout << "hull area:" << endl;
	cout << boost::geometry::area(hull) << endl;
	
	return 0;
}