The method I'm using seems super inefficient. What I did was define $4$ RVs namely $X_1,X_2,Y_1,Y_2$ and thus my two uniformly random points are $(X_1,Y_1),(X_2,Y_2)$ and hence the $y$ intercept is $-X_1(Y_2-Y_1)/(X_2-X_1)+Y_1$ but them I'll have to do Jacobian tranform multiple times.
Is there any simpler method I seem to be missing?

I'd approach this by first computing the conditional cdf of the slope given one point $(X_1,Y_1)$, hence the conditional cdf of the $y$-intercept and integrate to get the unconditioned cdf. It still requires some integration and separation into cases but it looks more approachable this way, because the conditional cdf is just really areas of triangles and similar shapes.