matlab图像算子变换,拉布拉斯算子锐化图像的空域和频域--matlab实现

本文将用matlab实现拉布拉斯算子提取边缘，锐化图像。分别采用空域和频域两种方法。

空域

步骤：

1.读取图像/(标定)

2.生成空间滤波器(fspecial('laplacian',0))

3.卷积(m = imfilter(f,w,'replicate'))

clear all;

f = imread('D:\moon.tif')

figure(1)

imshow(f,[]);

title('a')

% b

w = fspecial('laplacian',0); %生成空间滤波器

m = imfilter(f,w,'replicate'); %'replicate'通过复制外边界的值来扩展默认相关

figure(2)

imshow(m,[]);

title('b')

f2 = im2double(f); %标定，将负数变为0，有更好的锐化效果

n = imfilter(f2,w,'replicate');

figure(3)

imshow(n,[])

title('c')

m1 = f - m;

figure(4)

imshow(m1,[]);

title('d')

n1 = f2 - n;

figure(5)

imshow(n1,[]);

title('e')

频域

步骤：

1.读取图像

2.补0，使图像长、宽变为原来两倍

3.将图像移到补零后的中心

4.傅里叶变换

5.在频域进行运算

6.傅里叶逆变换取实部

7.把图像移回原来位置并取原来size的区域

%步骤1

clear all;

f = imread('D:\moon.tif')

f=double(f);

[m,n]=size(f);

figure(1)

imshow(f,[]);

title('origin')

%步骤2

p=2*m;

q=2*n;

A=zeros(p,q);

for i=1:m

for j=1:n

A(i,j)=f(i,j);

end

%步骤3

for i=1:p

for j=1:q

A(i,j)=A(i,j)*(-1)^(i+j);

end

%步骤4

I2=fft2(A);

%步骤5

for i=1:p

for j=1:q

D(i,j)=sqrt((i-p/2)^2+(j-q/2)^2);

end

X=zeros(p,q);

%pi=3.14

for i=1:p

for j=1:q

X(i,j)=-4*pi^2*D(i,j)^2;

Y(i,j)=X(i,j)*I2(i,j);

end

%步骤6

gp3=real(ifft2(Y));

%步骤7

for i=1:p

for j=1:q

gp3(i,j)=gp3(i,j).*(-1).^(i+j);

end

figure(3);

imshow(gp3,[]);

%归一化，目的是让后面gs的第一项和第二项同一个数量级

gp3max = max(gp3(:));

g3 = gp3(1:m,1:n);

k=10; %k自己取，用来改变锐化项的比例

gs = f-k*g3/gp3max;

figure(2)

imshow(gs,[]),title('拉普拉斯锐化图')

matlab图像算子变换,拉布拉斯算子锐化图像的空域和频域--matlab实现相关推荐

MATLAB中拉氏变换/拉氏反变换如何代入数值求解具体结果（简析syms函数和sym类型）
1.前言–函数介绍 1.1 laplace函数功能:拉氏变换示例: syms t f = t; laplace(f) 输出结果: 1.2 ilaplace函数功能:拉氏反变换示例: syms ...
图像清晰度识别之Laplacian算子
最近在做图像清晰度检测的时候,用到了Laplacian算子,所以系统整理下这方面的相关知识. Laplacce算子是一种各向同性算子,二阶微分算子,在只关心边缘的位置而不考虑其周围的象素灰度差值时比较 ...
matlab相位相关图像配准,数字图像处理，相位相关图像配准
前言图像配准是图像处理的基本任务之一,用于将不同时间.不同传感器.不同视角及不同拍摄条件下获取的关于同一目标或场景的两幅或多幅图像进行主要是几何意义上的匹配套和的过程.在对图像配准的研究过程中,大量 ...
OpenCV-Python图像形态变换概述及morphologyEx函数介绍
☞ ░ 前往老猿Python博客 https://blog.csdn.net/LaoYuanPython ░ 一.形态变换概念图像形态变换又称为形态学图像处理.图像形态学,它是基于数学形态学(Mat ...
数字图像处理——图像代数变换
图像代数变换以下是两个图像的像素矩阵: f = [f11,f12,f13; f21,f22,f23; f31,f32,f33] g = [g11,g12,g13; g21,g22,g23; g31, ...
matlab拉普拉斯变换锐化,图像的拉普拉斯锐化方法及讨论
图像的拉普拉斯锐化方法及讨论摘要:本文讲述了空域锐化中常用的二阶微分算法--拉普拉斯算子法.全文首先对拉普拉斯运算做了简单的描述,并简明地分析了其原理:通常是将原图像和对他实施拉式算子后的结果组合后 ...
第3章 Python 数字图像处理(DIP) - 灰度变换与空间滤波17 - 锐化高通滤波器 - 梯度图像（罗伯特，Sobel算子）
目录锐化(高通)空间滤波器使用一阶导数锐化图像-梯度锐化(高通)空间滤波器平滑通过称为低通滤波类似于积分运算锐化通常称为高通滤波微分运算高过(负责细节的)高频,衰减或抑制低频使用一阶 ...
利用matlab函数生成运动模糊图像,运动模糊图像，模糊图像，锐化图像
%用函数创建一个制定的滤波器模版,其语法格式为: %H=fspecial(type); %H=fspecial(type,parameters); %运动模糊图像,模糊图像,锐化图像 I = imre ...
GCN频域视角相关——傅里叶变换、拉普拉斯变换、拉普拉斯算子、拉普拉斯矩阵、卷积
试图通俗地捋清标题名词之间的关系 0. 前置知识 0.1 函数的正交 0.2 什么是卷积? 0.3 散度 0.4 欧拉公式 1. 卷积与傅里叶变换 1.1 傅里叶变换 1.2 时域的卷积等于频域的乘积 ...

matlab图像算子变换,拉布拉斯算子锐化图像的空域和频域--matlab实现

matlab图像算子变换,拉布拉斯算子锐化图像的空域和频域--matlab实现相关推荐

最新文章

热门文章