Compute the rotation matrix to transform a unit vector
into the local frame with axises defined by:
R[0] (first axis) = unit vector towards Zenith
R[1] (second axis) = unit vector towards East
R[2] (second axis) = unit vector towards North
(the projected x,y are along R[1] and R[2] axises)