Código fuente de cálculo de engranajes de Matlab

¿Yo? =?Estoy leído(' prueba 2 . jpg ');imshow(1);¿mantener? Abrir;

¿Yo? =?im2bw(一);

¿Yo? =?~Yo;

I1? =?imfill(I,'agujero');

se? =?strel('disco ', 50);

I2? =?imclose(I1, se);

[B, L]? =?Límite(1);

¿Estadísticas? =?regionprops(L,'centro de masa');

x0? =?ESTADÍSTICAS(1). Centroide(1);y0? =?ESTADÍSTICAS(1). Centroide(2);

Plot(x0, y0, 'r+')

¿Estadísticas? =?regionprops(L,' cuadro delimitador');

r? =?ESTADÍSTICAS(1). Cuadro delimitador(3)/2;

θ? =?0:0.01:2 * pi;

x? =?r * cos(θ)? +?x0;

y? =?r * pecado(θ)? +?y0;

plot(x, y, ' r ');

text(x0, y10, ['r= '?num2 str(r)]) ;