Changeset 15688 for trunk/ippScripts/scripts/magic_tree.pl

Timestamp:

Nov 26, 2007, 2:36:18 PM (18 years ago)

Author:

Paul Price

Message:

Fixing signs of rotation equation so that we can correctly parse the
example in WCS paper 2, example 1 (Tables 5 and 6).

File:

• trunk/ippScripts/scripts/magic_tree.pl (modified) (1 diff)

trunk/ippScripts/scripts/magic_tree.pl

@@ -164 +164 @@
     $crval2 = deg2rad($crval2);
     my $ra = $crval1 + atan2(cos($theta) * sin($phi),
-                             sin($theta) * cos($crval2) - cos($theta) * sin($crval2) * cos($phi));
-    my $dec = asin(sin($theta) * sin($crval2) + cos($theta) * cos($crval2) * cos($phi));
-    
+                             sin($theta) * cos($crval2) + cos($theta) * sin($crval2) * cos($phi));
+    my $dec = asin(sin($theta) * sin($crval2) - cos($theta) * cos($crval2) * cos($phi));
+
     # Rotate to boresight
     my $phi_new = atan2(cos($dec) * sin($ra - $ra0), 
-                        sin($dec) * cos($dec0) - cos($dec) * sin($dec0) * cos($ra - $ra0));
-    my $theta_new = asin(sin($dec) * sin($dec0) + cos($dec) * cos($dec0) * cos($ra - $ra0));
+                        sin($dec) * cos($dec0) + cos($dec) * sin($dec0) * cos($ra - $ra0));
+    my $theta_new = asin(sin($dec) * sin($dec0) - cos($dec) * cos($dec0) * cos($ra - $ra0));
 
     # Project