diff options
| author | hack <hack@stsci.edu> | 2010-03-05 16:16:30 -0500 |
|---|---|---|
| committer | hack <hack@stsci.edu> | 2010-03-05 16:16:30 -0500 |
| commit | 99f434a972afb4db89eafa7f134e9af752571b9c (patch) | |
| tree | cde701f322fe096fd695c5d866d96dfab2fa6bd6 /updatewcs/dgeo.py | |
| parent | fb55bf4f24017b1669f9da14d02a4b6fbca3ea76 (diff) | |
| download | stwcs_hcf-99f434a972afb4db89eafa7f134e9af752571b9c.tar.gz | |
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
Diffstat (limited to 'updatewcs/dgeo.py')
| -rw-r--r-- | updatewcs/dgeo.py | 6 |
1 files changed, 4 insertions, 2 deletions
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 |