From inside to outside: CT scanning as a tool to link internal knot structure and external branch diameter as a prerequisite for quality assessment

Branches are as essential for tree growth as knots are detrimental from the wood quality point of view. To bridge the gap between tree growth and the quality toward end-use, this study aims to establish a relationship between internal and external diameters of Douglas-fir whorl branches. The data comprised 102 trees of a wide age range (30–80 years old) from nine study sites in Southwest Germany. External branch measurements were performed in the field following an established protocol. Logs were scanned on a MiCROTEC CT.LOG, and knots were detected by applying an automated algorithm. Obvious detection artefacts by the CT algorithm were excluded to reveal the relationship between inner-outer branch diameters as clear as possible. Results showed a significant mean difference of 13.8 (± 10.0) mm between the methods (external diameter being larger), with a model indicating an offset of 9.75 mm and angular shift of 0.53 (RMSE = 7.12 mm; R2 = 0.57) between the methods. Separate calculations of sound and dead datasets did not reveal a statistically significant difference. By linking the internal knot structure to external branch measurements, the findings of this study constitute a first step toward the incorporation of CT data into growth models, providing a meaningful prediction of the maximum internal knot diameter at an early stage in the wood supply chain.


Introduction
A knot is a wood feature that can greatly affect the value recovery of logs, in particular, Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) logs, as the species is known for its high vigor and growth rate, which often translates into large knots (Hein et al. 2008b). The largest branches in Douglas-fir are usually observed in whorls, which is a group of knots that originates at approximately the same height in the tree, with each knot in a different angular position. From the wood production perspective, knots represent a discontinuity in the wood matrix (Grace et al. 2006), causing a reduction in mechanical properties (e.g., strength, stiffness) (Gartner 2005), as well as affecting the visual appearance of the final products (Nyrud et al. 2008; Barbour and Parry 2001;Gartner 2005).
In terms of the wood production chain, foresters are responsible for planting, managing, felling and, sometimes the bucking of trees, ultimately generating roundwood as their final product. By aiming to achieve higher log prices, foresters have, therefore, an interest in acquiring information on the effect of their stand treatment strategies and concepts on the inner quality of logs. To standardize the roundwood quality parameters for trading purposes, log sorting agreements have been established. Regarding knots, these agreements are based on the premise that branches visible from the outside of a log are related to the effect of their inner part (knot) on the sawn timber quality. The European Standard for log sorting (CEN 2008) considers the diameter of the largest knot as a quality criterion, which is measured after debranching at the log surface (surface knot diameter-SKD). The standard differentiates knot diameter thresholds between sound and dead (unsound, loose) knots. For log trades carried out in Germany, the RVR agreement (Anonymous 2015) is commonly used, which has a similar structure to the European standard but applies more restrictive thresholds with regards to allowed dead knot size.
Sawn timber quality is ultimately determined by tree internal characteristics, which cannot be exactly described from outside the log. Sawmillers, in turn, consider these rules when buying the roundwood to convert it to various sawn timber products (e.g., boards, blocks, veneers, etc.). Therefore, predicting internal log quality as early as possible within the forest-wood industry chain is of common interest to both log producers and log consumers, as the aim of both is primarily the same: to efficiently produce and sell good quality wood products. The quality aspect is key here, as it is one of the factors that will influence the product value in different stages of the production chain (Gartner 2005).
Growth models have covered part of the solution for an integration between forest and timber products, by exploring the causal effect of silvicultural treatments on external branch diameter (Väisänen et al. 1989;Hein et al. 2008b;Garber and Maguire 2005;Grace et al. 2015) or directly on saw timber properties (Rais et al. 2014;Högberg et al. 2010), as well as relating branch development to tree measurements (Duchateau et al. 2015) and stand structural data (Mäkinen and Colin 1998). Ultimately, the incorporation of such models into growth simulators (Yue et al. 2013;Dufour-Kowalski et al. 2012) led to its output being linked to product recovery models (Houllier et al. 1995). Tree and stand parameters have also been related to outputs from sawing simulations (Ikonen et al. 2009;Nordmark 2005;Weiskittel et al. 2006), aiming to analyze how silvicultural interventions might affect the final sawn timber quality under several breakdown scenarios. These relationships have mainly covered the effects of forest management on crown aspects, more precisely up to the branch diameter outside the stem (Hein et al. 2008a(Hein et al. , 2b, 2009. Other studies focused on knot size issues with regards to processing in sawmills (Fredriksson 2014), relating it to the end product, either through measurements (Krajnc et al. 2019) or simulated scenarios (Mäkelä et al. 2010). There is, however, a knowledge gap on the relationship between inner and outer branch diameter on the stem/roundwood/log level, due to the different methods required to obtain such measurement information. Such a relationship could pave the way for more robust predictions of internal sawn-timber-oriented parameters from tree or stand measurements.
The acquisition of reliable inner-log measurements has been made possible with advances in non-destructive technologies, of which computed tomography (CT) has been one of the most successful techniques to obtain information on the position and size of internal features of a log (Schad et al. 1996;Longuetaud et al. 2012). Several studies have initially applied the technology to extract internal features, including the knot geometry (Grundberg 1994;Funt and Bryant 1987), which evolved into more advanced algorithms capable of extracting knots automatically Johansson et al. 2013;Longuetaud et al. 2012;Roussel et al. 2014), ultimately leading to sawmill applications, more precisely regarding optimization in terms of log rotation before cutting Fredriksson 2014;Fredriksson et al. 2013;Johansson 2013;Belley et al. 2019), cutting patterns (Ursella et al. 2018) and control of the process (Grundberg 1999).
External branch, respectively, internal knot size as derived from CT scanning is a highly relevant quality indicator both for roundwood and sawn timber. However, CT scanning is primarily a stationary machine in the industry. Despite a study pilot-testing the idea of using the CT technique to decide on log grade for payment (Oja et al. 2010), it is mainly used only after logs have been bought and paid for. CT technology itself does not solve the problem of early prediction of internal features in the forest. In the present approach, CT scanning is used as a methodological tool to obtain internal knot diameter information prior to processing. This way, the internal knot measures are related to external branch size, a measure applied in the forest for log grading. In contrast to recent grading procedures, i.e., EN 1927-3 (CEN 2008 or RVR (Anonymous 2015), which consider only one (i.e., the largest) branch per log for assessment, in this study a higher measurement resolution was chosen, taking into account the largest branch of each whorl of a given log or stem. As a first step toward building a link between external and internal branch characteristics of Douglas-fir trees grown in Southwest Germany, this study aims to establish a quantitative relationship between early measurements of branches from outside the log (on felled trees, with possible application to standing trees in the future) and the inner knots in a log, which are relevant for the quality of processed sawn timber. In addition, this study aimed to analyze such a relationship further by investigating possible differences between sound and dead branches. Then, as a second step, based on such a relationship, it will be possible to predict the internal knot size based on parameters externally measured on whorl or tree level, which will be reported in a subsequent publication.

Sample trees
Data from 102 trees cross-cut to 295 logs, sampled in four earlier scientific projects on Douglas-fir, were used in this study. Within these projects, trees from nine study sites in Southwest Germany, including three age intervals, were sampled. Logs were obtained in 4-5 m sections of the tree merchantable length (up to a minimum log diameter of 7 cm). Detailed information on each project's motivation, goals and methodology can be found in the literature (Abetz 1971;Ehring 2006;Kenk and Hradetzky 1984;Kenk and Thren 1984;Kohnle et al. 2012;Šeho et al. 2013;Šeho and Kohnle 2014). The characteristics of these datasets relevant to this study are summarized in Table 1.

Data structure
This study considered different datasets (Fig. 1). Given that the largest branch per log is crucial for roundwood sorting rules (SKD as base measurement) and that branches usually present larger diameters in whorls than in internodes, this investigation was focused on the maximum branch diameter of the largest branch in a whorl. The CT data available was filtered first to identify only whorl knots, and in a second filtering step, the largest knot per whorl was selected, resulting in the dataset Raw. Subsequently, to ensure a reliable dataset for the aforementioned comparison, we visually inspected the CT outputs by comparing the CT images with the automated CT detected knot data. Obvious cases of misdetection between detected knots and the underlying CT image were eliminated, resulting in a subset of the previous dataset: Inspected. To evaluate possible differences in this relation due to knot status (sound and dead knots), the Inspected dataset was split according to this factor, resulting in the datasets Sound and Dead. Knot detection as mentioned in the previous procedure was completely automated. Therefore, it was validated against manual measurements on CT images, to provide criteria of precision and accuracy. The aim here was to verify how the knot detection algorithm performed in providing the maximum knot diameter in a given whorl. For this validation analysis, we selected only whorl knots from the CT data and compared them with the respective knot data from manual measurements on CT images, generating the Validation CT dataset (Table 2; Fig. 1).

Diameter of the largest branch per whorl (field data)
The technical preparation of the trees composing the material of this study followed an internal protocol, based on which they were felled and bucked into 3-5-m-length logs. The debranching of these trees was delayed until all measurements were finished. Groups of at least three branches that occurred at the same height in the stem were identified as whorls. Moreover, the distance between whorls and the counting of annual rings in both crosscut faces of each log were considered as a verification of the whorls. Description of the data structure. Colours represent the origin of the data being considered (red, blue and yellow depict, respectively, the field measurements, CT automated output and manual measurements on CT images; purple and green show, respectively, steps that considered the first two and the latter two data inputs. Bold terms indicate the name of the generated dataset. *Maximum knot diameter here refers to the largest diameter found in the radial direction, between stem pith and bark, for each largest knot in a whorl (colour figure online) The diameter of the largest branch per whorl was determined as the mean of two crossed diameter measurements perpendicular to the branch axis, taken with a caliper (± 0.5 mm). To avoid overestimation of this variable, the exact measurement position (axial distance from the stem surface) varied, since we tried to minimize the effect of the branch collar (Shigo 1985). The whorl longitudinal position along the stem was recorded as the distance to the ground of the lowest branch in this whorl, according to a specific protocol (Brüchert et al. 2017). Whorls that only consisted of (overgrown) branch scars, or in case the largest branch was broken, were counted but not measured.

Diameter of the largest knot per whorl (CT data)
Logs were transported to the Forest Research Institute of Baden-Württemberg (FVA), where they were scanned using the MiCROTEC CT.LOG scanner (Giudiceandrea et al. 2011). The scanning process generated sinograms that were converted into stacked grey-level images. Each image presented a cross-cut view with a pixel size of 1.107 × 1.107 mm representing 5 mm of the log length. By positioning a group of such consecutive images in a row, a 3D virtual log can be generated.
The internal knot structure was derived by processing these stacked images. This step was automatically carried out using the following established algorithms for the detection of: (1) the stem pith (Boukadida et al. 2012); (2) borders, such as heartwood-sapwood, wood-bark and outer border (Longuetaud et al. 2007;Baumgartner et al. 2010); and (3) knots . The algorithm configurations were set as defined by Longo et al. (2019a) in a study in which the authors presented knot detection filter configurations adjusted for Douglas-fir, as well as accuracy results.
The knot detection algorithm fits ellipses (Fig. 2B) to roughly round-shaped areas with distinctly higher values (i.e., lighter shades), which correspond to knots  Fig. 2A). The ellipse fitting process is performed radially in relation to the stem pith, in ten positions along the pith-bark gradient. Due to better image contrast between knots and the surrounding wood, the algorithm was designed to locate at least five positions within the heartwood area. Ultimately, based on ten ellipses detected along the stem radius, the algorithm fits a knot diameter model (Grönlund et al. 1995) with a unique parameterization for each knot (Eq. 1).
where φ(r) is the arc (in rad) corresponding to the knot aperture at a given r radial position from the stem pith (in mm); A and B are parameters generated uniquely for each knot based on the algorithm's knot detection. This model generates the arc corresponding to the knot diameter aperture at a given radial position. To acquire the linear knot diameter, i.e., the smaller distance between two points (both extremities of the knot aperture), a transformation of the variable was necessary (Eq. 2).
(1) where D(r) is the linear knot diameter at the r radial position (in mm) from the stem pith (in mm). The maximum diameter of sound knots extracted by the algorithm is at the outermost radial position within the stem, i.e., the position right under the bark. In this situation, the corresponding branches are still functional; thus, their size is expected to grow further. Dead knots, in contrast, show the peak of their diameter development inside the stem: their maximum diameter is considered within the algorithm as the point after which the knot size is expected either to directly decrease or to remain constant and then to decrease. In cases where dead knot diameters remain constant, the innermost maximum value is considered.

Validation of the CT knot detection algorithm
Although Longo et al. (2019a) have carried out a validation of the knot detection algorithm on Douglas-fir knots in general, a further validation was performed in this study focusing on the accuracy of the algorithm in detecting the maximum knot diameter of whorl knots (large or small). The material for this validation consisted of 22 logs selected randomly out of 295 logs available for this study. In the CT images of these logs, all whorls were visually identified, and the vertical maximum diameter of each knot belonging to a given whorl was manually measured. In total, the maximum diameter of 881 manually measured knots was paired with the CT automated output (Validation CT dataset), identified according to their longitudinal and azimuthal positions.

Knot-branch measurement pairing
The pairing of internal and external measurements was performed in two steps: (1) whorls were matched based on their longitudinal position (± 10 cm); (2) the largest knot detected per whorl was assigned to the largest branch measurement. Due to overestimation and underestimation from the knot detection algorithm (Longo et al. 2019a), the largest knot detected by the algorithm might not always be the largest knot in reality. Therefore, all whorls were verified in the CT images after the automatic knot detection to ensure that the CT data contained the largest knot (pairing step in Fig. 1). The verification was mainly visual, but in cases where two or more knots presented visually similar sizes, a measurement of each doubtful knot was performed (on the CT image) at its outermost position under bark. The largest knot was then manually assigned to that particular whorl (if it was not already the case). In cases where no matching between CT and field datasets was possible (due to, e.g., the largest knot being undetected or not measured in the field, height position differences higher than 10 cm), the whorl in question was ignored in the analysis. These three steps generated the first paired dataset (Raw dataset, Fig. 1).

Visual inspection of CT images
In a further step, the aim was to check for the presence of technical artefacts. Since the field data cannot be re-measured, a visual inspection was conducted only on CT data, aiming to check for detection inconsistencies (incomplete detection, and knot sizes and shapes beyond plausibility). As a result of this inspection, errors were identified as a consequence of natural occurrences (crack through the stem or through the knot, group of knots, subsequent knots, steep knots, extrapolated wood-bark border and large occluded knots). Implausible cases were identified when the detection was clearly not representing the respective knot (shape, size or orientation) at the peak of its internal development (i.e., maximum diameter position). The overall knot detection could have been very good, but if the detected knot clearly did not reflect the actual knot at the maximum diameter point, this knot was also removed during the inspection. Thus, only knots that were not affected by the aforementioned inconsistencies remained in the analysis, referred to as Inspected dataset, Fig. 1.

Data analysis
In the validation of the CT algorithm, a detection rate and the number of falsely detected instances (non-existent but detected as "knots") were derived by comparing manually measured knots with knots detected by the algorithm (Validation CT). The detection rate was calculated by dividing the number of correctly detected knots by the total number of knots. In addition, the number of falsely detected knots was recorded. The comparison between the measurements from both methods was described by the mean absolute error (MAE), standard deviation of the errors (SDE), the Pearson correlation coefficient (ρ) and the 95% limits of agreement (LAE, calculated as MAE ± 1.96*SDE) according to Bland and Altman (2003). The latter returns a probabilistic knot diameter error range. In addition, the concordance correlation coefficient (ρ c ) was calculated (Lin 1989(Lin , 2000, which provides a measure of how precise the data are compared to a perfect match scenario between the two methods. Lin's concordance correlation coefficient varies from − 1 to 1 and quantifies the distance from the reference line (x = y). According to the author, such a coefficient can be further explored by analyzing the scale (v) and location (u) shift parameters, which indicate the weight of each type of shift from the reference line, as they diverge from 1 and 0, respectively. In addition, the analysis of the C b parameter reveals how different the best-fit line is in comparison with the reference x = y line as it deviates from 1. The C b parameter is a measure of accuracy and is calculated as a correction factor between ρ c and ρ, i.e., the quotient between ρ c and ρ (Lin 1989(Lin , 2000. In the inner-outer knot/branch diameter analysis, a summary of the descriptive statistics was produced for each dataset (Raw, Inspected, Sound and Dead). To examine whether the knot diameter difference between field and CT dataset was statistically significant (D Field -D CT ), a paired t test with 0.05 of significance threshold was calculated for each of the four datasets. In addition, to establish how these two variables relate to each other a linear model was fitted (Eq. 3).
where D CT is the internal knot diameter automatically detected by the CT algorithm; D Field is the external branch diameter measured in the field; estimates b 0 and b 1 are the intercept and slope of the model, respectively; and ε i represents the information unexplained by the previous components, following a normal distribution with mean 0 and variance σ 2 .
For the comparison of the parametric mean differences between sound and dead knot groups, the Welch's two sample t test was applied. According to Delacre et al. (2017), the Welch t test should be preferred (over, e.g., the Student's t test) independently of the samples' variances assumption, for cases where the samples being tested follow a normal distribution but have a different number of observations (unpaired samples).
Data management and analyzes throughout this study were performed in R (R Core Team 2016), using the following additional packages: epiR (Stevenson et al. 2018

Validation of the CT algorithm
The CT automated algorithm detected 881 whorl knots, from a total of 1016 whorl knots that were identified in the CT images (Validation CT dataset). In addition, 33 falsely detected knots were observed in whorl areas. In this dataset, the manual measurement (on CT images) method had a mean knot diameter of 10.3 mm (standard deviation = 5.91 mm), while the automated CT method found a mean value of 10.8 mm (standard deviation = 5.6 mm). The mean absolute error between the two diameter measurements was 0.48 mm (standard deviation of 2.70 mm). The results presented good precision (ρ = 0.89), with slight shift from the reference line, as shown by both scale (v = 0.94) and location (u = 0.08) parameters. These parameters show that the best-fit line is unbiased and not shifted in comparison with the x = y reference line. This behavior can be seen in Fig. 3A, which illustrates how measurements from both methods compare. In addition, ρ c (0.89; with confidence interval of 0.87-0.90) and C b (0.99) corroborate the performance observed, as they approach 1. The limits of agreement were − 5.77 and 4.81 mm, as presented in Fig. 3B. This plot was recommended by Altman and Bland (1983) as a tool to avoid letting the range of measurements influence the appearance of graphic agreement, since it shows the difference between the methods against the mean between them.
The knot detection algorithm identified 86.7% of all whorl knots on CT images on the Validation CT dataset. In a previous study (Longo et al. 2019a), a validation of the same knot detection algorithm was performed on 15 Douglas-fir logs, in which several measurements along the radial gradient of whorl knots and internodes were accounted for. That study observed a higher detection rate of 93.9%. With regards to the algorithm's accuracy, the aforementioned study showed mean diameters of 22.4 and 23.7 mm, respectively, for manual and CT methods, resulting in a MAE of 1.34 mm (manual-CT), SDE of 7.72 mm and a ρ c of 0.82. The differences observed between the results of the current study and the ones presented in Longo et al. (2019a) are mainly due to the aim of each study and the data considered in each investigation. The aforementioned study used old Douglas-fir trees (approx. age of 78 years), covering knots in general, accounting for as much small knots as larger ones, in whorls and internodal ones, and with measurements taken along the knot length, providing an overview of the overall accuracy of the algorithm. In turn, the current study utilizes a broader range of tree age, focuses on whorl areas, and considers only the maximum diameter of such knots.

Comparison between external branch measurements and CT knot detection results
The analysis of the Raw dataset comprised 1198 matched pairs of internal knot and external branch diameter (Table 3). The descriptive statistics of the dataset revealed a statistically significant (p < 0.05, Table 4) mean difference of 23.7 mm (± 13.9 mm). The linear model showed a low capability to explain the variability present in this dataset (R 2 = 0.23). Figure 4A indicated possible agreement issues, in which points showed suspicious differences between the two diameters. A slight tendency could be seen in Fig. 4C, between the residuals to the reference line and the external branch diameter measurements. The Inspected dataset accounted for 288 knot pairs (from 84 trees and 156 logs) and was used to analyze the relation   . 4 Results of model fitting on the relation between CT and field measurements before visual inspection on CT images. Sample size: 1198 knots. Black lines represent reference lines (in A is the x = y relation, and in B is the zero deviation of the residuals). The blue line depicts the linear model fitted to the data and the shading area surrounding it represents the 95% confidence interval of the mean model. Regarding colors, the reader is referred to the digital version of this article (colour figure online) between knot and branch measurements (Fig. 5A). It showed a statistically significant difference between internal and external diameters according to the paired t test (p < 0.05), with no clear pattern in the residuals (Fig. 5B). Regarding the model depicted in Fig. 5A, the intercept value indicates that external measurements are 9.75 mm larger than the CT internal diameter (Table 5). The model's angular coefficient (b1 = 0.53) indicates that the relationship has a positive trend, but also that the difference between external branch diameter and the respective predicted internal diameter increases at a lower rate as branches increase in size.
The present results showed that the relation between the external branch diameter measurements performed in the field and the maximum internal knot diameter acquired by means of CT technology is linear. In Fig. 5A, the model line intercepts the reference (x = y) line at a branch diameter of 20.74 mm, indicating that below such a point, the model would estimate negative values of internal knot diameter. This was likely due to the fact that this range of data (small knots; below 20 mm of diameter) was not sampled in this study, as only the largest branch per whorl was considered. Pyörälä et al. (2018) compared external and internal knot diameter  1 3 measurements of Scots pine (Pinus sylvestris L.) using Terrestrial Laser Scanning (TLS) of standing trees and X-ray scanning of the corresponding logs, respectively. The authors found a significant mean difference of 6.5 mm for maximum knot diameter between the methods (TLS-X-ray). Such a difference was attributed to the internal knot diameter acquisition method (single-directional X-ray digital radiograph scanner) and to the estimation of the maximum knot diameter per whorl, which was based on "the length of the whorl along the stems longitudinal axis in the X-ray images, which is sensitive to noise and overlapping knots". Furthermore, the authors mentioned that the external measurements are affected by wind, which likely causes an increasing registration error on the TLS point cloud toward the treetops. In the present study, the observed differences between the measurement results might be explained by the distinct measurement positions along the axis of the branch and the branch bark thickness present in the external measurement. Differences are also found in the knot definition, i.e., the outlining rule that defines the boundaries of a knot, when examining knots internally and externally applying distinct methods. External branch measurements contribute to deviations in this inner-outer knot/branch diameter relation due to the presence of bark, characterizing the field measurement as a diameter with double bark thickness; branch collar (Shigo 1985), measurement errors, and considering that external branch diameter is inversely related to the position in branch length (Fernández and Norero 2006), to inaccuracies of measurement positions. In turn, internal knot measurements in the log periphery using CT techniques present poor knot-stem wood matrix contrast in the sapwood region (Funt and Bryant 1987;Longuetaud et al. 2012;Wei et al. 2009;Breinig et al. 2012;Fredriksson et al. 2017;Johansson et al. 2013). As a consequence, the increment in knot diameter in this region might be overlooked, due to smooth boundaries (especially in sound knot cases). Moreover, we observed that the CT method is influenced by technical aspects (image resolution, quality and contrast in different areas), as well as the biological configuration of the material (knot distribution and positioning). Furthermore, the visual identification of whorls in the field or when analyzing CT images can pose a challenge when analyzing paired data from different methods. In the field, the number of whorls can be roughly verified by comparing the annual rings at the top and bottom positions of the stem section and consequently comparing the difference between those to the number of whorls in that section. On the other hand, visually identifying whorls in large CT datasets is very onerous, time consuming and susceptible to visual errors. However, the density distribution along the stem height provided by analyzing CT images might reveal whorl footprints (Longo et al. 2019b;Longuetaud et al. 2005), which could be used for automated whorl detection, reducing errors related to whorl interpretation.
As the extraction of sound and dead knots' diameter in the CT detection follows different approaches, the Inspected dataset was split in two subgroups representing the different knot status. Figure 6 presents the visualization of the results for these two knot groups. The linear tendency was maintained, with field measurements presenting mean values of 14.2 (± 10.1) and 13.6 (± 10.0) mm higher than the CT knot diameter for sound and dead knots, respectively. Although the goodness of fit for sound knots is higher than for dead knots, the spread of data in both groups is very similar. According to the Welch test, the disparity of the mean difference between the groups is not statistically significant (t = 0.51; df = 203.6; p-value = 0.61). In addition, in Fig. 5A, the point in which the model line intercepts the reference line should split sound and dead knot groups. However, as also observed in Table 4 and Fig. 6, this is not the case. Negative minimum values in Table 4, as well as the region lower than the intercept between the lines in Figs. 5A, 6A1, A2 mean that we observed larger D CT than its respective external counterpart D Field at times. It was expected that inner-outer branch diameter differences for dead knots would gravitate around zero or even negative values (D Field -D CT ). This was based on the knot Fig. 6 Results of model fitting on the relation between CT and field measurements after visual inspection on CT images, detailing groups of sound (A) and dead (B) knots. Sample size: 102 (sound) and 186 (dead) knots. Black lines represent reference lines (in 1 is the x = y relation, in 2 is the zero deviation of the residuals). The blue line depicts the linear model fitted to the data, and the shaded area surrounding it represents the 95% confidence interval of the mean model. Regarding colors, the reader is referred to the digital version of this article (colour figure online) status definition within the applied CT algorithm , in which a knot is considered dead when it has reached its development peak within the stem, i.e., the maximum knot diameter found somewhere between pith and bark. Theoretically, after this point the knot stops growing, and its diameter might be constant for some time or directly decrease. However, this was not the main observed behavior in the present study. There were a few dead knots that would confirm this pattern, but the majority showed differences similar to a sound knot instance, i.e., further growth. This behavior could be due to thicker branch bark in dead branches, which was not noticed in the field. The most probable cause is a misclassification from the algorithm, identifying sound knots as dead ones, as found in a previous study (Longo et al. 2019a). Only considering the maximum knot diameter might not be sufficient to accurately identify the sound-dead border of a knot. In fact, evidence was found in a study on Douglas-fir branch radial growth that after the maximum diameter has been reached, the knot might still be sound for an average of 8 (± 4.7) years (Kershaw et al. 1990), presenting no differences in diameter. Therefore, we understand that there is an uncertainty regarding the sound-dead boundary of a knot when applying CT images to identify the status of a knot, since deviations in knot detection might lead to a false maximum diameter point. In this context, further studies on branch mortality and its identification on material grown and adapted in the region might lead to improvements of the sound-dead border detection within the algorithm.

Conclusion
This was, to the authors` knowledge, the first study on Douglas-fir logs to relate external whorl branch diameters from field measurements and internal knot diameters using a CT automated knot detection algorithm. A linear relation was evidenced after CT detection inconsistencies were excluded. No difference in the inner-outer branch diameter relation between sound and dead knot data was found to be significant. The findings of this study must be considered along with the limitations of the methods. Probable influencing factors identified were related to the automated CT detection (whorl identification, steep knots, sapwood effect, smooth knot boundaries in sound knots) and to the field data acquisition (the presence of bark, branch collar and measurement errors). Developments in the knot detection process (pith, borders and knot detection algorithms) might improve the accuracy of measurements acquired from CT data.
As the inner maximum knot diameter relates to the knot morphology more precisely than the external branch diameter, this information is also realistically closer to the branch diameters in the final timber products (Duchateau et al. 2013). Future work should concentrate on predicting the internal knot diameter based on external measurements obtained in the forest, which could be potentially either measured at the log surface (SKD) or remotely on standing trees (e.g., with photo-optical or laser sensing technology) in the course of inventory measurements in the field. As we understand the relationships between different stages of this process (remote sensing, felled tree measurements, internal measurements), we will likely unlock possibilities, such as: enabling a more reliable roundwood quality assessment (oriented toward the end-user) via model predictions and providing more accurate internal knot metrics when refining existing forest growth models to include a prediction of wood quality at an early point in the wood supply chain.