|
25 | 25 | hdu1 = astropy.io.fits.open(joinpath("data", "gc_2mass_k.fits"))[1]
|
26 | 26 | hdu2 = astropy.io.fits.open(joinpath("data", "gc_msx_e.fits"))[1]
|
27 | 27 |
|
28 |
| - @test isapprox(reproject(imgin, imgout, order = 0)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 0)[1], nans = true, rtol = 1e-4) |
29 |
| - @test isapprox(reproject(imgout, imgin, order = 0)[1]', rp.reproject_interp(hdu1, hdu2.header, order = 0)[1], nans = true, rtol = 1e-3) |
30 |
| - @test isapprox(reproject(imgin, imgout, order = 1)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 1)[1], nans = true, rtol = 1e-4) |
31 |
| - @test isapprox(reproject(imgin, imgout, order = 2)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 2)[1], nans = true, rtol = 1e-4) |
| 28 | + @test isapprox(reproject(imgin, imgout, order = 0)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 0)[1], nans = true, rtol = 1e-7) |
| 29 | + @test isapprox(reproject(imgout, imgin, order = 0)[1]', rp.reproject_interp(hdu1, hdu2.header, order = 0)[1], nans = true, rtol = 1e-6) |
| 30 | + @test isapprox(reproject(imgin, imgout, order = 1)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 1)[1], nans = true, rtol = 1e-7) |
| 31 | + @test isapprox(reproject(imgin, imgout, order = 2)[1]', rp.reproject_interp(hdu2, hdu1.header, order = 2)[1], nans = true, rtol = 6e-2) |
32 | 32 | @test isapprox(reproject(imgin[1], imgout[1], shape_out = (1000,1000))[1]',
|
33 |
| - rp.reproject_interp(hdu2, astropy.wcs.WCS(hdu1.header), shape_out = (1000,1000))[1], nans = true, atol = 1e-4) |
| 33 | + rp.reproject_interp(hdu2, astropy.wcs.WCS(hdu1.header), shape_out = (1000,1000))[1], nans = true, rtol = 1e-7) |
34 | 34 |
|
35 | 35 | wcs = WCSTransform(2; ctype = ["RA---AIR", "DEC--AIR"], radesys = "UNK")
|
36 | 36 | @test_throws ArgumentError reproject(imgin, wcs, shape_out = (100,100))
|
|
0 commit comments