xNDC = (xworld-Bx)*2/(Tx-Bx) -1 yNDC = (yworld-By)*2/(Ty-By) -1
xDC = (xNDC+1)(tx-bx)/2 + bx yDC = (-1*yNDC+1)(ty-by)/2 + by
| a b c d| |A| | e f g h| |B| | i j k l|x|C| = [Aa+Bb+cC+Dd, Ae+Bf+Cg+Dh, Ai+Bj+Ck+Dl, Am+Bn+Co+Dp] | m n o p| |D|
(xworld-Bx)*2/(Tx-Bx) -1
= 2xworld/(Tx-Bx) -2Bx/(Tx-Bx) -1
= 2/(Tx-Bx) xworld + (-2Bx - (Tx-Bx))/(Tx-Bx)
= 2/(Tx-Bx) xworld + (-Bx - Tx)/(Tx-Bx)
= 2/(Tx-Bx) xworld + -(Bx + Tx)/(Tx-Bx)
2 Bx+Tx
or a = ----- d = - -----
Tx-Bx Tx-Bx
e = 0, g = 0
2 By+Ty
f = ---- h = - -----
Ty-By Ty-By
| 2 Bx+Tx |
| ----- 0 0 - ----- |
| Tx-Bx Tx-Bx |
|
| 2 By+Ty |
T = | 0 ---- 0 - ----- |
| Ty-By Ty-By |
|
| 0 0 1 0 |
| 0 0 0 1 |