2011年08月23日

在 Python
中,比较常用的图像处理库是 PIL(Python Image Library),当前版本是 1.1.6 ,用起来非常方便。大家可以在
http://www.pythonware.com/products/pil/index.htm
下载和学习。

基本图像处理

使用 PIL 之前需要 import Image 模块:

import Image

然后你就可以使用Image.open(‘xx.bmp’)
来打开一个位图文件进行处理了。打开文件你不用担心格式,也不用了解格式,无论什么格式,都只要把文件名丢给 Image.open
就可以了。真所谓 bmp、jpg、png、gif……,一个都不能少。

img =
Image.open(‘origin.png’)   
# 得到一个图像的实例对象 img

图 1原图

图像处理中,最基本的就是色彩空间的转换。一般而言,我们的图像都是 RGB
色彩空间的,但在图像识别当中,我们可能需要转换图像到灰度图、二值图等不同的色彩空间。 PIL
在这方面也提供了极完备的支持,我们可以:

new_img = img.convert(‘L’)

把 img 转换为 256 级灰度图像, convert() 是图像实例对象的一个方法,接受一个 mode
参数,用以指定一种色彩模式,mode 的取值可以是如下几种:

· 1 (1-bit pixels, black and white, stored with one pixel per
byte)

· L (8-bit pixels, black and white)

· P (8-bit pixels, mapped to any other mode using a colour
palette)

· RGB (3x8-bit pixels, true colour)

· RGBA (4x8-bit pixels, true colour with transparency
mask)

· CMYK (4x8-bit pixels, colour separation)

· YCbCr (3x8-bit pixels, colour video format)

· I (32-bit signed integer pixels)

· F (32-bit floating point pixels)

怎么样,够丰富吧?其实如此之处,PIL 还有限制地支持以下几种比较少见的色彩模式:LA (L with alpha),
RGBX (true colour with padding) and RGBa (true colour with
premultiplied alpha)。

下面看一下 mode 为 ‘1’、’L’、’P’时转换出来的图像:

图 2 mode = ‘1‘

图 3 mode = ‘L‘

图 4 mode = ‘P‘

convert() 函数也接受另一个隐含参数 matrix,转换矩阵 matrix 是一个长度为4 或者16
tuple。下例是一个转换 RGB 空间到 CIE XYZ 空间的例子:

rgb2xyz = (

0.412453, 0.357580, 0.180423, 0,

0.212671, 0.715160, 0.072169, 0,

0.019334, 0.119193, 0.950227, 0 )

out =
im.convert("RGB", rgb2xyz)

除了完备的色彩空间转换能力外, PIL
还提供了resize()、rotate()等函数以获得改变大小,旋转图片等几何变换能力,在图像识别方面,图像实例提供了一个
histogram() 方法来计算直方图,非常方便实用。

图像增强

图像增强通常用以图像识别之前的预处理,适当的图像增强能够使得识别过程达到事半功倍的效果。 PIL 在这方面提供了一个名为
ImageEnhance 的模块,提供了几种常见的图像增强方案:

import ImageEnhance

enhancer =
ImageEnhance.Sharpness(image)

for i in range(8):

factor
= i / 4.0

enhancer.enhance(factor).show("Sharpness %f" % factor)

上面的代码即是一个典型的使用 ImageEnhance 模块的例子。 Sharpness 是 ImageEnhance
模块的一个类,用以锐化图片。这一模块主要包含如下几个类:Color、Brightness、Contrast和Sharpness。它们都有一个共同的接口
.enhance(factor) ,接受一个浮点参数 factor,标示增强的比例。下面看看这四个类在不同的 factor
下的效果

图 5 使用Color 进行色彩增强,factor 取值 [0, 4],步进 0.5

图 6 用 Birghtness 增强亮度,factor取值[0,4],步进0.5

图 7用 Contrast 增强对比度, factor 取值 [0,4],步进0.5

图 8用 Sharpness 锐化图像,factor取值 [0,4],步进0.5

