STL模型切片时出现异常伪影的技术求助

问题本质 #

你的伪影根源是扫描线填充算法在顶点交点、非轴对齐线段场景下的奇偶计数逻辑错误：

• 顶点处多线段交汇时，交点计数变为奇数，破坏"奇偶翻转填充"的规则；
• 非轴对齐线段的交点映射像素时，直接用floor导致同一像素内多个交点重叠，或线段覆盖的像素未被完整标记，进一步打乱填充状态。

现有代码的核心缺陷 #

1. 轮廓生成逻辑：直接将交点映射到单个像素并填充，未考虑同一行内多个交点的排序与成对关系，顶点处的重复交点会导致填充状态翻转异常；
2. 填充逻辑：逐像素判断非0值来翻转填充状态，无法区分"轮廓线像素"和"内部填充像素"，且未处理顶点处的特殊计数规则。

具体解决方案 #

1. 重构轮廓生成：收集并排序交点，成对处理 #

放弃直接填充单个交点像素的逻辑，改为先收集当前行的所有交点，排序后成对标记轮廓区间，从根源避免奇偶计数错误：

// 替换原轮廓生成代码的填充部分
std::vector<double> row_intersections;

const auto result = intersection(line_seg, segment->segment());
double length = line_seg.bbox().xmax() - line_seg.bbox().xmin();
if (result) {
  Point_2 point;
  if (assign(point, result)) {
    // 先收集所有交点的精确x坐标，暂不映射像素
    row_intersections.push_back(point.x());
  }
}

// 所有交点收集完成后，按x坐标排序
std::sort(row_intersections.begin(), row_intersections.end());

// 成对处理交点，标记轮廓区间（确保进入/退出成对）
for (size_t i = 0; i < row_intersections.size(); i += 2) {
  // 处理顶点导致的奇数个交点（后续优化顶点判断后可移除）
  if (i + 1 >= row_intersections.size()) break;

  double x_start = row_intersections[i];
  double x_end = row_intersections[i + 1];

  // 映射到像素列索引
  int col_start = static_cast<int>(floor(((x_start - x_min_) / length) * num_cols_));
  int col_end = static_cast<int>(floor(((x_end - x_min_) / length) * num_cols_));

  // 确保列索引在有效范围内
  col_start = std::max(0, std::min(col_start, num_cols_ - 1));
  col_end = std::max(0, std::min(col_end, num_cols_ - 1));

  // 填充轮廓区间内的像素（替代原单像素填充）
  for (int col = col_start; col <= col_end; ++col) {
    (*data).data[row * num_cols_ + col] = 255;
  }
}

2. 修复顶点交点的奇偶计数问题 #

当扫描线经过三角面顶点时，需判断顶点连接的两条边的y方向，避免错误计数：

• 若两条边分别在扫描线上下两侧（一条边的起点y > 扫描线y，终点y < 扫描线y；另一条相反），则计数1次（翻转填充状态）；
• 若两条边都在扫描线同侧（都在上方或都在下方），则计数2次（翻转两次，相当于不改变状态）。

示例判断逻辑（需结合三角面的顶点坐标）：

// 假设当前顶点为v，连接的两条边为v->v1和v->v2
bool edge1_above = (v1.y() > line_seg.y());
bool edge2_above = (v2.y() > line_seg.y());

// 只有当两条边分属扫描线两侧时，才将该顶点计入交点
if (edge1_above != edge2_above) {
  row_intersections.push_back(v.x());
}

3. 优化非轴对齐线段的像素覆盖 #

非轴对齐线段的伪影源于线段覆盖的像素未被完整标记，可改用Bresenham算法绘制轮廓线，确保线段经过的每个像素都被正确标记：

// Bresenham算法绘制线段（替代原交点填充）
void draw_line(int x0, int y0, int x1, int y1, unsigned char* data, int num_cols) {
  int dx = abs(x1 - x0);
  int dy = abs(y1 - y0);
  int sx = x0 < x1 ? 1 : -1;
  int sy = y0 < y1 ? 1 : -1;
  int err = dx - dy;

  while (true) {
    if (x0 >=0 && x0 < num_cols && y0 >=0) {
      data[y0 * num_cols + x0] = 255;
    }
    if (x0 == x1 && y0 == y1) break;
    int e2 = 2 * err;
    if (e2 > -dy) { err -= dy; x0 += sx; }
    if (e2 < dx) { err += dx; y0 += sy; }
  }
}

4. 简化填充逻辑：基于交点对直接填充 #

放弃逐像素判断的填充逻辑，直接用排序后的交点对填充区间，彻底避免状态翻转错误：

// 替换原填充代码
for (int y_index = 0; y_index < num_rows; ++y_index) {
  // 提前获取当前行的所有交点对（从轮廓生成阶段传递）
  std::vector<std::pair<double, double>>& row_pairs = get_row_intersection_pairs(y_index);
  
  for (auto& pair : row_pairs) {
    int col_start = floor(((pair.first - x_min_) / length) * num_cols_);
    int col_end = floor(((pair.second - x_min_) / length) * num_cols_);
    col_start = std::max(0, std::min(col_start, num_cols_ - 1));
    col_end = std::max(0, std::min(col_end, num_cols_ - 1));

    for (int x_index = col_start; x_index <= col_end; ++x_index) {
      grayscale_mat.at<uchar>(y_index + LENGTHOFFSET, x_index + WIDTHOFFSET) = 255;
    }
  }
}

总结 #

你的临时修改只是绕过了表面问题，核心要解决的是扫描线算法的交点计数规则和线段像素覆盖精度问题：通过收集排序交点、处理顶点特殊情况、用线段绘制替代单像素标记，就能彻底解决顶点和非轴对齐线段导致的伪影。

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