diff options
| author | dencheva <dencheva@stsci.edu> | 2009-02-11 12:40:42 -0500 |
|---|---|---|
| committer | dencheva <dencheva@stsci.edu> | 2009-02-11 12:40:42 -0500 |
| commit | 48391819c2284d83a2440b3ba7173a2e4c18354b (patch) | |
| tree | 5a51643851fe1c59921dc0c32c41c69150e24a94 | |
| parent | 03a2c5643a4c3c10ed7cb05a44facd150bdfc7a5 (diff) | |
| download | stwcs_hcf-48391819c2284d83a2440b3ba7173a2e4c18354b.tar.gz | |
Corrections to origin parameter in s2p
git-svn-id: http://svn.stsci.edu/svn/ssb/stsci_python/development/trunk/updatewcs@7540 fe389314-cf27-0410-b35b-8c050e845b92
| -rw-r--r-- | distortion/utils.py | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/distortion/utils.py b/distortion/utils.py index 26f8f6e..d2bb0a9 100644 --- a/distortion/utils.py +++ b/distortion/utils.py @@ -26,16 +26,16 @@ def output_wcs(list_of_wcsobj, ref_wcs=None, outwcs=None): 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 - tanpix = outwcs.wcs.s2p(fra_dec)['pixcrd'] + tanpix = outwcs.wcs.s2p(fra_dec, 1)['pixcrd'] outwcs.naxis1 = int(np.ceil(tanpix[:,0].max() - tanpix[:,0].min())) outwcs.naxis2 = int(np.ceil(tanpix[:,1].max() - tanpix[:,1].min())) crpix = np.array([outwcs.naxis1/2., outwcs.naxis2/2.]) outwcs.wcs.crpix = crpix - tanpix = outwcs.wcs.s2p(fra_dec)['pixcrd'] + tanpix = outwcs.wcs.s2p(fra_dec, 1)['pixcrd'] newcrpix = np.array([crpix[0]+np.ceil(tanpix[:,0].min()), crpix[1]+np.ceil(tanpix[:,1].min())]) - newcrval = outwcs.wcs.p2s([newcrpix])['world'][0] + newcrval = outwcs.wcs.p2s([newcrpix], 1)['world'][0] outwcs.wcs.crval = newcrval return outwcs |