def test_get_circle_cutout_pixels_iris_no_overlap():
header_filename = os.path.join(TESTDATA_DIR, IRIS_3D_CUBE_HEADER)
header = fits.Header.fromtextfile(header_filename)
circle = Circle(20.0, 20.0, 0.1)
test_subject = Transform()
try:
test_subject.get_circle_cutout_pixels(circle, header, 1, 2)
assert False, 'Should raise NoContentError.'
except NoContentError:
assert True
def test_get_circle_cutout_pixels_vlass():
header_filename = os.path.join(TESTDATA_DIR, VLASS_4D_CUBE_HEADER)
header = fits.Header.fromtextfile(header_filename)
circle = Circle(168.34719985367971, 76.18699791158396, 0.010)
test_subject = Transform()
pixels = test_subject.get_circle_cutout_pixels(circle, header, 1, 2)
assert pixels is not None
assert len(pixels) == 4
assert pixels[0] == 2940
assert pixels[1] == 3061
assert pixels[2] == 4193
assert pixels[3] == 4314
def test_get_circle_cutout_pixels_cgps_galactic():
header_filename = os.path.join(TESTDATA_DIR, CGPS_4D_CUBE_HEADER)
header = fits.Header.fromtextfile(header_filename)
circle = Circle(25.0, 60.0, 0.5)
test_subject = Transform()
pixels = test_subject.get_circle_cutout_pixels(circle, header, 1, 2)
# SODA returns [0][350:584,136:370]
assert pixels is not None
assert pixels[0] == 367
assert pixels[1] == 568
assert pixels[2] == 152
assert pixels[3] == 353
def test_circle():
header_filename = os.path.join(TESTDATA_DIR, QUERY_HEADER)
header = fits.Header.fromtextfile(header_filename)
circle = Circle(140.0, 0.0, 0.1)
test_subject = Transform()
pixels = test_subject.get_circle_cutout_pixels(circle, header, 1, 2)
# SODA returns [0][271:279,254:262,*]
assert pixels is not None
assert len(pixels) == 4
assert pixels[0] == 271
assert pixels[1] == 280
assert pixels[2] == 253
assert pixels[3] == 262
def test_get_circle_cutout_pixels_iris_all_overlap():
header_filename = os.path.join(TESTDATA_DIR, IRIS_3D_CUBE_HEADER)
header = fits.Header.fromtextfile(header_filename)
circle = Circle(140.0, 0.0, 10.0)
test_subject = Transform()
pixels = test_subject.get_circle_cutout_pixels(circle, header, 1, 2)
# cutout pixels: -125:676, -143:659
# cutout returning entire image returns empty list
assert pixels is not None
assert len(pixels) == 4
assert pixels[0] == 1
assert pixels[1] == 500
assert pixels[2] == 1
assert pixels[3] == 500
