Collects typical statistics on trees intree. Input trees can be organized as:

.single tree

.one group of trees: {tree1, tree2,... treen}

.many groups of trees: {{treei1,...,treein},{treej1,...,treejm},...}

With s a cell array of names defining the individual groups of trees can be passed on, option '-f' saves statistics to file with name name. option '-x' avoids time consuming statistics (use that!!) and option '-s' shows the result (see “dstats_tree”). This function will shortly be greatly enhanced.

Example:

>> dLPTCs = load_tree ('dLPTCs.mtr')

{1x15 cell} {1x10 cell} {1x10 cell} {1x10 cell} {1x10 cell}

>> stats_tree (dLPTCs, {'HSE', 'HSN', 'VS2', 'VS3', 'VS4'}, [], '-w -x')

  gstats: [1x5 struct]
  dstats: [1x5 struct]
  s: {'HSE' 'HSN' 'VS2' 'VS3' 'VS4'}
  dim: 3

>> stats.gstats (global stats)

1x5 struct array with fields:
  len,   total cable length
  max_plen,   maximum path length
  bpoints,   number of branch points
  mpeucl,   mean path/Euclidean distance
  maxbo,   maximum branch order
  mangleB,   mean branch angle
  mblen,   mean branch length
  mplen, vmean path length
  mbo,   mean branch order
  wh,   field height/width
  wz,   field depth/width
  chullx,   center of mass x
  chully,   center of mass y
  chullz,   center of mass z

>> stats.dstats (distributions)

1x5 struct array with fields:
  BO,   branch order distribution
  Plen,   path length distribution
  peucl,   path/Euclidean distance dist.
  angleB,   branch angle distribution
  blen,   branch lengths distribution

Creative Commons License
This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License