Image registration is one of today's challenging image processing problems, particularly in medical imaging. Since the problem is ill posed, one may like to add additional information about distortions. This applies, for example, to the registration of time series of contrast-enhanced images, where variations of substructures are not related to patient motion but to contrast uptake. Here, one may only be interested in registrations which do not alter the volume of any substructure. In this paper, we discuss image registration techniques with a focus on volume preserving constraints. These constraints can reduce the non-uniqueness of the registration problem significantly. Our implementation is based on a constrained optimization formulation. Upon discretization, we obtain a large, discrete, highly nonlinear optimization problem and the necessary conditions for the solution form a discretized nonlinear partial differential equation. To solve the problem we use a variant of the sequential quadratic programming method. Finally, we present results on synthetic as well as on real-life data.