Some minor fixes

CanlabCore/@region/isosurface.m

@@ -82,20 +82,20 @@
 
 cl = region2struct(r);
 
-% surface_handles = []; % This appears to be a bug, as it coerces the Patch object into a double. 05/23/2023 MS
-surface_handles = {};
+surface_handles = []; % This appears to be a bug, as it coerces the Patch object into a double. 05/23/2023 MS
+% surface_handles = {};
 
 for i = 1:k
 
     try
         out = imageCluster('cluster', cl(i), 'color', colors{i}, varargin{:});
 
-%         surface_handles(i) = out; % This appears to be a bug, as it
+        surface_handles(i) = out; % This appears to be a bug, as it
 %         coerces the Patch object into a double. 05/23/2023 MS
-        surface_handles{i} = out;
+%        surface_handles{i} = out;
     catch
         disp('Error imaging isosurface; too few voxels?')
-        %surface_handles(i) = [];
+        surface_handles(i) = [];
     end
 
 %     % Isocaps, if needed

CanlabCore/@region/region.m

@@ -384,9 +384,11 @@
                 obj(i).dim = mask.volInfo.dim;
                 obj(i).numVox = size(obj(i).XYZ, 2);
 
-                obj(i).center = center_of_mass(myXYZ, double(obj(i).Z));
-                obj(i).mm_center = center_of_mass(obj(i).XYZmm, double(obj(i).Z));
-
+                if ~isempty(myXYZ) % XYZ initialization may be skipped, so skip this if so MS 2/15/2023
+                    obj(i).center = center_of_mass(myXYZ, double(obj(i).Z));
+                    obj(i).mm_center = center_of_mass(obj(i).XYZmm, double(obj(i).Z));
+                end
+
                 if ~isempty(mask.volInfo.fname)
                     [~, ff, ee] = fileparts(mask.volInfo.fname);
                     obj(i).custom_info1 = [ff ee];