STL to binary mask conversion

Question

Nicolás 2024-11-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2169473-stl-to-binary-mask-conversion

回答： Hitesh 2024-11-29

I have an STL file and I would like to make a binary mask from it in MATLAB, that would be defined similarly to this:

stl_mask = zeros(size(stl_file));
stl_mask(stl_indices) = 1

I have searched around and there is a stlread function that lets me visualize the STL file in MATLAB, but I am a bit lost on how to work around it. This is the STL file that I'm using. For the moment, I have just this small code to visualize the STL file:

clearvars;
filename = 'example.STL';
[v, f, n, c, stltitle] = stlread(filename);
figure;
patch('Faces',f,'Vertices',v,'FaceVertexCData',c);
grid on; view(45, 25); daspect([1 1 1]); camlight(-30,-30);
figure; 
plot3(v(:,1),v(:,2),v(:,3),'.');    

I have also checked the stl_slice_and_plot function from https://es.mathworks.com/matlabcentral/fileexchange/62113-slice_stl_create_path-triangles-slice_height . It gives a nice plot of the contours, but again, I would need some help on how to convert the contour plot into (x y z) coordinates to generate a mask:

%this script shows how everything works together
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Sunil Bhandari 
%3/17/2017
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc; clear vars; close all;
triangles = read_binary_stl_file('example.STL');
original = triangles;
triangles = rotate_stl(triangles,'y',90);
slice_height = .1;
tic;[movelist, z_slices] = slice_stl_create_path(triangles, slice_height);toc;
%'plotting'
plot_slices(movelist,z_slices, 0)

Thank you

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Hitesh 2024-11-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2169473-stl-to-binary-mask-conversion#answer_1551518

在 MATLAB Online 中打开

Hi Nicolás,

You need to add "inpolyhedron" file exchange for checking whether the points are inside a 3D triangulated (faces/vertices) surface. Kindly refer to the below code ofr generating the binary mask.

filename = 'example.STL';
[TR,fileformat,attributes,solidID] = stlread(filename);
addpath('C:\Users\hiteshk\Downloads\inpolyhedron');
% Define the grid resolution and limits
gridSize = 100; % Adjust the resolution as needed
xRange = [min(TR.Points(:,1)), max(TR.Points(:,1))];
yRange = [min(TR.Points(:,2)), max(TR.Points(:,2))];
zRange = [min(TR.Points(:,3)), max(TR.Points(:,3))];
[xGrid, yGrid, zGrid] = ndgrid(...
    linspace(xRange(1), xRange(2), gridSize), ...
    linspace(yRange(1), yRange(2), gridSize), ...
    linspace(zRange(1), zRange(2), gridSize));
% Initialize the binary mask
binaryMask = false(size(xGrid));
% Use inpolyhedron or similar function to fill the binary mask
binaryMask = inpolyhedron(TR.ConnectivityList, TR.Points, [xGrid(:), yGrid(:), zGrid(:)]);
% Reshape the binary mask
binaryMask = reshape(binaryMask, size(xGrid))

You need to "trisurf" function to visualize the STL using the triangulation object. Kindly refer to the below code for visualazing the STL file.

% Create a figure to visualize the STL using the triangulation object
figure;
trisurf(TR, 'FaceColor', 'cyan', 'EdgeColor', 'none');
title('STL Visualization');
xlabel('X');
ylabel('Y');
zlabel('Z');
grid on;
axis equal;
view(45, 25);
camlight('headlight');
lighting gouraud;