Changeset 38599 for branches/eam_branches/ipp-20150625/Ohana/src/relastro/src/ImageOps.c

Timestamp:

Jul 17, 2015, 8:39:34 PM (11 years ago)

Author:

eugene

Message:

working on relastro

File:

• branches/eam_branches/ipp-20150625/Ohana/src/relastro/src/ImageOps.c (modified) (1 diff)

branches/eam_branches/ipp-20150625/Ohana/src/relastro/src/ImageOps.c

@@ -748 +748 @@
     n = measure[0].averef;
 
+    // apply proper-motion from average position to measure epoch:
+    float dTime = (measure[k].t - catalog[c].average[n].Tmean) / (86400*365.25) ; // time relative to Tmean in years
+
     /* apply the current image transformation or use the current value of R+dR, D+dD? */
     ref[i].R = catalog[c].average[n].R;
     ref[i].D = catalog[c].average[n].D;
+    
+    // XXX do this in a better way?
+    ref[i].R += dTime * catalog[c].average[n].uR / 3600.0 / cos(ref[i].D*RAD_DEG);
+    ref[i].D += dTime * catalog[c].average[n].uD / 3600.0;
 
     // if we are correcting for the Galaxy Motion Model, we assume the mean R,D is at the J2000 epoch position