解析二进制STL文件时三角形丢失问题排查求助

你的问题完全出在二进制STL的读取逻辑上，和软件导出时省略三角形无关，核心错误及修正如下：

1. 致命错误：文件打开模式错误 #

你用文本模式"r"打开二进制STL文件，在Windows系统下，文本模式会自动将文件中的0x0A（换行符）转换为0x0D0A（回车换行），导致文件读取时字节偏移错位。小文件（如12个三角形的立方体）因字节数少，刚好未触发错位，但大文件会出现后续数据全部读错，表现为“丢失三角形”。

修正：改用二进制模式打开文件：

fp = fopen(filename, "rb");

2. 读取整数/浮点数的不规范写法 #

二进制STL的结构定义中：

• 三角形数量是4字节无符号整数
• 法向量、顶点坐标都是4字节单精度浮点数
• 每个三角形末尾的属性字段是2字节无符号短整数

你的代码多次拆分读取单字节，虽然逻辑上能凑出数据，但易出错且效率低，还可能因类型对齐问题导致异常。

修正示例：

// 读取80字节头（一次性读取更高效）
char header[80];
fread(header, sizeof(char), 80, fp);

// 读取三角形数量（4字节无符号整数）
uint32_t tri_count;
fread(&tri_count, sizeof(uint32_t), 1, fp);
*numTriangles = tri_count;

// 读取单个三角形的法向量
fread(&triangles[i].normal.X, sizeof(float), 1, fp);
fread(&triangles[i].normal.Y, sizeof(float), 1, fp);
fread(&triangles[i].normal.Z, sizeof(float), 1, fp);

// 读取2字节属性字段
uint16_t throwaway;
fread(&throwaway, sizeof(uint16_t), 1, fp);

3. 额外注意：字节序问题 #

二进制STL采用小端字节序存储数据，多数现代CPU（x86/ARM）都是小端架构，无需转换；但如果运行在大端架构的系统上，需要对读取的整数、浮点数进行字节序转换。

修正后的完整读取函数示例 #

#include <stdint.h>
#include <stdio.h>
#include <stdlib.h>

// 假设你的结构体定义如下
struct point {
    float X, Y, Z;
};

struct tri {
    struct point normal;
    struct point p1, p2, p3;
};

struct point computeNormal(struct tri t); // 你的计算函数

// RETURNED POINTER MUST BE FREE'D
struct tri* readSTL(char* filename, int* numTriangles) {
    FILE* fp = fopen(filename, "rb");
    if (!fp) {
        *numTriangles = 0;
        return NULL;
    }

    // READ 80 BYTE HEADER
    char header[80];
    fread(header, sizeof(char), 80, fp);

    // READ 4 BYTE NUMBER OF TRIANGLES
    uint32_t tri_count;
    fread(&tri_count, sizeof(uint32_t), 1, fp);
    *numTriangles = tri_count;

    // READ ALL TRIANGLES
    struct tri* triangles = calloc(*numTriangles, sizeof(struct tri));
    if (!triangles) {
        fclose(fp);
        *numTriangles = 0;
        return NULL;
    }

    for (int i = 0; i < *numTriangles; i++) {
        // READ NORMAL VECTOR
        fread(&triangles[i].normal.X, sizeof(float), 1, fp);
        fread(&triangles[i].normal.Y, sizeof(float), 1, fp);
        fread(&triangles[i].normal.Z, sizeof(float), 1, fp);

        // READ VERTEX 1
        fread(&triangles[i].p1.X, sizeof(float), 1, fp);
        fread(&triangles[i].p1.Y, sizeof(float), 1, fp);
        fread(&triangles[i].p1.Z, sizeof(float), 1, fp);

        // READ VERTEX 2
        fread(&triangles[i].p2.X, sizeof(float), 1, fp);
        fread(&triangles[i].p2.Y, sizeof(float), 1, fp);
        fread(&triangles[i].p2.Z, sizeof(float), 1, fp);

        // READ VERTEX 3
        fread(&triangles[i].p3.X, sizeof(float), 1, fp);
        fread(&triangles[i].p3.Y, sizeof(float), 1, fp);
        fread(&triangles[i].p3.Z, sizeof(float), 1, fp);

        // READ ATTRIBUTE FIELD
        uint16_t attr;
        fread(&attr, sizeof(uint16_t), 1, fp);
        
        // 调试输出保留
        printf("FILE NORMAL: %f, %f, %f\n", triangles[i].normal.X,
                triangles[i].normal.Y, triangles[i].normal.Z);

        struct point Normal = computeNormal(triangles[i]);
        printf("Computed NORMAL: %f, %f, %f\n\n", Normal.X, Normal.Y, Normal.Z);
    }

    fclose(fp);
    return triangles;
}

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