手写代码实现卷积操作(Python)

手写代码实现卷积操作(Python)在卷积神经网络中,才用卷积技术实现对图片的降噪和特征提取。一般我们构建卷积神经网络都是使用成熟的框架,今天我就来自己实现一下卷积,并使用不同的卷积核来看看效果。卷积操作的原理可以由下图表示:一个3*3的卷积核,以滑动窗口的形式在图片上滑动,每滑动一次,就计算窗口中的数据的加权之和,权值就是卷积核的数据。通过这个过程将图片进行转化。准备图片数据:使用P…

大家好,又见面了,我是你们的朋友全栈君。

在卷积神经网络中,才用卷积技术实现对图片的降噪和特征提取。

一般我们构建卷积神经网络都是使用成熟的框架,今天我就来自己实现一下卷积,并使用不同的卷积核来看看效果。

 

卷积操作的原理可以由下图表示:

手写代码实现卷积操作(Python)

 一个3*3的卷积核,以滑动窗口的形式在图片上滑动,每滑动一次,就计算窗口中的数据的加权之和,权值就是卷积核的数据。

通过这个过程将图片进行转化。

 

准备图片数据:

手写代码实现卷积操作(Python)

 使用PIL库读取图片:

from PIL import Image
import matplotlib.pyplot as plt
import numpy as np

读取图片,并可视化:

img = Image.open('金鱼姬.jpg')
plt.axis("off")
plt.imshow(img)

手写代码实现卷积操作(Python)

 以灰度显示图片:

gray = img.convert('L')
plt.figure()
plt.imshow(gray, cmap='gray')
plt.axis('off')

 

手写代码实现卷积操作(Python)

 RGB三通道数据分离:

r, g, b = img.split()

        三色通道可视化:

plt.imshow(r, cmap='gray')
plt.axis('off')

手写代码实现卷积操作(Python)手写代码实现卷积操作(Python)手写代码实现卷积操作(Python)

 将图片数据转化为np矩阵:

np.array(img)
np.array(r)
np.array(g)
np.array(b)

单通道数据(r):

array([[ 52,  58,  55, ..., 139, 133,  91],
       [ 46,  53,  52, ..., 140, 138, 103],
       [ 45,  53,  52, ..., 142, 144, 119],
       ...,
       [ 40,  37,  37, ...,  60,  61,  47],
       [ 39,  36,  37, ...,  56,  58,  48],
       [ 35,  33,  35, ...,  54,  57,  48]], dtype=uint8)

定义卷积核:

#定义卷积核(3*3)
k = np.array([
    [0,1,2],
    [2,2,0],
    [0,1,2]
])

我们使用单通道图片数据进行卷积:

data = np.array(r)
n,m = data.shape

       可以看到,图片的大小为:720*1280

 

关键点来了,定义卷积函数(二维):

def convolution(k, data):
    n,m = data.shape
    img_new = []
    for i in range(n-3):
        line = []
        for j in range(m-3):
            a = data[i:i+3,j:j+3]
            line.append(np.sum(np.multiply(k, a)))
        img_new.append(line)
    return np.array(img_new)

采用双循环遍历图片数据,line存储一行卷积后的数据,img_new记录每一行卷积后的数据,形成新的图片

 

下面来看一看卷积结果:

img_new = convolution(k, data)#卷积过程

#卷积结果可视化
plt.imshow(img_new, cmap='gray')
plt.axis('off')

手写代码实现卷积操作(Python)

 因为卷积结果和卷积核有密切关系,虽然上图看上去没他大的区别,那是因为卷积核的关系。

 

我们来使用常见的几种卷积核来看看结果:

卷积核1:垂直边缘检测

k1 = np.array([
    [1,0,-1],
    [1,0,-1],
    [1,0,-1]
])

手写代码实现卷积操作(Python)

可以看到,图片中物体的纵向边界特征被提取出来,而其他的特征被弱化了

 

卷积核2:水平边缘检测

k2 = np.array([
    [1,1,1],
    [0,0,0],
    [-1,-1,-1]
])

手写代码实现卷积操作(Python)

可以看出,横向特征被凸显。

我们换一个,纵向特征比较明显的图片来试一下:

手写代码实现卷积操作(Python)

 例如该图中有很多竖向的扶手,我们看看他的卷积后的效果:

