Bitterlich sampling

Revision as of 20:31, 5 March 2011

Introduction to angle count sampling

Figure 1 Illustration of selection proportional to size (basal area) in Bitterlich sampling (Kleinn 2007^[1]).

Angle count sampling was developed by Walter Bitterlich (1948), an Austrian forester. It is some times also referred to point sampling, horizontal point sampling, variable plot sampling, angle count technique, prism cruising, angle gauge sampling, and simply Bitterlich sampling.

The idea to employ nested sub-plots was introduced because we wished to have a balanced number of trees in all dimension classes; that is, we wanted to assign a higher probability of selection to the larger trees of which there are usually less in a stand (see also sampling with unequal selection probabilities).

For nested sub-plots, we defined a number of fixed plot areas in which the different dimension classes are observed. One may now further pursue that idea and develop a plot design in which the inclusion probability is strictly proportional to dimension. That is each tree has its particular and own plot size (and therefore also its particular inclusion zone). This is exactly what Bitterlich sampling does: from a selected sample point, the neighboring trees are selected strictly proportional to their basal area.

While this sounds complicated, the technique itself is very simple; the only device one needs is one that produced a defined opening angle. That can be a dendrometer, a relascope or simply your own thumb. While standing on the sample point and aiming over e.g. the thumb at the stretched out arm to the dbhs of the surrounding trees; you sweep around 360° and count all trees that appear larger than your thumb. It is obvious then, that larger trees have a larger probability of being taken as sample tree. From this counting alone you can produce then an estimate of basal area per hectare.

The only additional information that you need is the “calibration factor” of your measurement device, as, obviously, the number of trees counted depends on the opening angle which is produced by the instrument (in this case by your arm length and thumb width).

This calibration factor is also called basal area factor. Figure 1 illustrates that approach. Around each tree the individual inclusion zone is drawn. Sample points that fall into this inclusion zone will lead to an inclusion of this particular tree as sample tree. We may also imagine that these specific circular inclusion zones are concentrically placed around the sample point; then we may view them as a set of nested sample sub-plots, one for each sample tree.

The principle of Bitterlich sampling

In a first step of developing this technique we need to calculate the size of these tree‑specific virtual sample plots. For each tree it holds the statement that it is only counted if it appears at 1.3 m larger than the opening angle of the measurement device. For each tree there are exactly three different situations as depicted in Figure 2: (1) either the tree is larger than the opening angle; then the tree is inside the virtual nested sub-plot and is counted; or (2) the tree is smaller than the opening angle (that is: completely covered by the thumb); that is a non-tally tree which is outside the virtual plot; or (3) the tree is exactly covered by the opening angle; then, the tree is exactly on the perimeter line of the virtual circular plot and the distance from the sample point to the center of that tree is (for practical purposes sufficiently good approximation to) the radius of the virtual circular plot. We need to find this radius, because the virtual plot will allow us to calculate the tree-specific expansion factor that we need to expand the per-plot observation to the per-hectare observation.

Inclusion probabilities and estimators for Bitterlich sampling

For the inclusion zone approach where for each tree an inclusion zone is defined, see the corresponding article Infinite population approach. If used,the inclusion probability is then proportional to the size of this inclusion zone – which actually defines the probability that the correspondent tree is included in a sample.

We saw, for example, that angle count sampling (Bitterlich sampling) selects the trees with a probability proportional to their basal area and we emphasized that this fact makes Bitterlich sampling so efficient for basal area estimation. In contrast, point to tree distance sampling, or k-tree sampling, has inclusion zones that do not depend on any individual tree characteristic but only on the spatial arrangement of the neighboring trees; therefore, point-to tree distance sampling is not particularly precise for any tree characteristic.

In Bitterlich sampling, the selection probability of a particular tree i results from the inclusion zone F_i and the size of the reference area, for example the hectare

\[\pi_i = \frac {F_i}{10000}\]

with the Horvitz-Thompson estimator, we have the total

\[\hat \tau = \sum_{i=1}^m \frac {y_i}{\pi_i}\]

for any tree attribute \(y_i\). Applied to estimating basal area \(y_i = g_i = \frac {\pi}{4} d_i^2\) and its per hectare estimation, we have

\[\hat \tau = \sum_{i=1}^m \frac {y_i}{\pi_i} = \sum_{i=1}^m \cfrac {\cfrac {\pi}{4} d_i^2}{\cfrac {F_i}{10000}}\]

and with \( F_i = \pi r_i^2 = \pi c^2 \, d_i^2\), we have the same as Bitterlich sampling

\[\hat \tau = \sum_{i=1}^m \frac {y_i}{\pi_i} = \sum_{i=1}^m \cfrac {\cfrac {\pi}{4} d_i^2}{\cfrac {\pi c^2 \, d_i^2}{10000}} = \frac {2500 \pi}{\pi c^2} \sum_{i=1}^m \frac {d_i^2}{d_i^2} = \frac {2500}{c^2} m\]

which is the estimated basal area per hectare from one sample point where m trees were tallied. The factor 2500/c² is the basal area factor, for details see Bitterlich sampling.

References

1. ↑ Kleinn, C. 2007. Lecture Notes for the Teaching Module Forest Inventory. Department of Forest Inventory and Remote Sensing. Faculty of Forest Science and Forest Ecology, Georg-August-Universität Göttingen. 164 S.