Examples

Example histogram files used below:

Isotropic system

Running in ipython:

run histograms_brain.py isotropy_histogram_10x10.txt --skip-rows 1 --histogram --histogram-fraction --card-orien --card-vectors --threshold 0.5 --prefix iso --isotropy -s

Results:

___Threshold - evaluation of cardinal directions___
n_max 50
___Evaluation of isotropy___
chi_square 0.0
distribution is isotropic, \alpha =  0.0
fraction ditribution: mean \pm SD 0.01 \pm 6.93889390391e-18

shape, vector, vector_means
shape = (h_azim, v_elev)
vector = [azim1, azim2, elev1, elev2]
vector_mean = [azim_mean, elev_mean]

elevs: ...
azims: ...
vector: ...
vector_mean: ...

___Cardinal directions___
elevs: ...
azims: ...
vector: ...
vector_mean: ...
Length density Cardinal directions
_images/iso_histogram_mean_fraction_010_010.png
_images/iso_orientation.png

Single angle system

Running in ipython:

run histograms_brain.py one_angle_histogram_10x10.txt --skip-rows 1     --histogram --histogram-fraction --card-orien --card-vectors --threshold 0.5 --prefix one_angle --isotropy -s

Results:

___Threshold - evaluation of cardinal directions___
n_max 1
___Evaluation of isotropy___
chi_square 99.0
distribution is anisotropic and there are preffered directions, \alpha =  1.0
fraction ditribution: mean \pm SD 0.01 \pm 0.0994987437107

shape, vector, vector_means
shape = (h_azim, v_elev)
vector = [azim1, azim2, elev1, elev2]
vector_mean = [azim_mean, elev_mean]

elevs: ...
azims: ...
vector: ...
vector_mean: ...

___Cardinal directions___
elevs: ...
azims: ...
vector: ...
vector_mean: ...
Length density Cardinal directions
_images/one_angle_histogram_mean_fraction_010_010.png
_images/one_angle_orientation.png

Azimuthal oriented system

Running in ipython:

run histograms_brain.py azimuth_histogram_10x10.txt --skip-rows 1 --histogram --histogram-fraction --card-orien --card-vectors --threshold 0.5 --prefix azimuth --isotropy -s

Results:

___Threshold - evaluation of cardinal directions___
n_max 31
___Evaluation of isotropy___
chi_square 0.272727272727
distribution is anisotropic, \alpha =  0.398491865559
fraction ditribution: mean \pm SD 0.01 \pm 0.00522232967867

shape, vector, vector_means
shape = (h_azim, v_elev)
vector = [azim1, azim2, elev1, elev2]
vector_mean = [azim_mean, elev_mean]

elevs: ...
azims: ...
vector: ...
vector_mean: ...

___Cardinal directions___
elevs: ...
azims: ...
vector: ...
vector_mean: ...
Length density Cardinal directions
_images/azimuth_histogram_mean_fraction_010_010.png
_images/azimuth_orientation.png

Elevation oriented system

Running in ipython:

run histograms_brain.py elevation_histogram_10x10.txt --skip-rows 1 --histogram --histogram-fraction --card-orien --card-vectors --threshold 0.5 --prefix elevation --isotropy -s

Results:

___Threshold - evaluation of cardinal directions___
n_max 31
___Evaluation of isotropy___
chi_square 0.272727272727
distribution is anisotropic, \alpha =  0.398491865559
fraction ditribution: mean \pm SD 0.01 \pm 0.00522232967867

shape, vector, vector_means
shape = (h_azim, v_elev)
vector = [azim1, azim2, elev1, elev2]
vector_mean = [azim_mean, elev_mean]

elevs: ...
azims: ...
vector: ...
vector_mean: ...

___Cardinal directions___
elevs: ...
azims: ...
vector: ...
vector_mean: ...
Length density Cardinal directions
_images/elevation_histogram_mean_fraction_010_010.png
_images/elevation_orientation.png

Anisotropic system with required threshold

Runing in ipython:

run histograms_brain.py anisotropy_histogram_10x10.txt --skip-rows 1 --histogram --histogram-fraction --card-orien --card-vectors --threshold 0.5 --prefix anisotropy --isotropy -s

Results:

___Threshold - evaluation of cardinal directions___
n_max 2
___Evaluation of isotropy___
chi_square 24.5555555556
distribution is anisotropic and there are preffered directions, \alpha =  0.99999927802116295
fraction ditribution: mean \pm SD 0.01 \pm 0.0495535624911

shape, vector, vector_means
shape = (h_azim, v_elev)
vector = [azim1, azim2, elev1, elev2]
vector_mean = [azim_mean, elev_mean]

elevs: ...
azims: ...
vector: ...
vector_mean: ...

___Cardinal directions___
elevs: ...
azims: ...
vector: ...
vector_mean: ...
Length density Cardinal directions
_images/anisotropy_histogram_mean_fraction_010_010.png
_images/anisotropy_orientation.png

Anisotropic system with given number of required directions

Runing in ipython:

run histograms_brain.py anisotropy_histogram_10x10.txt --skip-rows 1 --histogram --histogram-fraction --card-orien --card-vectors --n-max 3 --prefix anisotropy_n_max_3 --isotropy -s

Results:

___Evaluation of isotropy___
chi_square 24.5555555556
distribution is anisotropic and there are preffered directions, \alpha =  0.99999927802116295
fraction ditribution: mean \pm SD 0.01 \pm 0.0495535624911

shape, vector, vector_means
shape = (h_azim, v_elev)
vector = [azim1, azim2, elev1, elev2]
vector_mean = [azim_mean, elev_mean]

elevs: ...
azims: ...
vector: ...
vector_mean: ...

___Cardinal directions___
elevs: ...
azims: ...
vector: ...
vector_mean: ...
Length density Cardinal directions
_images/anisotropy_n_max_3_histogram_mean_fraction_010_010.png
_images/anisotropy_n_max_3_orientation.png

Random oriented system

Run in ipython:

run histograms_brain.py histogram_10x10.txt  --histogram --histogram-fraction --card-orien --card-vectors --threshold 0.5 --isotropy -s

Results:

___Threshold - evaluation of cardinal directions___
n_max 45
___Evaluation of isotropy___
chi_square 0.0131227089718
distribution is isotropic, \alpha =  0.0912016692997
fraction ditribution: mean \pm SD 0.01 \pm 0.00114554393071

shape, vector, vector_means
shape = (h_azim, v_elev)
vector = [azim1, azim2, elev1, elev2]
vector_mean = [azim_mean, elev_mean]

elevs: ...
azims: ...
vector: ...
vector_mean: ...

___Cardinal directions___
elevs: ...
azims: ...
vector: ...
vector_mean: ...
Length density Cardinal directions
_images/histogram_mean_fraction_010_010.png
_images/orientation.png