Changeset 5353 for trunk/archive/modules/src/pmAstrom.c

Timestamp:

Oct 17, 2005, 9:59:28 AM (21 years ago)

Author:

eugene

Message:

added pmAstromMatchedListFit

File:

• trunk/archive/modules/src/pmAstrom.c (modified) (1 diff)

trunk/archive/modules/src/pmAstrom.c

@@ -223 +223 @@
     return (match);
 }
+
+pmAstromMatchedListFit (psArray *st1, psArray *st2, psArray *match, pmAstromOpt *opt) {
+
+  psVector *X = psVectorAlloc (match->n);
+  psVector *Y = psVectorAlloc (match->n);
+  psVector *x = psVectorAlloc (match->n);
+  psVector *y = psVectorAlloc (match->n);
+  
+  // take the matched stars, first fit 
+  for (i = 0; i < match->n; i++) {
+
+    pair = match->data[i];
+    ob1 = st1->data[pair->i1];
+    ob2 = st2->data[pair->i2];
+
+    X->data.F64[i] = ob1->P;
+    Y->data.F64[i] = ob1->Q;
+
+    x->data.F64[i] = ob2->P;
+    y->data.F64[i] = ob2->Q;
+
+    // use one or the other...
+    psPolynomial2D *xt = psVectorFitPolynomial2D (NULL, stats, NULL, 0, X, NULL, x, y)
+    psPolynomial2D *xt = psVectorClipFitPolynomial2D (NULL, stats, NULL, 0, X, NULL, x, y)
+
+    psPolynomial2D *yt = psVectorFitPolynomial2D (NULL, stats, NULL, 0, Y, NULL, x, y)
+    psPolynomial2D *yt = psVectorClipFitPolynomial2D (NULL, stats, NULL, 0, Y, NULL, x, y)