diff options
Diffstat (limited to 'updatewcs')
| -rw-r--r-- | updatewcs/corrections.py | 4 | ||||
| -rw-r--r-- | updatewcs/makewcs.py | 8 |
2 files changed, 6 insertions, 6 deletions
diff --git a/updatewcs/corrections.py b/updatewcs/corrections.py index fa0a290..3d504ca 100644 --- a/updatewcs/corrections.py +++ b/updatewcs/corrections.py @@ -168,8 +168,8 @@ class CompSIP(object): if n >= m and n>=2: idcval = np.array([[cx[n,m]],[cy[n,m]]]) sipval = np.dot(imatr, idcval) - akeys1[m,n-m] = sipval[0] - bkeys1[m,n-m] = sipval[1] + akeys1[m,n-m] = sipval[0] * ext_wcs.binned + bkeys1[m,n-m] = sipval[1] * ext_wcs.binned Akey="A_%d_%d" % (m,n-m) Bkey="B_%d_%d" % (m,n-m) kw2update[Akey] = sipval[0,0] diff --git a/updatewcs/makewcs.py b/updatewcs/makewcs.py index b455640..77c8aa9 100644 --- a/updatewcs/makewcs.py +++ b/updatewcs/makewcs.py @@ -82,9 +82,9 @@ class MakeWCS(object): off = sqrt((v2-v2ref)**2 + (v3-v3ref)**2)/(R_scale*3600.0) if v3 == v3ref: - theta=0.0 + theta=0.0 else: - theta = atan2(ext_wcs.parity[0][0]*(v2-v2ref), ext_wcs.parity[1][1]*(v3-v3ref)) + theta = atan2(ext_wcs.parity[0][0]*(v2-v2ref), ext_wcs.parity[1][1]*(v3-v3ref)) if ref_wcs.idcmodel.refpix['THETA']: theta += ref_wcs.idcmodel.refpix['THETA']*pi/180.0 @@ -106,9 +106,9 @@ class MakeWCS(object): # Account for subarray offset # Angle of chip relative to chip if ext_wcs.idcmodel.refpix['THETA']: - dtheta = ext_wcs.idcmodel.refpix['THETA'] - ref_wcs.idcmodel.refpix['THETA'] + dtheta = ext_wcs.idcmodel.refpix['THETA'] - ref_wcs.idcmodel.refpix['THETA'] else: - dtheta = 0.0 + dtheta = 0.0 rrmat = fileutil.buildRotMatrix(dtheta) # Rotate the vectors |