You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
{{ message }}
Repository navigation
New WCS TPV support doesn't work with DECAM images #3904
Sorry to say this, @mdboom - but I just attempted the procedure I described in my comment in #3845, and it still doesn't seem to work even after #3878 . I've copied my procedure below, although not a critical change: it seem to work with neither all_pix2worldnorwcs_pix2world now...
Hit "search", and then find the first result the search returns that has "Product"="image". Hit "retreive" on that row and it will download a DECam image.
You should see this with that image:
>>> f = fits.open('c4d_130831_235507_ooi_g_d3.fits.fz')
>>> w = wcs.WCS(f[1].header)
>>> ra, dec = w.all_pix2world([0,10, 15],[0,10, 20],0)
>>> ra
array([ 301.2065433, 301.2065433, 301.2065433])
>>> dec
array([-50.93399553, -50.93399553, -50.93399553])
which is clearly wrong: they should be different ra/dec.
4. Also, I see this:
>>> ra, dec = w.wcs_pix2world([0,10, 15],[0,10, 20],0)
>>> ra
array([ 301.2065433, 301.2065433, 301.2065433])
>>> dec
array([-50.93399553, -50.93399553, -50.93399553])
which is the same incorrect answer, as opposed to when I tried this in #3845, where wcs_pix2worlddid work.
I think this is (another) instance of wcscopy being broken. The polynomial distortion coefficients are somehow ending up as all zeros (which is why the results are always the same). Looking into it now.
Activity
mdboom commentedon Jun 30, 2015
I think this is (another) instance of wcscopy being broken. The polynomial distortion coefficients are somehow ending up as all zeros (which is why the results are always the same). Looking into it now.
mdboom commentedon Jun 30, 2015
It is an issue with copying the wcs struct, but it appears to have already been fixed in wcslib 5.7. I have a PR almost ready to go on that.
Add test for astropy#3904
eteq commentedon Jul 1, 2015
Looks like #3905 fixes this (see my comment there)
Add test for astropy#3904
bsipocz commentedon Jun 28, 2017
#3905 missed to auto-close this, closing now.