USAGE: shapecompare [-m metric] [-t type] [-T thresh] [-e edgemode] [-i invert] infile1 infile2
USAGE: shapecompare [-h or -help]

-m ... metric ..... compare difference metric; choices are: sad, said, sand, rssd, all;
................... default=sad
-t ... type ....... type of image to convert input images for compare processing;
................... choices are: grayscale (g), binary (b) or outline (o);
................... default=grayscale
-T ... thresh ..... threshold for converting to binary; 0<=thresh<=100 (percent);
................... default=50
-e ... edgemode ... mode of morphologic edge processing; choices are: edge, edgin or
................... edgeout; default=edge
-i ... invert ..... invert (negate) colors in input image; choices are: yes (y) or
................... no (n); default=no

PURPOSE: To compares two images' shapes using metrics computed from Hu image moments.

DESCRIPTION: SHAPECOMPARE compares two images' shapes using metrics computed from Hu Image Moments. Several metrics are available. The image will first be converted to grayscale and can then be converted to binary or an edge outline. The script is designed to work with primarily binary or outline shapes on black backgrounds. However, grayscale shapes on black background or antialiased black/white shapes on black backgrounds should work, also. My analysis (see chart) indicates that outline mode does not work as well as the other two modes.

ARGUMENTS:

-m metric ... METRIC is the compare metric to compute the difference between the two images. Smaller is better than larger. The choices are: sad (sum absolute difference), said (sum absolute inverse difference), sand (sum absolute normalized difference), rssd (root sum squared difference) or all (report all of them). The default=sad.

The equations for the metrics are:
SAD = sum( abs(Hu1 - Hu2) )
SAID = sum( abs(1/Hu1 - 1/Hu2) )
SAND = sum( abs(Hu1 - Hu2) / sqrt( abs(Hu1*Hu2) ) )
RSSD = sqrt( sum( (Hu1 - Hu2)^2 ) )
where there are 7 different (log) Hu image moments that are summed

-t type ... TYPE of image to convert the input images for processing. The choices are: grayscale (g), binary (b) or outline (o). The default=grayscale.

-T thresh ... THRESH is the threshold for converting to binary. Values are: 0<=thresh<=100 (percent). The default=50.

-e edgemode. EDGEMODE is the mode of morphologic edge processing. The choices are: edge, edgin or edgeout. Edge will be 2 pixels thick. The other two will be only 1 pixel thick. All 3 use a diamond shaped kernel. The edgein will be on the white side and the edgeout will be on the black side. The default=edge.

-i invert ... INVERT (negate) the colors in input image as preprocessing. Choices are: yes (y) or no (n). The default=no

REFERENCE: https://www.learnopencv.com/shape-matching-using-hu-moments-c-python/

CAVEAT: No guarantee that this script will work on all platforms, nor that trapping of inconsistent parameters is complete and foolproof. Use At Your Own Risk.