From c5a9470d7db659918d4fbd8530b8981eddad04e7 Mon Sep 17 00:00:00 2001 From: dencheva Date: Wed, 10 Aug 2011 16:15:01 +0000 Subject: (past) last minute check in - corrected aproblem when wcsname was overwritten by an old value git-svn-id: http://svn.stsci.edu/svn/ssb/stsci_python/stsci_python/trunk/stwcs@13565 fe389314-cf27-0410-b35b-8c050e845b92 --- lib/stwcs/wcsutil/altwcs.py | 11 +++++++---- 1 file changed, 7 insertions(+), 4 deletions(-) (limited to 'lib/stwcs/wcsutil/altwcs.py') diff --git a/lib/stwcs/wcsutil/altwcs.py b/lib/stwcs/wcsutil/altwcs.py index dd5f1e9..1f6dc01 100644 --- a/lib/stwcs/wcsutil/altwcs.py +++ b/lib/stwcs/wcsutil/altwcs.py @@ -102,7 +102,7 @@ def archiveWCS(fname, ext=None, wcskey=" ", wcsname=" ", reusekey=False): else: wkey = wcskey wname = wcsname - + for e in ext: try: w = pywcs.WCS(f[e].header, fobj=f) @@ -113,14 +113,17 @@ def archiveWCS(fname, ext=None, wcskey=" ", wcsname=" ", reusekey=False): hwcs = w.to_header() wcsnamekey = 'WCSNAME' + wkey f[e].header.update(key=wcsnamekey, value=wname) + if w.wcs.has_cd(): pc2cd(hwcs) + try: + old_wcsname=hwcs.pop('WCSNAME') + except: + pass for k in hwcs.keys(): + key = k[:7] + wkey f[e].header.update(key=key, value=hwcs[k]) - #norient = np.rad2deg(np.arctan2(hwcs['CD1_2'],hwcs['CD2_2'])) - #okey = 'ORIENT%s' % wkey - #f[e].header.update(key=okey, value=norient) closefobj(fname, f) def restoreWCS(f, ext=None, fromext=None, toext=None, wcskey=" ", wcsname=" "): -- cgit