harris位_【论文阅读】Harris角点算法

#coding=utf-8

from PIL import Image

import numpy as np

from scipy.ndimage import filters

import matplotlib.pyplot as plt

import scipy.signal

def mean_filter2d(arr):

n = 3

# 3*3 滤波器, 每个系数都是 1/9

w = np.ones((n, n)) / n ** 2

# 使用滤波器卷积图像

# mode = same 表示输出尺寸等于输入尺寸

# boundary 表示采用对称边界条件处理图像边缘

s = scipy.signal.convolve2d(arr, w, mode='same', boundary='symm')

return s

def harris_response1(im,sigma=1.1):

"""计算图像的harris响应函数"""

#计算导数

imx = np.zeros(im.shape)

imx = scipy.ndimage.sobel(im,axis=0,mode='reflect')

# imx = filters.gaussia