图像 Filter

PIL 在 Filter
方面的支持是非常完备的,除常见的模糊、浮雕、轮廓、边缘增强和平滑,还有中值滤波、ModeFilter等,简直方便到可以做自己做一个Photoshop。这些
Filter 都放置在 ImageFilter 模块中,ImageFilter主要包括两部分内容,一是内置的 Filter,如
BLUR、DETAIL等,另一部分是 Filter 函数,可以指定不同的参数获得不同的效果。示例如下:

import ImageFilter

im1 =
im.filter(ImageFilter.BLUR)

im2 =
im.filter(ImageFilter.MinFilter(3))

im3 =
im.filter(ImageFilter.MinFilter()) # same as MinFilter(3)

可以看到 ImageFilter 模块的使用非常简单,每一个 Filter
都只需要一行代码就可调用,开发效率非常高。

图 9使用 BLUR

图 10使用 CONTOUR

图 11使用 DETAIL

图 12使用 EMBOSS

图 13使用 EDGE_ENHANCE

图 14使用 EDGE_ENHANCE_MORE

图 15使用 FIND_EDGES

图 16使用 SHARPEN

图 17使用 SMOOTH

图 18使用 SMOOTH_MORE

以上是几种内置的 Filter 的效果图,除此之外, ImageFilter 还提供了一些 Filter
函数,下面我们来看看这些可以通过参数改变行为的 Filter 的效果:

图 19使用 Kernel(),参数:size = (3, 3), kernel = (0.5,
0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5)

图 20使用 MaxFilter,默认参数

图 21使用 MinFilter,默认参数

图 22使用 MedianFilter,默认参数

图 23使用 ModeFilter,参数 size = 3

图 24使用 RankFilter,参数 size = 3, rank = 3

小结

到此,对 PIL 的介绍就告一段落了。总的来说,对于图像处理和识别,PIL 内建了强大的支持,从各种增强算法到 Filter
,都让人无法怀疑使用 Python 的可行性。
Python唯一的劣势在于执行时间过慢,特别是当实现一些计算量大的算法时候,需要极强的耐心。我曾用 Hough
Transform(霍夫变换)来查找图像中的直线,纯 Python 的实现处理一个 340 * 100 的图片也要花去数秒时间(P4
3.0G + 1G memory)。但使用 PIL 无需关注图像格式、内建的图像增强算法和 Filter 算法,这些优点使
Python 适合用于构造原型和进行实验,在这两方面Python 比 matlab 更加方便。商业的图像识别产品开发,可以考虑已经被
boost accepted的来自 adobe 的开源 C++ 库 gil,可以兼顾执行性能和开发效率。

转自:http://blog.csdn.net/lanphaday/article/details/1852726

时间: 2024-10-13 17:31:02

2011年08月23日的相关文章

2014年4月23日 10:22:08

step 1 : 做tcp网络编程,要解析一批批的数据,可是数据是通过Socket连接的InputStream一次次读取的,读取到的不是需要转换的对象,而是要直接根据字节流和协议来生成自己的数据对象. 按照之前的编程思维,总是请求然后响应,当然Socket也是请求和响应,不过与单纯的请求响应是不同的. 这里Socket连接往往是要保持住的,也就是长连接,然后设置一个缓冲区,网络流不断的追加到缓冲区.然后后台去解析缓冲区的字节流. http://cuisuqiang.iteye.com/blog/

8月23日28家中国域名商六类国际域名注册保有量统计

IDC评述网(idcps.com)08月27日报道:根据RegistrarStats公布的实时数据显示,截止至2015年8月23日,国内域名注册保有量二十八强排名情况,与上期对比,发生变化.环球互易实现超越,升至第24名,被其取代的广州名扬信息(琥珀网)则跌至第25名,排名互换.接下来,请看IDC评述网整理的具体数据分析. (图1)中国域名注册商(国际域名)保有量TOP15分布图 如图1所示,截止至8月23日,国内域名注册保有量十五强分别是阿里云(万网) .新网.易名中国.商务中国.西部数码.爱

