求助MATLAB多切片合并为完整图像网格的实现代码

根据你提供的数据首尾行（首行1,1,1,1、末行2,108,160,16），推测数据为四列结构：[数值, X坐标, Y坐标, Z切片索引]，其中X范围1-108、Y范围1-160、Z范围1-16，每个Z值对应一个2D切片。以下是合并切片并可视化的具体代码：

1. 数据读取与3D网格重构 #

假设数据存储在文本文件data.txt中，执行以下代码将零散数据整理为3D数组：

% 读取原始数据
data = readmatrix('data.txt');

% 拆分数值与坐标列
values = data(:, 1);
x_coords = data(:, 2);
y_coords = data(:, 3);
z_slices = data(:, 4);

% 获取坐标的唯一值与维度
x_unique = unique(x_coords);
y_unique = unique(y_coords);
z_unique = unique(z_slices);
nx = length(x_unique);
ny = length(y_unique);
nz = length(z_unique);

% 初始化3D网格并填充数据
grid3D = zeros(nx, ny, nz);
for k = 1:nz
    % 筛选当前Z切片的所有数据
    slice_mask = z_slices == z_unique(k);
    slice_vals = values(slice_mask);
    % 重构为2D切片并填充到3D网格
    grid3D(:, :, k) = reshape(slice_vals, nx, ny);
end

如果你的切片是单独存储为.mat文件（如slice_1.mat至slice_16.mat），可直接用循环读取拼接：

% 预设网格维度（匹配你的数据范围）
grid3D = zeros(108, 160, 16);
for k = 1:16
    load(['slice_', num2str(k), '.mat']); % 假设每个文件中切片变量名为slice_data
    grid3D(:, :, k) = slice_data;
end

2. 完整3D网格可视化 #

方式1：多切片截面展示（匹配目标效果） #

用slice函数绘制正交截面，还原类似目标图的展示效果：

figure;
% 生成坐标网格
[X, Y, Z] = meshgrid(x_unique, y_unique, z_unique);
% 绘制指定Z切片的截面
slice(X, Y, Z, grid3D, [], [], z_unique([1, 8, 16])); % 可自行调整展示的切片位置

xlabel('X'); ylabel('Y'); zlabel('Z');
title('完整3D网格切片可视化');
colormap jet;
colorbar;
axis tight;
grid on;

方式2：等值面空间分布展示 #

如果需要直观呈现数值的空间分布，使用isosurface绘制等值面：

figure;
isosurface(X, Y, Z, grid3D, 1.5); % 1.5为等值面阈值，可根据你的数据调整

axis equal;
xlabel('X'); ylabel('Y'); zlabel('Z');
title('3D网格等值面可视化');
colormap jet;
colorbar;
camlight;
lighting gouraud;

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