def test_interval_polygon():
header_filename = os.path.join(TESTDATA_DIR, QUERY_HEADER)
header = fits.Header.fromtextfile(header_filename)
# CW winding direction
polygon = Polygon([[139.9, -0.1], [140.1, -0.1],
[140.1, 0.1], [139.9, 0.1]])
test_subject = Transform()
try:
test_subject.get_polygon_cutout_pixels(polygon, header, 1, 2)
assert False
except ValueError:
assert True
def test_get_polygon_cutout_pixels_vlass():
header_filename = os.path.join(TESTDATA_DIR, VLASS_4D_CUBE_HEADER)
header = fits.Header.fromtextfile(header_filename)
polygon = Polygon([[168.34, 76.18], [168.34, 76.19], [168.35, 76.19]])
test_subject = Transform()
pixels = test_subject.get_polygon_cutout_pixels(polygon, header, 1, 2)
# SODA returns cutout=[0][2997:3011,4211:4272,*,*]
assert pixels is not None
assert len(pixels) == 4
assert pixels[0] == 2996
assert pixels[1] == 3012
assert pixels[2] == 4211
assert pixels[3] == 4272
def test_polygon():
header_filename = os.path.join(TESTDATA_DIR, QUERY_HEADER)
header = fits.Header.fromtextfile(header_filename)
# CCW winding direction
polygon = Polygon([[139.9, 0.1], [140.1, 0.1],
[140.1, -0.1], [139.9, -0.1]])
test_subject = Transform()
pixels = test_subject.get_polygon_cutout_pixels(polygon, 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] == 263
