| 查看: 965 | 回复: 6 | ||
| 当前主题已经存档。 | ||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 sacooler 的 18 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会 | ||
[交流]
【求助】求matlab边缘检测图像的源代码
|
||
|
小虫是学化学的，对matlab不了解，最近课题碰到难题，需要对图像进行边缘检测，然后求出图像中边缘区域内的面积。可是还有3个月就要毕业了，重头开始学已经来不及，求高手帮助，不胜感激。 |
回复此楼
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
小木虫(金币+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 |
2楼2010-03-11 17:05:56
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
4楼2010-03-16 15:44:53
5楼2010-03-19 08:56:02
6楼2010-03-20 08:46:38
7楼2010-03-30 20:41:11