Index: /branches/eam_rel8_b0/psLib/src/astro/psCoord.c
===================================================================
--- /branches/eam_rel8_b0/psLib/src/astro/psCoord.c	(revision 5546)
+++ /branches/eam_rel8_b0/psLib/src/astro/psCoord.c	(revision 5547)
@@ -10,6 +10,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.88 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-10-12 21:02:20 $
+*  @version $Revision: 1.88.4.1 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-11-18 20:22:35 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -215,6 +215,6 @@
 
     psPlaneTransform *pt = psAlloc(sizeof(psPlaneTransform));
-    pt->x = psPolynomial2DAlloc(n1-1, n2-1, PS_POLYNOMIAL_ORD);
-    pt->y = psPolynomial2DAlloc(n1-1, n2-1, PS_POLYNOMIAL_ORD);
+    pt->x = psPolynomial2DAlloc(n1, n2, PS_POLYNOMIAL_ORD);
+    pt->y = psPolynomial2DAlloc(n1, n2, PS_POLYNOMIAL_ORD);
 
     psMemSetDeallocator(pt, (psFreeFunc) planeTransformFree);
@@ -286,6 +286,6 @@
 
     psPlaneDistort *pt = psAlloc(sizeof(psPlaneDistort));
-    pt->x = psPolynomial4DAlloc(n1-1, n2-1, n3-1, n4-1, PS_POLYNOMIAL_ORD);
-    pt->y = psPolynomial4DAlloc(n1-1, n2-1, n3-1, n4-1, PS_POLYNOMIAL_ORD);
+    pt->x = psPolynomial4DAlloc(n1, n2, n3, n4, PS_POLYNOMIAL_ORD);
+    pt->y = psPolynomial4DAlloc(n1, n2, n3, n4, PS_POLYNOMIAL_ORD);
 
     psMemSetDeallocator(pt, (psFreeFunc) planeDistortFree);
@@ -428,6 +428,6 @@
 
     // Apply plate scales
-    out->x *= projection->Xs;
-    out->y *= projection->Ys;
+    out->x /= projection->Xs;
+    out->y /= projection->Ys;
 
     // Return output
@@ -448,6 +448,6 @@
 
     // Remove plate scales
-    psF64  x = coord->x/projection->Xs;
-    psF64  y = coord->y/projection->Ys;
+    psF64  x = coord->x*projection->Xs;
+    psF64  y = coord->y*projection->Ys;
 
     // Perform inverse projection