Changeset 7543

Timestamp:

Jun 13, 2006, 12:56:52 PM (20 years ago)

Author:

Paul Price

Message:

Adding check on the gradient.

File:

• trunk/stac/src/stacRejection.c (modified) (1 diff)

trunk/stac/src/stacRejection.c

@@ -144 +144 @@
                                     // Calculate the gradient
                                     grads->data.F32[j] = stacGradient(inputs->data[j], xPix, yPix);
-                                    gradsMask->data.U8[j] = 0;
-                                    numGrads++;
+                                    if (finite(grads->data.F32[j])) {
+                                        gradsMask->data.U8[j] = 0;
+                                        numGrads++;
+                                    } else {
+                                        gradsMask->data.U8[j] = 1;
+                                    }
                                 } else {
                                     gradsMask->data.U8[j] = 1; // Mask this one