代码之家 › 专栏 › 技术社区 › daiana05

在倍频程中实现Sobel梯度

sobel octave image-processing matlab

daiana05 · 技术社区 · 8 年前

因此,我应该在倍频程的图像上应用Sobel滤波器,但我不允许使用图像包中的任何函数。我写了代码,但我的输出只是一个黑色图像。以下是我目前掌握的情况:

%Sobel Gradient

kx= [1 ,0 ,-1; 2,0,-2; 1, 0 ,-1];
ky= [1,2,1; 0,0, 0; -1, -2 ,-1];

H = conv2(kx,im2double(my_img),'same');
V = conv2(ky,im2double(my_img),'same');
E = sqrtm(H.*H + V.*V);
figure 4
imshow(E, [])

谢谢你的帮助!

1 回复 | 直到 8 年前

Ander Biguri 8 年前

很少有小错误,您的问题已经解决:

1) 您通过内核来卷积图像,而不是相反。主要是因为您正在使用 'same' 这将输出第一个输入的大小3x3。

2) 你不想 matrix square root ,只是正常的元素平方根。

my_img=imread('cameraman.tif') % available in MATLAB, not sure octave
kx= [1 ,0 ,-1; 2,0,-2; 1, 0 ,-1];
ky= [1,2,1; 0,0, 0; -1, -2 ,-1];

H = conv2(im2double(my_img),kx,'same');
V = conv2(im2double(my_img),ky,'same');
E = sqrt(H.*H + V.*V); % or sqrt(H.^2+V.^2)
figure
imshow(E, [])

推荐文章

SH_IQ · 在MATLAB条形图中为特定条形图颜色添加其他图例

1 年前

Tom Mozdzen · 如何使用来自向量B的条件和来自向量a值的条件对Matlab向量a进行操作

1 年前

Vlad Vadean · Matlab数组乘法

1 年前

Catalin Baba · 如何在MATLAB中连接两个向量

1 年前

Costinean Sebastian · 如何优化使用for循环处理大型数组的MATLAB脚本的执行时间?"

1 年前

Cole Pratt · 使用MATLAB Actxserver时,避免在Excel公式中添加“隐式交集”符号

1 年前

Leander Heine · 用python从matlab文件中读取字符串数据集

2 年前

servoz · 在matlab和python中创建类似的矩阵对象

2 年前

CircAnalyzer · MATLAB中字符串的十进制数列表

2 年前

I Like Algebra · 在Julia中交换行的最简单方法

2 年前