{"id":4854,"date":"2013-10-04T09:00:22","date_gmt":"2013-10-04T13:00:22","guid":{"rendered":"https:\/\/blogs.mathworks.com\/pick\/?p=4854"},"modified":"2013-10-03T23:00:50","modified_gmt":"2013-10-04T03:00:50","slug":"easy-visualization-of-volumetric-data","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/pick\/2013\/10\/04\/easy-visualization-of-volumetric-data\/","title":{"rendered":"Easy visualization of volumetric data"},"content":{"rendered":"\r\n

Brett<\/a>'s Pick this week is \"vol3d v2\"<\/a>, by Oliver Woodford<\/a>.<\/p>

Today's Pick is brought to you by Oliver Woodford, one of the most popular contributors to the File Exchange, and author of export_fig<\/a>, the<\/i> long-time \"most-downloaded-file\" award winner--itself the recipient of a coveted :) Pick of the Week<\/a> banner.<\/p>

When I'm demonstrating the use of MATLAB for image processing, I often show how one can use the output of the isosurface<\/tt> command as an input<\/i> to the patch<\/tt> command to create easily a rotatable model of a 3D image stack.<\/p>

>> whos data\r\n<\/pre>
Name       Size<\/span>                   Bytes<\/span>    Class<\/span>    Attributes<\/span>\r\ndata       60x256x256<\/span>             3932160<\/span>  uint8<\/span>\r\n<\/pre>
white_vol = isosurface(data,2.5);\r\ngray_vol  = isosurface(data,1.5);\r\nfigure('renderer'<\/span>,'openGL'<\/span>,'windowstyle'<\/span>,'normal'<\/span>,'color'<\/span>,'w'<\/span>);\r\npatch(white_vol,'FaceColor'<\/span>,'b'<\/span>,'EdgeColor'<\/span>,'none'<\/span>);\r\npatch(gray_vol,'FaceColor'<\/span>,'y'<\/span> ,'EdgeColor'<\/span>,'none'<\/span>, 'FaceAlpha'<\/span>,0.3);\r\n<\/pre>
\"\" <\/p>
I recognize that that visualization could be improved--there are all kinds of options for controlling the rendering, the lighting, the smoothing, the coloration, etc. But the ease with which MATLAB can create such visualizations impresses me.<\/p>
Oliver's \"vol_3d v2\" provides yet another approach to visualization that I will likely be incorporating into my \"3D visualization story.\" With a simple interface modified from the (quite old!) code of fellow MathWorker Joe Conti, Oliver allows us to visualize data, and manipulate its display quite easily:<\/p>
vol3d('cdata'<\/span>,data,'texture'<\/span>,'3D'<\/span>);\r\nview(45,15); daspect(1.\/voxel_size2); axis tight<\/span>;axis off<\/span>;\r\ncamlight; camlight(-90,-10); camlight(180,-10);lighting phong<\/span>;\r\nalphamap('rampup'<\/span>);\r\nalphamap(0.05 .* alphamap);\r\n<\/pre>
\"\" <\/p>
Oliver writes really nice code, and has the downloads and reviews to prove it. Thanks, Oliver, for this cool submission! (And thank you, Joe, for being amenable to having your code modified and improved by others in the user community. I love to see the File Exchange utilized in this manner!)<\/p>
As always, I welcome your thoughts and comments<\/a>. Or leave feedback for Oliver here<\/a>.<\/p>