### python

#### Camera rotation matrix

I am trying to raytrace a video. For that matter, I need the camera roation matrix for every frame, in world coordinates. The camera is at the origin. There are no translation. I have the trajectory of the camera as changes in rotation for every frame. So for every frame, I have three values (roll, yaw, pitch) that describe how much the camera should rotate from this frame to the next. These rotations are to be understand in the camera coordinate system. How do I compute the world coordinate rotation matrices for the frames? What I tried: def rot_x(angle): cosa = np.cos(angle) sina = np.sin(angle) return np.array([[1,0,0], [0, cosa, -sina], [0, sina, cosa]]) def rot_y(angle): cosa = np.cos(angle) sina = np.sin(angle) return np.array([[cosa, 0, sina], [0,1,0], [-sina, 0, cosa]]) def rot_z(angle): cosa = np.cos(angle) sina = np.sin(angle) return np.array([[cosa, -sina, 0], [sina, cosa, 0], [0,0,1]]) matrices = [initial_rot] for pitch, yaw, roll in frames_data: rx = rot_x(pitch) ry = rot_y(yaw) rz = rot_z(roll) last_matrix = matrices[-1] matrices.append(last_matrix.T.dot(rx).dot(ry).dot(rz)) (Since last_matrix should be orthonormal, so the inverse of it should be the transpose). However, something is horribly wrong, the rendered video just flickers in the y dimension. I am sure there is something wrong with the math here..

The order of the matrix multiplications matters. Applying another rotation should be done by left-multiplying (assuming a standard convention). Since this is simply composition multiple rotations, there should not be any need to inverse the last rotation. The complete rotation that should be computed for frame N is: R_n = R(yaw_n, pitch_n, roll_n) R_{n - 1} R_{n - 2} ... R_1 R_0 with: R_0: the initial rotation (i.e. initial_rot) R_n: the complete rotation for the frame N R(yaw_n, pitch_n, roll_n): the rotation derived from the yaw / pitch / roll values applied between frame N - 1 and N (i.e. rx.dot(ry).dot(rz)) Hence the last line in the code excerpt should instead be: rotation = rx.dot(ry).dot(rz) matrices.append(rotation.dot(last_matrix))

### Related Links

python foreign character in csv

How to pack python files and its dependencies in a single executable file?

Printing Variable names and contents as debugging tool; looking for emacs/Python shortcut

Cheking added file to upload python, pylons?

How to refer to the local module in Python?

Is close() necessary when using iterator on a Python file object [duplicate]

Django Admin “Edit Selection” Action?

How to change firefox proxy from webdriver?

Is it possible to hook up a more robust HTML parser to Python mechanize?

Enable Unicode “globally” in Python

Dynamically import a callable given the full module path?

python chaining

py2app and xml.etree.ElementTree

What is the difference between isinstance('aaa', basestring) and isinstance('aaa', str)?

Is this essential functional programming feature missing from python?

Hooking into a wave-out on different platforms