写一段matlab代码，用来提取等倾干涉同心圆干涉条纹的间距，需要包括转为灰度图像、转二值、高通滤波去噪、、骨化、去毛刺这些步骤，最后给出条纹间距长度

以下是Matlab代码，用于提取等倾干涉同心圆干涉条纹的间距：

% 读入图像
img = imread('your_image_path.jpg');
% 转为灰度图像
gray_img = rgb2gray(img);
% 转二值（使用Otsu自适应阈值法）
level = graythresh(gray_img);
bw_img = im2bw(gray_img, level);
% 高通滤波去噪
g = fspecial('gaussian', [3 3], 0.5); % 高斯滤波
bw_img = imfilter(bw_img, g, 'same');
h = fspecial('laplacian'); % 拉普拉斯滤波
bw_img = imfilter(bw_img, h, 'replicate');
% 骨化
skel = bwskel(bw_img);
% 去毛刺（去除较小的干扰区域）
final_img = bwareaopen(skel, 10);
% 利用骨架图像，提取等倾干涉同心圆干涉条纹的间距
[H, theta, rho] = hough(final_img);
P = houghpeaks(H, 5); % 查找峰值点
lines = houghlines(final_img, theta, rho, P);
distances = zeros(length(lines)-1,1);
for i = 1:length(lines)-1
    distances(i) = abs(lines(i+1).rho - lines(i).rho);
end
avg_distance = mean(distances);
disp(['The average distance of the interferogram fringe lines is: ' num2str(avg_distance) 'pixels']);

希望这个代码可以帮到你。