MRC格式3D图像栈的读取与逐图处理技术问询

一、当前导入方案是否最优？ #

你当前的文件流读取方式并非最优，原因有两点：

1. 代码仅能处理单个MRC文件，无法满足读取文件夹下多个MRC文件组成图像栈的需求；
2. 手动操作文件流需要自行解析MRC头信息、字节序等细节，容易出错，DM Script本身提供了更简便的批量/栈导入API，能自动处理格式解析。

针对你的需求，推荐两种更优的导入方式：

1. 读取文件夹下多个MRC文件组成图像栈 #

// 选择目标文件夹
string folder = GetDirectoryDialog("选择存放MRC文件的文件夹", "")
if (folder == "") exit()

// 获取文件夹内所有MRC文件列表
string fileList[]
FindFiles(folder, "*.mrc", fileList, 0)
number fileCount = len(fileList)
if (fileCount == 0) {
    Result("文件夹中无MRC文件\n")
    exit()
}

// 获取第一张图的尺寸，创建对应大小的3D图像栈
Image firstImg := ImageImport(fileList[0], 0)
number sizeX = firstImg.ImageGetDimensionSize(0)
number sizeY = firstImg.ImageGetDimensionSize(1)
ImageStack stack := IntegerImage("MRC图像栈", 1, 0, sizeX, sizeY, fileCount)

// 批量导入文件到栈中（保持原文件顺序）
for (number i = 0; i < fileCount; i++) {
    Image tempImg := ImageImport(fileList[i], 0)
    stack.ImageCopySlice(tempImg, i)
    tempImg.Close()
}
stack.ShowImage()

如果需要严格按文件名数字排序（如file_001.mrc、file_002.mrc），可以对fileList做排序处理，避免系统默认排序混乱。

2. 读取单个MRC文件内部的3D栈 #

若你的需求是读取单个MRC文件自带的3D图像栈，直接用ImageImport即可：

string filePath = OpenFileDialog("选择单个MRC栈文件", "", "*.mrc")
if (filePath == "") exit()
ImageStack stack := ImageImport(filePath, 0)
stack.ShowImage()

二、图像栈的逐图处理 #

先定义你的单图处理函数（示例为阈值处理，替换为你的实际逻辑），再遍历栈的每个切片完成处理：

方式1：原地修改原栈 #

// 自定义处理函数
void ProcessSingleImage(Image img) {
    // 示例：二值化处理，替换为你的实际操作（如滤波、增强等）
    img = img > 100 ? 255 : 0
}

// 遍历所有切片处理
number sliceCount = stack.ImageGetDimensionSize(2)
for (number i = 0; i < sliceCount; i++) {
    Image slice := stack.ImageGetSlice(i)
    ProcessSingleImage(slice)
}

方式2：创建新栈存储处理结果（保留原栈） #

number sizeX = stack.ImageGetDimensionSize(0)
number sizeY = stack.ImageGetDimensionSize(1)
number sliceCount = stack.ImageGetDimensionSize(2)
ImageStack processedStack := IntegerImage("处理后图像栈", 1, 0, sizeX, sizeY, sliceCount)

for (number i = 0; i < sliceCount; i++) {
    Image slice := stack.ImageGetSlice(i).Clone() // 复制原切片避免修改原数据
    ProcessSingleImage(slice)
    processedStack.ImageCopySlice(slice, i)
    slice.Close()
}
processedStack.ShowImage()

三、按原顺序输出处理后的图像 #

根据需求选择输出方式：

1. 保存为单个3D MRC栈文件 #

string savePath = SaveFileDialog("保存处理后的MRC栈", "", "*.mrc")
if (savePath != "") {
    processedStack.ImageSave(savePath)
}

2. 按原文件名顺序保存为单个MRC文件 #

string saveFolder = GetDirectoryDialog("选择保存单个文件的文件夹", "")
if (saveFolder == "") exit()

for (number i = 0; i < sliceCount; i++) {
    string originalName = GetFileName(fileList[i])
    string saveName = saveFolder + "/" + ReplaceString(originalName, ".mrc", "_processed.mrc")
    Image slice := processedStack.ImageGetSlice(i)
    slice.ImageSave(saveName)
    slice.Close()
}
Result("所有文件处理并保存完成\n")

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