基于二维连续小波变换条纹图相位提取方法附matlab代码

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法       神经网络预测       雷达通信       无线传感器        电力系统

信号处理              图像处理               路径规划       元胞自动机        无人机 

⛄ 内容介绍

提出了一种基于二维连续小波变换的电子散斑干涉(ESPI)条纹图相位提取方法.通过检测二维小波脊确定条纹相位,并引入条纹频率作为向导,有效地避免了相位解调过程中的符号奇异性问题,从而使该算法既能处理开条纹图也能处理闭条纹图,且对散斑噪声具有较强的抑制能力.数值模拟和实验结果表明,该方法在抑制散斑噪声的同时能够有效地提取出条纹相位,对开条纹图和闭条纹图都能处理.

⛄ 部分代码

close all;clc;clear;

%% 调制光栅进行  不适用于横条纹

% img=(imread('face 4.jpg'));  %读入图片

% img=(imread('5D.bmp'));  %读入图片

% img=(imread('2B.jpg'));  %读入图片

 img=(imread('4P.jpg'));  %读入图片

% img=(imread('2C.bmp'));  %读入图片

if numel( size( img ) ) > 2

    img = rgb2gray( img ); % 3通道图像转换为单通道   

end

global img_height;

global img_width;

[img_height, img_width ] = size( img ); 

% sigma=25;%添加噪声的标准差

% imgnoise=imnoise(img,'gaussian',0,(sigma/255)^2); %添加均值为0，标准差为sigma高斯噪声

% imshow(imgnoise);

% img = imgnoise;

%% 二维连续小波变换

% ：paul

% ：fan、gabmexh、cauchy

%  dog

% ：sinc endstop1 gaus3 gaus2 gaus escauchy esmexh esmorl

% ：rmorl endstop2 isomorl dog2 isodog dogpow pethat wheel   

cwtstruct = cwtft2(img,'wavelet','fan','scales',1:0.5:20,'angles',0,'plot'); axis equal;

global Len_scales;

global Len_angles; % 设为全局变量，供部分子函数调用

Len_scales = size(cwtstruct.cfs, 4);  % 获得尺度的个数

Len_angles = size(cwtstruct.cfs, 5); % 获得角度变化的个数

%% 分别对每行平移量利用代价函数提取小波脊

[opt_Scale, opt_Angle] = phase_extract_cost( cwtstruct );

⛄ 运行结果

⛄ 参考文献

[1] 韩民, 成旭. 基于小波变换的多尺度锥束CT图像快速三维重建方法:, CN106228601A[P]. 2016.

[2] 王翠珍, 陈振学, 王燕. 基于小波和双边滤波的活体牙CT图像去噪及三维重建[J].  2016.

[3] 韩民, 成旭, 李登旺. 基于小波变换的多分辨率锥束CT图像快速三维重建算法[J]. 电子与信息学报, 2017, 39(10):5.

[4] 孔令凯, 温俊海. 基于二维小波变换的SPECT图像重建算法研究[C]// 2007中国生物医学工程联合学术年会. 0.

⛳️ 代码获取关注我

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料