Mexican Hat Filter - Laplacian of Gaussian

Filtro de sobrero mexicano, es un filtro Laplaciano de un Gaussiano

 

 

Les comparto el código en Python2 para que dibujen.

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import axes3d

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
dx = 0.1
dy = 0.1
x = np.arange(-20, 20, dx)
y = np.arange(-20, 20, dy)
x2d, y2d = np.meshgrid(x,y)
sigma = 3
kernel_2d = (x2d**2 + y2d**2 - (2*sigma**2))/sigma**4
kernel_2d = kernel_2d * np.exp(-(x2d**2 + y2d**2)/(2*sigma**2))
ax.plot_surface(x2d, y2d, -kernel_2d,cmap='jet')
plt.show() 

Una imagen del código

Filtro de sombrero mexicano
Mexican Hat Filter