Changeset 3610 for trunk/stac/src/stacTransform.c

Timestamp:

Mar 31, 2005, 5:07:12 PM (21 years ago)

Author:

Paul Price

Message:

Working very nicely now.

File:

• trunk/stac/src/stacTransform.c (modified) (4 diffs)

trunk/stac/src/stacTransform.c

@@ -113 +113 @@
                    const psArray *masks, // Masks of input images
                    const psImage *region, // Region of interest for transformation
+                   const psVector *scales, // Relative scales
+                   const psVector *offsets, // Relative offsets
                    const stacConfig *config // Configuration
     )
@@ -123 +125 @@
     assert(images->n == maps->n);
     assert(!errors || (images->n == errors->n));
+    assert(!scales || scales->n == images->n);
+    assert(!offsets || offsets->n == images->n);
+    assert(!scales || scales->type.type == PS_TYPE_F32);
+    assert(!offsets || offsets->type.type == PS_TYPE_F32);
 
     // Allocate the output images if required, otherwise check the number
@@ -170 +176 @@
         psImage *outImage = (*outputs)->data[n]; // The output image
         psImage *outError = (*outErrors)->data[n]; // The output error image
+        float offset = 0.0;             // Relative offset
+        float scale = 1.0;              // Relative scale
+        if (offsets) {
+            offset = offsets->data.F32[n];
+        }
+        if (scales) {
+            scale = scales->data.F32[n];
+        }
 
 #if 0
@@ -206 +220 @@
                                                                                             detector->y,
                                                                                             mask, 1, 0.0);
+
+                    outImage->data.F32[y][x] = (outImage->data.F32[y][x] - offset) / scale;
+                    outImage->data.F32[y][x] = outImage->data.F32[y][x] / SQUARE(scale);
+
                 } // Pixels of interest