CGAL users discussion list

[VSav <>]

Hi 
Thanks a ton for the above sugestion. I tried using the ipelets function.
My code is as follows.

typedef CGAL::Exact_predicates_inexact_constructions_kernel Kernel_Exact;
typedef CGAL::Apollonius_graph_traits_2<Kernel_Exact>                         
    
AGT2_K;
typedef CGAL::Apollonius_graph_2<AGT2_K>   AG2;
typedef CGAL::Apollonius_site_2<Kernel_Exact>                   Site_2_Apo;
typedef Site_2_Apo::Point_2                                    
Site_2_Point_2;
typedef typename Kernel_Exact::Segment_2                                      
 
Segment_2;
typedef typename Kernel_Exact::Ray_2                                          
 
Ray_2;
typedef typename Kernel_Exact::Line_2                                         
 
Line_2;
typedef typename Kernel_Exact::Iso_rectangle_2                                
 
Iso_rectangle_2;




struct Voronoi_from_tri{  //Class using stream to get the voronoi diagram
      std::list<Ray_2> ray_list;
      std::list<Line_2> line_list;
      std::list<Segment_2> seg_list;

      void operator<<(const Ray_2& p){ray_list.push_back(p);}
      void operator<<(const Line_2& p){line_list.push_back(p);}
      void operator<<(const Segment_2& p){seg_list.push_back(p);}

    };


  template <class T,class output_iterator>
    bool
    cast_into_seg(const T& obj,const Iso_rectangle_2& bbox,output_iterator
out_it) {
      CGAL::Object obj_cgal = CGAL::intersection(obj,bbox);
      Segment_2 s;
      bool ret=CGAL::assign(s, obj_cgal);
     
      if (ret){
       *out_it++=s;
     }
      return ret;
    }

    //Convert infinite objects into drawable segments
    template<class iterator,class output_iterator>
    void
    cast_into_seg(const iterator first,const iterator end,
                  const Iso_rectangle_2& bbox, output_iterator out_it) 
    {
      for (iterator it=first;it!=end;++it)
      {
        cast_into_seg(*it,bbox,out_it);
      }
    }


 void draw_dual_(Voronoi_from_tri& v_recup,const Iso_rectangle_2& bbox,bool
makegrp) 
    {
      std::vector<Segment_2> seg_cont;
      //filter degenerate segments
      for(typename std::list<Segment_2>::iterator iteS =
v_recup.seg_list.begin();iteS!=v_recup.seg_list.end();){
        typename std::list<Segment_2>::iterator itc=iteS++;
        if (itc->is_degenerate())
          v_recup.seg_list.erase(itc);
      }
      
     
cast_into_seg(v_recup.ray_list.begin(),v_recup.ray_list.end(),bbox,std::back_inserter(seg_cont));//cast
rays into segments in bbox
     
cast_into_seg(v_recup.line_list.begin(),v_recup.line_list.end(),bbox,std::back_inserter(seg_cont));//cast
lines into segments in bbox
     
cast_into_seg(v_recup.seg_list.begin(),v_recup.seg_list.end(),bbox,std::back_inserter(seg_cont));//cast
lines into segments in bbox
      //draw_in_ipe(seg_cont.begin(), seg_cont.end(),makegrp);
    }



using namespace std;
int main(){
  std::ifstream ifs("sites.cin");
  assert( ifs );
  AG2 ag;         
  Site_2_Apo  site;
  while ( ifs >> site )
  {
    ag.insert(site);
    
  }

 Voronoi_from_tri v_recup;
  ag.draw_dual(v_recup);
  
  Iso_rectangle_2 bbox(-2000,-2000,2000,2000);
  draw_dual_(v_recup,bbox,true);
  return 0;
} 


Now I am not able to print the cropped apollonius diagram.
If i am able to print the diagram I will be able to construct the bounding
box around each cell of apollonius diagram.
How can I print the cropped diagram?

Thanks,
Virti




--
Sent from: http://cgal-discuss.949826.n4.nabble.com/