diff --git a/package/lib/fits.py b/package/lib/fits.py
index 67a62d9..cbacc85 100755
--- a/package/lib/fits.py
+++ b/package/lib/fits.py
@@ -16,7 +16,7 @@ from astropy.io import fits
 from astropy.wcs import WCS
 
 from .convex_hull import clean_ROI
-from .utils import wcs_PA
+from .utils import wcs_PA, princ_angle
 
 
 def get_obs_data(infiles, data_folder="", compute_flux=False):
diff --git a/package/lib/utils.py b/package/lib/utils.py
index ee9fc5b..9e36dbc 100755
--- a/package/lib/utils.py
+++ b/package/lib/utils.py
@@ -177,4 +177,4 @@ def wcs_PA(PC21, PC22):
         orient = np.arccos(PC22) * 180.0 / np.pi
     elif (abs(PC21) < abs(PC22)) and (PC22 < 0):
         orient = -np.arccos(PC22) * 180.0 / np.pi
-    return orient
+    return princ_angle(orient)