diff --git a/CEP/PyBDSM/src/python/functions.py b/CEP/PyBDSM/src/python/functions.py
index 61d6fa6417bf32f8e581ed17a77d479213f92d0d..9b707b2d9098cf120e85dba729921fe6ed4ec35a 100755
--- a/CEP/PyBDSM/src/python/functions.py
+++ b/CEP/PyBDSM/src/python/functions.py
@@ -1120,12 +1120,11 @@ def read_image_from_file(filename, img, indir, quiet=False):
         key_val_raw = hdr['CTYPE' + str(i+1)]
         key_val = key_val_raw.split('-')[0]
         ctype_in.append(key_val.strip())
-    if img.use_io == 'fits':
-        ctype_in.reverse() # Need to reverse order, as pyfits does this
+    ctype_in.reverse() # Need to reverse order, as pyfits does this
 
     if 'RA' not in ctype_in or 'DEC' not in ctype_in:
         raise RuntimeError("Image data not found")
-    if 'FREQ' not in ctype_in:
+    if len(ctype_in) > 2 and 'FREQ' not in ctype_in:
         from pywcs import WCS
         t = WCS(hdr)
         t.wcs.fix()
diff --git a/CEP/PyBDSM/src/python/islands.py b/CEP/PyBDSM/src/python/islands.py
index e42b8fb4acb8fe8fc70c2248e220cae8dde524f4..16734659555dcb26e4199d1d110bb2e62936e633 100644
--- a/CEP/PyBDSM/src/python/islands.py
+++ b/CEP/PyBDSM/src/python/islands.py
@@ -73,7 +73,7 @@ class Op_islands(Op):
             img.island_labels = det_img.island_labels
             corr_islands = []
             for isl in det_img.islands:
-                corr_islands.append(isl.copy(img))
+                corr_islands.append(isl.copy(img, img.pixel_beamarea))
             img.islands = corr_islands
             img.nisl = len(img.islands)
             img.pyrank = det_img.pyrank