We have developed a method for analyzing motion at skeletal joints based on the 3D reconstruction of magnetic resonance (MR) image data. Since the information about each voxel in MR images includes its location in the scanner, it follows that information is available for each organ whose 3D surface is computed from a series of MR slices. In addition, there is information on the shape and orientation of each organ, and the contact areas of adjacent bones. By collecting image data in different positions we can calculate the motion of the individual bones. We have used this method to study human foot bones, in order to understand normal and abnormal foot function. It has been used to evaluate patients with tarsal coalitions, various forms of pes planus, ankle sprains, and several other conditions. A newly described feature of this system is the ability to visualize the contact area at a joint, as determined by the region of minimum distance. The display of contact area helps understand abnormal joint function. Also, the use of 3D imaging reveals motions in joints which cannot otherwise be visualized, such as the subtalar joint, for more accurate diagnosis of joint injury.