💥💥💞💞欢迎来到本博客❤️❤️💥💥
🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。
⛳️座右铭:行百里者,半于九十。
📋📋📋本文目录如下:🎁🎁🎁
目录
💥1 概述
在频域中滤除噪声和插值显微镜图像可以通过以下步骤实现:
-
将显微镜图像转换为频域:使用快速傅里叶变换(FFT)将显微镜图像从空域转换为频域。可以使用Matlab中的fft2函数来实现。
-
滤除噪声:在频域中,噪声通常表现为高频成分。可以通过在频域中应用低通滤波器来滤除噪声。常用的滤波器包括高斯滤波器、中值滤波器等。选择适当的滤波器并在频域中应用它,可以减少图像中的噪声。
-
插值:如果显微镜图像中存在缺失的像素或者需要增加图像的分辨率,可以在频域中进行插值操作。可以使用插值算法(如双线性插值、双三次插值等)来填充缺失的像素或增加图像的分辨率。
-
将图像转换回空域:使用逆傅里叶变换(IFFT)将经过滤波和插值处理的图像从频域转换回空域。可以使用Matlab中的ifft2函数来实现。
显微镜(任何具有有限孔径的成像系统)记录的强度受到衍射限制 - 记录的特征不能比系统的光学截止更精细。滤除超出光学截止的空间频率提供了简单而有效的降低噪声的方法。由于显微镜数据是有频带限制的,因此频域填充可提供精确的插值。
📚2 运行结果
部分代码:
% Profiles from airy pattern.
AiryProfiles=cellfun(@(img) xprof(img),...
{AiryPattern,AiryFiltI1,AiryFiltI2,AiryFiltI3},...
'UniformOutput',false);
AiryAxis=cellfun(@(img) xaxis(img,xrange),...
{AiryPattern,AiryFiltI1,AiryFiltI2,AiryFiltI3},...
'UniformOutput',false);
% Profiles from noisy pattern.
NoisyProfiles=cellfun(@(img) xprof(img),...
{NoisyPattern,NoisyFiltI1,NoisyFiltI2,NoisyFiltI3},...
'UniformOutput',false);
NoisyAxis=cellfun(@(img) xaxis(img,xrange),...
{NoisyPattern,NoisyFiltI1,NoisyFiltI2,NoisyFiltI3},...
'UniformOutput',false);
figure(2); clf; set(2,'defaultaxesfontsize',16);
plot(AiryAxis{1},AiryProfiles{1},AiryAxis{2},AiryProfiles{2},...
AiryAxis{3},AiryProfiles{3},AiryAxis{4},AiryProfiles{4},...
'LineWidth',2);
legend('Airy','AiryFiltInterpolation1','AiryFiltInterpolation2','AiryFiltInterpolation3');
snapnow;
figure(3); clf; set(3,'defaultaxesfontsize',16);
plot(NoisyAxis{1},NoisyProfiles{1},NoisyAxis{2},NoisyProfiles{2},...
NoisyAxis{3},NoisyProfiles{3},NoisyAxis{4},NoisyProfiles{4},...
'LineWidth',2);
legend('Noisy','NoisyFiltInterpolation1','NoisyFiltInterpolation2','NoisyFiltInterpolation3');
snapnow;
🎉3 参考文献
部分理论来源于网络,如有侵权请联系删除。
-
Gonzalez, R.C., Woods, R.E., & Eddins, S.L. (2009). Digital Image Processing Using MATLAB. Gatesmark Publishing.
-
Pratt, W.K. (2007). Digital Image Processing: PIKS Inside. John Wiley & Sons.
-
Jain, A.K. (1989). Fundamentals of Digital Image Processing. Prentice-Hall.