/*Calculates the quadratic average of the matching cardiac phase ensemble of a 1000 frames acquisition with an estimated 75 frames/cardiac cycle*/

//Crop non-essential parts of the echocardiographic image. Adjust crop-size to actual acquisition.
makeRectangle(285, 120, 576, 600);
run("Crop");

//Covert from 8-bit RGB to 8-bit grayscale and then to 32-bit.
run("8-bit");
run("32-bit");

//Measure quadratic average of all frames at one cardiac phase at a time. Slices in substack are defined by local maxima in correlation value over frame curve.

//Phase 1

rename("X");
run("Make Substack...", "  slices=70,103,139,207,276,345,413,482,551,584,620,688,757,826,895,963");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=16");
run("Square Root");
rename("1");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 2
selectWindow("X");
run("Make Substack...", "  slices=2,70,139,208,277,345,414,483,552,620,689,757,826,895,964");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("2");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 3
selectWindow("X");
run("Make Substack...", "  slices=3,71,140,209,277,346,415,483,552,621,690,759,827,896,965");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("3");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 4
selectWindow("X");
run("Make Substack...", "  slices=4,73,141,210,279,348,416,485,516,554,622,691,723,760,828,897,966");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=17");
run("Square Root");
rename("4");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 5
selectWindow("X");
run("Make Substack...", "  slices=5,74,142,211,280,349,417,486,554,623,692,722,761,829,898,967");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=16");
run("Square Root");
rename("5");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 6
selectWindow("X");
run("Make Substack...", "  slices=6,75,144,212,281,350,419,487,556,625,693,722,762,831,900,968");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=16");
run("Square Root");
rename("6");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 7
selectWindow("X");
run("Make Substack...", "  slices=7,76,103,144,213,282,350,419,488,556,625,694,723,763,832,900,969");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=17");
run("Square Root");
rename("7");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 8
selectWindow("X");
run("Make Substack...", "  slices=8,76,103,144,214,282,351,420,488,557,626,695,723,763,791,832,901,970");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=18");
run("Square Root");
rename("8");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 9
selectWindow("X");
run("Make Substack...", "  slices=9,79,103,147,216,284,354,423,491,560,584,628,698,722,765,834,903,972");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=18");
run("Square Root");
rename("9");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 10
selectWindow("X");
run("Make Substack...", "  slices=10,79,148,216,285,354,423,491,560,629,698,766,835,903,972");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("10");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 11
selectWindow("X");
run("Make Substack...", "  slices=11,79,148,216,285,354,423,492,560,629,698,767,835,903,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("11");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 12
selectWindow("X");
run("Make Substack...", "  slices=12,79,120,149,216,259,285,354,406,424,492,532,561,600,630,699,767,791,807,837,876,903,943,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=24");
run("Square Root");
rename("12");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 13
selectWindow("X");
run("Make Substack...", "  slices=13,79,149,214,285,354,379,423,492,516,563,600,629,699,722,767,791,837,903,974");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=20");
run("Square Root");
rename("13");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 14
selectWindow("X");
run("Make Substack...", "  slices=14,53,83,131,149,190,220,289,326,358,395,423,447,463,492,532,563,584,600,629,701,722,767,791,807,837,876,907,944");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=30");
run("Square Root");
rename("14");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 15
selectWindow("X");
run("Make Substack...", "  slices=15,62,84,103,149,221,290,326,359,406,427,463,496,532,563,600,634,702,767,808,840,876,908,944,977");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=25");
run("Square Root");
rename("15");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 16
selectWindow("X");
run("Make Substack...", "  slices=16,51,86,103,131,154,173,199,222,258,291,326,360,379,395,423,463,492,532,566,600,635,654,670,705,723,739,772,807");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=35");
run("Square Root");
rename("16");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 17
selectWindow("X");
run("Make Substack...", "  slices=17,35,53,86,120,155,183,223,258,292,310,326,360,379,406,423,463,492,532,567,600,636,670,705,740,773,807,841,876");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=32");
run("Square Root");
rename("17");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 18
selectWindow("X");
run("Make Substack...", "  slices=18,51,87,120,156,190,224,258,293,326,362,395,430,463,500,532,567,600,636,670,705,740,775,807,843,876,912,943,980");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=29");
run("Square Root");
rename("18");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 19
selectWindow("X");
run("Make Substack...", "  slices=19,53,87,120,156,190,225,293,326,364,395,432,463,500,532,569,600,637,670,706,775,807,844,876,912,981");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=26");
run("Square Root");
rename("19");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 20
selectWindow("X");
run("Make Substack...", "  slices=20,51,89,120,158,183,226,295,364,395,433,500,532,569,600,637,670,706,775,845,876,913,981");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=23");
run("Square Root");
rename("20");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 21
selectWindow("X");
run("Make Substack...", "  slices=21,51,90,120,160,190,228,258,296,326,365,395,433,463,500,532,569,600,640,670,706,740,776,808,846,876,914,983");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=28");
run("Square Root");
rename("21");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 22
selectWindow("X");
run("Make Substack...", "  slices=22,53,92,120,160,190,229,258,297,326,365,395,436,463,504,532,571,600,640,670,710,739,775,847,876,914,983");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=27");
run("Square Root");
rename("22");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 23
selectWindow("X");
run("Make Substack...", "  slices=23,51,92,120,160,229,258,298,326,367,395,436,463,504,532,571,600,641,670,710,779,808,848,876,918,983");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=26");
run("Square Root");
rename("23");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 24
selectWindow("X");
run("Make Substack...", "  slices=24,53,91,120,160,230,298,368,395,436,505,532,572,600,642,670,710,779,848,876,918,983");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=22");
run("Square Root");
rename("24");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 25
selectWindow("X");
run("Make Substack...", "  slices=25,51,94,120,162,231,301,326,368,395,436,505,532,575,600,643,670,712,780,848,876,918,988");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=23");
run("Square Root");
rename("25");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 26
selectWindow("X");
run("Make Substack...", "  slices=26,51,95,120,164,232,258,301,326,368,395,438,507,532,575,600,643,670,712,780,808,851,876,918,988");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=25");
run("Square Root");
rename("26");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 27
selectWindow("X");
run("Make Substack...", "  slices=27,97,120,164,199,233,302,370,439,508,578,645,714,785,808,851,920,989");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=18");
run("Square Root");
rename("27");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 28
selectWindow("X");
run("Make Substack...", "  slices=28,98,120,166,233,259,302,371,439,508,578,646,714,785,808,851,876,922,989");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=19");
run("Square Root");
rename("28");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 29
selectWindow("X");
run("Make Substack...", "  slices=29,51,98,120,167,190,237,259,304,373,395,441,511,532,579,600,648,670,717,740,785,808,853,876,923,991");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=26");
run("Square Root");
rename("29");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 30
selectWindow("X");
run("Make Substack...", "  slices=11,30,51,79,98,120,167,237,259,305,354,373,395,423,442,511,532,579,600,629,648,670,698,717,785,855,876,923,991");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=29");
run("Square Root");
rename("30");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 31
selectWindow("X");
run("Make Substack...", "  slices=11,31,53,98,120,168,237,305,354,374,395,423,442,492,512,580,600,629,649,718,787,835,855,876,923,972,993");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=27");
run("Square Root");
rename("31");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 32
selectWindow("X");
run("Make Substack...", "  slices=11,32,101,120,149,169,201,238,307,375,423,444,492,513,532,582,600,629,650,698,719,767,788,856,925,973,994");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=27");
run("Square Root");
rename("32");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 33
selectWindow("X");
run("Make Substack...", "  slices=11,33,79,102,120,149,171,216,239,258,308,354,376,423,445,492,516,532,560,583,600,629,651,698,720,767,789,857,876");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=31");
run("Square Root");
rename("33");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 34
selectWindow("X");
run("Make Substack...", "  slices=11,34,79,103,148,172,216,240,309,354,379,423,447,492,515,563,584,629,654,698,721,767,791,835,859,928,973,996");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=28");
run("Square Root");
rename("34");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 35
selectWindow("X");
run("Make Substack...", "  slices=35,53,79,104,148,172,216,241,285,310,354,379,423,447,492,516,560,584,629,653,698,722,767,791,808,835,859,903,928");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=31");
run("Square Root");
rename("35");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 36
selectWindow("X");
run("Make Substack...", "  slices=36,79,104,149,173,216,242,285,311,354,379,423,448,492,517,563,586,629,654,698,723,767,792,808,837,860,903,929,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=30");
run("Square Root");
rename("36");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 37
selectWindow("X");
run("Make Substack...", "  slices=37,79,106,149,174,216,243,285,312,354,380,423,449,492,518,560,587,629,655,698,724,767,793,835,861,903,930,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=28");
run("Square Root");
rename("37");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 38
selectWindow("X");
run("Make Substack...", "  slices=11,38,79,107,148,175,216,244,285,313,354,381,423,450,492,519,560,587,629,656,698,725,767,794,835,862,903,931,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=30");
run("Square Root");
rename("38");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 39
selectWindow("X");
run("Make Substack...", "  slices=39,79,108,148,176,216,245,285,314,354,383,423,452,492,520,560,600,629,658,698,727,767,795,835,864,903,933,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=28");
run("Square Root");
rename("39");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 40
selectWindow("X");
run("Make Substack...", "  slices=11,40,79,109,148,177,216,246,285,315,354,383,423,452,492,521,560,600,629,658,698,727,767,796,835,864,903,933,972");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=29");
run("Square Root");
rename("40");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 41
selectWindow("X");
run("Make Substack...", "  slices=11,41,79,110,148,178,216,247,285,316,354,384,423,453,492,522,560,600,629,659,698,728,767,796,835,865,903,934,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=30");
run("Square Root");
rename("41");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 42
selectWindow("X");
run("Make Substack...", "  slices=42,79,120,148,183,216,251,285,326,354,395,423,463,492,532,560,600,629,663,698,740,767,807,835,869,903,944,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=28");
run("Square Root");
rename("42");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 43
selectWindow("X");
run("Make Substack...", "  slices=11,43,79,120,148,183,216,251,285,326,354,395,423,463,492,532,560,600,629,663,698,740,767,807,835,869,903,944,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=29");
run("Square Root");
rename("43");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 44
selectWindow("X");
run("Make Substack...", "  slices=11,44,79,120,148,183,216,251,285,326,354,395,423,463,492,532,560,600,629,663,698,740,767,807,835,869,903,944,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=29");
run("Square Root");
rename("44");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 45
selectWindow("X");
run("Make Substack...", "  slices=11,45,79,120,148,183,251,326,354,395,423,463,492,532,560,600,629,663,698,739,767,807,835,869,944,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=26");
run("Square Root");
rename("45");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 46
selectWindow("X");
run("Make Substack...", "  slices=11,46,79,120,148,183,251,326,354,395,423,463,492,532,560,600,629,663,698,740,767,807,835,869,944,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=26");
run("Square Root");
rename("46");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 47
selectWindow("X");
run("Make Substack...", "  slices=11,47,120,148,183,259,326,395,463,492,532,560,600,629,670,740,767,807,835,876,943,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=22");
run("Square Root");
rename("47");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 48
selectWindow("X");
run("Make Substack...", "  slices=11,48,79,120,148,190,258,326,354,395,463,492,532,560,600,629,670,698,740,767,807,835,876,943,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=25");
run("Square Root");
rename("48");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 49
selectWindow("X");
run("Make Substack...", "  slices=11,49,120,190,258,326,395,423,463,492,532,560,600,629,670,698,740,767,807,835,876,943,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=23");
run("Square Root");
rename("49");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 50
selectWindow("X");
run("Make Substack...", "  slices=11,50,79,120,148,190,258,326,354,395,423,463,492,532,560,600,629,670,698,739,767,807,835,876,944,972");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=26");
run("Square Root");
rename("50");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 51
selectWindow("X");
run("Make Substack...", "  slices=11,51,120,190,258,326,395,423,463,492,532,560,600,629,670,698,739,767,807,876,944,973");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=22");
run("Square Root");
rename("51");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 52
selectWindow("X");
run("Make Substack...", "  slices=11,52,121,190,258,327,395,464,492,532,600,670,698,739,808,877,945");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=17");
run("Square Root");
rename("52");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 53
selectWindow("X");
run("Make Substack...", "  slices=11,53,122,190,259,328,396,465,492,534,600,672,740,808,877,946");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=16");
run("Square Root");
rename("53");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 54
selectWindow("X");
run("Make Substack...", "  slices=11,54,122,191,259,328,397,466,492,534,603,672,740,809,877,946");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=16");
run("Square Root");
rename("54");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 55
selectWindow("X");
run("Make Substack...", "  slices=8,55,123,193,261,329,399,467,535,604,672,740,811,880,946");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("55");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 56
selectWindow("X");
run("Make Substack...", "  slices=11,56,123,193,262,330,398,467,537,605,674,743,811,880,949");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("56");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 57
selectWindow("X");
run("Make Substack...", "  slices=8,57,126,194,263,332,400,469,538,607,674,744,813,881,949");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("57");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 58
selectWindow("X");
run("Make Substack...", "  slices=8,58,126,195,263,332,401,470,539,607,676,745,814,882,951");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("58");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 59
selectWindow("X");
run("Make Substack...", "  slices=8,59,127,199,265,334,402,471,540,608,677,746,814,883,951");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("59");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 60
selectWindow("X");
run("Make Substack...", "  slices=6,60,131,199,266,335,403,476,543,610,683,749,816,888,953");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("60");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 61
selectWindow("X");
run("Make Substack...", "  slices=8,61,131,199,270,335,406,476,543,613,679,749,816,888,957");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("61");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 62
selectWindow("X");
run("Make Substack...", "  slices=8,62,131,200,270,338,406,476,543,613,682,749,821,888,957");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("62");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 63
selectWindow("X");
run("Make Substack...", "  slices=2,63,131,201,270,338,406,476,544,613,682,750,821,888,957");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("63");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 64
selectWindow("X");
run("Make Substack...", "  slices=2,64,132,201,270,339,407,476,545,613,682,751,821,888,957");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("64");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 65
selectWindow("X");
run("Make Substack...", "  slices=8,65,134,202,271,340,408,478,546,615,683,752,821,889,959");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("65");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 66
selectWindow("X");
run("Make Substack...", "  slices=8,66,134,203,272,341,409,478,547,615,684,753,822,890,959");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("66");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 67
selectWindow("X");
run("Make Substack...", "  slices=2,67,136,205,273,342,411,479,548,617,685,754,823,891,960");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("67");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 68
selectWindow("X");
run("Make Substack...", "  slices=1,68,137,205,274,343,411,480,549,617,686,756,823,892,961");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("68");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 69
selectWindow("X");
run("Make Substack...", "  slices=1,69,137,206,275,344,412,481,550,618,687,756,825,894,962");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("69");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 70
selectWindow("X");
run("Make Substack...", "  slices=2,70,139,208,276,345,414,482,551,620,688,757,826,895,963");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=15");
run("Square Root");
rename("70");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 71
selectWindow("X");
run("Make Substack...", "  slices=2,71,140,208,277,346,414,483,552,621,689,723,758,827,895,964");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=16");
run("Square Root");
rename("71");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 72
selectWindow("X");
run("Make Substack...", "  slices=4,72,141,210,279,347,416,485,553,622,691,722,759,828,897,927,966");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=17");
run("Square Root");
rename("72");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 73
selectWindow("X");
run("Make Substack...", "  slices=4,35,73,103,142,210,279,348,416,485,554,622,692,760,829,898,966");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=17");
run("Square Root");
rename("73");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 74
selectWindow("X");
run("Make Substack...", "  slices=5,74,143,211,280,349,417,486,516,555,624,692,761,830,898,927,967");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=17");
run("Square Root");
rename("74");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();

//Phase 75
selectWindow("X");
run("Make Substack...", "  slices=6,75,104,144,212,281,350,419,487,556,625,694,762,791,831,900,968");
rename("Y");
run("Z Project...", "projection=[Average Intensity]");
imageCalculator("Subtract create stack", "Y","AVG_Y");
selectWindow("Result of Y");
run("Square", "stack");
run("Z Project...", "projection=[Sum Slices]");
run("Divide...", "value=17");
run("Square Root");
rename("75");
selectWindow("Result of Y");
close();
selectWindow("AVG_Y");
close();
selectWindow("Y");
close();
selectWindow("X");
close();
run("Images to Stack", "name=Stack title=[] use");