七周二次课(1月23日) 10.6 监控io性能 10.7 free命令 10.8 ps命令 10.9 查看网络状态 10.10 linux下抓包

七周二次课(1月23日)10.6 监控io性能10.7 free命令10.8 ps命令10.9 查看网络状态10.10 linux下抓包 =====================================================================================================================================================================================

十三周五次课(3月23日)

十三周五次课(3月23日)13.4 mysql用户管理 grant 授权  all  所有的权限  *.* 所有的库和表 @ 限定ip %通配所有的ip identified by"123456a" 密码 查看指定用户的授权 13.5 常用sql语句 查 insert 插入 update 更改 truncate 清空内容,结构还在 drop 删 13.6 mysql数据库备份恢复 备份恢复库 备份恢复表 备份所有库 只备份表结构 扩展SQL语句教程  http://www.runoob

软考中高项学员:2016年3月23日作业

软考中高项学员:2016年3月23日作业 一.项目进度管理1.进度管理包括哪六个过程?2.什么是滚动式规划?3.什么是控制账户.规划组合?4.请说明FS.FF.SS.SF的含义.5.虚活动的含义?6.三种依赖关系是哪三种?7.活动资源估算的方法.工具和技术?(记)8.活动历时估算的方法.工具和技术?(记)9.何时可以用类比估算?10.参数估算的含义?11.三点估算公式?标准差计算公式?正负一倍的标准差情况下,概率是多少?2倍的呢?3倍的呢?12.后备分析中,应急时间=时间储备=缓冲时间吗?13.

阿西莫夫 - 神们自己(2015年6月23日)

<神们自己> 作 者:阿西莫夫译 者:崔正男系 列:出 版:江苏凤凰文艺出版社字 数:230千字阅读完成:2015年6月23日

第17期中国智能家居主题沙龙将于5月23日在京举行

中国智能家居主题沙龙已经办16期,为中国智能家居的发展贡献了自己的力量,在这轮智能家居热潮中,沙龙将继续为政府.投资机投.智能家居产品提供商.智能家居方案提供商,智能家居集成商,智能家居渠道提供多方位的合作,促进行业的快速爆发. 活动宗旨: 搭建行业用户与企业沟通的桥梁,以交流.合作.服务为目的,探讨行业现状及走势,推广普及新理念.建立智能家居沟通平台以宣传企业品牌,推广新技术,展示新产品.根据智能家居行业各产业链上下游人士的需求,有针对性地提供服务,促进行业发展. 活动背景: 2013年初<国

高级:2015年10月23日作业答案

高级学员:2015年10月23日作业一.业务流程管理与重组1.流程管理的本质与核心是什么: 流程管理的核心是流程,流程管理的本质就是构造卓越的业务流程.2.流程管理与BPR的区别: 流程管理是以一种规范化的构造端到端的卓越业务流程为中心,以持续的提高组织业务绩效为目的的系统化方法.与业务流程重组的定义相似.流程管理的思想包含了BPR,但比BPR的概念更广泛,更适合现实需要.3.业务流程改进的过程和方法中,分析问题用哪2个图:寻找原因用什么图. 分析问题采用头脑风暴法.矩阵图法,寻找原因用鱼骨图法

高级学员:2015年10月23日作业

高级学员:2015年10月23日作业一.业务流程管理与重组1.流程管理的本质与核心是什么:2.流程管理与BPR的区别:3.业务流程改进的过程和方法中,分析问题用哪2个图:寻找原因用什么图.4.流程的分析与设计方法是哪六条(只答小标题):5.BPR的首要问题是什么:6.BPR的实施会引起企业哪3个方面的变化.二.知识管理1.对于项目组织来说,什么是知识.2.隐性知识的共享方法有哪4条:三.项目整体绩效评估1.什么是3E审计:2.霍尔从哪3个维度考察系统工程.3.用什么分别解决:确定性.可测定的不确