Changeset 23634 for trunk/ippTools/src/difftool.c

Timestamp:

Mar 31, 2009, 12:51:29 PM (17 years ago)

Author:

Paul Price

Message:

Adding limit on great circle distance, using Haversine formula

File:

• trunk/ippTools/src/difftool.c (modified) (1 diff)

trunk/ippTools/src/difftool.c

@@ -1274 +1274 @@
                    "ABS(TIME_TO_SEC(TIMEDIFF(inputRawExp.dateobs, templateRawExp.dateobs)))", "<=");
     PXOPT_COPY_F32(config->args, selectWhere, "-rotdiff", "ABS(inputRawExp.posang - templateRawExp.posang)", "<=");
+    // Haversine formula for great circle distance
+    PXOPT_COPY_F32(config->args, selectWhere, "-distance",
+                   "DEGREES(2*ASIN(SQRT(POW(SIN(inputRawExp.decl - templateRawExp.decl),2) + COS(inputRawExp.decl)*COS(templateRawExp.decl)*POW(SIN(inputRawExp.ra - templateRawExp.ra),2))))", "<=");
 
     // Restrictions for inserting skycells