function  [Objects] = f_imageAnalysis(Label, chThioS, chSynucleinP, chSynuclein, chNuclei, PreviewPath, MetaTable)
    
    %vol(chHoechst, 0, 10000)% Hoechst
    %vol(chMap2, 0, 5000)
    %vol(chSynucleinP, 0, 2000)
    %vol(chSynuclein, 0,1000)

    %% Segment Organoid
    chOrganoid = imlincomb(0.25, chNuclei, 0.25, chThioS, 0.25, chSynuclein, 0.25, chSynucleinP); %vol(chOrganoid, 0, 4000)
    OrganoidMask = chOrganoid > 175;
    %OrganoidMask = Iris_KeepbiggestObct3D(OrganoidMask);

    
    OrganoidMask = bwareaopen(OrganoidMask, 1e2);
    OrganoidMask = imclose(OrganoidMask, strel('disk', 50));
    OrganoidMask = imclose(OrganoidMask, strel('sphere', 5));
    OrganoidMask = logical(imfill(OrganoidMask, 'holes')); % vol(OrganoidMask,0,1)
    
    OrganoidMask = bwareaopen(OrganoidMask, 1e6);
    
    OrganoidMaskEroded = imerode(OrganoidMask, strel('disk', 250));
    OrganoidMaskOuterRingCandidates = OrganoidMask & ~OrganoidMaskEroded; % imtool(max(OrganoidMaskOuterRingCandidates,[],3))
    OrganoidMaskOuterRingObjects = regionprops('table', bwlabeln(OrganoidMaskOuterRingCandidates), {'Area', 'PixelIdxList'});
    OrganoidMaskOuterRingObjects = OrganoidMaskOuterRingObjects(OrganoidMaskOuterRingObjects.Area > 5e6, :); % todo  > 1e7
    OrganoidMaskOuterRing = f_Create_Mask_from_ObjectList_Pixel_IDX(OrganoidMaskOuterRingObjects, 'PixelIdxList', OrganoidMask); % imtool(max(OrganoidMaskOuterRing,[],3))
    OrganoidMaskCurated = OrganoidMask & ~OrganoidMaskOuterRing; % vol(OrganoidMaskCurated,0,1)
    
    %% Segment Thioflavin S structures
    ThioSMask = chThioS > 2500;
    ThioSMask = bwareaopen(ThioSMask, 50);
    ThioSMask = ThioSMask & OrganoidMaskCurated; % vol(NucMask, 0,1)

      %% Segment SynucleinP
    % Please exclude autofluoresc. and check if spots are synapses (i.e. flip 2ndary ABs)
    
    SNCAMaskP = chSynucleinP > 6500; % vol(SNCAMaskP)
    SNCAMaskP = bwareaopen(SNCAMaskP, 30);

    SNCAMaskP = SNCAMaskP & OrganoidMaskCurated;
    
    %% Segment Synuclein
    % Please exclude autofluoresc. 
    
    SNCAMask = chSynuclein > 6000; % vol(SNCAMask)
    SNCAMask = bwareaopen(SNCAMask, 50);
   
    SNCAMask = SNCAMask & OrganoidMaskCurated;
    
    %% Segment SYTOX Nuclei
    NucleiMask = medfilt3(chNuclei) > 500;
    NucleiMask = bwareaopen(NucleiMask, 50);
    ThioSMask = ThioSMask & OrganoidMaskCurated; % vol(NucMask, 0,1)
 
    %% Feature extraction
    Objects = table();
    Objects.OrganoidIdx = MetaTable.Idx;
    Objects.Well = MetaTable.Well;
    Objects.AreaName = MetaTable.AreaName;
    Objects.OrganoidID = MetaTable.OrganoidID;
    
    % Volumes in voxels
    Objects.ThioSArea = sum(ThioSMask(:));
    Objects.NucleiArea = sum(NucleiMask(:));
    Objects.SNCAArea = sum(SNCAMask(:));
    Objects.SNCAPArea = sum(SNCAMaskP(:));
    Objects.OrganoidVolume = sum(OrganoidMask(:));
    
    % Normalized voxel/voxel
    Objects.OrganoidVolumeCurated = sum(OrganoidMaskCurated(:));
    Objects.ThioSVoxelsByNucVoxels = Objects.ThioSArea / Objects.NucleiArea;
    Objects.SNCAVoxelsByNucVoxels = Objects.SNCAArea / Objects.NucleiArea;
    Objects.SNCAPVoxelsByNucVoxels = Objects.SNCAPArea / Objects.NucleiArea;
    Objects.SNCAVoxelsByThioSVoxels = Objects.SNCAArea / Objects.ThioSArea;
    Objects.SNCAPVoxelsByThioSVoxels = Objects.SNCAPArea / Objects.ThioSArea;
    
    % means in ROIs
    %Objects.MeanSNCAIntensityInOrganoidVoxels = mean(chSynuclein(OrganoidMask));
    Objects.MeanSNCAIntensityInNucVoxels = mean(chSynuclein(NucleiMask));
    Objects.MeanSNCAIntensityInThioSVoxels = mean(chSynuclein(ThioSMask));
    %Objects.MeanSNCAPIntensityInOrganoidVoxels = mean(chSynucleinP(OrganoidMask));
    Objects.MeanSNCAPIntensityInNucVoxels = mean(chSynucleinP(NucleiMask));
    Objects.MeanSNCAPIntensityInThioSVoxels = mean(chSynucleinP(ThioSMask));

  
    %% Previews 
    % Scalebar
    imSize = [size(chNuclei, 1), size(chNuclei, 2)];
    [BarMask, BarCenter] = f_barMask(200, 0.32393102760889064, imSize, imSize(1)-200, 200, 25);
    %it(BarMask)

    RGB = cat(3, imadjust(max(chSynucleinP, [], 3), [0 0.09], [0 1]), imadjust(max(chNuclei, [], 3), [0 0.1], [0 1]), imadjust(max(chSynuclein, [], 3), [0 0.01], [0 1]));
    RGB = imoverlay2(RGB, BarMask, [1 1 1]);
    %imtool(RGB)
     
    %PreviewThioS = imoverlay2(imadjust(max(chThioS,[],3), [0 0.5], [0 1]), bwperim(max(ThioSMask,[],3)), [1 0 0]);
    PreviewThioS = f_imoverlayAlpha(imadjust(max(chThioS,[],3), [0 0.05], [0 1]), max(ThioSMask,[],3), [0 1 0], 0.1);
    PreviewThioS = imoverlay2(PreviewThioS, BarMask, [1 1 1]);
    %imtool(PreviewThioS,[])
    
    PreviewNuclei = f_imoverlayAlpha(imadjust(max(chNuclei,[],3), [0 0.3], [0 1]), max(NucleiMask,[],3), [0 0 1], 0.1);
    PreviewNuclei = imoverlay2(PreviewNuclei, BarMask, [1 1 1]);
    %imtool(PreviewNuclei,[])

    PreviewSNCA = f_imoverlayAlpha(imadjust(max(chSynuclein,[],3), [0 0.4], [0 1]), max(SNCAMask,[],3), [1 1 0], 0.1);
    PreviewSNCA = imoverlay2(PreviewSNCA, BarMask, [1 1 1]);
    %imtool(PreviewSNCA,[])
    
    PreviewSNCAP = f_imoverlayAlpha(imadjust(max(chSynucleinP,[],3), [0 0.4], [0 1]), max(SNCAMaskP,[],3), [0 1 1], 0.1);
    PreviewSNCAP = imoverlay2(PreviewSNCAP, BarMask, [1 1 1]);
    %imtool(PreviewSNCAP,[])

    %PreviewOrganoid = f_imoverlayAlpha(imadjust(max(chNuclei,[],3), [0 0.5], [0 1]), max(OrganoidMask,[],3), [1 0 1], 0.05);
    %PreviewOrganoid = f_imoverlayAlpha(PreviewOrganoid, max(OrganoidMaskCurated,[],3), [0 1 0], 0.05);
    %PreviewOrganoid = imoverlay2(PreviewOrganoid, BarMask, [1 1 1]);
    %imtool(PreviewOrganoid,[])

    %% Write Previews 
    IdentityString = ['PreviewOrganoid', '_', MetaTable.OrganoidID{:}];
    imwrite(RGB, [PreviewPath, filesep, IdentityString, '_RGB.png'])
    imwrite(PreviewThioS, [PreviewPath, filesep, IdentityString, '_ThioS.png'])
    imwrite(PreviewNuclei, [PreviewPath, filesep, IdentityString, '_Nuclei.png'])
    imwrite(PreviewSNCA, [PreviewPath, filesep, IdentityString, '_SNCA.png'])
    imwrite(PreviewSNCAP, [PreviewPath, filesep, IdentityString, '_SNCAP.png'])
    %imwrite(PreviewOrganoid, [PreviewPath, filesep, IdentityString, '_Organoid.png'])

end
