6 图像处理实现螺纹识别案例（matlab程序）-Toy模板网

这篇具有很好参考价值的文章主要介绍了6 图像处理实现螺纹识别案例（matlab程序）。希望对大家有所帮助。如果存在错误或未考虑完全的地方，请大家不吝赐教，您也可以点击"举报违法"按钮提交疑问。

学习目的：学习识别案例掌握识别方法

文章来源地址https://www.toymoban.com/news/detail-526065.html

2.代码

clear;clc;close all
I=imread('luowen1.bmp'); %读取螺纹图片
try
I=rgb2gray(I); %如果是RGB图，则转换成灰度图
catch
end
figure
imshow(I)
title('原图（半边螺纹）')
for K=1:15
I=wiener2(I,[5 5]); %通过滤波15次把多余点滤除
end
[m,n]=size(I);

I=edge(I,'canny'); %通过边缘检测得到螺纹的波形
[m,n]=size(I); %计算图像的尺寸

I=I(20:m-20,20:n-20); %把图像的边角去掉，只留下有用部分
[m,n]=size(I); %计算去掉边角之后的图像尺寸
figure
imshow(I) %显示去掉边角之后的图像
title('螺纹波形')

%% 把白像素点的位置得到（像素为1的点的坐标）
N=1; %计数器
for i=1:m
for j=1:n
if I(i,j)==1
x(N)=i; %保存白色像素的横坐标x
y(N)=j; %保存白色像素的纵坐标y
N=N+1; %计数器+1
end
end
end

[x,IX]=sort(x); %将x按升序排列
y=y(IX); %对应的y顺序依次和x对应
figure
plot(x,y) %绘制出螺纹线
xlabel('横轴x') %标注横轴x
ylabel('纵轴y') %标注纵轴y
title('计算结果')
view([90 90])
hold on
axis equal
axis([1 m 1 n])

M=10; %设置波峰和波谷的个数
lengthxy=length(x); %计算x和y数据点的个数
dlength=floor(lengthxy/M); %将x和y等分成M段，计算每一段的长度

for K=1:M
xx{K}=x((K-1)*dlength+1:K*dlength); %保存每一段x的值
yy{K}=y((K-1)*dlength+1:K*dlength); %保存每一段y的值
end

for K=1:M
[bofengy(K),index1]=max(yy{K}); %寻找每一段y里的最大值
bofengx(K)=xx{K}(index1); %得到与波峰对应的横坐标x
[boguy(K),index2]=min(yy{K}); %寻找每一段y里的最小值
bogux(K)=xx{K}(index2); %得到与波谷对应的横坐标x
end
scatter(bofengx,bofengy) %绘制出波峰的点
scatter(bogux,boguy,'k') %绘制出波谷的点

%% 对所有的波峰进行直线拟合，计算出直线斜率k1和截距b1
A1=[bofengx',ones(length(bofengx),1)];
kb1=A1\bofengy';
k1=kb1(1);
b1=kb1(2);
%% 绘制出波峰拟合得到的直线
x0=[1 m];
Y1=k1*x0+b1;
plot(x0,Y1,'m');

%% 对所有的波谷进行直线拟合，计算出直线的斜率k2和截距b2
A2=[bogux',ones(length(bogux),1)];
kb2=A2\boguy';
k2=kb2(1);
b2=kb2(2);
%% 绘制出波谷拟合直线
x0=[1 m];
Y2=k2*x0+b2;
plot(x0,Y2,'r');

3.运行结果

6 图像处理实现螺纹识别案例（matlab程序）,matlab,图像处理,数学建模,学习