theano中对图像进行convolution 运算

2024-04-10 08:18

本文主要是介绍theano中对图像进行convolution 运算,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

(1) 定义计算过程中需要的symbolic expression

 1 """
 2 定义相关的symbolic experssion
 3 """
 4 # convolution layer的输入,根据theano,它应该是一个4d tensor
 5 input = T.tensor4(name='input')
 6 # 共享权值W,它的shape为2,3,9,9
 7 w_shp = (2,3,9,9);w_bound = numpy.sqrt(3*9*9)
 8 W = theano.shared(numpy.asarray(rng.uniform(low= -1.0/w_bound, high = 1.0/w_bound,size=w_shp),dtype=input.dtype),name='W')
 9 # 利用卷积核W对input进行卷积运算
10 conv_out = conv.conv2d(input,W)
11 # 偏执向量b
12 b_shp = (2,)  # b是一个只有1个元素2的tuple
13 b = theano.shared(numpy.asarray(rng.uniform(low= -.5, high = .5,size=b_shp),dtype=input.dtype),name='b')
14 # 计算sigmoid函数
15 output = T.nnet.sigmoid(conv_out+b.dimshuffle('x',0,'x','x'))
16 # 输入输出function
17 f = theano.function([input],output)

 

 

(2)利用真实数据计算

 1 """
 2 开始使用具体数值
 3 """
 4 # 读入图像
 5 img = Image.open('3wolfmoon.jpg', mode='r')
 6 # 将输入图像存入在array中
 7 img = numpy.array(img,dtype='float64')/256
 8 # 对输入图像进行reshape
 9 img_=img.transpose(2,0,1).reshape(1,3,639,516)
10 # 利用convolution kernel对输入图像进行卷积运算
11 filtered_img=f(img_)

(3)绘制需要显示的图像

 1 """
 2 绘制图像
 3 """
 4 # 显示原始图像
 5 pylab.subplot(1,3,1);pylab.axis('off');pylab.imshow(img);pylab.gray()
 6 # 显示filter后的图像的channel1
 7 pylab.subplot(1,3,2);pylab.axis('off');pylab.imshow(filtered_img[0,0,:,:])
 8 # 显示filter后的图像的channel2
 9 pylab.subplot(1,3,3);pylab.axis('off');pylab.imshow(filtered_img[0,1,:,:])
10 # 显示
11 pylab.show()

 

整个代码段

 1 # -*- coding: utf-8 -*-
 2 
 3 # 导入相关的模块
 4 import theano
 5 from theano import tensor as T
 6 from theano.tensor.nnet import conv
 7 import numpy
 8 import pylab
 9 from PIL import Image
10 
11 
12 # 产生随机数的种子
13 rng = numpy.random.RandomState(23455)
14 
15 """
16 定义相关的symbolic experssion
17 """
18 # convolution layer的输入,根据theano,它应该是一个4d tensor
19 input = T.tensor4(name='input')
20 # 共享权值W,它的shape为2,3,9,9
21 w_shp = (2,3,9,9);w_bound = numpy.sqrt(3*9*9)
22 W = theano.shared(numpy.asarray(rng.uniform(low= -1.0/w_bound, high = 1.0/w_bound,size=w_shp),dtype=input.dtype),name='W')
23 # 利用卷积核W对input进行卷积运算
24 conv_out = conv.conv2d(input,W)
25 # 偏执向量b
26 b_shp = (2,)  # b是一个只有1个元素2的tuple
27 b = theano.shared(numpy.asarray(rng.uniform(low= -.5, high = .5,size=b_shp),dtype=input.dtype),name='b')
28 # 计算sigmoid函数
29 output = T.nnet.sigmoid(conv_out+b.dimshuffle('x',0,'x','x'))
30 # 输入输出function
31 f = theano.function([input],output)
32 
33 """
34 开始使用具体数值
35 """
36 # 读入图像
37 img = Image.open('3wolfmoon.jpg', mode='r')
38 # 将输入图像存入在array中
39 img = numpy.array(img,dtype='float64')/256
40 # 对输入图像进行reshape
41 img_=img.transpose(2,0,1).reshape(1,3,639,516)
42 # 利用convolution kernel对输入图像进行卷积运算
43 filtered_img=f(img_)
44 
45 """
46 绘制图像
47 """
48 # 显示原始图像
49 pylab.subplot(1,3,1);pylab.axis('off');pylab.imshow(img);pylab.gray()
50 # 显示filter后的图像的channel1
51 pylab.subplot(1,3,2);pylab.axis('off');pylab.imshow(filtered_img[0,0,:,:])
52 # 显示filter后的图像的channel2
53 pylab.subplot(1,3,3);pylab.axis('off');pylab.imshow(filtered_img[0,1,:,:])
54 # 显示
55 pylab.show()
View Code

 

