In this paper, we present anew deformable volume rendering algorithm. The volume deformation is modeled by a landmark- based volume morphing method using Hardy's scattered data interpolation. The algorithm is able to directly render the deformed volume without going through the expensive volume construction process. Piecewise linear approximation of the deformation function by adaptive space subdivision and template-based block projection are used to speed up the rendering process. Our algorithm can render the morphing of a 2563 volume in seconds, instead of minutes and hours with the traditional morphing-rendering pipeline.