IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 10586


Ignore:
Timestamp:
Dec 8, 2006, 2:30:57 PM (20 years ago)
Author:
eugene
Message:

fixed errors in psEllipse functions

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/doc/pslib/psLibADD.tex

    r9785 r10586  
    1 %%% $Id: psLibADD.tex,v 1.94 2006-10-30 21:45:46 eugene Exp $
     1%%% $Id: psLibADD.tex,v 1.95 2006-12-09 00:30:57 eugene Exp $
    22\documentclass[panstarrs]{panstarrs}
    33
     
    11461146\begin{equation}
    11471147\left( \begin{array}{c} x^\prime \\ y^\prime \end{array} \right) =
    1148 \left| \begin{array}{cc} \cos \theta & \sin \theta \\
    1149                         -\sin \theta & \cos \theta
     1148\left| \begin{array}{cc} \cos \theta & -\sin \theta \\
     1149                         \sin \theta & \cos \theta
    11501150\end{array} \right|
    11511151\left( \begin{array}{c} x \\ y \end{array} \right)
     
    11541154(aligned) ellipse.  Applying this rotation to (\ref{aligned-ellipse}) yields:
    11551155\begin{equation}
    1156 z = \frac{x^2 \cos^2 \theta + y^2 \sin^2 \theta + 2 x y \sin \theta \cos \theta}{2\sigma_a^2} +
    1157     \frac{x^2 \sin^2 \theta + y^2 \cos^2 \theta - 2 x y \sin \theta \cos \theta}{2\sigma_b^2}
     1156z = \frac{x^2 \cos^2 \theta + y^2 \sin^2 \theta - 2 x y \sin \theta \cos \theta}{2\sigma_a^2} +
     1157    \frac{x^2 \sin^2 \theta + y^2 \cos^2 \theta + 2 x y \sin \theta \cos \theta}{2\sigma_b^2}
    11581158\end{equation}
    11591159Grouping these terms together, we find:
    11601160\begin{equation}
    11611161z = \frac{x^2}{2}(\sigma_a^{-2} \cos^2 \theta + \sigma_b^{-2}\sin^2 \theta) +
    1162     \frac{y^2}{2}(\sigma_b^{-2} \cos^2 \theta + \sigma_a^{-2}\sin^2 \theta) -
     1162    \frac{y^2}{2}(\sigma_b^{-2} \cos^2 \theta + \sigma_a^{-2}\sin^2 \theta) +
    11631163    \frac{xy}{2} \sin (2 \theta) (\sigma_b^{-2} - \sigma_a^{-2})
    11641164\end{equation}
     
    11921192From these, we may derive the equations for $\sigma_a$, $\sigma_b$, and $\theta$:
    11931193\begin{eqnarray}
    1194 \theta & = & \frac{1}{2} \arg (-2 \sigma_{xy}, f_2) \\
     1194\theta & = & \frac{1}{2} \arg (2 \sigma_{xy}, f_2) \\
    11951195\sigma_a & = & \sqrt{\frac{2}{f_1 - f_3}} \\
    11961196\sigma_b & = & \sqrt{\frac{2}{f_1 + f_3}}
     
    12101210\left|
    12111211\begin{array}{cc}
    1212 +\cos \theta & +\sin \theta \\
    1213 -\sin \theta & +\cos \theta \\
     1212+\cos \theta & -\sin \theta \\
     1213+\sin \theta & +\cos \theta \\
    12141214\end{array} \right|
    12151215\left|
     
    12201220\left|
    12211221\begin{array}{cc}
    1222 +\cos \theta & -\sin \theta \\
    1223 +\sin \theta & +\cos \theta \\
     1222+\cos \theta & +\sin \theta \\
     1223-\sin \theta & +\cos \theta \\
    12241224\end{array} \right|
    12251225\end{equation}
     
    12281228m_{x,x} & = & \sigma_a^{2} \cos^2 \theta + \sigma_b^{2}\sin^2 \theta \\
    12291229m_{y,y} & = & \sigma_b^{2} \cos^2 \theta + \sigma_a^{2}\sin^2 \theta \\
    1230 m_{x,y} & = & -\frac{1}{2} \sin (2 \theta) (\sigma_a^2 - \sigma_b^2)
     1230m_{x,y} & = & \frac{1}{2} \sin (2 \theta) (\sigma_a^2 - \sigma_b^2)
    12311231\end{eqnarray}
    12321232Using the double-angle relationships, these become:
     
    12341234m_{x,x} & = & \frac{1}{2}(\sigma_a^{2} + \sigma_b^{2}) + \frac{1}{2}(\sigma_a^{2} - \sigma_b^{2}) \cos (2 \theta) \\
    12351235m_{y,y} & = & \frac{1}{2}(\sigma_a^{2} + \sigma_b^{2}) - \frac{1}{2}(\sigma_a^{2} - \sigma_b^{2}) \cos (2 \theta) \\
    1236 m_{x,y} & = & -\frac{1}{2} \sin (2 \theta) (\sigma_a^{-2} - \sigma_b^{-2})
     1236m_{x,y} & = & \frac{1}{2} \sin (2 \theta) (\sigma_a^{2} - \sigma_b^{2})
    12371237\end{eqnarray}
    12381238These three formulae define the second moments in terms of $\sigma_a$, $\sigma_b$, and $\theta$.
     
    12421242g_1 = m_{x,x} + m_{y,y}          & = & \sigma_a^{2} + \sigma_b^{2} \\
    12431243g_2 = m_{x,x} - m_{y,y}          & = & (\sigma_a^{2} - \sigma_b^{2}) \cos (2 \theta) \\
    1244 g_3 = \sqrt{f_2^2 + 4 m_{x,y}^2} & = & \sigma_a^{2} - \sigma_b^{-2}
     1244g_3 = \sqrt{g_2^2 + 4 m_{x,y}^2} & = & \sigma_a^{2} - \sigma_b^{2}
    12451245\end{eqnarray}
    12461246From these, we may derive the equations for $\sigma_a$, $\sigma_b$, and $\theta$:
    12471247\begin{eqnarray}
    1248 \theta   & = & \frac{1}{2} \arg (-2 m_{x,y}, g_2) \\
    1249 \sigma_a & = & \sqrt{\frac{g_1 - g_3}{2}} \\
    1250 \sigma_b & = & \sqrt{\frac{g_1 + g_3}{2}}
     1248\theta   & = & \frac{1}{2} \arg (2 m_{x,y}, g_2) \\
     1249\sigma_a & = & \sqrt{\frac{g_1 + g_3}{2}} \\
     1250\sigma_b & = & \sqrt{\frac{g_1 - g_3}{2}}
    12511251\end{eqnarray}
    12521252
Note: See TracChangeset for help on using the changeset viewer.