这篇关于theano中对图像进行convolution 运算的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!


原文地址:
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.chinasem.cn/article/890528

相关文章

Nginx中配置使用非默认80端口进行服务的完整指南

《Nginx中配置使用非默认80端口进行服务的完整指南》在实际生产环境中,我们经常需要将Nginx配置在其他端口上运行,本文将详细介绍如何在Nginx中配置使用非默认端口进行服务,希望对大家有所帮助... 目录一、为什么需要使用非默认端口二、配置Nginx使用非默认端口的基本方法2.1 修改listen指令

MySQL按时间维度对亿级数据表进行平滑分表

《MySQL按时间维度对亿级数据表进行平滑分表》本文将以一个真实的4亿数据表分表案例为基础,详细介绍如何在不影响线上业务的情况下,完成按时间维度分表的完整过程,感兴趣的小伙伴可以了解一下... 目录引言一、为什么我们需要分表1.1 单表数据量过大的问题1.2 分表方案选型二、分表前的准备工作2.1 数据评估

MySQL进行分片合并的实现步骤

《MySQL进行分片合并的实现步骤》分片合并是指在分布式数据库系统中,将不同分片上的查询结果进行整合,以获得完整的查询结果,下面就来具体介绍一下,感兴趣的可以了解一下... 目录环境准备项目依赖数据源配置分片上下文分片查询和合并代码实现1. 查询单条记录2. 跨分片查询和合并测试结论分片合并(Shardin

SpringBoot结合Knife4j进行API分组授权管理配置详解

《SpringBoot结合Knife4j进行API分组授权管理配置详解》在现代的微服务架构中,API文档和授权管理是不可或缺的一部分,本文将介绍如何在SpringBoot应用中集成Knife4j,并进... 目录环境准备配置 Swagger配置 Swagger OpenAPI自定义 Swagger UI 底

基于Python Playwright进行前端性能测试的脚本实现

《基于PythonPlaywright进行前端性能测试的脚本实现》在当今Web应用开发中,性能优化是提升用户体验的关键因素之一,本文将介绍如何使用Playwright构建一个自动化性能测试工具,希望... 目录引言工具概述整体架构核心实现解析1. 浏览器初始化2. 性能数据收集3. 资源分析4. 关键性能指

Nginx进行平滑升级的实战指南(不中断服务版本更新)

《Nginx进行平滑升级的实战指南(不中断服务版本更新)》Nginx的平滑升级(也称为热升级)是一种在不停止服务的情况下更新Nginx版本或添加模块的方法,这种升级方式确保了服务的高可用性,避免了因升... 目录一.下载并编译新版Nginx1.下载解压2.编译二.替换可执行文件,并平滑升级1.替换可执行文件

Python进行JSON和Excel文件转换处理指南

《Python进行JSON和Excel文件转换处理指南》在数据交换与系统集成中,JSON与Excel是两种极为常见的数据格式,本文将介绍如何使用Python实现将JSON转换为格式化的Excel文件,... 目录将 jsON 导入为格式化 Excel将 Excel 导出为结构化 JSON处理嵌套 JSON:

一文解密Python进行监控进程的黑科技

《一文解密Python进行监控进程的黑科技》在计算机系统管理和应用性能优化中,监控进程的CPU、内存和IO使用率是非常重要的任务,下面我们就来讲讲如何Python写一个简单使用的监控进程的工具吧... 目录准备工作监控CPU使用率监控内存使用率监控IO使用率小工具代码整合在计算机系统管理和应用性能优化中,监

如何使用Lombok进行spring 注入

《如何使用Lombok进行spring注入》本文介绍如何用Lombok简化Spring注入,推荐优先使用setter注入,通过注解自动生成getter/setter及构造器,减少冗余代码,提升开发效... Lombok为了开发环境简化代码,好处不用多说。spring 注入方式为2种,构造器注入和setter

基于Python开发一个图像水印批量添加工具

《基于Python开发一个图像水印批量添加工具》在当今数字化内容爆炸式增长的时代,图像版权保护已成为创作者和企业的核心需求,本方案将详细介绍一个基于PythonPIL库的工业级图像水印解决方案,有需要... 目录一、系统架构设计1.1 整体处理流程1.2 类结构设计(扩展版本)二、核心算法深入解析2.1 自