垂直边缘检测:

手写代码实现卷积操作(Python)

水平边缘检测:

手写代码实现卷积操作(Python)

这样看,特征就很明显了。

通过卷积后的图片,大小变化了,可以使用零填充技术使得输出图片和原图片大小一致,这里就不实现了,本文主要是实现卷积过程。

还有其他常用的卷积核,以下分别展示一下各自的效果:

 

sobel滤波(纵向,横向):

手写代码实现卷积操作(Python)手写代码实现卷积操作(Python)

scharr滤波(纵向,横向):

手写代码实现卷积操作(Python)手写代码实现卷积操作(Python)

 大家猜猜,最后我使用的是出自哪里的图片吧!

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/140232.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • postgresal去重_postgresql数据库去重方法

    postgresal去重_postgresql数据库去重方法数据库去重有很多方法,下面列出目前理解与使用的方法第一种通过groupby分组,然后将分组后的数据写入临时表然后再写入另外的表,对于没有出现再groupby后面的field可以用函数max,min提取,效率较高–适合情况:这种情况适合重复率非常高的情况,一般来说重复率超过5成则可以考虑用这个方法–优点:对于重复率高的数据集的去重,十分推荐用这种方法–缺点:uuid不能用max或min提取…

    2022年10月1日
    0
  • 如何卸载赛门铁克(Symantec)企业防病毒客户端软件SEP(Symantec Endpoint Protection)?[通俗易懂]

    如何卸载赛门铁克(Symantec)企业防病毒客户端软件SEP(Symantec Endpoint Protection)?[通俗易懂]本文分三个部分介绍如何卸载赛门铁克(Symantec)企业防病毒客户端软件SEP的Windows版本,Linux版本和MacOS版本。1卸载SEPforWindows版本SEPforWindows版本卸载使用传统办法–通过Windows控制面板卸载是比较容易出问题的,常见问题如卸载过程发生回滚而导致卸载失败,卸载后存在有文件残余导致无法重装SEP或者第三方厂商防病毒软件等问题。最佳的解决办法是使用官方的卸载工具CleanWipe进行卸载。由于CSDN上CleanWipe工具

    2022年6月5日
    334
  • 浅析Anycast技术[通俗易懂]

    浅析Anycast技术[通俗易懂]什么是AS号码AS号码即自治系统号码,是用来标识独立的自治系统的,在同一个自治系统内,使用相同内部路由协议,自治系统间使用外部路由协议(通常是BGP协议)。申请AS号码的单位需要与两家以上(包括两家)、有不同AS号码的网络接入商进行网络互联,并计划三个月内与他们同时运行BGP协议进行外部路由。什么是BGPAnyCast?BGPanycast就是利用一个(多个)as号码在不同的地区广播相同的一个ip段。利用bgp的寻路原则,短的aspath会选成最优路径(bgp寻路原则之n),从.

    2022年5月10日
    58
  • QThread源码浅析[通俗易懂]

    QThread源码浅析[通俗易懂]Qt版本Qt5.6.0,下面以Windows平台为例简单研究下QThread源码实现。1.仅研究下QThread::start()函数,其他细节在次不涉及:src\qtbase\src\corelib\thread\qthread_win.cppvoidQThread::start(Prioritypriority){Q_D(QThread);QMutexLocker…

    2022年5月28日
    73
  • 5.1活动规划_上海市城市规划条例

    5.1活动规划_上海市城市规划条例规划范围管理是为记录如何定义,确认和控制项目范围及产品范围,而创建范围管理计划的过程。本过程的主要作用是:在整个项目期间对如何管理项目提供指南和方向。范围管理计划是项目或项目管理计划的组成部分,描述将如何定义、制定、监督、控制和确认项目范围。 …

    2022年9月15日
    0
  • python安装步骤(pycharm运行python)

    文章目录一。pycharm下载安装二。python下载安装三.pycharm上配置python一。pycharm下载安装pycharm下载地址:http://www.jetbrains.com/pycharm/download/#section=windows下载详细步骤:1-2-3-4-5-67-8-直接finish二。python下载安装9-python官网:https://www.python.org/进去网址后点击:1011-下载好后12

    2022年4月10日
    155

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号