Replaced `is` with `==` for comparison
Test script
Replace `is` with `==` for comparison
Use np.subtract for subtraction to avoid ubyte overflow
imshow() to display the output
Change for 1 channel test image
Convert input array to `int` type to avoid ubyte subtraction overflow
Changes to include the BAND points as well
Smaller test image