So am I going to try my best to answer this, don't be too harsh on me if it is all wrong!
I think you will want to look at the Cauchy-Riemann equations in polar coordinates, namely:
$$\frac{\partial u}{\partial r} = \frac{1}{r} \frac{\partial v}{\partial \theta}$$ and $$\frac{ \partial v}{\partial r} = - \frac{1}{r} \frac{\partial u}{\partial \theta}$$
The reason I say this is because we are going to be integrating using polar coordinates. What we are going to do is given our harmonic function $u(r,\theta)$ we will take partial derivatives with respect to $r$ and $\theta$ giving us $u_r$ and $u_\theta$ respectively. We will then solve for $v$ using the Cauchy Riemann equations and integrating.
To this end, set
$$v(r,\theta) = \int_a^r -\frac{1}{t} u_\theta(t,0) dt + \int_0^{\theta} r u_r(r,\psi) d\psi$$
We can then check that $u$ and $v$ satisfy the Cauchy-Riemann equations by taking partial derivatives.
We get:
$$\frac{\partial v}{\partial \theta} = r u_r = r \frac{\partial u}{\partial r}$$
by the Fundamental Theorem of calculus and
$$\frac{\partial v}{\partial r} = - \frac{1}{r}u_\theta(r, 0) + \int_0^\theta \frac{\partial}{\partial r} r u_r (r, \psi) d\psi$$
(we can differentiate under the integral sign because $u$ is $C^2$ and the region is compact)
$$= \int_0^\theta u_r(r,\psi)d\psi + \int_0^\theta r \frac{\partial^2 u}{\partial r^2} (r, \psi) d\psi$$
Now, $u$ is harmonic, so the Laplace equation in polar coordinates gives us
$$\frac{\partial^2 u}{\partial r^2} = -\frac{1}{r}\frac{\partial u}{\partial r} -\frac{1}{r^2}\frac{\partial^2 u}{\partial \theta^2}$$
Hence, the above is
$$= \int_0^\theta u_r(r,\psi)d\psi + \int_0^\theta r \left( -\frac{1}{r}u_r(r,\psi) -\frac{1}{r^2}u_{\theta\theta}(r,\psi) \right ) d\psi$$
$$=\int_0^\theta u_r(r,\psi)d\psi - \int_0^\theta u_r(r,\psi)d\psi + \int_0^\theta -\frac{1}{r}u_{\theta\theta}(r,\psi) d\psi= -\frac{1}{r}u_\theta$$
Hence, $$\frac{\partial v}{\partial r} = -\frac{1}{r}\frac{\partial u}{\partial \theta}$$
So $u,v$ satisfy the Cauchy Riemann equations.