// Morphological2Derivative3.txt
// Morphological2Derivative1-Morphological2Derivative2
// by G. Landini at bham. ac. uk
// 13/Nov/03

a=getNumber("Radius",1);
run("Duplicate...", "title=original");
run("Duplicate...", "title=dilation");
run("Duplicate...", "title=smoothed");
run("Duplicate...", "title=original2");
selectWindow("dilation");
run("Maximum...", "radius="+a);
selectWindow("smoothed");
run("Minimum...", "radius="+a);
run("Image Calculator...", "image1=smoothed operation=Average image2=dilation");
selectWindow("dilation");
run("Close");

selectWindow("original");
run("Divide...", "value=2");
run("Add...", "value=127");
selectWindow("smoothed");
run("Divide...", "value=2");
run("Image Calculator...", "image1=original operation=Subtract image2=smoothed");
selectWindow("smoothed");
run("Close");

selectWindow("original2");
run("Duplicate...", "title=dilation");
run("Duplicate...", "title=smoothed");

selectWindow("dilation");
run("Maximum...", "radius="+a);
run("Minimum...", "radius="+a);
selectWindow("smoothed");
run("Minimum...", "radius="+a);
run("Maximum...", "radius="+a);
run("Image Calculator...", "image1=smoothed operation=Average image2=dilation");
selectWindow("dilation");
run("Close");

selectWindow("original2");
run("Divide...", "value=2");
run("Add...", "value=127");
selectWindow("smoothed");
run("Divide...", "value=2");
run("Image Calculator...", "image1=original2 operation=Subtract image2=smoothed");
selectWindow("smoothed");
run("Close");
selectWindow("original");
run("Add...", "value=127");
run("Image Calculator...", "image1=original operation=Subtract image2=original2");

selectWindow("original2");
run("Close");