24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 sacooler 的 18 个金币，回帖就立即获得 1 个金币，每人有 1 次机会

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

[交流] 【求助】求matlab边缘检测图像的源代码

小虫是学化学的，对matlab不了解，最近课题碰到难题，需要对图像进行边缘检测，然后求出图像中边缘区域内的面积。可是还有3个月就要毕业了，重头开始学已经来不及，求高手帮助，不胜感激。

回复此楼

» 猜你喜欢

2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
孩子确诊有中度注意力缺陷已经有6人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

1楼2010-03-11 16:58:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wujixian

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 2003
帖子: 423
在线: 250.1小时
虫号: 126308
注册: 2005-12-07
性别: GG
专业: 化学反应工程

sacooler(金币+5): 2010-03-12 14:39

f=imread('0255.jpg');
>> g=edge(f,'Sobel');
>> figure,imshow(g)

%上面是自动阈值的，如果要取特定阈值，可以这样：g=edge(f,'Sobel'，0.075);剩下的几个算法是把Sobel改为prewitt,roberts,rob,log,Marr,Canny。“Marr”算法太老了，matlab会自动用“laplace”代替，汗...
%加入噪声的语句：
f1 = imnoise(f,'salt & pepper', 0.02);
figure,imshow(f1);
%这是加椒盐噪声“salt & pepper”的，浓度为0.02

赞一下

回复此楼

3楼2010-03-11 17:08:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 7 个回答

wujixian

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 2003
帖子: 423
在线: 250.1小时
虫号: 126308
注册: 2005-12-07
性别: GG
专业: 化学反应工程

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
sacooler(金币+5):谢谢您的帮助，我试试看。希望有机会具体向你讨教 2010-03-12 14:35
adu886886(金币+2):谢谢交流，欢迎常来仿真模拟板块 2010-03-16 18:57

clear all;
load wbarb;
I = ind2gray(X,map);imshow(I);
I1 = imadjust(I,stretchlim(I),[0,1]);figure;imshow(I1);
[N,M] = size(I);
h = [0.125,0.375,0.375,0.125];
g = [0.5,-0.5];
delta = [1,0,0];
J = 3;
a(1:N,1:M,1,1:J+1) = 0;
dx(1:N,1:M,1,1:J+1) = 0;
dy(1:N,1:M,1,1:J+1) = 0;
d(1:N,1:M,1,1:J+1) = 0;
a(:,:,1,1) = conv2(h,h,I,'same');
dx(:,:,1,1) = conv2(delta,g,I,'same');
dy(:,:,1,1) = conv2(g,delta,I,'same');
x = dx(:,:,1,1);
y = dy(:,:,1,1);
d(:,:,1,1) = sqrt(x.^2+y.^2);
I1 = imadjust(d(:,:,1,1),stretchlim(d(:,:,1,1)),[0 1]);figure;imshow(I1);
lh = length(h);
lg = length(g);
for j = 1:J+1
lhj = 2^j*(lh-1)+1;
lgj = 2^j*(lg-1)+1;
hj(1:lhj)=0;
gj(1:lgj)=0;
for n = 1:lh
hj(2^j*(n-1)+1)=h(n);
end
for n = 1:lg
gj(2^j*(n-1)+1)=g(n);
end
a(:,:,1,j+1) = conv2(hj,hj,a(:,:,1,j),'same');
dx(:,:,1,j+1) = conv2(delta,gj,a(:,:,1,j),'same');
dy(:,:,1,j+1) = conv2(gj,delta,a(:,:,1,j),'same');
x = dx(:,:,1,j+1);
y = dy(:,:,1,j+1);
dj(:,:,1,j+1) = sqrt(x.^2+y.^2);
I1 = imadjust(dj(:,:,1,j+1),stretchlim(dj(:,:,1,j+1)),[0 1]);figure;imshow(I1);
End

赞一下(3人)

回复此楼

2楼2010-03-11 17:05:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