#統(tǒng)計直方圖，用load()載入圖片的像素pix，再分別讀取每個像素點的R\G\B值進行統(tǒng)計（分別為0-255）
#將256個顏色值的統(tǒng)計情況投影到32個，返回R\G\B投影后的統(tǒng)計值數(shù)組，共32*3=96個元素
def calc_Hist(img):
  '''
  #120張圖片，4.43s
  w,h = img.size
  pix = img.load() #載入圖片，pix存的是像素
  calcR = [0 for i in range(0,32)]
  calcG = [0 for i in range(0,32)]
  calcB = [0 for i in range(0,32)]
  for i in range(0,w):
    for j in range(0,h):
      (r,g,b) = pix[i,j]
      #print (r,g,b)
      calcR[r/8] += 1
      calcG[g/8] += 1
      calcB[b/8] += 1
  calcG.extend(calcB)
  calcR.extend(calcG)

  return calcR
  '''
  #120張圖，3.49s

  w,h = img.size
  pix = img.load() #載入圖片，pix存的是像素
  calcR = [0 for i in range(0,256)]
  calcG = [0 for i in range(0,256)]
  calcB = [0 for i in range(0,256)]
  for i in range(0,w):
    for j in range(0,h):
      (r,g,b) = pix[i,j]
      #print (r,g,b)
      calcR[r] += 1
      calcG[g] += 1
      calcB[b] += 1
  calcG.extend(calcB)
  calcR.extend(calcG) #256*3

  #calc存放最終結(jié)果，32*3
  calc = [0 for i in range(0,96)]
  step = 0 #calc的下標，0~95
  start = 0 #每次統(tǒng)計的開始位置
  while step < 96:
    for i in range(start,start+8): #8個值為1組，統(tǒng)計值相加，eg：色彩值為0~7的統(tǒng)計值全部轉(zhuǎn)換為色彩值為0的統(tǒng)計值
      calc[step] += calcR[i]
    start = start+8
    step += 1
  #print calc 
  return calc 

#該函數(shù)用于統(tǒng)一圖片大小為256*256，并且分割為16個塊，返回值是16個局部圖像句柄的數(shù)組
def split_Img(img, size = (64,64)):
  img = img.resize((256,256)).convert('RGB')
  w,h = img.size
  sw,sh = size
  return [img.crop((i,j,i+sw,j+sh)).copy() for i in xrange(0,w,sw) for j in xrange(0,h,sh)]

#計算兩個直方圖之間的相似度，h1和h2為直方圖，zip表示同步遍歷
def calc_Similar(h1,h2):
  return sum(1 - (0 if g==s else float(abs(g-s))/max(g,s)) for g,s in zip(h1,h2)) / len(h1)

results = {} #記錄結(jié)果
reverse = True #correlation/intersection方法reverse為true，另外兩種為false

imgCV = cv2.imread(self.testImg.encode('utf-8'))
#self.testImg為待匹配圖片
testHist = cv2.calcHist([imgCV],[0,1,2],None,[8,8,8],[0,256,0,256,0,256])
#提取直方圖
testHist = cv2.normalize(testHist,testHist,0,255,cv2.NORM_MINMAX).flatten()
#均衡化

#計算self.testImg和其他圖片的直方圖差異，INTERSECTION方法效果比較好
for (k, hist) in self.index_cv.items(): 
#self.index_cv保存的是圖片庫中圖片的直方圖信息
  d = cv2.compareHist(testHist,hist, cv2.cv.CV_COMP_INTERSECT)
  results[k] = d
  #對結(jié)果排序，以v即上面的d作為關(guān)鍵字
  results = sorted([(v, k) for (k, v) in results.items()], reverse = reverse) 
  end = time.time()
  print 'OpenCV Time:'
  print end-start     
self.DisplayTotalPics(results)

#獲取排序時的關(guān)鍵值（即漢明距離）    
def getKey(x): 
  return int(x[1])

#由灰度圖得到2值“指紋”，從而計算漢明距離
def getCode(img):
  w,h = img.size
  pixel = []
  for i in range(0,w):
    for j in range(0,h):
      pixel_value = img.getpixel((i,j))
      pixel.append(pixel_value) #加入pixel數(shù)組
  avg = sum(pixel)/len(pixel) #計算像素平均值

  cp = [] #二值數(shù)組
  for px in pixel:
    if px > avg:
      cp.append(1)
    else:
      cp.append(0)
  return cp

#計算兩個編碼之間的漢明距離
def compCode(code1,code2):
  num = 0
  for index in range(0,len(code1)):
    if code1[index] != code2[index]:
      num+=1
  #print num
  #print '\n'
  return num