matlab源文件 gear.m

matlab原文件 版本为2019b
MATLAB [1] 是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。
MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域

bw1=imread（‘gear.bmp‘）; %读取图片
subplot（221）;
imshow（bw1）;%展示读取到的图片
bw2=im2bw（bw1）;%将图片二值化
subplot（222）;
imshow（bw2）;%展示二值化后的图片
bw3=edge（bw2‘canny‘）;%对二值化的图片边缘提取
subplot（223）;
imshow（bw3）;%展示边缘提取后的图片
bw4 = bwmorph（bw3‘remove‘）;%对边缘提取后的图片进行细化
subplot（224）;
imshow（bw4）;%展示细化提取后的图片

figure;%另开一个图像窗口验证边缘追踪的结果是否正确
imshow（bw4）;
x1 = 74; y1 =390;%由之前运行结果可以查bw4的数值，该点应该在齿轮边缘上，由该点开始边缘追踪
gearline = bwtraceboundary（bw4[x1 y1]‘E‘）;%使用边缘追踪函数获得齿轮部分的坐标
hold on
plot（gearline（:2）gearline（:1）‘g‘‘LineWidth‘2）;%在图像中展示出来，验证是否提取成功
x2=250; y2=315;
circle = bwtraceboundary（bw4[x2 y2]‘NE‘）;%使用边缘追踪函数获得内孔部分的坐标
hold on
plot（circle（:2）circle（:1）‘r‘‘LineWidth‘2）;%在图像中展示出来，验证是否提取成功

%分别求取齿轮边缘线和内孔边缘线的中心
%齿轮边缘线部分
[n~]=size（gearline）;
centerx1=0; centery1=0;
for i=1:n