Apparatus for Determining a disease parameter (κ) describing a type of cardiac disease, the apparatus comprising the following units: - a unit for reading in a set of input parameters (f), wherein the input parameters (f) represent a field of velocity vectors measuring 3D myocardial velocity fields with voxel wise 3D velocity information, each velocity vector representing a velocity of motion of a voxel of a region of the heart at a predefined time point or time interval and wherein a velocity vector depicts the voxel wise three-dimensional velocity of the myocardium with different magnitude; - a unit for applying a processing algorithm on the set of input parameters (f) in order to obtain at least a rotation parameter (R) and/or at least a strain-rate-tensor (SRT), the rotation parameter (R) representing the local curl of the set of input parameters (f) and the strain-rate-tensor (SRT) representing a local rate of strain at a certain spatial point at a certain time point or time interval, wherein the strain-rate-tensor (SRT) is calculated by the Jacobian matrix (Df) of the considered velocity field (f), using the formula SRT= 1⁄2(Df+Df T ) wherein Df T denotes a transposed Jacobian matrix and wherein the rotation parameter (R) is derived from the formula curl = ∇ x f wherein ∇ represents the Nabla-operator and x denotes the cross product or wherein the rotation parameter (R) is calculated using the formula RRT= 1⁄2(Df-Df T ); and - a unit for calculating at least the disease parameter (κ) on the basis of the rotation parameter (R) and the strain-rate-tensor (SRT), wherein the unit for calculating is configured for calculating the disease parameter (κ) by dividing an absolute value of the rotation parameter (R) by at least one Eigenvalue, the at least one Eigenvalue resulting from solving the eigenvector problem SRT f = λ f, wherein λ denotes the eigenvalue and wherein the disease parameter (κ) reflects the proportion to overall deformation of local layer shearing with regard to local pure deformation. - a unit for reading in a set of input parameters (f), wherein the input parameters (f) represent a field of velocity vectors measuring 3D myocardial velocity fields with voxel wise 3D velocity information, each velocity vector representing a velocity of motion of a voxel of a region of the heart at a predefined time point or time interval and wherein a velocity vector depicts the voxel wise three-dimensional velocity of the myocardium with different magnitude; - a unit for applying a processing algorithm on the set of input parameters (f) in order to obtain at least a rotation parameter (R) and/or at least a strain-rate-tensor (SRT), the rotation parameter (R) representing the local curl of the set of input parameters (f) and the strain-rate-tensor (SRT) representing a local rate of strain at a certain spatial point at a certain time point or time interval, wherein the strain-rate-tensor (SRT) is calculated by the Jacobian matrix (Df) of the considered velocity field (f), using the formula SRT= 1⁄2(Df+Df T ) wherein Df T denotes a transposed Jacobian matrix and wherein the rotation parameter (R) is derived from the formula curl = ∇ x f wherein ∇ represents the Nabla-operator and x denotes the cross product or wherein the rotation parameter (R) is calculated using the formula RRT= 1⁄2(Df-Df T ); and - a unit for calculating at least the disease parameter (κ) on the basis of the rotation parameter (R) and the strain-rate-tensor (SRT), wherein the unit for calculating is configured for calculating the disease parameter (κ) by dividing an absolute value of the rotation parameter (R) by at least one Eigenvalue, the at least one Eigenvalue resulting from solving the eigenvector problem SRT f = λ f, wherein λ denotes the eigenvalue and wherein the disease parameter (κ) reflects the proportion to overall deformation of local layer shearing with regard to local pure deformation. <EOS>