diff options
Diffstat (limited to 'lib/stwcs/distortion')
| -rw-r--r-- | lib/stwcs/distortion/mutil.py | 2 | ||||
| -rw-r--r-- | lib/stwcs/distortion/utils.py | 8 |
2 files changed, 5 insertions, 5 deletions
diff --git a/lib/stwcs/distortion/mutil.py b/lib/stwcs/distortion/mutil.py index 8f80b67..4965bd1 100644 --- a/lib/stwcs/distortion/mutil.py +++ b/lib/stwcs/distortion/mutil.py @@ -424,7 +424,7 @@ def readWCSCoeffs(header): refpix['PSCALE'] = _scale refpix['V2REF'] = 0. refpix['V3REF'] = 0. - refpix['THETA'] = RADTODEG(_theta) + refpix['THETA'] = np.rad2deg(_theta) refpix['XDELTA'] = 0.0 refpix['YDELTA'] = 0.0 refpix['DEFAULT_SCALE'] = yes diff --git a/lib/stwcs/distortion/utils.py b/lib/stwcs/distortion/utils.py index 7f937f7..1d395cb 100644 --- a/lib/stwcs/distortion/utils.py +++ b/lib/stwcs/distortion/utils.py @@ -72,15 +72,15 @@ def computeFootprintCenter(edges): This algorithm should be more robust against discontinuities at the poles. """ - alpha = fileutil.DEGTORAD(edges[:,0]) - dec = fileutil.DEGTORAD(edges[:,1]) + alpha = np.deg2rad(edges[:,0]) + dec = np.deg2rad(edges[:,1]) xmean = np.mean(np.cos(dec)*np.cos(alpha)) ymean = np.mean(np.cos(dec)*np.sin(alpha)) zmean = np.mean(np.sin(dec)) - crval1 = fileutil.RADTODEG(np.arctan2(ymean,xmean))%360.0 - crval2 = fileutil.RADTODEG(np.arctan2(zmean,np.sqrt(xmean*xmean+ymean*ymean))) + crval1 = np.rad2deg(np.arctan2(ymean,xmean))%360.0 + crval2 = np.rad2deg(np.arctan2(zmean,np.sqrt(xmean*xmean+ymean*ymean))) return crval1,crval2 |