3D骨架化体素分支长度估计方法咨询及代码报错排查

你的代码存在关键问题，直接导致了NaN和Inf值的出现：

1. bwdistgeodesic参数使用错误：该函数第一个输入应为距离计算的掩码（即骨架本身），而非分支点；且必须限定在骨架体素范围内计算测地线距离。
2. 未区分独立分支：直接计算分支点到端点的距离会混淆不同分支，无法得到单个分支的准确长度。

修正后的实现步骤 #

• 第一步：确认骨架与关键点提取正确
  先保证bwskel和bwmorph3输出正确的3D二值数组：

  skeleton = bwskel(logical(V)); % V为原始3D体素数据
  branch_points = bwmorph3(skeleton, 'branchpoints');
  end_points = bwmorph3(skeleton, 'endpoints');
  

• 第二步：标记所有关键点
  将分支点和端点合并为带唯一ID的标记数组，方便后续区分路径：

  [r, c, z] = ind2sub(size(skeleton), find(branch_points | end_points));
  markers = zeros(size(skeleton));
  for i = 1:length(r)
      markers(r(i), c(i), z(i)) = i;
  end
  

• 第三步：计算骨架上的测地线距离
  指定骨架为掩码、标记点为起点，仅在骨架体素上计算距离：

  dist_matrix = bwdistgeodesic(skeleton, markers);
  

• 第四步：提取单个分支长度
  遍历每个分支点，找到其连通的端点，通过距离差得到分支长度：

  branch_indices = find(branch_points);
  branch_lengths = [];
  for idx = branch_indices
      bp_dist = dist_matrix(idx);
      % 筛选与当前分支点连通的端点
      conn_endpoints = find(end_points & bwdistgeodesic(skeleton, idx) ~= Inf);
      for ep_idx = conn_endpoints
          branch_len = dist_matrix(ep_idx) - bp_dist;
          branch_lengths = [branch_lengths, branch_len];
      end
  end
  

原代码出现NaN/Inf的原因 #

• Inf：分支点与端点不在同一连通分支，或未限定在骨架上计算，导致距离无法获取。
• NaN：若branch_points/end_points为空，或bwdistgeodesic输入参数维度不匹配，会返回NaN。

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