$\dfrac{z-i}{-iz+1} = \dfrac{x+iy-i}{-i(x+iy)+1} = \dfrac{x+i(y-1)}{1+y-ix} = \dfrac{2x+[x^2+y^2-1]i}{(y+1)^2+x^2}$
so $G(x,y;0,1) = C \log \dfrac{4x^2 + (x^2+y^2-1)^2}{((y+1)^2+x^2)^2}$
$=C \log \dfrac{(x^2+1)^2 + 2y^2(x^2+1) + y^4 - 4y^2}{((y+1)^2+x^2)^2} = C\log \dfrac{(x^2+y^2+2y+1)(x^2+y^2-2y+1)}{((y+1)^2+x^2)^2} = C \log \dfrac{x^2+(y-1)^2}{x^2+(y+1)^2}$
omg @Semiclassical I derived your answer
now $G(x,y;\xi,\eta) = G\left(\frac{x-\xi}\eta,\frac y\eta;0,1\right) = C \log \dfrac{(x-\xi)^2 + (y-\eta)^2}{(x-\xi)^2 + (y+\eta)^2}$
yes indeed it matches with your answer 100%
now to test what $C$ is, we use a test function $\varphi$ right
we want $\displaystyle \int_{\Bbb R \times (\Bbb R_{\ge0})} \varphi(\xi,\eta) \nabla^2 G(x,y;\xi,\eta) \ \mathrm d\xi \ \mathrm d\eta = \varphi(x,y)$ formally
using Green's identity $\displaystyle \int_M u \nabla^2 v - v \nabla^2 u = \int_{\partial M} uNv-vNu$
the LHS becomes $\displaystyle \int_{\Bbb H} G \nabla^2 \varphi + \int_{-\infty}^\infty \varphi(x,0) G_y(x,0;\xi,\eta) - \varphi_y(x,0) G(x,0;\xi,\eta) \ \mathrm dx$
since $G$ vanishes on the x-axis, the third term becomes zero
well $G_y = C \dfrac{2(y-\eta)}{(x-\xi)^2 + (y-\eta)^2} - C \dfrac{2(y+\eta)}{(x-\xi)^2 + (y+\eta)^2}$
if we let $\varphi = 1$ then $\nabla^2 \varphi = 0$, so the equation becomes
$\displaystyle \int_{-\infty}^\infty G_y(x,0;\xi,\eta) \ \mathrm d\xi = 1$
$2C \left[ \arctan\left(\dfrac{\xi-x}{-\eta}\right) - \arctan\left(\dfrac{\xi-x}{\eta}\right) \right]_{-\infty}^\infty = 1$
$2C\left[\left(-\dfrac\pi2 - \dfrac\pi2\right) - \left( \dfrac\pi2 + \dfrac\pi2 \right)\right] = 1$
so in conclusion, $G(x,y;\xi,\eta) = \dfrac1{4\pi} \log \dfrac{(x-\xi)^2 + (y+\eta)^2}{(x-\xi)^2 + (y-\eta)^2}$