# Chapter 21

2D Intersection of Curves

*Baruch Zukerman, Ron Wein, and Efi Fogel*

## Table of Contents

## 21.1 Introduction

Let $$* = {C*_{1}, C_{2}, ..., C_{n}} be a set of curves.
We wish to compute all
intersection
points between
two curves in the set in an output-sensitive manner, without having to
go over all $$*O(n*^{2}) curve pairs. To this end, we sweep an imaginary line $$*l*
from $$*x = - * to $$*x = * over the plane. While sweeping
the plane, we keep track of the order of curves intersecting it.
This order changes at a finite number of *event points*, such that
we only have to calculate the
intersection
points
between two curves when they become contiguous. For more details on the
*sweep-line algorithm* see, for example, [dBvKOS00, Chapter 2].

This chapter describes three functions implemented using the sweep-line
algorithm: given a collection of input curves, compute all
intersection
points,
compute the set of subcurves that are pairwise interior-disjoint induced by them,
and checking whether there is at least one pair of curves among them that intersect in their interior.

The implementation is robust. It supports general
curves and handles all degenerate cases, including overlapping curves,
vertical segments, and tangency between curves. The robustness of the
algorithm is guaranteed if the functions are instantiated with a traits
class that employs certified computations. This traits class must be a model
of the *ArrangementTraits_2* concept - see the
Chapter 20 for more details.

The complexity of the sweep-line algorithm is $$*O((n + k)*log$$*n)* where $$*n*
is the number of the input curves and $$*k* is the number of
intersection
points induced by these curves.

## 21.2 Example

The simple program listed below computes
intersection
points induced by
a set of four input segments illustrated in figure 21.1.

**Figure 21.1: **Four input segments

**File: **examples/Arrangement_2/sweep_line.cpp

#include <CGAL/Cartesian.h>
#include <CGAL/MP_Float.h>
#include <CGAL/Quotient.h>
#include <CGAL/Sweep_line_2_algorithms.h>
#include <list>
typedef CGAL::Quotient<CGAL::MP_Float> NT;
typedef CGAL::Cartesian<NT> Kernel;
typedef Kernel::Point_2 Point_2;
typedef Kernel::Segment_2 Segment_2;
int main()
{
// Construct the input segments.
Segment_2 segments[] = {Segment_2 (Point_2 (1, 5), Point_2 (8, 5)),
Segment_2 (Point_2 (1, 1), Point_2 (8, 8)),
Segment_2 (Point_2 (3, 1), Point_2 (3, 8)),
Segment_2 (Point_2 (8, 5), Point_2 (8, 8))};
// Compute all intersection points.
std::list<Point_2> pts;
CGAL::get_intersection_points (segments, segments + 4,
std::back_inserter (pts));
// Print the result.
std::copy (pts.begin(), pts.end(),
std::ostream_iterator<Point_2>(std::cout, "\n"));
return (0);
}

## 21.3 Design and Implementation History

The current version of the sweep-line algorithm was written by
Baruch Zukerman, based on previous implementations by Ester Ezra and
Tali Zvi.