Skip to main content

graphics3d - How to Discretize a Tube


In Is there a way to outline text?, I was amazed that BoundaryDiscretizeGraphics can discretize a Text primitive. So in a similar way, I wanted to apply DiscretizeGraphics to a Tube like this:


DiscretizeGraphics @ Graphics3D[Tube[{{0, 0, 0}, {1, 1, 1}}, .1]]

or


tube = Tube[{{0, 0, 0}, {1, 1, 1}}, .1];
BoundaryDiscretizeGraphics[tube, #] & /@ {_Tube, _Polygon, _Line, _Point}


But I failed to do it. Is it a bug in Mathematica, or is it some mistake in my code?



Answer



It took quite a while, but I've finally come up with a way to generate discretized tubes. This again is based on work in this previous answer (from the thread mentioned earlier by Michael). In the interest of keeping things short, I will not be repeating the definitions of orthogonalDirections[], extend[], and crossSection[] from that answer. Here, then, is the tube generator:


makeCap[type : ("Butt" | "Round" | "Square"), s : (-1 | 1), path_?MatrixQ, tube_?ArrayQ,
        r_?NumericQ, h_?NumericQ] := 
    Module[{d = Take[path, 2 s], cs, p, t0, t1},
           cs = tube[[s]];
           t0 = h; t1 = 1 - h Boole[type =!= "Square"];
           If[s == -1, {t0, t1} = {t1, t0}];
           If[type === "Butt",

              {d[[s]], 
               Table[ScalingTransform[{t, t, t}, d[[s]]] @ cs, {t, t0, t1, s h}]},
              p = (s r/(EuclideanDistance @@ d)) ({1, -1}.d);
              {d[[s]] + p, Switch[type,
               "Round",
               Table[Composition[TranslationTransform[p Cos[Ï€ t/2]], 
                                 ScalingTransform[{1, 1, 1} Sin[Ï€ t/2], d[[s]]]][cs],
                     {t, t0, t1, s h}],
               "Square",
               Table[Composition[TranslationTransform[p],

                                 ScalingTransform[{t, t, t}, d[[s]]]][cs],
                     {t, t0, t1, s h}]]}]]

Options[TubeMesh] = {"CapForm" -> None, "CirclePoints" -> Automatic, 
                     "MeshType" -> Automatic, Tolerance -> Automatic};

TubeMesh[path_?MatrixQ, r_?NumericQ, opts : OptionsPattern[{TubeMesh, MeshRegion}]] := 
    Module[{c0, c1, cf, mt, dims, h, idx, m, n, p0, p1, t0, t1, tol, tube},
           cf = OptionValue["CapForm"]; mt = OptionValue["MeshType"]; 
           If[mt === Automatic, 

              mt = If[MatchQ[cf, "Butt" | "Round" | "Square"], 
                      BoundaryMeshRegion, MeshRegion]];
           tol = OptionValue[Tolerance] /. Automatic -> 0.0015;
           n = OptionValue["CirclePoints"]; 
           If[n === Automatic, n = Round[17 tol^(-1/3)/5 - 57 tol^(1/3)/59]; 
              n += Boole[OddQ[n]]]; h = 2/n;
           tube = FoldList[Function[{p, t},
                                    extend[p, t[[2]], t[[2]] - t[[1]],
                                           orthogonalDirections[t]]], 
                           crossSection[path, r, n], Partition[path, 3, 1, {1, 2}, {}]];

           If[MatchQ[cf, "Butt" | "Round" | "Square"],
              {p0, c0} = makeCap[cf, 1, path, tube, r, h];
              {p1, c1} = makeCap[cf, -1, path, tube, r, h];
              tube = Join[c0, tube, c1]];
           dims = Most[Dimensions[tube]]; tube = Apply[Join, tube];
           m = Times @@ dims; idx = Partition[Range[m], Last[dims]];  t0 = t1 = {};
           If[MatchQ[cf, "Butt" | "Round" | "Square"],
              PrependTo[tube, p0]; AppendTo[tube, p1]; idx += 1;
              t0 = PadLeft[Partition[First[idx], 2, 1], {Automatic, 3}, 1];
              t1 = PadRight[Reverse /@ Partition[Last[idx], 2, 1],

                            {Automatic, 3}, m + 2]];
           mt[tube, Triangle[Join[t0,
                                  Flatten[Apply[{Append[Reverse[#1], Last[#2]], 
                                                 Prepend[#2, First[#1]]} &,
                                  Partition[idx, {2, 2}, {1, 1}], {2}], 2], t1]],
              FilterRules[{opts}, Options[MeshRegion]]]]

As designed, TubeMesh[] will return a MeshRegion[] object for the setting "CapForm" -> None, and a BoundaryMeshRegion[] for the other "CapForm"s.


Using once more the example path in the previous answer:


path = First @ Cases[ParametricPlot3D[

                     BSplineFunction[{{0, 0, 0}, {1, 1, 1}, {2, -1, -1},
                                      {3, 0, 1}, {4, 1, -1}}][u] // Evaluate, {u, 0, 1}, 
                     MaxRecursion -> 1], Line[l_] :> l, ∞];

TubeMesh[path, 1/5, "CapForm" -> "Round", PlotTheme -> "SmoothShading"]

some tube converted into a region


Compute the volume:


Volume[%]
   0.727654081738915

Comments

Post a Comment

Popular posts from this blog

functions - Get leading series expansion term?

Given a function f[x] , I would like to have a function leadingSeries that returns just the leading term in the series around x=0 . For example: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x)] x and leadingSeries[(1/x + 2 + (1 - 1/x^3)/4)/(4 + x)] -(1/(16 x^3)) Is there such a function in Mathematica? Or maybe one can implement it efficiently? EDIT I finally went with the following implementation, based on Carl Woll 's answer: lds[ex_,x_]:=( (ex/.x->(x+O[x]^2))/.SeriesData[U_,Z_,L_List,Mi_,Ma_,De_]:>SeriesData[U,Z,{L[[1]]},Mi,Mi+1,De]//Quiet//Normal) The advantage is, that this one also properly works with functions whose leading term is a constant: lds[Exp[x],x] 1 Answer Update 1 Updated to eliminate SeriesData and to not return additional terms Perhaps you could use: leadingSeries[expr_, x_] := Normal[expr /. x->(x+O[x]^2) /. a_List :> Take[a, 1]] Then for your examples: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x), x] leadingSeries[Exp[x], x] leadingSeries[(1/x + 2 + (1 - 1/x...
Post a Comment
Read more

mathematical optimization - Minimizing using indices, error: Part::pkspec1: The expression cannot be used as a part specification

I want to use Minimize where the variables to minimize are indices pointing into an array. Here a MWE that hopefully shows what my problem is. vars = u@# & /@ Range[3]; cons = Flatten@ { Table[(u[j] != #) & /@ vars[[j + 1 ;; -1]], {j, 1, 3 - 1}], 1 vec1 = {1, 2, 3}; vec2 = {1, 2, 3}; Minimize[{Total@((vec1[[#]] - vec2[[u[#]]])^2 & /@ Range[1, 3]), cons}, vars, Integers] The error I get: Part::pkspec1: The expression u[1] cannot be used as a part specification. >> Answer Ok, it seems that one can get around Mathematica trying to evaluate vec2[[u[1]]] too early by using the function Indexed[vec2,u[1]] . The working MWE would then look like the following: vars = u@# & /@ Range[3]; cons = Flatten@{ Table[(u[j] != #) & /@ vars[[j + 1 ;; -1]], {j, 1, 3 - 1}], 1 vec1 = {1, 2, 3}; vec2 = {1, 2, 3}; NMinimize[ {Total@((vec1[[#]] - Indexed[vec2, u[#]])^2 & /@ R...
Post a Comment
Read more

plotting - Plot 4D data with color as 4th dimension

I have a list of 4D data (x position, y position, amplitude, wavelength). I want to plot x, y, and amplitude on a 3D plot and have the color of the points correspond to the wavelength. I have seen many examples using functions to define color but my wavelength cannot be expressed by an analytic function. Is there a simple way to do this? Answer Here a another possible way to visualize 4D data: data = Flatten[Table[{x, y, x^2 + y^2, Sin[x - y]}, {x, -Pi, Pi,Pi/10}, {y,-Pi,Pi, Pi/10}], 1]; You can use the function Point along with VertexColors . Now the points are places using the first three elements and the color is determined by the fourth. In this case I used Hue, but you can use whatever you prefer. Graphics3D[ Point[data[[All, 1 ;; 3]], VertexColors -> Hue /@ data[[All, 4]]], Axes -> True, BoxRatios -> {1, 1, 1/GoldenRatio}]
Post a Comment
Read more