图像导数

最新推荐文章于 2023-09-16 13:36:14 发布

TurnaroundBin

最新推荐文章于 2023-09-16 13:36:14 发布

阅读量749

点赞数

分类专栏：计算机视觉文章标签： python 计算机视觉

本文链接：https://blog.csdn.net/qq_31714485/article/details/70918512

版权

计算机视觉专栏收录该内容

0 篇文章 0 订阅

订阅专栏

图像导数

/*
python计算机视觉——图像导数
2017.04.28
*/

from PIL import Image
from numpy import *
from scipy.ndimage import filters


imo = Image.open(r'C:/Users/BInZhuang/Desktop/1.jpg')
im = array(Image.open(r'C:/Users/BInZhuang/Desktop/1.jpg').convert('L'))

sigma = 5

imx = zeros(im.shape)
filters.gaussian_filter(im,(sigma,sigma),(0,1),imx)

imy = zeros(im.shape)
filters.gaussian_filter(im,(sigma,sigma),(1,0),imy)

magnitude = sqrt(imx**2 + imy**2)

theta = (cos(imy/imx)/sin(imy/imx))

imsave('im.jpg',im)
imsave('imx.jpg',imx)
imsave('imy.jpg',imy)
imsave('imxy2.jpg',magnitude)
imsave('imxy.jpg',theta)


figure()
gray()

subplot(1,5,1)
imshow(im)
axis('equal')
axis('off')

subplot(1,5,2)
imshow(imx)
axis('equal')
axis('off')

subplot(1,5,3)
imshow(imy)
axis('equal')
axis('off')

subplot(1,5,4)
imshow(magnitude)
axis('equal')
axis('off')

subplot(1,5,5)
imshow(theta)
axis('equal')
axis('off')
show()