It appears to me that the notation represents taking an arbitrary 3-vector, applying a rotation matrix to the 3-vector based on the bearing 3-vector, and then removing one of the dimensions of the result. If the arbitrary 3-vector and the bearing 3-vector are the same, then the result of the rotation should be [ 1 0 0 ], and the first coordinate would be removed, since the projection of that vector onto the plane orthogonal to it is zero. For an arbitrary 3-vector, you get the components of the projection onto a plane that is normal to the bearing 3-vector.
As to how to compute the rotation matrix based on the bearing 3-vector, I would defer to the following:
https://math.stackexchange.com/questions/114512/how-to-find-the-orthonormal-transformation-that-will-rotate-a-vector-to-the-x-ax