This function polygon2function will convert a Triangulated Mesh into a Voxel Volume which will contain the discretized mesh.
Discretization of a polygon is done by splitting the polygon in smaller polygons, until the longest edge is smaller than 0.5 voxel, then the voxel beneath the vertice coordinates is set to one.
The function is fast because the main function is implemented in c-code /mex file (takes seconds on a core 2 duo with function compiled by Microsoft Visual Studio 2008.).
Try the example.
Bugs, successes and other comments are welcome!
As Stefan Roth pointed out, there is a problem on UNIX platforms. I had the same problem when trying to compile under ubuntu 11.10. Adding the macros indeed helps.
Great job thanks!
On Mac OS compyling the mex file leads to undefined symbol errors.
Adding the macros helps:
#define min(X,Y) ((X) < (Y) ? (X) : (Y))
#define max(X,Y) ((X) > (Y) ? (X) : (Y))
Only integer volume sizes are supported
please check some kind of limitation.
I test a FV data,the data is from real CT images,so the range of object's x,y,z coordinate is about: 140:230 ,90:170, 725:820 , but in your code, I can not set the real range,only use [100 100 100], after test 'none','auto','center'.
only in 'auto', I can see the medical structure in the top left corner, but the structure is scaled down and original details can not be restored.
Thanks you for your comment. The code is basic, thus not based on any publications.
Great Job, Exactly what I needed. Is you algorithm based on any publication(s)?
Great tool. Just solve my need. Well documented.
I have updated the code today with an new version which outputs a logical volume.
THis is indeed fantastic, however I would suggest to change the output Volume from datatype double to uint8 in c file and eventually to logical in m file; this will save lots of memory and make it faster, especially when volumeSize is large;
Incredible performance! Very good
very helpful. Good job.
Function now also working without compiling c-code
Added boundary options
Now first index volume is 1,1,1 instead of 0,0,0, and YXZ dimensions as in Matlab convention (Thanks Yuanming SUO).
Changed output to Logicals instead of Double datatype
Small Bug fixed, for the case if all vertices are outside the volume. Only not draw the face if they are outside at the same side of the volume.
Linux Ubuntu Tested
Added center and resize options
Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.