Python + OpenCV

import cv2
import numpy as np
import math
import matplotlib.pyplot as plt
import imageio
 
#Загружаю изображения
from random import random
image1 = cv2.imread('D:\ph4.png')
plt.imshow(cv2.cvtColor(image1, cv2.COLOR_BGR2RGB))
print(image1.shape)
image2 = cv2.imread('D:\ph5.png')
plt.imshow(cv2.cvtColor(image2, cv2.COLOR_BGR2RGB))
print(image2.shape)
 
#Нахожу угловые точки на первом
gray = cv2.cvtColor(image1, cv2.COLOR_BGR2GRAY)
gray = np.float32(gray)
dst = cv2.cornerHarris(gray,5,3,0.04)
ret, dst = cv2.threshold(dst,0.1*dst.max(),255,0)
dst = np.uint8(dst)
ret, labels, stats, centroids = cv2.connectedComponentsWithStats(dst)
criteria = (cv2.TERM_CRITERIA_EPS + cv2.TERM_CRITERIA_MAX_ITER, 100, 0.001)
corners1 = cv2.cornerSubPix(gray,np.float32(centroids),(5,5),(-1,-1),criteria)
for i in range(1, len(corners1)):
    print(corners1[i])
    
image1[dst>0.1*dst.max()]=[0,0,255]
plt.imshow(cv2.cvtColor(image1, cv2.COLOR_BGR2RGB))
 
#Нахожу угловые точки на втором
gray = cv2.cvtColor(image2, cv2.COLOR_BGR2GRAY)
gray = np.float32(gray)
dst = cv2.cornerHarris(gray,5,3,0.04)
ret, dst = cv2.threshold(dst,0.1*dst.max(),255,0)
dst = np.uint8(dst)
ret, labels, stats, centroids = cv2.connectedComponentsWithStats(dst)
criteria = (cv2.TERM_CRITERIA_EPS + cv2.TERM_CRITERIA_MAX_ITER, 100, 0.001)
corners2 = cv2.cornerSubPix(gray,np.float32(centroids),(5,5),(-1,-1),criteria)
for i in range(1, len(corners2)):
    print(corners2[i])
    
image2[dst>0.1*dst.max()]=[0,0,255]
plt.imshow(cv2.cvtColor(image2, cv2.COLOR_BGR2RGB))
 
#Создаю HOG описания для каждой угловой точки и вывожу их
hog = cv2.HOGDescriptor()
gray = cv2.cvtColor(image1, cv2.COLOR_BGR2GRAY)
winStride=(8, 8)
padding=(4, 4)
locations = arr1[1]
fv1 = np.zeros((arr1.shape[0], 3780, 1), 'float')
print('fv1')
for i in range (0, len(corners1)):
    fv1[i, :, :] = hog.compute(gray, winStride, padding, (arr1[i],))
    print(fv1[i])
gray = cv2.cvtColor(image2, cv2.COLOR_BGR2GRAY)
fv2 = np.zeros((arr1.shape[0], 3780, 1), 'float')
print('fv2')
for i in range (0, len(corners2)):
    fv2[i, :, :] = hog.compute(gray, winStride, padding, (arr2[i],))
    print(fv2[i])