Seems like this works
IMlist={ROI1,ROI2,ROI3};
AXlist={ax1,ax2,ax3};
FWHMX=[]; FWHMY=[]; FWHM2D=[]; CenInt=[];
for idx = 1:length(IMlist)
IM = IMlist{idx};
AX = AXlist{idx};
% IM=ROI1; AX=ax1;
[x1, y1]=myAnalysisFunctions.FindSpotsGaussian(IM,fwhm,thresh);
app.x1=x1; app.y1=y1;
F=isnan(x1); % 1 means its a nan, i.e no spots found
if F==0
hold(AX,'on'); plot(AX,x1,y1,'r.');
[mask,medMaxI]=myAnalysisFunctions.createMaskFromCentroids(x1,y1,IM,4);
[fwhmBlockX,fwhmBlockY,fwhm2D] = getFHWMbyACF(app,AX,IM,x1,y1,nl, linespacing, maxFwhm); drawnow
FWHMX = [FWHMX;fwhmBlockX]; FWHMY = [FWHMY;fwhmBlockY]; FWHM2D = [FWHM2D;fwhm2D];
CenInt = [CenInt; medMaxI];
end
end
And for the Table
params = ["ROI_L";"ROI_C";"ROI_R"];
% FWHMX = [fwhmBlockX1;fwhmBlockX2;fwhmBlockX3];
% FWHMY = [fwhmBlockY1;fwhmBlockY2;fwhmBlockY3];
% FWHM2D = [fwhm2D1; fwhm2D2;fwhm2D3];
% CenInt = [medMaxI1; medMaxI2;medMaxI3];
format bank
tbl = table(params,FWHMX,FWHMY,FWHM2D,CenInt);