import numpy as np
import cv2 as cv


img=cv.imread('N.ppm',0)
img1=cv.imread('N45.ppm',0)
cv.imshow('image',img)
cv.imshow('image 45',img1)
t,trsh=cv.threshold(img,127,255,cv.THRESH_BINARY)
cv.imshow('Thresholded image',trsh)
                 
raw_m=cv.moments(img,True)
raw_m45=cv.moments(img1,True)
#trsh_m=cv.moments(trsh,True)
hu_m=cv.HuMoments(raw_m)
hu_m45=cv.HuMoments(raw_m45)
print(hu_m)
print(hu_m45)
#print(np.sign(hu_m)*np.log10(abs(hu_m)) )
#print(np.sign(hu_m45)*np.log10(abs(hu_m45)))
#for x in raw_m:
#    print(x,'=',np.sign(raw_m[x])*np.log10(abs(raw_m[x])),"\t: ",np.sign(trsh_m[x])*np.log10(abs(trsh_m[x])))
cv.waitKey(0)
cv.destroyAllWindows()