/usr/lib/python2.7/dist-packages/dipy/core/subdivide_octahedron.py is in python-dipy 0.10.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 | """Create a unit sphere by subdividing all triangles of an octahedron
recursively.
The unit sphere has a radius of 1, which also means that all points in this
sphere (assumed to have centre at [0, 0, 0]) have an absolute value (modulus)
of 1. Another feature of the unit sphere is that the unit normals of this
sphere are exactly the same as the vertices.
This recursive method will avoid the common problem of the polar singularity,
produced by 2d (lon-lat) parameterization methods.
"""
from .sphere import unit_octahedron, HemiSphere
def create_unit_sphere( recursion_level=2 ):
""" Creates a unit sphere by subdividing a unit octahedron.
Starts with a unit octahedron and subdivides the faces, projecting the
resulting points onto the surface of a unit sphere.
Parameters
------------
recursion_level : int
Level of subdivision, recursion_level=1 will return an octahedron,
anything bigger will return a more subdivided sphere. The sphere will
have $4^recursion_level+2$ vertices.
Returns
---------
Sphere :
The unit sphere.
See Also
----------
create_unit_hemisphere, Sphere
"""
if recursion_level > 7 or recursion_level < 1:
raise ValueError("recursion_level must be between 1 and 7")
return unit_octahedron.subdivide(recursion_level - 1)
def create_unit_hemisphere( recursion_level=2 ):
"""Creates a unit sphere by subdividing a unit octahedron, returns half
the sphere.
Parameters
-------------
recursion_level : int
Level of subdivision, recursion_level=1 will return an octahedron,
anything bigger will return a more subdivided sphere. The sphere will
have $(4^recursion_level+2)/2$ vertices.
Returns
---------
HemiSphere :
Half of a unit sphere.
See Also
----------
create_unit_sphere, Sphere, HemiSphere
"""
sphere = create_unit_sphere( recursion_level )
return HemiSphere.from_sphere(sphere)
|