遇到一个三维的二进制文件，通过FILE指针进行IO操作，并存入三维Mat中，通过二维Mat进行分层显示

一、打开文件

fopen_s(&file, "XXX.uint16_scv", "rb"，通过file文件指针打开XXX.uint16_scv文件，rb以二进制可读的形式
fopen_s：打开文件成功返回0，失败返回非0。

FILE* file;
if (fopen_s(&file, "XXX.uint16_scv", "rb") != 0)
{std::cerr << "Error opening file: " << std::endl;return false;
}
else
{std::cerr << "Opening file is successful " << std::endl;
}

二、判断FILE是否为空指针

if (file == nullptr)
{std::cerr << "Error opening file, file is nullptr: " << std::endl;return false;
}
else
{std::cerr << "File is no empty,Okk! " << std::endl;
}

三、文件中前1024个字节是注释信息无效，需要通过前1024给字节

int fseek(FILE *stream, long int offset, int whence)
参数一：指向 FILE 对象的指针
参数二：相对 whence 的偏移量，以字节为单位
参数三：开始添加偏移 offset 的位置，其中SEEK_SET为文件开头、SEEK_CUR为文件当前位置、SEEK_END为文件末尾

fseek(file, 1024, SEEK_SET)
SEEK_SET从文件头开始，偏移1024个字节，相当于跳过前1024个字节

fseek：操作成功返回0，失败返回非0。

    // 跳过前1024个字节，如果成功，则该函数返回零，否则返回非零值。if (fseek(file, 1024, SEEK_SET) != 0){std::cerr << "Error seeking file\n";fclose(file);return false;}else{std::cerr << "Now skip 1024 byte is successful " << std::endl;}

四、通过三维mat接收字节流

已知三维数据的各个维度的大小xs、ys和zs
size_t count{ static_cast<size_t>(xs) * static_cast<size_t>(ys) * static_cast<size_t>(zs) }; count为数据的总大小
mat = Mat(3, sizes, CV_16UC1);定义一个三维的mat对象，用于存储数据流信息
fread(mat.data, sizeof(uint16_t), count, file)从file中读取count字符到mat.data中且每个数据信息为uint16_t大小

    int xs{ 1032 }, ys{ 1760 }, zs{ 1213 };size_t count{ static_cast<size_t>(xs) * static_cast<size_t>(ys) * static_cast<size_t>(zs) };const int sizes[3] = { zs, ys, xs };mat = Mat(3, sizes, CV_16UC1);if (fread(mat.data, sizeof(uint16_t), count, file) != count){std::cerr << "Error reading file(fread)" << std::endl;return false;}else{std::cerr << "read file is successful, count is :  " << count << std::endl;}

此时三维的mat已经存放了二进制文件中存放的数据流信息

五、以Y轴为例，计算每层的数据之和

mat.ptr<uint16_t>(z_, y_)[x_]获取三维mat的xyz信息

    unsigned long arr[kys]{};//用于存储Y轴的所有平面数据之和for (int y_ = 0; y_ < ys; y_++)//y{unsigned long value = 0;for (int x_ = 0; x_ < xs; x_++)//x{for (int z_ = 0; z_ < zs; z_++)//z{value += unsigned long(mat.ptr<uint16_t>(z_, y_)[x_]);}}arr[y_] = value;}

六、将指定的Y轴某一层进行赋值给二维Mat

因为这里给到的数据XYZ是相反的，故mat.ptr<uint16_t>(z_,300)[x_]表示Y轴的300位置的切片
*mat1.ptr<uint16_t>(z_, x_) = mat.ptr<uint16_t>(z_,300)[x_] / 255;将这层切片赋值给二维的mat1，记得维度信息得对应且每层的大小也要清楚，Y轴方向切得到的图片大小为Z × X

Mat mat1 = Mat::zeros(zs,xs, CV_8UC1);
std::cout << "mat1.size() is : " << mat1.size() << std::endl;
for (int x_ = 0; x_ < xs; x_++)//x
{for (int z_ = 0; z_ < zs; z_++)//z{*mat1.ptr<uint16_t>(z_, x_) = mat.ptr<uint16_t>(z_,300)[x_] / 255;}
}

七、显示二维Mat

传入mat1对象即可

namedWindow("mat1_y", WINDOW_FREERATIO);
imshow("mat1_y", mat1);

其他维度信息也都类似，因为项目需要，故进行了记录一下，方便后续回顾