File Exchange

image thumbnail

Converting a 3D logical array into an STL surface mesh

version 1.6 (74.1 KB) by

Conversion of a voxelised object contained within a 3D logical array into an STL surface mesh.

25 Downloads

Updated

View License

Converting a 3D logical array into an STL surface mesh
======================================================
 
Adam H. Aitkenhead
adam.aitkenhead@christie.nhs.uk
The Christie NHS Foundation Trust
24th May 2010
 
 
USAGE
=====
 
>> CONVERT_voxels_to_stl(STLname,gridDATA,gridX,gridY,gridZ,'ascii')
..or..
>> [faces,vertices] = CONVERT_voxels_to_stl(STLname,gridDATA,gridX,gridY,gridZ,STLformat)

converts the voxelised object contained within the 3D logical array <gridINPUT> into an STL surface mesh, which is saved to the ascii file <STLin>. The x,y,z axes for <gridINPUT> are defined by <gridX>, <gridY> and <gridZ>. The (faces,vertices) data are optional outputs.
 
 
INPUTS
======
 
STLname - string - Filename of the STL file.
gridINPUT - 3D logical array of size (P,Q,R) - The voxelised object (1 => Inside the object, 0 => Outside the object)
gridX - A 1xP array - List of the X axis coordinates.
gridY - A 1xQ array - List of the Y axis coordinates.
gridZ - A 1xR array - List of the Z axis coordinates.
STLformat - string (optional) - STL file format: 'binary' or 'ascii'.
 
 
OUTPUTS
=======
 
faces - Nx3 array - A list of the vertices used in each facet of the mesh, identified using the row number in the array vertices.
vertices - Nx3 array - A list of the x,y,z coordinates of each vertex in the mesh.
 
 
EXAMPLE
=======
 
For an example, run the following script:
>> CONVERT_voxels_to_stl_example
 
 
NOTES
=====
 
- This code does not apply any smoothing. The stl mesh will be exactly the same geometry as the original voxelised object.
 
 

Comments and Ratings (7)

I used this function to convert a stack of masked MRI images into a number of very useful 3D STL models, and I have to say that it works like a charm! Many thanks for sharing this with the community, it will continue to be a useful tool in the future.

Hi Adam,
can you explain me how to compute "GridInput"?
I mean, I did a segmentation and now i have three coordinates vectors X Y Z. how can y use this program to have a Mesh?

Thank you!

Erez

Erez (view profile)

I have a logical array that is 2500 x 2500 x 2400. Will this run or am I wasting my time trying it on something so big?

Rik

Rik (view profile)

Nice program!
I am currently working on a 3d printer and I am looking for a function which can do the opposit, so converting a .STL or any other 3d cad model to a 3d logical matrix.
Does anyone know if there already exists such a function?

Kind regards,

Rik

Nathan

Nathan (view profile)

Hi Adam,
Great program! Very helpful, and quite fast considering the computations involved.

Adam A

Adam A (view profile)

Hi Johannes. If you're able to send an example of an array and the resulting mesh I will see if I can duplicate and fix any problems.

Johannes Korsawe

Really good and fast.
Sometimes there seem to be some artifacts (missing faces) in the resulting mesh. I have not already discovered where, but suspect a missing line somewhere.

Thanks for sharing!

Updates

1.6

Now optionally output the mesh (Faces,Vertices) data

1.5

Allow the STL to be written in either binary or ascii format.

1.4

Improved documentation

1.3

Major speed improvement

1.2

Improved speed.

1.1

Improved memory usage.

MATLAB Release
MATLAB 7.10 (R2010a)
Acknowledgements

Inspired: Export Voxel Data

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video