Without loss of generality you may assume the corners to have the following coordinates:
$$(a, 0)\qquad (b, 0)\qquad (c, h)\qquad (d, h)$$
Then the points you mention have coordinates
$$
\begin{array}{ccc}
\left(\frac{a+b}2, 0\right) &\qquad&
\left(\frac{c+d}2, h\right) \\
\left(\frac{ad-bc}{a-b-c+d}, \frac{(a-b)h}{a-b-c+d}\right)
&\qquad&
\left(\frac{ac-bd}{a-b+c-d}, \frac{(a-b)h}{a-b+c-d}\right)
\end{array}
$$
All of them are incident to the line
$$2h\,x + (a+b-c-d)y=(a+b)h$$
I actually computed these coordinates and the corresponding line using projective geometry, where both operations (joining points and intersecting lines) can be expressed easily using the cross product.