From 99f434a972afb4db89eafa7f134e9af752571b9c Mon Sep 17 00:00:00 2001 From: hack Date: Fri, 5 Mar 2010 21:16:30 +0000 Subject: Modified 'updatewcs' in 'stwcs' to report the CRPIX value corrected for the additional pixel in the NAXIS. This correction was tested using a the artificial DGEOFILE generation code added to 'reftools.test_small_dgeo'. WJH git-svn-id: http://svn.stsci.edu/svn/ssb/stsci_python/stsci_python/trunk/stwcs@8834 fe389314-cf27-0410-b35b-8c050e845b92 --- updatewcs/dgeo.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) (limited to 'updatewcs/dgeo.py') diff --git a/updatewcs/dgeo.py b/updatewcs/dgeo.py index b430fd6..63d5f52 100644 --- a/updatewcs/dgeo.py +++ b/updatewcs/dgeo.py @@ -251,8 +251,10 @@ class DGEOCorr(object): for i in range(1, naxis+1): si = str(i) kw_val1['NAXIS'+si] = npol_header.get('NAXIS'+si) - kw_val1['CRPIX'+si] = (kw_val1['NAXIS'+si]-1)/2. - kw_val1['CDELT'+si] = float(npol_header.get('ONAXIS'+si))/ ((kw_val1['NAXIS'+si]-1) * binned) + cdelt = (float(npol_header.get('ONAXIS'+si))/ ((kw_val1['NAXIS'+si]-1) * binned)) + kw_val1['CDELT'+si] = cdelt + if cdelt == 0: cdelt = 1.0 + kw_val1['CRPIX'+si] = (kw_val1['NAXIS'+si]-1)/2. + 1/cdelt kw_val1['CRVAL'+si] = (npol_header.get('ONAXIS'+si)/2. + \ sciheader.get('LTV'+si, 0.)) / binned -- cgit