Create plot_region.py

The original code I used to plot leonidas' region files over the fits files (before I got ds9 to work). Saving this version here.