######################## Data Description This folder contains 10 groups of data. In each group, there are 12 randomly generated images with bifurcation number increasing from 1 to 56, with step = 5. All images are saved in MetaImage format. For details of MetaImage header, please see: http://www.itk.org/Wiki/MetaIO/Documentation#Quick_Start The node information of each image is saved with a structure array in a MAT-file under the same sub-folder, with the following 5 fields: type -- if it's a root, terminal or bifurcation node coord -- coordinate of the node in the volume image parent -- index of parent node of the current node; if it's a root node, then the field is empty children -- indices of two child nodes of the current node; if it's root or terminal, then the field is empty parent_radius -- radius of the parent branch; if it's a root node, then the field is empty ######################## License This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License. For details of this license, please see: http://creativecommons.org/licenses/by-nc-sa/3.0/legalcode ######################## Contact and References VascuSynth Homepage: vascusynth.cs.sfu.ca Mengliu Zhao: mengliuz@sfu.ca Ghassan Hamarneh: hamarneh@sfu.ca Medical Image Analysis Lab, School of Computing Science, Simon Fraser University, Canada If you're using this dataset, please cite the following two papers: Preet Jassi and Ghassan Hamarneh. VascuSynth: Vascular Tree Synthesis Software. Insight Journal, January-June:1-12, 2011 Ghassan Hamarneh and Preet Jassi. VascuSynth: Simulating Vascular Trees for Generating Volumetric Image data with Ground Truth Segmentation and Tree Analysis. Computerized Medical Imaging and Graphics, 34(8):605-616, 2010. You can find the corresponding bibtex entries at: http://vascusynth.cs.sfu.ca/Publications.html ###################### Sample MATLAB code to read the data % read volume image into MATLAB: % use following toolbox: % http://www.mathworks.com/matlabcentral/fileexchange/27983-3d-slicer cd('Group2/data3') info = metaImageInfo('testVascuSynth3_101_101_101_uchar.mhd'); V = metaImageRead(info); % read in tree structure information load treeStructure_3.mat % locate root, bifurcation and terminal indices b = struct2cell(node); root_index = find(sum(squeeze(strcmp(b,'root')))); bif_index = find(sum(squeeze(strcmp(b,'bif')))); term_index = find(sum(squeeze(strcmp(b,'terminal')))); % visualize volume image isosurface(V, 20) axis([0 101 0 101 0 101]) alpha(.4) % visualize bifurcation locations coord = [node(bif_index(:)).coord]; coord = reshape(coord,3, length(coord)/3); x = coord(1,:); y = coord(2,:); z = coord(3,:); hold on scatter3(x, y, z, 80, 'r', 'filled')