Changeset 35541

Timestamp:

May 9, 2013, 10:54:51 AM (13 years ago)

Author:

eugene

Message:

in multi-input mode, some output readouts may not have valid inputs on the first pass

File:

• branches/eam_branches/ipp-20130419/pswarp/src/pswarpUpdateMetadata.c (modified) (1 diff)

branches/eam_branches/ipp-20130419/pswarp/src/pswarpUpdateMetadata.c

@@ -48 +48 @@
                     psAssert(covariances, "Should be there");
                     psArray *covars = psListToArray(covariances); // Array of covariance matrices
-                    psKernel *covar = psImageCovarianceAverage(covars);
-                    psFree(covars);
-                    psMetadataRemoveKey(readout->analysis, PSWARP_ANALYSIS_COVARIANCES);
+                    if (covars->n) {
+                      psKernel *covar = psImageCovarianceAverage(covars);
+                      psMetadataRemoveKey(readout->analysis, PSWARP_ANALYSIS_COVARIANCES);
 
-                    // Correct covariance matrix scale for the mean (square root of the) Jacobian
-                    double jacobian = psMetadataLookupF64(NULL, readout->analysis, PSWARP_ANALYSIS_JACOBIAN); // Jacobian
-                    int goodPixels = psMetadataLookupS32(NULL, readout->analysis, PSWARP_ANALYSIS_GOODPIX);   // Good pixels
-                    jacobian /= goodPixels;
-                    readout->covariance = psImageCovarianceScale(covar, jacobian);
-                    psFree(covar);
+                      // Correct covariance matrix scale for the mean (square root of the) Jacobian
+                      double jacobian = psMetadataLookupF64(NULL, readout->analysis, PSWARP_ANALYSIS_JACOBIAN); // Jacobian
+                      int goodPixels = psMetadataLookupS32(NULL, readout->analysis, PSWARP_ANALYSIS_GOODPIX);   // Good pixels
+                      jacobian /= goodPixels;
+                      readout->covariance = psImageCovarianceScale(covar, jacobian);
+                      psFree(covar);
 
-                    if (readout->variance) {
+                      if (readout->variance) {
                         psImageCovarianceTransfer(readout->variance, readout->covariance);
+                      }
+
+                      psFree(covars);
                     }
                 }