读取条纹图像 img_stripes = imread7bmp; 转换为灰度图像 img_stripes = rgb2grayimg_stripes; 计算每个条纹的曲率 stats = regionpropsimg_stripes PixelList Perimeter; for i = 1lengthstats 拟合曲线 x = statsiPixelList1; y = stat

% 读取条纹图像 img_stripes = imread('7.bmp'); % 转换为灰度图像 img_stripes = rgb2gray(img_stripes); % 计算每个条纹的曲率 stats = regionprops(img_stripes, 'PixelList', 'Perimeter');

for i = 1:length(stats) % 拟合曲线 x = stats(i).PixelList(:,1); y = stats(i).PixelList(:,2); if length(x) <= 2 continue; end % 计算曲线上每个点的位置 xfit = linspace(min(x), max(x), 100); % 在最小值和最大值之间生成100个等间距的点 p = polyfit(x, y, 2); yfit = polyval(p, xfit);

% 计算相邻点之间的距离 
dx = diff(xfit); 
dy = diff(yfit); 
d = sqrt(dx.^2 + dy.^2); 

% 计算条纹的弯曲程度 
curvature = sum(d) / stats(i).Perimeter; 

% 显示条纹和弯曲程度 
subplot(length(stats), 2, 2*i-1), imshow(img_stripes); 
hold on; 
plot(stats(i).PixelList(:,1), stats(i).PixelList(:,2), 'r.'); 
plot(xfit, yfit, 'g-'); 
hold off; 
title(['Stripe ', num2str(i)]); 

subplot(length(stats), 2, 2*i), text(0.5, 0.5, ['Curvature = ', num2str(curvature)]); 
axis off; 

end