Changeset 10586
- Timestamp:
- Dec 8, 2006, 2:30:57 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/doc/pslib/psLibADD.tex (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/doc/pslib/psLibADD.tex
r9785 r10586 1 %%% $Id: psLibADD.tex,v 1.9 4 2006-10-30 21:45:46eugene Exp $1 %%% $Id: psLibADD.tex,v 1.95 2006-12-09 00:30:57 eugene Exp $ 2 2 \documentclass[panstarrs]{panstarrs} 3 3 … … 1146 1146 \begin{equation} 1147 1147 \left( \begin{array}{c} x^\prime \\ y^\prime \end{array} \right) = 1148 \left| \begin{array}{cc} \cos \theta & \sin \theta \\1149 -\sin \theta & \cos \theta1148 \left| \begin{array}{cc} \cos \theta & -\sin \theta \\ 1149 \sin \theta & \cos \theta 1150 1150 \end{array} \right| 1151 1151 \left( \begin{array}{c} x \\ y \end{array} \right) … … 1154 1154 (aligned) ellipse. Applying this rotation to (\ref{aligned-ellipse}) yields: 1155 1155 \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}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} 1158 1158 \end{equation} 1159 1159 Grouping these terms together, we find: 1160 1160 \begin{equation} 1161 1161 z = \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) + 1163 1163 \frac{xy}{2} \sin (2 \theta) (\sigma_b^{-2} - \sigma_a^{-2}) 1164 1164 \end{equation} … … 1192 1192 From these, we may derive the equations for $\sigma_a$, $\sigma_b$, and $\theta$: 1193 1193 \begin{eqnarray} 1194 \theta & = & \frac{1}{2} \arg ( -2 \sigma_{xy}, f_2) \\1194 \theta & = & \frac{1}{2} \arg (2 \sigma_{xy}, f_2) \\ 1195 1195 \sigma_a & = & \sqrt{\frac{2}{f_1 - f_3}} \\ 1196 1196 \sigma_b & = & \sqrt{\frac{2}{f_1 + f_3}} … … 1210 1210 \left| 1211 1211 \begin{array}{cc} 1212 +\cos \theta & +\sin \theta \\1213 -\sin \theta & +\cos \theta \\1212 +\cos \theta & -\sin \theta \\ 1213 +\sin \theta & +\cos \theta \\ 1214 1214 \end{array} \right| 1215 1215 \left| … … 1220 1220 \left| 1221 1221 \begin{array}{cc} 1222 +\cos \theta & -\sin \theta \\1223 +\sin \theta & +\cos \theta \\1222 +\cos \theta & +\sin \theta \\ 1223 -\sin \theta & +\cos \theta \\ 1224 1224 \end{array} \right| 1225 1225 \end{equation} … … 1228 1228 m_{x,x} & = & \sigma_a^{2} \cos^2 \theta + \sigma_b^{2}\sin^2 \theta \\ 1229 1229 m_{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)1230 m_{x,y} & = & \frac{1}{2} \sin (2 \theta) (\sigma_a^2 - \sigma_b^2) 1231 1231 \end{eqnarray} 1232 1232 Using the double-angle relationships, these become: … … 1234 1234 m_{x,x} & = & \frac{1}{2}(\sigma_a^{2} + \sigma_b^{2}) + \frac{1}{2}(\sigma_a^{2} - \sigma_b^{2}) \cos (2 \theta) \\ 1235 1235 m_{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})1236 m_{x,y} & = & \frac{1}{2} \sin (2 \theta) (\sigma_a^{2} - \sigma_b^{2}) 1237 1237 \end{eqnarray} 1238 1238 These three formulae define the second moments in terms of $\sigma_a$, $\sigma_b$, and $\theta$. … … 1242 1242 g_1 = m_{x,x} + m_{y,y} & = & \sigma_a^{2} + \sigma_b^{2} \\ 1243 1243 g_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}1244 g_3 = \sqrt{g_2^2 + 4 m_{x,y}^2} & = & \sigma_a^{2} - \sigma_b^{2} 1245 1245 \end{eqnarray} 1246 1246 From these, we may derive the equations for $\sigma_a$, $\sigma_b$, and $\theta$: 1247 1247 \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}} 1251 1251 \end{eqnarray} 1252 1252
Note:
See TracChangeset
for help on using the changeset viewer.
