| Scripting Resources for DigitalMicrograph™ |
Calculate Areal Percentages in Binary Images |
|
Function |
Sums the front-most image to compute the areal percentages of black and white pixels in a binary image. |
Version |
version:20131110, v1.0 |
Author |
D. R. G. Mitchell |
Acknowledgements |
- |
Comments |
The image data type does not need to be binary for this script to work, since greyscale images can also contain values of only 0 and 1. However, the script does check the minimum and maximum values and only computes percentages if these are 0 and 1 respectively. The script will work on both 1D and 2D images. |
System Requirements |
Should be compatible with all recent versions of DigitalMicrograph. |
Known Issues |
- |
Supported |
Yes |
Included Files |
Main script file. |
Source Code |
// Simple script to compute the areal percentage // of black and white pixels in a binary imagee
// D. R. G. Mitchell, adminnospam@dmscripting.com (remove the nospam to make this work) // version:20131110, v1.0, Nov. 2013, www.dmscripting.com
// Source the front-most image and get some information from it
number nodocs=countdocumentwindowsoftype(5) if(nodocs<1) { showalert("Ensure an image containing binary values (0 and 1) is displayed front-most.",2) exit(0) }
image front:=getfrontimage() string imgname=getname(front) number xsize,ysize
getsize(front, xsize, ysize) number maxvalue, minvalue front.minmax(minvalue, maxvalue)
// Check if the image contains binary values (0 and 1)
result("\n\n\nMeasurement of Image: "+imgname+"\n")
if(minvalue!=0 || maxvalue!=1) { result("\nWARNING: This script expects the values to be binary - values of either 0 or 1.") result("\nThe image values are not compatible with this script:\n\nMinimum value ="+minvalue+"\nMaximum value = "+maxvalue+"\n") exit(0) }
// Sum the white pixels and compute the fractions
number sum=sum(front) number whitepct=(sum/(xsize*ysize))*100 number blackpct=100-whitepct
result("\nFraction of white pixels = "+format(whitepct, "%4.2f")+" %") result("\nFraction of black pixels = "+format(blackpct, "%4.2f")+" %")
|