diff --git a/package/lib/fits.py b/package/lib/fits.py
index feb828a..4ea9a47 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):
@@ -72,13 +72,13 @@ def get_obs_data(infiles, data_folder="", compute_flux=False):
                 for key in keys:
                     header.remove(key, ignore_missing=True)
                 new_cdelt = np.linalg.eigvals(wcs.wcs.cd)
-                new_cdelt.sort()
+                # new_cdelt.sort()
             new_wcs.wcs.pc = wcs.wcs.cd.dot(np.diag(1.0 / new_cdelt))
             new_wcs.wcs.cdelt = new_cdelt
             for key, val in new_wcs.to_header().items():
                 header[key] = val
         try:
-            _ = header["ORIENTAT"]
+            header["ORIENTAT"] = princ_angle(float(header["ORIENTAT"]))
         except KeyError:
             header["ORIENTAT"] = wcs_PA(new_wcs.wcs.pc[1, 0], np.diag(new_wcs.wcs.pc).mean())
 
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)