az_L = 118.38333333333334
D_L = 0.002386
ua = 1.495978707061E8
R_L = 1737.4
α_L: atan(2 * R_L / (D_L * ua)) * 180 / π
A = (210.4479658469958, 389.0306870855365)
B = (386.80621413722935, 390.4196002860302)
c_L = Circle(A, B)
r_L = Radius(c_L)
dppx: α_L / (r_L * 2)
P_B = (518.1495614985844, 342.35248576053334)
P_H = (520.7753216694111, 526.1972221894678)
l_H: y = y(P_H)
l_B: y = y(P_B)
l = Distance(l_H, l_B)
β: l * dppx
h_{col} = 6.2
d_T: h_{col} / tan(β * π / 180)
C = (25.837407842274814, 29.52365661240635)
D = (131.84984976172146, 29.99692644240388)
gmap: 200 / (x(D) - x(C))
d_{pix}: d_T / gmap
T = (759.1905483328298, 149.903478409643)
a: -az_L + 180 + 90
c_O = Circle(T, d_T / gmap)
l_V = Ray(T, Vector[(cos((-az_L + 180 + 90) * π / 180), sin((-az_L + 180 + 90) * π / 180))])
P = Intersect(c_O, l_V)