Changes default displacement values for excluded image-frames

was 0, now -1 if image-frame excluded or no overlap with masked region

also cleans up displacement related debug logging statements

packages/segmentation/applications/reconstructionCardiac.cc

@@ -1334,25 +1334,24 @@ int main(int argc, char **argv)
 
       // Calculate Displacements Relative to Alignment
       mean_displacement.push_back(reconstruction.CalculateDisplacement(transformation_recon_to_ref));
-      if (debug)
-        cout<<"Mean Displacement (iter "<<iter<<") = "<<mean_displacement[iter]<<" mm."<<endl;
       mean_weighted_displacement.push_back(reconstruction.CalculateWeightedDisplacement(transformation_recon_to_ref));
-      if (debug)
-        cout<<"Mean Weighted Displacement (iter "<<iter<<") = "<<mean_weighted_displacement[iter]<<" mm."<<endl;
 
     } 
     else {
 
       // Calculate Displacement
       mean_displacement.push_back(reconstruction.CalculateDisplacement());
-      if (debug)
-        cout<<"Mean Displacement (iter "<<iter<<") = "<<mean_displacement[iter]<<" mm."<<endl;
       mean_weighted_displacement.push_back(reconstruction.CalculateWeightedDisplacement());
-      if (debug)
-        cout<<"Mean Weighted Displacement (iter "<<iter<<") = "<<mean_weighted_displacement[iter]<<" mm."<<endl;
 
     }
 
+    // Display Displacements and TRE
+    if (debug) {
+      cout<<"Mean Displacement (iter "<<iter<<") = "<<mean_displacement[iter]<<" mm."<<endl;
+      cout<<"Mean Weighted Displacement (iter "<<iter<<") = "<<mean_weighted_displacement[iter]<<" mm."<<endl;
+    }
+
+    // Save Info for Iteration
     if(debug)
     {
       cout<<"SlicesInfoCardiac4D"<<endl;

packages/segmentation/src/irtkReconstructionCardiac4D.cc

@@ -2031,7 +2031,7 @@ double irtkReconstructionCardiac4D::CalculateDisplacement()
 
     disp_sum_slice = 0;
     num_voxel_slice = 0;
-    slice_disp = 0;
+    slice_disp = -1;
     tx_slice = 0;
     ty_slice = 0;
     tz_slice = 0;
@@ -2077,7 +2077,7 @@ double irtkReconstructionCardiac4D::CalculateDisplacement()
   if (num_voxel_total>0)
     mean_disp = disp_sum_total / num_voxel_total;
   else
-    mean_disp = 0;
+    mean_disp = -1;
 
 	return mean_disp;
 
@@ -2137,7 +2137,7 @@ double irtkReconstructionCardiac4D::CalculateWeightedDisplacement()
 
     disp_sum_slice = 0;
     weight_slice = 0;
-    slice_disp = 0;
+    slice_disp = -1;
 
     if (_slice_excluded[inputIndex]==0) {    
   		for (int i = 0; i < _slices[inputIndex].GetX(); i++) {
@@ -2168,7 +2168,7 @@ double irtkReconstructionCardiac4D::CalculateWeightedDisplacement()
   if (weight_total>0)
     mean_disp = disp_sum_total / weight_total;
   else
-    mean_disp = 0;
+    mean_disp = -1;
 
 	return mean_disp;