A preliminary comment : It is much less usual to have curves defined by equations $\phi=F(r)$ than curves of the form $r=G(\phi)$.
When $F$ has a known explicit inverse formula, we can revert to the second form. But here, this is not the case. Therefore we have to use another "trick".
- Your last formula is correct due to
$$\phi'(r)=\arctan(r)+\frac{1}{1+r^2}$$
which is positive for all $r \ge 0$. Therefore function $\phi$ is increasing on any interval $[0,r_0]$.
- The fact that $\phi$ is an increasing function is the key for finding the bounds. Let us recall the following property : a continuous increasing function on a closed bounded interval $[a,b]$ is a bijection from $[a,b]$ onto $[f(a),f(b)]$.
Here, it is sufficient to verify that
$$[\phi(a),\phi(b)]=[\phi(0),\phi(\sqrt{3})].$$
Indeed, one can check that
Therefore, using (implicitly) the reciprocal function $\phi^{-1}$ the bounds for variable $r$ are $r_1=0$ and $r_2=\sqrt{3}.$

Straight line segment has length $\sqrt{3}$ and polar angle $\tfrac{\pi}{\sqrt{3}} \approx \tfrac{4\pi}{7}.$
Matlab program for this figure :
clear all;close all;hold on;axis equal
k=0;
for r=0:0.01:(pi/sqrt(3))
k=k+1;phi=r*atan(r);
T(1,k)=r*cos(phi);T(2,k)=r*sin(phi);
end;
plot(T(1,:),T(2,:),'r'); % curve
plot([0,T(1,end)],[0,T(2,end)],'r'); % straight line