0-Ph82ddlZddlZddlmZddlmZ dddddddd d Zd Zd Z ej ddgddgddgddgddgddgddgddgddgddgddgddgg d Z ej ddgddgddgddgddgddgddgddgddgddgddgddgg d Z ej ddgddgddgddgddgddgddgddgddgddgddgddgg d Z dZdZdS)N)_marching_cubes_lewiner_luts)_marching_cubes_lewiner_cy)?rrdescentTlewiner)spacinggradient_direction step_sizeallow_degeneratemethodmaskc pd}|dkrd}n|dkrtdt||||||||S)aMarching cubes algorithm to find surfaces in 3d volumetric data. In contrast with Lorensen et al. approach [2]_, Lewiner et al. algorithm is faster, resolves ambiguities, and guarantees topologically correct results. Therefore, this algorithm generally a better choice. Parameters ---------- volume : (M, N, P) ndarray Input data volume to find isosurfaces. Will internally be converted to float32 if necessary. level : float, optional Contour value to search for isosurfaces in `volume`. If not given or None, the average of the min and max of vol is used. spacing : length-3 tuple of floats, optional Voxel spacing in spatial dimensions corresponding to numpy array indexing dimensions (M, N, P) as in `volume`. gradient_direction : string, optional Controls if the mesh was generated from an isosurface with gradient descent toward objects of interest (the default), or the opposite, considering the *left-hand* rule. The two options are: * descent : Object was greater than exterior * ascent : Exterior was greater than object step_size : int, optional Step size in voxels. Default 1. Larger steps yield faster but coarser results. The result will always be topologically correct though. allow_degenerate : bool, optional Whether to allow degenerate (i.e. zero-area) triangles in the end-result. Default True. If False, degenerate triangles are removed, at the cost of making the algorithm slower. method: {'lewiner', 'lorensen'}, optional Whether the method of Lewiner et al. or Lorensen et al. will be used. mask : (M, N, P) array, optional Boolean array. The marching cube algorithm will be computed only on True elements. This will save computational time when interfaces are located within certain region of the volume M, N, P-e.g. the top half of the cube-and also allow to compute finite surfaces-i.e. open surfaces that do not end at the border of the cube. Returns ------- verts : (V, 3) array Spatial coordinates for V unique mesh vertices. Coordinate order matches input `volume` (M, N, P). If ``allow_degenerate`` is set to True, then the presence of degenerate triangles in the mesh can make this array have duplicate vertices. faces : (F, 3) array Define triangular faces via referencing vertex indices from ``verts``. This algorithm specifically outputs triangles, so each face has exactly three indices. normals : (V, 3) array The normal direction at each vertex, as calculated from the data. values : (V,) array Gives a measure for the maximum value of the data in the local region near each vertex. This can be used by visualization tools to apply a colormap to the mesh. See Also -------- skimage.measure.mesh_surface_area skimage.measure.find_contours Notes ----- The algorithm [1]_ is an improved version of Chernyaev's Marching Cubes 33 algorithm. It is an efficient algorithm that relies on heavy use of lookup tables to handle the many different cases, keeping the algorithm relatively easy. This implementation is written in Cython, ported from Lewiner's C++ implementation. To quantify the area of an isosurface generated by this algorithm, pass verts and faces to `skimage.measure.mesh_surface_area`. Regarding visualization of algorithm output, to contour a volume named `myvolume` about the level 0.0, using the ``mayavi`` package:: >>> >> from mayavi import mlab >> verts, faces, _, _ = marching_cubes(myvolume, 0.0) >> mlab.triangular_mesh([vert[0] for vert in verts], [vert[1] for vert in verts], [vert[2] for vert in verts], faces) >> mlab.show() Similarly using the ``visvis`` package:: >>> >> import visvis as vv >> verts, faces, normals, values = marching_cubes(myvolume, 0.0) >> vv.mesh(np.fliplr(verts), faces, normals, values) >> vv.use().Run() To reduce the number of triangles in the mesh for better performance, see this `example `_ using the ``mayavi`` package. References ---------- .. [1] Thomas Lewiner, Helio Lopes, Antonio Wilson Vieira and Geovan Tavares. Efficient implementation of Marching Cubes' cases with topological guarantees. Journal of Graphics Tools 8(2) pp. 1-15 (december 2003). :DOI:`10.1080/10867651.2003.10487582` .. [2] Lorensen, William and Harvey E. Cline. Marching Cubes: A High Resolution 3D Surface Construction Algorithm. Computer Graphics (SIGGRAPH 87 Proceedings) 21(4) July 1987, p. 163-170). :DOI:`10.1145/37401.37422` FlorensenTrz/method should be either 'lewiner' or 'lorensen') use_classicr) ValueError_marching_cubes_lewiner) volumelevelr r r r r rrs g/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/skimage/measure/_marching_cubes_lewiner.pymarching_cubesr sezK  9  JKKK "      cDt|tjr |jdkrt d|jddks"|jddks|jddkrt dtj|tj}|-d|| zz}nNt|}||ks|| krt d t|dkrt d t|}|dkrt d t|}t}||j|jkst d tj} | ||||||\} } } } t| st#d tj| } tj| } d| _|dkrtj| } n|dkst d|dtj|ds| tj|z} |r| | | | fStj}|| tj| | | S)zdLewiner et al. algorithm for marching cubes. See marching_cubes_lewiner for documentation. z(Input volume should be a 3D numpy array.rrz#Input array must be at least 2x2x2.N?z/Surface level must be within volume data range.z'`spacing` must consist of three floats.zstep_size must be at least one.z)volume and mask must have the same shape.z(No surface found at the given iso value.)rrascentzIncorrect input z( in `gradient_direction`, see docstring.)rrr) isinstancenpndarrayndimrshapeascontiguousarrayfloat32minmaxfloatlenintbool _get_mc_lutsrr RuntimeErrorfliplr array_equalr_remove_degenerate_facesastype)rrr r r r rrLfuncverticesfacesnormalsvaluesfuns rrrs fbj ) )EfkQ.>.>CDDD |Afl1o11V\!_q5H5H>???  !&"* 5 5F }vzz||fjjll23e  6::<<  56::<<#7#7NOO O 7||qBCCCII1}}:;;;{##K A zV\))HII I & 4D'+tq)[$(($HeWf x==GEFFFy""Hi  GEKY&& %   8 + + 1       >'9 - --beGn,H//(@s8??2:..wGGGrc|\}}tj|d}tj|d}||_|S)Nzutf-8int8)dtype)base64 decodebytesencoder frombufferr#)argsr#textbytsars r _to_arrayrEsGKE4  dkk'22 3 3D t6 * * *BBH Irr;c ttdstjgtt t ttjttj ttj ttj ttj ttj ttjttjttjttjttjttjttjttjttjttjttjttjttjttjttjttjttjttjttj ttj!ttj"ttj#ttj$ttj%ttj&ttj'ttj(ttj)ttj*ttj+ttj,ttj-ttj.ttj/ttj0ttj1ttj2ttj3ttj4ttj5ttj6ttj7Rt_8tj8S)z#Kind of lazy obtaining of the luts.THE_LUTS)9hasattrmclutsr LutProviderEDGETORELATIVEPOSXEDGETORELATIVEPOSYEDGETORELATIVEPOSZrE CASESCLASSICCASESTILING1TILING2 TILING3_1 TILING3_2 TILING4_1 TILING4_2TILING5 TILING6_1_1 TILING6_1_2 TILING6_2 TILING7_1 TILING7_2 TILING7_3 TILING7_4_1 TILING7_4_2TILING8TILING9 TILING10_1_1 TILING10_1_1_ TILING10_1_2 TILING10_2 TILING10_2_TILING11 TILING12_1_1 TILING12_1_1_ TILING12_1_2 TILING12_2 TILING12_2_ TILING13_1 TILING13_1_ TILING13_2 TILING13_2_ TILING13_3 TILING13_3_ TILING13_4 TILING13_5_1 TILING13_5_2TILING14TEST3TEST4TEST6TEST7TEST10TEST12TEST13 SUBCONFIG13rGrrr,r,s[ 6: & &5 4@4 4 4  4  f) * * 4 fl # # 4 fn % % 4  fn % %4  f& ' '4  f& ' '4  f& ' '4  f& ' '4  fn % %4  f( ) )4  f( ) )4  f& ' '4 f& ' '!4 " f& ' '#4 $ f& ' '%4 & f( ) )'4 ( f( ) ))4 * fn % %+4 , fn % %-4 . f) * */4 0 f* + +14 2 f) * *34 4 f' ( (54 6 f( ) )74 8 fo & &94 : f) * *;4 < f* + +=4 > f) * *?4 @ f' ( (A4 B f( ) )C4 D f' ( (E4 F f( ) )G4 H f' ( (I4 J f( ) )K4 L f' ( (M4 N f( ) )O4 P f' ( (Q4 R f) * *S4 T f) * *U4 V fo & &W4 X fl # #Y4 Z fl # #[4 \ fl # #]4 ^ fl # #_4 ` fm $ $a4 b fm $ $c4 d fm $ $e4 f f( ) )g4 4 4 l ?rc||}|dddddf|dddddfz }|dddddf|dddddfz }~tj||dzddzdz S)a=Compute surface area, given vertices and triangular faces. Parameters ---------- verts : (V, 3) array of floats Array containing coordinates for V unique mesh vertices. faces : (F, 3) array of ints List of length-3 lists of integers, referencing vertex coordinates as provided in `verts`. Returns ------- area : float Surface area of mesh. Units now [coordinate units] ** 2. Notes ----- The arguments expected by this function are the first two outputs from `skimage.measure.marching_cubes`. For unit correct output, ensure correct `spacing` was passed to `skimage.measure.marching_cubes`. This algorithm works properly only if the ``faces`` provided are all triangles. See Also -------- skimage.measure.marching_cubes Nrrr)axisrg@)r crosssum)vertsr6 actual_vertsabs rmesh_surface_arear;s>rs 444444(((((( K    KKKKK\OHOHOHd RX1qeQqE1Q%!A!uaUAa5STUVRWYZ[\X]_`ab^cefghdikmsttRX1qeQqE1Q%!A!uaUAa5STUVRWYZ[\X]_`ab^cefghdikmsttRX1qeQqE1Q%!A!uaUAa5STUVRWYZ[\X]_`ab^cefghdikmstt999x%B%B%B%B%Br