Matlab中2D平面图转3D遇异常：疑似roicolor函数出错求助

你提到的问题确实和roicolor的使用逻辑有关，而且现有代码还缺少绕Z轴旋转的关键步骤，我来一步步帮你梳理和修正：

先分析现有代码的问题 #

1. roicolor的行为：你说的没错，roicolor(i,0,1)会生成一个二值掩码，其中灰度在0-1范围内的区域（也就是你的黑色线条）设为1，其余区域设为0。但你后续执行i = 10 + peak，会把原来的非黑色区域变成10，黑色区域变成11——两者Z轴高度差只有1，视觉上几乎看不出“延长”效果。
2. 缺少绕Z轴旋转的操作：当前代码只是直接用surf(i)绘制了高度图，完全没有实现绕Z轴旋转的逻辑。
3. imagesc叠加的问题：imagesc(i)会把2D图像叠加在3D表面上，这会掩盖3D效果，反而让画面混乱。

修正后的实现方案 #

下面是调整后的代码，我会逐行解释逻辑：

% 1. 读取并预处理图像
im = imread('F:\floorplan.jpg'); 
gray_im = rgb2gray(im); 

% 2. 提取黑色线条的二值掩码（这里用更直观的阈值化，和roicolor效果一致但更清晰）
black_mask = gray_im <= 1; % 灰度0-1的区域设为1，其余为0

% 3. 构建3D高度矩阵：黑色线条Z轴高度设为50（可调整），其余区域设为0
z_height = zeros(size(gray_im));
z_height(black_mask) = 50; % 让黑色线条明显高于其他区域

% 4. 创建网格坐标（用于3D绘图和旋转）
[X, Y] = meshgrid(1:size(gray_im,2), 1:size(gray_im,1));

% 5. 绕Z轴旋转（这里旋转45度，可根据需求调整角度）
theta = 45; % 旋转角度（单位：度）
rot_matrix = [cosd(theta) -sind(theta) 0; sind(theta) cosd(theta) 0; 0 0 1];
% 把X,Y,Z坐标合并成矩阵，执行旋转
coords = [X(:) Y(:) z_height(:)] * rot_matrix;
X_rot = reshape(coords(:,1), size(X));
Y_rot = reshape(coords(:,2), size(Y));
Z_rot = reshape(coords(:,3), size(z_height));

% 6. 绘制3D效果
figure;
surf(X_rot, Y_rot, Z_rot, 'EdgeColor', 'black', 'FaceColor', 'white');
axis equal; % 保持比例，避免拉伸
view(3); % 切换到3D视角
grid on;

关键调整点说明 #

• 高度区分：把黑色线条的Z值设为远高于其他区域（比如50），这样“延长”的视觉效果会很明显。
• 绕Z轴旋转：通过旋转矩阵手动处理X、Y坐标，实现绕Z轴的任意角度旋转，比直接用rotate函数更灵活可控。
• 绘图优化：去掉了imagesc的叠加，改用surf的颜色参数设置线条和表面颜色，让3D效果更清晰。

如果需要调整旋转角度或者线条高度，直接修改theta和z_height(black_mask)的数值即可。

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