The VARIOGRAM Procedure

Pair Formation

The basic starting point in computing the empirical semivariance is the enumeration of pairs of points for the spatial data. Figure 102.19 shows the spatial domain D and the set of n measurements , $i=1,\dots ,n$ , that have been sampled at the indicated locations in D. Two data points and , with coordinates $\bm {s}_1=(x_1,y_1)$ and $\bm {s}_2=(x_2,y_2)$ , respectively, are selected for illustration.

A vector, or directed line segment, is drawn between these points. If the length

$\mid P_ iP_ j \mid = \mid \bm {s}_2 - \bm {s}_1 \mid = (x_2-x_1)^2 + (y_2-y_1)^2$

of this vector is smaller than the specified DEPSILON= value, then the pair is excluded from the continuity measure calculations because the two points and are considered to be at zero distance apart (or collocated). Spatial collocation might appear due to different scales in sampling, observations made at the same spatial location at different time instances, and errors in the data sets. PROC VARIOGRAM excludes such pairs from the pairwise distance and semivariance computations because they can cause numeric problems in spatial analysis.

If this pair is not discarded on the basis of collocation, it is then classified—first by orientation of the directed line segment $\bm {s}_2-\bm {s}_1$ , and then by its length $\mid P_ iP_ j \mid$ . For example, it is unlikely for actual data that the distance $\mid P_ iP_ j \mid$ between any pair of data points and located at $\bm {s}_ i$ and $\bm {s}_ j$ , respectively, would exactly satisfy $\mid P_ iP_ j \mid \ = \ \mid \bm {h} \mid \ = h$ in the preceding computation of $\hat{\gamma }_ z(\bm {h})$ . A similar argument can be made for the orientation of the segment $\bm {s}_2-\bm {s}_1$ . Consequently, the pair is placed into an angle and distance class.

The following subsections give more details about the nature of these classifications. You can also find extensive discussions about the size and the number of classes to consider for the computation of the empirical semivariogram.

Figure 102.19: Selection of Points and in Spatial Domain D

Selection of Points P1 and P2 in Spatial Domain D