在Jupyter Notebook上使用Python实现下述代码的边缘提取、图像滤波功能,这个过程中实现某些功能处理出来的图像可能会有点粗糙。关于opencv库的安装可以参考:Python下opencv库的安装过程与一些问题汇总。
一、边缘提取
import cv2#导入opencv库
#读取原灰度图片
image=cv2.imread("test.bmp")
cv2.imshow("image", image)#将原图片命名为“image”显示出来
#图像的阈值分割处理,即将图像处理成非黑即白的二值图像
ret,image1 = cv2.threshold(image, 80, 255, cv2.THRESH_BINARY) # binary(黑白二值),ret代表阈值,80是低阈值,255是高阈值
cv2.imshow('image1', image1)#将阈值分割处理后的图片命名为“image1”显示出来
#二值图像的反色处理,将图片像素取反
height,width = image1.shape #返回图片大小
image2 = image1.copy()
for i in range(height):
for j in range(width):
image2[i,j] = (255-image1[i,j])
cv2.imshow('image2', image2)#将反色处理后的图片命名为“image2”显示出来
#边缘提取,使用Canny函数
image2_3 = cv2.Canny(image2,80,255) #设置80为低阈值,255为高阈值
cv2.imshow('image2_3', image2_3)#将边缘提取后的图片命名为“image2_3”显示出来
#再次对图像进行反色处理使提取的边缘为黑色,其余部分为白色,方法同image2
height1,width1 = image2_3.shape
image3 = image2_3.copy()
for i in range(height1):
for j in range(width1):
image3[i,j] = (255-image2_3[i,j])
cv2.imshow('image3', image3)#将边缘提取后反色处理的图片命名为“image3”显示出来
cv2.waitKey(0)# 等待键盘输入,不输入则无限等待
cv2.destroyAllWindows() #销毁所有窗口
运行结果:
二、图像滤波
import cv2
import numpy
import matplotlib.pyplot as plt
from PIL import Image, ImageDraw, ImageFont
#用于给图片添加中文字符
def ImgText_CN(img, text, left, top, textColor=(0, 255, 0), textSize=20):
if (isinstance(img, numpy.ndarray)): #判断是否为OpenCV图片类型
img = Image.fromarray(cv2.cvtColor(img, cv2.COLOR_BGR2RGB))
draw = ImageDraw.Draw(img)
fontText = ImageFont.truetype("font/simhei.ttf", textSize, encoding="utf-8")
draw.text((left, top), text, textColor, font=fontText)
return cv2.cvtColor(numpy.asarray(img), cv2.COLOR_RGB2BGR)
#读取原灰度图片
src=cv2.imread("test2.bmp")
cv2.imshow("src", src)#将原图片命名为“src”显示出来
src_cn=ImgText_CN(src, "原图像", 5, 5, (255, 0, 0), 20)
cv2.imshow("src_cn", src_cn)#将添加文字的src命名为“src_cn”并显示出来
# 均值滤波
des = cv2.blur(src, (5,5))
cv2.imshow("des", des)#将均值滤波处理后的src命名为“des”并显示出来
des_cn=ImgText_CN(des, "均值滤波", 5, 5, (255, 0, 0), 20)
cv2.imshow("des_cn", des_cn)#将添加文字的des命名为“des_cn”并显示出来
# 中值滤波
med = cv2.medianBlur(src, 5)
cv2.imshow("med", med)#将中值滤波处理后的src命名为“med”并显示出来
med_cn=ImgText_CN(med, "中值滤波", 5, 5, (255, 0, 0), 20)
cv2.imshow("med_cn", med_cn)#将添加文字的med命名为“med_cn”并显示出来
# 高斯滤波
gauss = cv2.GaussianBlur(src,(5,5),0)
cv2.imshow("gauss", gauss)#将高斯滤波处理后的src命名为“gauss”并显示出来
gauss_cn=ImgText_CN(gauss, "高斯滤波", 5, 5, (255, 0, 0), 20)
cv2.imshow("gauss_cn", gauss_cn)#将添加文字的gauss命名为“gauss_cn”并显示出来
# 高斯边缘检测
gaussedge = cv2.Canny(gauss,0,50)
cv2.imshow("gaussedge", gaussedge)#将边缘检测处理后的gauss命名为“gaussedge”并显示出来
gaussedge_cn=ImgText_CN(gaussedge, "高斯边缘检测", 5, 5, (0, 255, 0), 20)
cv2.imshow("gaussedge_cn", gaussedge_cn)#将添加文字的gaussedge命名为“gaussedge_cn”并显示出来
cv2.waitKey(0)
运行结果: