时域的信号你可以理解为 1111=(1+1+1++1)
傅立叶你可以理解为 1111=1000+100+10+1 就是可以由不同频率的信号表示
傅立叶是将一维的信号,变换成为 1维的另一种信号
小波你可以理解为 8421=(4+4)1000+(2+2)100+(1+1)10+1
小波是将一维的信号变换成为 几个 一维的信号。变成2维了
他们的本质都是将信号分解为不同的部分,这些部分也可以组成原本的信号。
只是坐标图的时候 时域信号和傅立叶的图是平面图,小波是 三维图
一、确定方法:
小波分析把图像分解为两部分:低频信息+高频信息低频信息是变化缓慢的部分,是图像的框架,也是轮廓,占全部信息的大部分;高频信息是变化迅速的部分(如从黑色跳变到白色),它反映的是图像的细节信息,占全部信息的小部分以上是第一层分解
在第一层的基础上把高频信息部分再分解为两部分:低频+高频第三层是把第二层分解出来的高频信息分解为低频+高频依次类推
以上是单尺度分解,低频部分不进行再分解若是多尺度分解,会把低频部分像高频部分那样一层一层分解
二、小波分解的概念:
[c,l] = wavedec(s,3,'db1');
l是length的意思,记录的是由高到低各级的长度。
s代表进行分解的变量;
3代表分解层数
对1张图象进行小波分解,可以在MATLAB中实现。在COMMAND WINDOWS窗口中直接输入wavedemo进入说明,wavemenu进使用程序,也可以直接编程。程序在wavedemo里面自带。
小波变换的二层分解是小波包分解;
小波分析的二维分解是平面分析,二元的分解而已,有很多区别,给你我做过的代码:
那么下面一段 Matlab程序便是对于小波熵的求法:
这是我本科毕业时候写的小波包的例子
%求小波包的信息熵,中南大学数学学院信科0302班张扬编写
%ECG为信号源,n为分解的层数,wpname为小波名字,我常用的是db3,haar,db4
function y=waveletentropy(ECG,n,wpname)
wpt1=wpdec(ECG,n,wpname); %对数据进行小波包分解
for i=1:2n %wpcoef(wpt1,[n,i-1])是求第n层第i个节点的系数
disp('每个节点的能量E(i)');
E(i)=norm(wpcoef(wpt1,[n,i-1]),2)norm(wpcoef(wpt1,[n,i-1]),2)
%求第i个节点的范数平方,其实也就是平方和
end
disp('小波包分解总能量E_total');
E_total=sum(E) %求总能量
disp('以下是每个节点的概率P');
for i=1:2n
p(i)= E(i)/E_total %求每个节点的概率
end
%以下计算小波熵,即-sum(pjlnpj),
disp('小波熵的值S_wt');
for i=1:2n
m(i)=p(i)log(p(i));
end
S_wt=sum(m)(-1)
disp('小波包分解图形');
subplot(1,3,1);
text(10,10,'原始图象') ;
plot(ECG);
subplot(1,3,2);
text(10,200,'概率分布图象');
plot(p);
subplot(1,3,3);
plot(E);
至于小波分析图像的二维分解:
我给你说个大体思想,你首先对行向量每个行向量元素进行小波变换,然后再去变换后的新矩阵的列向量元素进行小波变换就可以了,具体的书可以参考《matlab小波分析》,在迅雷上面有的下
wavedec2函数用于对图像进行二维小波分解,其函数调用格式如下:
[c,l]=wavedec2(X,n,’wname’);
其中,X表示原始图像,n表示分解层数,wname表示小波函数,c表示各层系数,l表示各层系数对应的长度
ddencmp用于得到全局阀值,其调用格式如下:
[thr,sorh,keepapp]=ddencmp(‘cmp’,’wp’,X);
[thr,sorh,keepapp]=ddencmp(‘cmp’,’wv’,X);
其中cmp表示压缩,wp表示小波包,wv表示小波,X表示原始信号,thr表示阀值,sorh表示阀值类型,s表示软阀值,h表示硬阀值,keepapp=1表示保持近似系数不变
wdencmp用于对数据或图像进行阀值去噪或压缩,其调用格式如下:
[xcomp,c1,l1,perf0,perfl2]=wdencmp(‘gbl’,c,l,’wname’,n,thr,sorh,keepapp);
glb表示利用全局阀值,perf0表示恢复比,perfl2表示压缩比
示例:利用二维小波对图像进行压缩
编写对应的m文件如下:
clc;
load woman;
subplot(1,2,1)
imshow(X,map);
title('原始图像');
[c,l]=wavedec2(X,3,'sym4');
%%获取全局阀值%%
[thr,sorh,keepapp]=ddencmp('cmp','wp',X);
[xcmp,c1,l1,perf0,perfl2]=wdencmp('gbl',c,l,'sym4',3,thr,sorh,keepapp);
subplot(1,2,2)
imshow(xcmp,map);
title('压缩后');
程序运行结果如下图:
小波变换用与图像去噪,噪声会影响图像处理的输入,采集,处理的各个环节及输出结果等全过程,因此对于图像的噪声处理是一个不可忽略的重要的问题,去噪已经成为图像处理中不可或缺的一部分
示例:对图像进行二维小波去噪
编写对应的m文件如下:
load julia;
%%产生噪声信号%%
init=3718025452;
rand('seed',init);
xnoise=X+8rand(size(X));
colormap(map);
subplot(1,3,1)
imshow(X,map);
title('原始信号')
subplot(1,3,2)
imshow(xnoise,map);
title('含有噪声的信号');
%%获取全局阀值%%
[thr,sorh,keepapp]=ddencmp('den','wp',xnoise);
[xden,c1,l1]=wdencmp('gbl',xnoise,'sym4',3,thr,sorh,keepapp);
subplot(1,3,3)
imshow(xden,map);
title('去除噪声后信号');
程序运行结果如下图:
小波分析用于图像增强,图像增强是对图像进行一定处理,使图像比原图更加清晰,视觉效果更好。
示例:利用小波分析对图像进行增强
编写对应的m文件如下:
clc;
load facets;
subplot(1,2,1)
imshow(X,map);
title('原始信号');
[c,l]=wavedec2(X,3,'sym4');
sizec=size(c);
for i=1:sizec(2)
if(c(i)>250)
c(i)=2c(i);
else
c(i)=05c(i);
end
end
y=waverec2(c,l,'sym4');
subplot(1,2,2)
imshow(y,map);
title('增强图像');
程序运行结果如下图:
图像钝化
图像的钝化可以在时域中,也可以在频域中,在时域中处理较为简单,只需要加一个平滑滤波器,使图像中每个点与其邻点做平滑处理即可,在此主要说明图像钝化在频域中的处理。图像钝化是为了突出低频信息,弱化高频信息。
示例:对图像进行频域钝化处理,
编写对应的m文件如下:
load chess;
subplot(1,2,1)
imshow(X,map);
title('原始图像');
[c,l]=wavedec2(X,3,'db4');
sizec=size(c);
for i=1:sizec(2)
if(c(i)>280)
c(i)=c(i)2;
else
c(i)=c(i)05;
end
end
y=waverec2(c,l,'db4');
subplot(1,2,2)
imshow(y,map);
title('采用小波方法钝化图像');
程序运行结果如下图:
图像锐化,与图像钝化刚好相反,是为了突出高频信息,弱化低频信息,从快速变化的成分中分离出系统边界成分,以便进一步识别或者分割等 *** 作。
示例:对图像进行锐化处理
编写对应的m文件如下:
load chess;
subplot(1,2,1)
imshow(X,map);
title('原始图像');
[c,l]=wavedec2(X,3,'db5');
sizec=size(c);
%%突出高频信息,弱化低频信息%%
for i=1:sizec(2)
if(abs(c(i))<280)
c(i)=c(i)2;
else
c(i)=c(i)05;
end
end
y=waverec2(c,l,'db5');
subplot(1,2,2)
imshow(y,map);
title('采用小波方法锐化图像');
程序运行结果如下图:
小波分析用于图像融合
图像融合是将同一图像的两个部分或者不同图像合成一张图,以便合成之后的图形比原来更容易理解。
示例:利用二维小波变换将两幅图像融合在一起
编写对应的m文件如下:
clear all;
load bust;
X1=X;
map1=map;
load woman;
X2=X;
map2=map;
subplot(1,3,1)
imshow(X1,map1);
title('第一幅图像');
subplot(1,3,2)
imshow(X2,map2);
title('第二幅图像');
%%对第二幅图形低频部分和高频部分进行处理%%
for i=1:256
for j=1:256
if (X2(i,j)>120)
X2(i,j)=X2(i,j)2;
else
X2(i,j)=X2(i,j)05;
end
end
end
[c1,l1]=wavedec2(X1,3,'sym4');
[c2,l2]=wavedec2(X2,3,'sym4');
%%对图像进行融合%%
c=c1+c2;
%%减少图像的亮度%%
c=c05;
y=waverec2(c,l1,'sym4');
subplot(1,3,3)
imshow(y,map2);
title('融合后图像');
程序运行结果如下图:
小波分析用于图像分解
对图像分解的目地在于可以更好的观察图像的细节,对图像做出更好的判断,swt2函数用于对图像进行分解,其调用格式如下:
[sa,sh,sv,sd]=swt2(X,N,’wname’);
其中sa,sh,sv,sd分别表示近似系数,水平系数,竖直系数,对角系数,x分解图像,N分解的层数,wname表示小波基名称
示例:对图像进行分解
编写对应的m文件如下:
clear all;
load woman;
[sa,sh,sv,sd]=swt2(X,3,'db3');
s=1;
for i=1:3
subplot(3,4,s)
image(wcodemat(sa(:,:,i),192));
title(['第',num2str(i),'层近似系数']);
subplot(3,4,s+1)
image(wcodemat(sh(:,:,i),192));
title(['第',num2str(i),'层水平系数']);
subplot(3,4,s+2)
image(wcodemat(sv(:,:,i),192));
title(['第',num2str(i),'层竖直系数']);
subplot(3,4,s+3)
image(wcodemat(sd(:,:,i),192));
title(['第',num2str(i),'层对角系数']);
s=s+4;
end
程序运行结果如下图:
机器视觉是通过计算机来模拟人类视觉功能,以让机器获得相关视觉信息和加以理解。可分为“视”和“觉”两部分原理。
“视”是将外界信息通过成像来显示成数字信号反馈给计算机,需要依靠一整套的硬件解决方案,包括光源、相机、图像采集卡、视觉传感器等。“觉”则是计算机对数字信号进行处理和分析,主要是软件算法。
机器视觉在工业上应用领域广阔,核心功能包括:测量、检测、识别、定位等。
产业链可以分为上游部件级市场、中游系统集成/整机装备市场和下游应用市场。
机器视觉上游有光源、镜头、工业相机、图像采集卡、图像处理软件等软硬件提供商,中游有集成和整机设备提供商,行业下游应用较广,主要下游市场包括电子制造行业、 汽车 、印刷包装、烟草、农业、医药、纺织和交通等领域。
机器视觉全球市场主要分布在北美、欧洲、日本、中国等地区,根据统计数据,2014年,全球机器视觉系统及部件市场规模是 367 亿美元,2015年全球机器视觉系统及部件市场规模是42亿美元,2016年全球机器视觉系统及部件市场规模是62亿美元,2002-2016年市场年均复合增长率为12%左右。而机器视觉系统集成,根据北美市场数据估算,大约是视觉系统及部件市场的6倍。
中国机器视觉起步于80年代的技术引进,随着98年半导体工厂的整线引进,也带入机器视觉系统,06年以前国内机器视觉产品主要集中在外资制造企业,规模都较小,06年开始,工业机器视觉应用的客户群开始扩大到印刷、食品等检测领域,2011年市场开始高速增长,随着人工成本的增加和制造业的升级需求,加上计算机视觉技术的快速发展,越来越多机器视觉方案渗透到各领域,到2016年我国机器视觉市场规模已达近70亿元。
机器视觉中,缺陷检测功能,是机器视觉应用得最多的功能之一,主要检测产品表面的各种信息。在现代工业自动化生产中,连续大批量生产中每个制程都有一定的次品率,单独看虽然比率很小,但相乘后却成为企业难以提高良率的瓶颈,并且在经过完整制程后再剔除次品成本会高很多(例如,如果锡膏印刷工序存在定位偏差,且该问题直到芯片贴装后的在线测试才被发现,那么返修的成本将会是原成本的100倍以上),因此及时检测及次品剔除对质量控制和成本控制是非常重要的,也是制造业进一步升级的重要基石。
1在检测行业,与人类视觉相比,机器视觉优势明显
1)精确度高:人类视觉是64灰度级,且对微小目标分辨力弱;机器视觉可显著提高灰度级,同时可观测微米级的目标;
2)速度快:人类是无法看清快速运动的目标的,机器快门时间则可达微秒级别;
3)稳定性高:机器视觉解决了人类一个非常严重的问题,不稳定,人工目检是劳动非常枯燥和辛苦的行业,无论你设计怎样的奖惩制度,都会发生比较高的漏检率。但是机器视觉检测设备则没有疲劳问题,没有情绪波动,只要是你在算法中写好的东西,每一次都会认真执行。在质控中大大提升效果可控性。
4)信息的集成与留存:机器视觉获得的信息量是全面且可追溯的,相关信息可以很方便的集成和留存。
2机器视觉技术近年发展迅速
1)图像采集技术发展迅猛
CCD、CMOS等固件越来越成熟,图像敏感器件尺寸不断缩小,像元数量和数据率不断提高,分辨率和帧率的提升速度可以说日新月异,产品系列也越来越丰富,在增益、快门和信噪比等参数上不断优化,通过核心测试指标(MTF、畸变、信噪比、光源亮度、均匀性、色温、系统成像能力综合评估等)来对光源、镜头和相机进行综合选择,使得很多以前成像上的难点问题得以不断突破。
2)图像处理和模式识别发展迅速
图像处理上,随着图像高精度的边缘信息的提取,很多原本混合在背景噪声中难以直接检测的低对比度瑕疵开始得到分辨。
模式识别上,本身可以看作一个标记过程,在一定量度或观测的基础上,把待识模式划分到各自的模式中去。图像识别中运用得较多的主要是决策理论和结构方法。决策理论方法的基础是决策函数,利用它对模式向量进行分类识别,是以定时描述(如统计纹理)为基础的;结构方法的核心是将物体分解成了模式或模式基元,而不同的物体结构有不同的基元串(或称字符串),通过对未知物体利用给定的模式基元求出编码边界,得到字符串,再根据字符串判断它的属类。在特征生成上,很多新算法不断出现,包括基于小波、小波包、分形的特征,以及独二分量分析;还有关子支持向量机,变形模板匹配,线性以及非线性分类器的设计等都在不断延展。
3)深度学习带来的突破
传统的机器学习在特征提取上主要依靠人来分析和建立逻辑,而深度学习则通过多层感知机模拟大脑工作,构建深度神经网络(如卷积神经网络等)来学习简单特征、建立复杂特征、学习映射并输出,训练过程中所有层级都会被不断优化。在具体的应用上,例如自动ROI区域分割;标点定位(通过防真视觉可灵活检测未知瑕疵);从重噪声图像重检测无法描述或量化的瑕疵如橘皮瑕疵;分辨玻璃盖板检测中的真假瑕疵等。随着越来越多的基于深度学习的机器视觉软件推向市场(包括瑞士的vidi,韩国的SUALAB,香港的应科院等),深度学习给机器视觉的赋能会越来越明显。
4)3d视觉的发展
3D视觉还处于起步阶段,许多应用程序都在使用3D表面重构,包括导航、工业检测、逆向工程、测绘、物体识别、测量与分级等,但精度问题限制了3D视觉在很多场景的应用,目前工程上最先铺开的应用是物流里的标准件体积测量,相信未来这块潜力巨大。
3要全面替代人工目检,机器视觉还有诸多难点有待攻破
1)光源与成像:机器视觉中优质的成像是第一步,由于不同材料物体表面反光、折射等问题都会影响被测物体特征的提取,因此光源与成像可以说是机器视觉检测要攻克的第一个难关。比如现在玻璃、反光表面的划痕检测等,很多时候问题都卡在不同缺陷的集成成像上。
2)重噪音中低对比度图像中的特征提取:在重噪音环境下,真假瑕疵的鉴别很多时候较难,这也是很多场景始终存在一定误检率的原因,但这块通过成像和边缘特征提取的快速发展,已经在不断取得各种突破。
3)对非预期缺陷的识别:在应用中,往往是给定一些具体的缺陷模式,使用机器视觉来识别它们到底有没有发生。但经常遇到的情况是,许多明显的缺陷,因为之前没有发生过,或者发生的模式过分多样,而被漏检。如果换做是人,虽然在 *** 作流程文件中没让他去检测这个缺陷,但是他会注意到,从而有较大几率抓住它,而机器视觉在这点上的“智慧”目前还较难突破。
4机器视觉产业链情况
1)上游部件级市场
主要包括光源、镜头、工业相机、图像采集卡、图像处理软件等提供商,近几年智能相机、工业相机、光源和板卡都保持了不低于20%的增速。根据中国机器视觉产业联盟(CMVU)调查统计,现在已进入中国的国际机器视觉品牌已近200多家(如康耐视、达尔萨、堡盟等为代表的核心部件制造商,以基恩士、欧姆龙、松下、邦纳、NI等为代表的则同时涉足机器视觉核心部件和系统集成),中国自有的机器视觉品牌也已有100多家(如海康、华睿、盟拓光电、神州视觉、深圳灿锐、上海方诚、上海波创电气等),机器视觉各类产品代理商超过300家(如深圳鸿富视觉、微视新纪元、三宝兴业、凌云光、阳光视觉等)。很多国内机器视觉的部件市场都是从代理国外品牌开始,很多企业均与国外的同行有较好的合作,且这种合作具有一定的排他性,这给潜在进入者带来了一定的门槛,因此优质产品的代理商也都有不错的市场竞争力和利润表现。同时,以海康、华睿为代表的国产工业视觉核心部件正在快速崛起。
2)中游系统集成和整机装备市场
国内中游的系统集成和整机装备商有100多家,他们可以给各行业自动化公司提供综合的机器视觉方案,如凌云光、微视新纪元、嘉恒、凌华、阳光视觉、鼎信、大恒图像等。由于国内产品与国际依然有不小差距,很多中游系统集成商和整机装备商又是从核心零部件的贸易做起来的,因此很多在视觉产品的选择方面,依然更为青睐国外品牌。国内品牌为推广自己的软硬件产品,往往需要发展自己的方案集成能力,才能更好的面对市场竞争。
3)下游应用市场
机器视觉下游,主要是给终端用户提供非标自动化综合解决方案的公司,行业属性非常强,核心竞争力是对行业和生产的综合理解和多类技术整合。由于行业自动化的更迭有一定周期性,深受行业整体升级速度、出货量和利润状况影响,因此近两年来看,拉动机器视觉应用普及最主要的还是在电子制造业,其次是 汽车 和制药。
i 半导体和电子生产行业:从国内机器视觉工业上的应用分布来看,46%都集中在电子及半导体制造行业,包括晶圆加工制造的分类切割、PCB检测(底片、内/外层板、成品外观终检等)、SMT贴装检测、LCD全流程的AOI缺陷检测、各种3c组件的表面缺陷检测、3c产品外观检测等
ii 汽车 :车身装配检测、零件的几何尺寸和误差测量、表面和内部缺陷检测、间隙检测等
iii 印刷、包装检测:烟草外壳印刷、食品的包装和印刷、药品的铝塑板包装和印刷等
iv 农业:对农产品的分级、检验和分类
v 纺织:对异纤、云织、经疵、纬疵等瑕疵检测、织物表面绒毛鉴定、纱线结构分析等等。
5机器视觉系统未来发展趋势
1)嵌入式解决方案发展迅猛,智能相机性能与成本优势突出,嵌入式PC会越来越强大
2)模块化的通用型软件平台和人工智能软件平台将降低开发人员技术要求和缩短开发周期
3)3d视觉将走向更多应用场景
以上就是关于请问一下高手们,小波分解到底是怎么一回事,分解后得到的是图形还是表达式呢全部的内容,包括:请问一下高手们,小波分解到底是怎么一回事,分解后得到的是图形还是表达式呢、小波分解层数如何确定、小波变换二层分解与小波分析的二维分解的异同比较等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)