永发信息网

matlab程序输入后没有显示错误信息

答案:2  悬赏:60  手机版
解决时间 2021-02-14 02:56
I=imread('F:\ 20140422-1\42.bmp');
%RGB图像转为灰度图像
I1=0.3*I(:,:,1)+0.59*I(:,:,2)+0.11*I(:,:,3);
%维纳滤波
J5=wiener2(I1,[5 5]);
%归一化处理
A=J5;
M=mean2(A);
V=std2(A)^2;
[m,n]=size(A);
A1=double(A);
A2=ones(size(A1));
for i=1:m
for j=1:n
if A1(i,j)>M
A2(i,j)=80+sqrt(255*((A1(i,j)-M)^2)/V);
else
A2(i,j)=80-sqrt(255*((A1(i,j)-M)^2)/V);
end
end
end
A2=uint8(A2);
%直方图均衡
B=histeq(A2);
%采用反锐化掩膜法进行锐化处理
H=fspecial('gaussian');
I2=imfilter(B,H);
I3=B-I2;
I4=B+9*I3;
%阈值分割
C=I4;
C1=C;
C2=ones(size(C1));[m,n]=size(C2);
for i=41:(m-40)
for j=41:(n-40)
im=C1((i-40):(i+40),(j-40):(j+40));%81*81块
level=graythresh(im);%otsu法求阈值
level2=256*level;
if C1(i,j) C2(i,j)=0;
end
end
end
%将小于25000像素的单元去除
C3=~C2;
result=bwareaopen(C3,25000);
%中值滤波
Z=medfilt2(result,[10 10]);
%再次去噪声斑块
Z1=bwareaopen(Z,500);
Z2=~Z1;Z3=bwareaopen(Z2,500);Z4=~Z3;
%细化
Z5=bwmorph(Z4,'thin',Inf);
%细化后去除毛刺
Z6=bwmorph(Z5,'spur',30);
最佳答案
我改下就有结果了:
%I=imread('F:\ 20140422-1\42.bmp');
I=imread('齿轮.jpg');
figure,imshow(I),title('原图')
%RGB图像转为灰度图像
I1=0.3*I(:,:,1)+0.59*I(:,:,2)+0.11*I(:,:,3);
figure,imshow(I1),title('灰度图')
%维纳滤波
J5=wiener2(I1,[5 5]);
figure,imshow(J5),title('维纳滤波灰度图')
%归一化处理
A=J5;
M=mean2(A);
V=std2(A)^2;
[m,n]=size(A);
A1=double(A);
A2=ones(size(A1));
for i=1:m
for j=1:n
if A1(i,j)>M
A2(i,j)=80+sqrt(255*((A1(i,j)-M)^2)/V);
else
A2(i,j)=80-sqrt(255*((A1(i,j)-M)^2)/V);
end
end
end
A2=uint8(A2);
figure,imshow(A2),title('归一处理灰度图')
%直方图均衡
B=histeq(A2);
figure,imshow(B),title('直方图均衡灰度图')
%采用反锐化掩膜法进行锐化处理
H=fspecial('gaussian');
I2=imfilter(B,H);
I3=B-I2;
I4=B+9*I3;
figure,imshow(I4),title('锐化处理灰度图')
%阈值分割
C=I4;
C1=C;
C2=ones(size(C1));
[m,n]=size(C2);
for i=41:(m-40)
for j=41:(n-40)
im=C1((i-40):(i+40),(j-40):(j+40));%81*81块
level=graythresh(im);%otsu法求阈值
level2=256*level;
if C1(i,j) C2(i,j)=0;
end
end
end
figure,imshow(C2),title('阈值分割二值图')
%将小于25000像素的单元去除
C3=~C2;
result=bwareaopen(C3,25000);
figure,imshow(result),title('去除小于25000像素二值图')
%中值滤波
Z=medfilt2(result,[10 10]);
figure,imshow(Z),title('中值滤波二值图')
%再次去噪声斑块
Z1=bwareaopen(Z,500);
Z2=~Z1;
Z3=bwareaopen(Z2,500);
Z4=~Z3;
figure,imshow(Z4),title('再次去噪声斑块图')
%细化
Z5=bwmorph(Z4,'thin',Inf);
figure,imshow(Z5),title('细化图')
%细化后去除毛刺
Z6=bwmorph(Z5,'spur',30);
figure,imshow(Z6),title('细化后去除毛刺图')
全部回答
你的程序只是计算,没有给图,你可以用imshow(x,[])的方式绘制你的计算结果
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
亲戚借去我家的房产证复印件 明天还要我妈去
中国移动(前甸镇靠山村移动营业厅)地址好找么
伊洁雅净衣馆巫山巫峡店这个地址在什么地方,
htc好吗,说说
某种商品按成本价的25%为利润定价,然后为吸
企业期末结转前“本年利润”的借方余额为7万
基金理财和P2P理财哪个风险更加大?
此字下一口字读什么
南极人只有吊牌没有商标 是真的假的
中国移动(河东装饰城侧面)地址在什么地方,想
湖南农业大学动物疾病防治研究所地址在什么地
战国时的“宾馆”被称作什么?
关于营养性缺铁性贫血的诊断下列哪项是错误的
小米note2曲面屏有功能么 小米Note2曲面屏
平顶山宝莲升汽车销售服务有限公司地址有知道
推荐资讯
【实变函数与泛函分析】实变函数基础知识在看
邂逅花屋地址在什么地方,我要处理点事
中国联通(天和通讯)地址在什么地方,想过去办
厦门永成宾馆地址在哪,我要去那里办事
梦见自己吃肉圆子 m.sosuo.name
疑犯追踪2 中 卡拉·史丹顿
2011年3月11日13时46分.日本发生9.0级地震.地
什么季节吃香蕉
兴艺网吧地址好找么,我有些事要过去
银行零售业务代表怎么样?是不是要有广泛的交
不参加房交会的楼盘是不是都是小盘不出名咯?
蜜途风格酒店地址在什么地方,想过去办事
正方形一边上任一点到这个正方形两条对角线的
阴历怎么看 ?