// run("Close All");

rename("canal verde 100")
run("Subtract Background...", "rolling=10 sliding stack");
run("Green")
setSlice(15);
run("Duplicate...", "title=[canal rojo 100] duplicate");
run("Red")
setSlice(15);

selectWindow("canal rojo 100");
run("Duplicate...", "title=[canal rojo 80] duplicate");
run("Multiply...", "value=.80 stack");
setSlice(15);

selectWindow("canal rojo 100");
run("Duplicate...", "title=[canal rojo 40] duplicate");
run("Multiply...", "value=.40 stack");
setSlice(15);

selectWindow("canal rojo 100");
run("Duplicate...", "title=[canal rojo 100 noise 15] duplicate");
run("Add Specified Noise...", "stack standard=15");
setSlice(15);

selectWindow("canal rojo 100");
run("Duplicate...", "title=[canal rojo 100 noise 25] duplicate");
run("Add Specified Noise...", "stack standard=25");
setSlice(15);

selectWindow("canal rojo 100");
run("Duplicate...", "title=[canal rojo 100 noise 50] duplicate");
run("Add Specified Noise...", "stack standard=50");
setSlice(15);

selectWindow("canal rojo 100");
run("Duplicate...", "title=[canal rojo 100 shift 02] duplicate");
run("Convolve...", "text1=[1 0 0 0 0 0 0\n] normalize stack");
setSlice(15);

selectWindow("canal rojo 100");
run("Duplicate...", "title=[canal rojo 100 shift 07] duplicate");
run("Convolve...", "text1=[1 0 0 0 0 0 0 0 0 0 0 0 0 0 0\n] normalize stack");
setSlice(15);


// Solo para visualizacion
/*
run("Merge Channels...", "c1=[canal rojo 100] c2=[canal verde 100] create keep");
setSlice(30);
run("Merge Channels...", "c1=[canal rojo 80] c2=[canal verde 100] create keep");
setSlice(30);
run("Merge Channels...", "c1=[canal rojo 40] c2=[canal verde 100] create keep");
setSlice(30);
run("Merge Channels...", "c1=[canal rojo 100 noise 15] c2=[canal verde 100] create keep");
setSlice(30);
run("Merge Channels...", "c1=[canal rojo 100 noise 25] c2=[canal verde 100] create keep");
setSlice(30);
run("Merge Channels...", "c1=[canal rojo 100 noise 50] c2=[canal verde 100] create keep");
setSlice(30);
run("Merge Channels...", "c1=[canal rojo 100 shift 07] c2=[canal verde 100] create keep");
setSlice(30);
run("Merge Channels...", "c1=[canal rojo 100 shift 02] c2=[canal verde 100] create keep");
setSlice(30);
*/

run("Tile");

/*
run("Coloc 2", "channel_1=[canal rojo 100] channel_2=[canal verde 100] roi_or_mask=<None> threshold_regression=Costes display_images_in_result li_histogram_channel_1 li_histogram_channel_2 li_icq manders'_correlation 2d_instensity_histogram psf=3 costes_randomisations=10");
run("Coloc 2", "channel_1=[canal rojo 80] channel_2=[canal verde 100] roi_or_mask=<None> threshold_regression=Costes display_images_in_result li_histogram_channel_1 li_histogram_channel_2 li_icq manders'_correlation 2d_instensity_histogram psf=3 costes_randomisations=10");
run("Coloc 2", "channel_1=[canal rojo 40] channel_2=[canal verde 100] roi_or_mask=<None> threshold_regression=Costes display_images_in_result li_histogram_channel_1 li_histogram_channel_2 li_icq manders'_correlation 2d_instensity_histogram psf=3 costes_randomisations=10");
run("Coloc 2", "channel_1=[canal rojo 100 noise 15] channel_2=[canal verde 100] roi_or_mask=<None> threshold_regression=Costes display_images_in_result li_histogram_channel_1 li_histogram_channel_2 li_icq manders'_correlation 2d_instensity_histogram psf=3 costes_randomisations=10");
run("Coloc 2", "channel_1=[canal rojo 100 noise 25] channel_2=[canal verde 100] roi_or_mask=<None> threshold_regression=Costes display_images_in_result li_histogram_channel_1 li_histogram_channel_2 li_icq manders'_correlation 2d_instensity_histogram psf=3 costes_randomisations=10");
run("Coloc 2", "channel_1=[canal rojo 100 noise 50] channel_2=[canal verde 100] roi_or_mask=<None> threshold_regression=Costes display_images_in_result li_histogram_channel_1 li_histogram_channel_2 li_icq manders'_correlation 2d_instensity_histogram psf=3 costes_randomisations=10");
run("Coloc 2", "channel_1=[canal rojo 100 shift 02] channel_2=[canal verde 100] roi_or_mask=<None> threshold_regression=Costes display_images_in_result li_histogram_channel_1 li_histogram_channel_2 li_icq manders'_correlation 2d_instensity_histogram psf=3 costes_randomisations=10");
run("Coloc 2", "channel_1=[canal rojo 100 shift 07] channel_2=[canal verde 100] roi_or_mask=<None> threshold_regression=Costes display_images_in_result li_histogram_channel_1 li_histogram_channel_2 li_icq manders'_correlation 2d_instensity_histogram psf=3 costes_randomisations=10");
*/