This paper presents a boundary integral equation method for conformal mapping of unbounded multiply connected regions onto circular slit regions. Three linear boundary integral equations are constructed from a boundary relationship satisfied by an analytic function on an unbounded multiply connected region. The integral equations are uniquely solvable. The kernels involved in these integral equations are the classical and the adjoint generalized Neumann kernels. Several numerical examples are presented.