Changeset 21282 for branches/pap_branch_20090128/pswarp/src/pswarpLoop.c

Timestamp:

Feb 4, 2009, 10:24:33 AM (17 years ago)

Author:

Paul Price

Message:

Calculate covariance matrix for output.

File:

• branches/pap_branch_20090128/pswarp/src/pswarpLoop.c (modified) (1 diff)

branches/pap_branch_20090128/pswarp/src/pswarpLoop.c

@@ -273 +273 @@
     }
 
+    // Set covariance matrix for output
+    {
+        psList *covariances = psMetadataLookupPtr(&mdok, output->analysis,
+                                                  PSWARP_ANALYSIS_COVARIANCES); // Covariance matrices
+        psAssert(covariances, "Should be there");
+        psArray *covars = psListToArray(covariances); // Array of covariance matrices
+        output->covariance = psImageCovarianceAverage(covars);
+        psFree(covars);
+        psMetadataRemoveKey(output->analysis, PSWARP_ANALYSIS_COVARIANCES);
+    }
+
     if (!pmConceptsAverageCells(outCell, cells, NULL, NULL, false)) {
         psError(PS_ERR_UNKNOWN, false, "Unable to average cell concepts.");