correct second moment transformation

Author: PFLeget

python/lsst/pipe/tasks/finalizeCharacterization.py

@@ -637,12 +637,12 @@ def _compute_sky_moments(self, wcs, x, y, ixx, iyy, ixy):
             cd22 = cd_matrix[1, 1]
 
             # Transform moments: M_sky = CD * M_pixel * CD^T
-            # Iuu = CD11*(Ixx*CD11 + Ixy*CD21) + CD12*(Ixy*CD11 + Iyy*CD21)
-            iuu[i] = (cd11 * (ixx[i] * cd11 + ixy[i] * cd21)
-                      + cd12 * (ixy[i] * cd11 + iyy[i] * cd21))
-            # Ivv = CD21*(Ixx*CD12 + Ixy*CD22) + CD22*(Ixy*CD12 + Iyy*CD22)
-            ivv[i] = (cd21 * (ixx[i] * cd12 + ixy[i] * cd22)
-                      + cd22 * (ixy[i] * cd12 + iyy[i] * cd22))
+            # Iuu = CD11*(Ixx*CD11 + Ixy*CD12) + CD12*(Ixy*CD11 + Iyy*CD12)
+            iuu[i] = (cd11 * (ixx[i] * cd11 + ixy[i] * cd12)
+                      + cd12 * (ixy[i] * cd11 + iyy[i] * cd12))
+            # Ivv = CD21*(Ixx*CD21 + Ixy*CD22) + CD22*(Ixy*CD21 + Iyy*CD22)
+            ivv[i] = (cd21 * (ixx[i] * cd21 + ixy[i] * cd22)
+                      + cd22 * (ixy[i] * cd21 + iyy[i] * cd22))
             # Iuv = (CD11*Ixx + CD12*Ixy)*CD21 + (CD11*Ixy + CD12*Iyy)*CD22
             iuv[i] = ((cd11 * ixx[i] + cd12 * ixy[i]) * cd21
                       + (cd11 * ixy[i] + cd12 * iyy[i]) * cd22)