summaryrefslogtreecommitdiff
path: root/distortion/utils.py
diff options
context:
space:
mode:
authordencheva <dencheva@stsci.edu>2011-02-04 09:12:23 -0500
committerdencheva <dencheva@stsci.edu>2011-02-04 09:12:23 -0500
commit3c67ebb351a3695fc5207835423b366109cd2894 (patch)
tree5cc6968a363b9123f159b2102de962e058b751ab /distortion/utils.py
parentd6db300e8b08a606d821eaaf16bb7d589e9c16d5 (diff)
downloadstwcs_hcf-3c67ebb351a3695fc5207835423b366109cd2894.tar.gz
Make a copy of the output WCS before using it.
git-svn-id: http://svn.stsci.edu/svn/ssb/stsci_python/stsci_python/trunk/stwcs@11870 fe389314-cf27-0410-b35b-8c050e845b92
Diffstat (limited to 'distortion/utils.py')
-rw-r--r--distortion/utils.py5
1 files changed, 3 insertions, 2 deletions
diff --git a/distortion/utils.py b/distortion/utils.py
index d132eda..ad58f95 100644
--- a/distortion/utils.py
+++ b/distortion/utils.py
@@ -7,7 +7,7 @@ from stwcs import wcsutil
from numpy import sqrt, arctan2
from pytools import fileutil
-def output_wcs(list_of_wcsobj, ref_wcs=None, outwcs=None, undistort=True):
+def output_wcs(list_of_wcsobj, ref_wcs=None, owcs=None, undistort=True):
"""
Create an output WCS.
@@ -31,7 +31,7 @@ def output_wcs(list_of_wcsobj, ref_wcs=None, outwcs=None, undistort=True):
crval1,crval2 = computeFootprintCenter(fra_dec)
crval = np.array([crval1,crval2], dtype=np.float64)
- if outwcs is None:
+ if owcs is None:
if ref_wcs == None:
ref_wcs = list_of_wcsobj[0].deepcopy()
if undistort:
@@ -43,6 +43,7 @@ def output_wcs(list_of_wcsobj, ref_wcs=None, outwcs=None, undistort=True):
outwcs.pscale = sqrt(outwcs.wcs.cd[0,0]**2 + outwcs.wcs.cd[1,0]**2)*3600.
outwcs.orientat = arctan2(outwcs.wcs.cd[0,1],outwcs.wcs.cd[1,1]) * 180./np.pi
else:
+ outwcs = owcs.deepcopy()
outwcs.pscale = sqrt(outwcs.wcs.cd[0,0]**2 + outwcs.wcs.cd[1,0]**2)*3600.
outwcs.orientat = arctan2(outwcs.wcs.cd[0,1],outwcs.wcs.cd[1,1]) * 180./np.pi