This article addresses visualization of deformation or shape differences between bones while conventional visualization techniques are often about a single bone such as its 3D reconstruction. Given a pair of bones with a set of corresponding anatomical landmarks, we compute displacement vectors describing the deformation from one bone to the other at the landmark points on one of the bones. Out of these prescribed ones, a displacement vector at each vertex on the bone surface is derived using a multi-level approximation technique of scattered data. Based on the value of inner product between a displacement vector and a surface normal at each vertex, color is mapped. Considering error-prone estimation of the landmark location in real applications, approximations at different levels are realized instead of exact interpolation as usually done in elastic image registration. Experimental results on a pair of femoral bones are presented.