MATLAB中2D荧光显微图像血管直径提取的问题咨询

任务背景 #

处理存储为HDF5（.h5）格式的脑部血管荧光显微图像时间序列，目标是从每一帧中提取血管直径，并分析其随时间的变化情况。

当前MATLAB处理流程 #

• 加载.h5文件并提取单帧图像
• 应用高斯滤波和对比度增强，提升血管可见度
• 执行Canny边缘检测和骨架化操作，提取血管中心线
• 利用基于骨架的欧氏距离变换计算血管直径

现存问题 #

• 边缘检测效果不稳定，常捕获血管内部的伪影
• 提取的骨架有时会延伸至血管外部，或出现不连续的情况
• 计算得到的直径数值看似恒定，无法反映血管实际存在的尺寸变化

当前MATLAB代码实现 #

% Load a single frame from HDF5 file
frame = double(h5read('file.h5', '/recording', [1,1,1,1000], [512,512,1,1]));

% Preprocessing: Smoothing & contrast enhancement
frame = imgaussfilt(frame, 2);
frame = imadjust(mat2gray(frame), stretchlim(mat2gray(frame), [0.01, 0.99]));

% Edge detection & binarization
bw = imbinarize(frame, 'adaptive', 'Sensitivity', 0.5);
bw = bwareaopen(bw, 50);

% Skeletonization & distance transform
skeleton = bwskel(bw);
distance_map = bwdist(~bw);
diameter_map = distance_map * 2;

% Extract diameters along the skeleton
[y, x] = find(skeleton);
diameters = arrayfun(@(i) diameter_map(y(i), x(i)), 1:length(x));

% Plot results
figure;
imshow(frame, []); hold on;
scatter(x, y, 15, diameters, 'filled'); % Color vessels by width
colormap(jet); colorbar;
title('Vessel Skeleton & Diameters');

% Diameter profile plot
figure;
plot(diameters, 'b.-');
xlabel('Skeleton Point Index');
ylabel('Vessel Diameter (pixels)');
title('Vessel Diameter Profile');

附：HDF5单帧图像示例、当前边缘与血管直径检测结果图

核心问题 #

如何实现准确的边缘检测以捕获真实血管边界，并正确计算其直径？

内容的提问来源于stack exchange，提问作者vesseling_stuff