A method of deformable image registration for thoracic 4-D computed tomography (CT) images includes: receiving, by a processing device, a set of thoracic 4-D CT images; and iteratively solving, by the processing device, an energy function applied to subsequent images of the set of thoracic 4-D CT images to transform the subsequent images into respective optical flow fields between the subsequent images, the energy function enforcing the following constraints on the subsequent images: intensity constancy; mass conservation; gradient constancy; and spatio-temporal smoothness. A method of determining regional lung function includes: receiving, by a processing device, a set of thoracic 4-D CT images of a lung; transforming, by the processing device, the set of 4-D CT images of the lung into respective spatial voxel-wise deformation maps; and transforming, by the processing device, the spatial voxel-wise deformation maps into respective spatial voxel-wise strain maps of the lung indicating regional mechanics of the lung.