Python边缘检测之prewitt,sobel和laplace算子详解
作者：微小冷 发布时间：2023-09-30 06:01:56
标签：Python,边缘,检测
滤波算子简介
ndimage中提供了卷积算法，并且建立在卷积之上，提供了三种边缘检测的滤波方案：prewitt, sobel以及laplace。
在convolve中列举了一个用于边缘检测的滤波算子，统一维度后，其x xx和y yy向的梯度算子分别写为
此即prewitt算子。
Sobel算子为Prewitt增添了中心值的权重，记为
这两种边缘检测算子，均适用于某一个方向，ndimage还提供了lapace算子，其本质是二阶微分算子，其3×3卷积模板可表示为
具体实现
ndimage封装的这三种卷积滤波算法，定义如下
prewitt(input, axis=-1, output=None, mode='reflect', cval=0.0)
sobel(input, axis=-1, output=None, mode='reflect', cval=0.0)
laplace(input, output=None, mode='reflect', cval=0.0)
其中，mode表示卷积过程中对边缘效应的弥补方案，设待滤波数组为a b c d，则在不同的模式下，对边缘进行如下填充
| 左侧填充 | 数据 | 右侧填充 | |
|---|---|---|---|
| reflect | d c b a | a b c d | d c b a |
| constant | k k k k | a b c d | k k k k |
| nearest | a a a a | a b c d | d d d d |
| mirror | d c b | a b c d | c b a |
| wrap | a b c d | a b c d | a b c d |
测试
接下来测试一下
from scipy.ndimage import prewitt, sobel, laplace
from scipy.misc import ascent
import matplotlib.pyplot as plt
img = ascent()
dct = {
"origin" : lambda img:img,
"prewitt" : prewitt,
"sobel" : sobel,
"laplace" : lambda img : abs(laplace(img))
}
fig = plt.figure()
for i,key in enumerate(dct):
ax = fig.add_subplot(2,2,i+1)
ax.imshow(dct[key](img), cmap=plt.cm.gray)
plt.ylabel(key)
plt.show()
为了看上去更加简洁，代码中将原图、prewitt滤波、sobel滤波以及laplace滤波封装在了一个字典中。其中origin表示原始图像，对应的函数是一个lambda表达式。
在绘图时，通过将cmap映射到plt.cm.gray，使得绘图之后表现为灰度图像。
效果如下
来源：https://blog.csdn.net/m0_37816922/article/details/130218971
0
投稿猜你喜欢
- HttpRequest.FILES表单上传的文件对象存储在类字典对象request.FILES中，表单格式需为multipart/form-
- 本文实例为大家分享了Python基于Socket实现简单聊天室，供大家参考，具体内容如下服务端#!/usr/bin/env python#
前言在设计爬虫项目的时候，首先要在脑内明确人工浏览页面获得图片时的步骤一般地，我们去网上批量打开壁纸的时候一般操作如下：1、打开壁纸网页2、- 安装cuda更新nvidia驱动打开GeForce Game Ready Driver或在GeForce Experience中下载符合自己
- Doug Bowman，Google的Visual Design Lead离职了，一封带有感 * 彩的离职信惹发了大家不少的讨论。甚至还有人用
- 大家先看一下Python os模块中的部分函数python 路径相关的函数 os.listdir(dirname)：列出dirname下的目
- 先说下自己之前的环境（都是Linux系统，差别不大）：Centos7.6NVIDIA Driver Version 440.33.01（等会
- 为了检验自己前期对机器学习中线性回归部分的掌握程度并找出自己在学习中存在的问题，我使用C语言简单实现了单变量简单线性回归。本文对自己使用C语
- Python的matplotlib包可以轻松的将数据可视化，博主最近遇到了一个问题，博主想同时在两个窗口展示两张图，但是代码运行结果总是显示
- php mysql获取表字段名称和字段信息的三种方法先给出本实例中使用的表的信息：使用desc获取表字段信息php代码如下：<?php
- 在IE浏览器调试代码，我们可以选择使用 IE WebDeveloper但是我个人用惯了ff浏览器下的firebug，所以在网上搜了一下，如果
- <html><head><title>遍历表格</title><script lang
- 给静态网页加密的方法有很多，有的简单有的复杂。前两天看见有人问静态网页加密问题，就写了这个代码思路：加密时：先把用户的密钥A用md5加密为B
- python的应用实践zipkin，需要py_zipkin，使用pip进行安装py_zipkin的时候出现问题, 根据stackoverfl
- 在Qtdesigner中新建一个主界面如下所示：ctrl+R 预览从预览图中可以看出这时的界面不支持伸缩，拖动过小的话会导致部分界面遮住不可
- 引言所有的层都具有的参数，如name, type, bottom, top和transform_param请参看我的前一篇文章：Caffe卷
- 由于不同的项目需要用不同的python版本，于是使用Anaconda来进行版本管理，现记录一下经验：在官网下载并安装好Anaconda以后（
- 1.安装Pillowpip install Pillow2.安装tesseract-ocrgithub地址: https://gi
- 本文介绍用python实现的搜索本地文本文件内容的小程序。从而学习Python I/O方面的知识。代码如下：import os#根据文件扩展
- 前言：有读者可能会一脸懵？啥是索引潜水？你给起的名字的吗？有没有索引蛙泳？这个名字还真不是我起的，今天要讲的知识点就叫索引潜水（Index
