import numpy as np
from astropy.io import fits
from astropy.utils.data import get_pkg_data_filename
from matplotlib import pyplot as plt
from regions import Regions

image_file = get_pkg_data_filename('tutorials/FITS-images/HorseHead.fits')
data = fits.getdata(image_file, ext=0, memmap=False)
data = np.rot90(data, k=-1)

fig, ax = plt.subplots(figsize=(8, 8))
ax.imshow(data, cmap='gray', origin='lower')

region_file = get_pkg_data_filename('data/plot_image.reg',
                                    package='regions.io.ds9.tests')
regions = Regions.read(region_file, format='ds9')
for i, region in enumerate(regions):
    region.plot(ax=ax)