Corrects cylinder radii
Examples
## TreeQSM Processing Chain
file <- system.file("extdata/QSM.mat", package = "rTwig")
qsm <- import_qsm(file)
#> Importing TreeQSM
cylinder <- qsm$cylinder
cylinder <- update_cylinders(cylinder)
#> Updating Cylinder Ordering
#> Calculating Total Children
#> Building Cylinder Network
#> Caching Network
#> Verifying Topology
#> Calculating Growth Length
#> Calculating Reverse Branch Order
#> Calculating Branch Segments
#> Calculating Alternate Branch Numbers
#> Calculating Path Metrics
cylinder <- correct_radii(cylinder, twig_radius = 4.23)
#> Generating Paths
#> Filtering Paths
#> Modeling Paths
#> Updating Radii
str(cylinder)
#> Classes ‘tidytable’, ‘tbl’, ‘data.table’ and 'data.frame': 1136 obs. of 34 variables:
#> $ length : num 0.09392 0.07216 0.06654 0.00938 0.06795 ...
#> $ start.x : num 0.768 0.768 0.768 0.769 0.769 ...
#> $ start.y : num -16.4 -16.4 -16.4 -16.3 -16.3 ...
#> $ start.z : num 254 254 254 254 254 ...
#> $ axis.x : num 0.00995 -0.0111 0.01364 0.01571 0.01449 ...
#> $ axis.y : num 0.0912 0.0391 0.0367 0.0271 0.0267 ...
#> $ axis.z : num 0.996 0.999 0.999 1 1 ...
#> $ end.x : num 0.769 0.768 0.769 0.769 0.77 ...
#> $ end.y : num -16.4 -16.4 -16.3 -16.3 -16.3 ...
#> $ end.z : num 254 254 254 254 254 ...
#> $ added : num 0 0 0 0 0 0 0 0 0 0 ...
#> $ BranchOrder : int 0 0 0 0 0 0 0 0 0 0 ...
#> $ reverseBranchOrder : num 20 20 20 20 20 20 20 20 20 20 ...
#> $ PositionInBranch : int 1 2 3 4 5 6 7 8 9 10 ...
#> $ segment : int 1 1 1 1 1 1 1 1 1 1 ...
#> $ parentSegment : int 0 0 0 0 0 0 0 0 0 0 ...
#> $ mad : num 0.00072 0.000538 0.000523 0.000335 0.000438 ...
#> $ SurfCov : num 0.875 1 1 1 1 1 1 1 1 1 ...
#> $ UnmodRadius : num 0.0465 0.0454 0.0442 0.0437 0.0429 ...
#> $ OldRadius : num 0.0465 0.0454 0.0442 0.0437 0.0429 ...
#> $ growthLength : num 32.8 32.7 32.6 32.6 32.6 ...
#> $ branch : int 1 1 1 1 1 1 1 1 1 1 ...
#> $ branch_alt : int 0 0 0 0 0 0 0 0 0 0 ...
#> $ parent : num 0 1 2 3 4 5 6 7 8 9 ...
#> $ extension : int 1 2 3 4 5 6 7 8 9 10 ...
#> $ totalChildren : int 1 1 1 1 1 1 1 1 1 1 ...
#> $ index : int 1 2 3 4 5 6 7 8 9 10 ...
#> $ distanceFromBase : num 0 0.0722 0.1387 0.1481 0.216 ...
#> $ distanceToTwig : num 3.4 3.31 3.24 3.17 3.16 ...
#> $ reversePipeAreaBranchorder : int 82 82 82 82 82 82 82 82 82 82 ...
#> $ reversePipeRadiusBranchorder: num 9.06 9.06 9.06 9.06 9.06 ...
#> $ vesselVolume : num 279 271 266 260 259 ...
#> $ radius : num 0.0465 0.0454 0.0442 0.0437 0.0429 ...
#> $ modified : num 0 0 0 0 0 0 0 0 0 0 ...
#> - attr(*, ".internal.selfref")=<externalptr>